Skip to main content

Characterization and clinical significance of right ventricular mechanics in pulmonary hypertension evaluated with cardiovascular magnetic resonance feature tracking



Prognosis in pulmonary hypertension (PH) is related to right ventricular (RV) function. Quantification of RV mechanics may offer additive value. The objective of our study is to determine the feasibility and clinical and prognostic value of RV strain analysis by cardiovascular magnetic resonance (CMR) based feature tracking (FT) in PH.


We retrospectively enrolled 116 patients (age 52.2 ± 12 years, 73.6 % women) referred to CMR for PH evaluation who underwent right heart catheterization within 1 month. Using dedicated FT software, peak global longitudinal and circumferential RV strain and strain rates (GLS, GCS, GLSR, and GCSR, respectively) were quantified from standard cine images. Using multivariate regression analysis, we evaluated the associations of strain with a composite endpoint of death, lung transplantation, or functional class deterioration.


RV strain analysis was feasible in 110 (95 %) patients. Patients were classified into: Group A (no PH, normal right ventricular ejection fraction [RVEF]; n = 17), Group B (PH, normal RVEF; n = 26), or Group C (PH, abnormal RVEF; n = 67). All strain and strain rate values were reduced in Group C. Furthermore, GCSR was significantly reduced in Group B (-0.92 [-1.0/-0.7]; p < 0.001) compared to Group A (-1.12 [-1.3/-0.9]; p < 0.001). After adjustment for six clinically meaningful covariates, GLS (hazard ratio 1.06; p = 0.026), GLSR (hazard ratio 2.52; p = 0.04), and GCSR (hazard ratio 4.5; p = 0.01) were independently associated with the composite endpoint. GCSR successfully discriminated patients with and without events (p = 0.01).


Quantification of RV strain with CMR-FT is feasible in the majority of patients, correlates with disease severity, and is independently associated with poor outcomes in PH.


Pulmonary hypertension (PH) is a condition characterized by increased pressure and resistance in the pulmonary vasculature and is associated with high mortality. The prognosis of this disorder is directly related to right ventricular (RV) function, one of the main predictors of long-term outcome in patients with PH irrespective of its etiology [13]. Thus, accurate, reproducible, and clinically meaningful noninvasive methods for the detection and quantification of global or regional RV systolic function are important in the assessment of PH.

Echocardiography is the most widely available imaging modality for the evaluation of RV structure and function; however, image quality of the RV is often inadequate, and quantification can be subjective and limited by the complex geometry of this chamber [4]. Cardiovascular magnetic resonance (CMR) has emerged as the gold standard for the quantification of RV volume and ejection fraction (RVEF), particularly using steady-state free precession (SSFP) cine imaging [5]; however, other metrics of RV function such as quantification of myocardial deformation may allow detection of early abnormalities and provide independent prognostic information, as demonstrated in echocardiographic studies [68].

Recently, a novel method of “feature tracking” (FT) that allows quantification of myocardial deformation from CMR cine images without the need for additional imaging or lengthy analysis has been developed [9, 10]. Of note, CMR-FT has been validated against myocardial tagging for left ventricular (LV) strain analysis [11, 12] and recently against RV speckle tracking imaging for RV longitudinal strain evaluation in tetralogy of Fallot patients [13].

To the best of our knowledge, the significance of quantifying RV myocardial deformation in patients with PH using CMR-FT has not been investigated. Therefore, the aims of this study were to 1) evaluate the feasibility of performing routine RV strain analysis using CMR-FT, 2) characterize the pattern and severity of RV strain abnormalities in patients with PH, and 3) establish the prognostic significance of RV strain measurements in PH.


Patient population

We retrospectively evaluated patients referred for CMR evaluation of known or suspected chronic PH who also underwent right heart catheterization (RHC) within 1 month of the CMR. The presence of PH was defined as a mean pulmonary artery pressure >25 mm Hg at RHC [14]. Patients without PH and with RV dysfunction (RVEF <50 %) [15] suggestive of underlying myocardial disease, patients with cardiac shunts, or those with PH Groups 2-4 were excluded, resulting in 116 patients from Group 1 (pulmonary arterial hypertension) and Group 5 (unclear multifactorial mechanisms) of the PH Nice Classification [2]. After exclusion of 6 patients with inappropriate image quality due to arrhythmia, 110 patients were included in the present analysis. Of those, 17 patients who had no PH on RHC and a normal echocardiogram and CMR were used as controls.

Medical records were reviewed for clinical, hemodynamic, and CMR data. In addition, a composite endpoint of clinical worsening was recorded, defined as 1) all-cause mortality, 2) lung transplantation, or 3) worsening of New York Heart Association (NYHA) functional class. Only the most severe endpoint (death > transplant > worsening NYHA class) was used for analysis if more than one outcome occurred in the same patient. All deaths were confirmed by the Social Security Death Index. The Mount Sinai Institutional Review Board (New York, NY, USA) approved the study with a waiver of informed consent.

CMR acquisition

CMR studies were performed on a 1.5 Tesla (MAGNETOM Sonata or MAGNETOM Avanto, Siemens Medical Solutions, Erlanger, Germany) or 3.0 Tesla (Ingenia Philips Healthcare, Best, The Netherlands) clinical magnets using 12- or 32-channel phased-array surface coils as receivers. Images were acquired during end-expiratory breath holds with retrospective electrocardiographic or pulse gating. Standard long-axis 4-chamber cine images were obtained using SSFP imaging. In addition, contiguous cine short-axis slices covering both ventricles from base to apex were also acquired with cine SSFP (typical acquisition parameters: repetition time/echo time 3.2–3.9/0.6–2 ms, flip angle 45–90°, slice thickness 6 mm, in-plane spatial resolution 1.5–2 mm, temporal resolution 33–45 ms, 25–30 reconstructed cardiac phases). RV and LV end-diastolic and end-systolic volumes and ejection fractions were obtained according to the Simpson method using specialized software (Argus, Siemens Medical Solutions or Extended MR WorkSpace, Philips Healthcare). To calculate RV mass, RV free wall epicardial and endocardial borders were traced on each end-diastolic short-axis. Volumes and mass were indexed to body surface area. Right atrial area was planimetered in the 4-chamber view in the phase showing maximal atrial dimension [16].

Late gadolinium enhancement (LGE) short and long axis images were obtained approximately 10 min after infusion of 0.2 mmol/kg of gadolinium-diethylenetriamine pentaacetic acid (Magnevist, Berlex Laboratories, Montville, New Jersey) using a T1-weighted, 2-dimensional, fast gradient echo sequence. The presence of LGE in the interventricular septum and/or RV insertion points was assessed visually as previously described [17]. No patient had LGE in other portions of the left ventricle.

CMR-FT analysis

Strain imaging measures the percentage change in myocardial deformation, while its derivative, strain rate, measures the rate of myocardium deformation over time (s-1) [18]. Myocardial fiber lengthening (expansion), by convention, is represented as a positive value for strain, while shortening (compression) is represented by a negative value. Myocardial fibers can deform in 3 spatial directions or planes (x, y, and z axis) (Fig. 1), respectively measured as circumferential, longitudinal, and radial strain.

Fig. 1
figure 1

Schematic representation of myocardial fiber orientation and deformation in 3 orthogonal axis

CMR-FT Diogenes software (2D CPA MR, version; TomTec Imaging Systems GmbH, Unterschleissheim, Germany) was used for deformation analysis on 3 short-axis (basal, mid, and apical levels) and a 4-chamber view. To ensure a standardized analysis for each patient, the basal slice in short axis view was defined as the first slice below the atrioventricular level showing circumferential LV myocardium, the mid-ventricular slice was localized at the level of both papillary muscles, and the apical slice at an apical location with the same distance to the mid-ventricular level as the basal slice. We used end diastolic phase to start the analysis. As described previously [19], RV endocardial contours were traced manually on one frame and then the software automatically propagated the contour throughout the remainder of the cardiac cycle. The contours were checked and manually adjusted if needed. RV global longitudinal strain (GLS) and strain rate (GLSR) were calculated as the average of peak systolic values from 2 basal, 2 mid, and 2 apical segments obtained from the 4-chamber view (Fig. 2). RV global circumferential strain (GCS) and strain rate (GCSR) were calculated as the average of 4 basal, 4 mid, and 4 apical segmental peak systolic values obtained from the short-axis views. In this study we decided not to evaluate radial deformation because it is our experience, as well as others’ [20], that this parameter is less reliably quantified using CMR-FT. The entire analytic process required approximately 4 min.

Fig. 2
figure 2

RV longitudinal strain analysis obtained from a four chamber view. Examples of longitudinal strain measurements in a patient without (bottom) and with pulmonary hypertension (top). The dotted line represents the average strain of the four segments

A single investigator (M.E.M.S) reviewed all CMR studies and performed the image analyses as detailed above. To assess intra-observer variability, the same investigator re-analyzed 20 studies 2 weeks after the first analysis. A second investigator (K.M) also evaluated 20 studies for the assessment of inter-observer reproducibility.

RHC protocol

RHC was performed under fluoroscopic guidance using a Swan-Ganz catheter. Zero-pressure calibration was performed at the level of the mid-axillary line with the patient in the supine position. Baseline measurements included mean right atrial pressure, mean and systolic pulmonary artery pressures, pulmonary artery wedge pressure, cardiac index (obtained by thermodilution), pulmonary oxygen saturation, and pulmonary vascular resistance index. Systolic aortic pressure was measured by sphygmomanometry.

We classified patients according to the results of the RHC and CMR into one of three groups: A (control group; no PH and normal RVEF), B (PH and normal RVEF), or C (PH and decreased RVEF [<50 %]).

Statistical analysis

Continuous variables were expressed as mean ± SD or median [interquartile range] depending on their distribution, and categorical variables were described as total number (percentage). When comparing patients with and without PH, the χ 2 test and Fisher exact test were used for categorical variables where appropriate. Differences in normally and non-normally distributed continuous variables were established using a 2-tailed unpaired Student t test and Mann-Whitney U test, respectively, or in cases involving multiple groups, the ANOVA test and Kruskal Wallis test, respectively. Correlations between hemodynamic variables and those derived from CMR-FT strain were assessed by Pearson r or Spearman rho correlation coefficients, as appropriate. To test for intra- and inter-observer variability in strain measurements, we used Bland-Altman plots and intraclass correlation coefficient with a 2-way random model of absolute agreement.

CMR global strain parameters as well as other clinical, morphological, and hemodynamic variables significantly associated with the combined endpoint were identified. After discarding variables that showed collinearity, multivariate models were created using a selection of 6 clinically relevant variables in order to avoid overfitting. Then individual strain parameters were separately added into the models, and a Cox regression model was derived with a backward stepwise method for each strain/strain rate. Hereof, variables independently associated with the endpoint and predictive models were obtained. Results were presented as hazard ratios with 95 % confidence intervals. Receiver operating characteristic (ROC) curves were used to determine the accuracy of the global strain parameters in predicting the primary combined endpoint (death, transplant, or worsening of NYHA functional class). In addition, associations between the strain parameters and time to the primary endpoint were evaluated with adjusted survival Cox analysis using the best cut-off value derived from the ROC curves. Results were considered statistically significant when the 2-tailed p value was <0.05. Analyses were performed using SPSS 18.0 (IBM, Armonk, NY, USA).


Patient characteristics

Demographic, clinical, hemodynamic, and CMR-derived parameters for the whole sample divided according to the presence or absence of PH and RV dysfunction are shown in Table 1. Among 110 patients, PH was absent in 17 (15.5 %) and present in 93 (84.5 %). There were 70 patients (75 %) with pulmonary arterial hypertension in PH Group 1 and 23 patients (25 %) in PH Group 5. The etiologic disease responsible for the placement in Group 1 was connective tissue disease in 25 patients, idiopathic PH in 23, portopulmonary syndrome in 11, human immunodeficiency virus infection in 10, and anorexigen abuse in 1. Among the PH Group 5 patients, sarcoidosis was the cause in 23 and sickle cell disease was the cause in 2. Diseases underlying the presumed diagnosis of PH in the 17 control subjects (Group A) included scleroderma in 5, sarcoidosis in 3, hepatitis in 2, and no disease in 7. Among those with PH, 26 patients had normal RVEF and 67 had decreased RVEF (comprising Groups B and C, respectively).

Table 1 Demographic, clinical, hemodynamic and cardiac magnetic resonance data according to the presence of pulmonary hypertension and right ventricular ejection fraction

As shown in Table 1, there were no differences among Groups A, B, and C with respect to age, sex, body surface area, or cardiovascular risk factors. Patients with PH were more likely to be symptomatic (NYHA functional class ≥2) and to use diuretics, phosphodiesterase inhibitors, and prostanoids. Those with preserved RVEF used calcium channel blockers more often, while those with RV dysfunction were more frequently treated with endothelin receptor antagonists and digoxin. As expected, mean pulmonary artery pressure and pulmonary vascular resistance index increased progressively from Group A to Group C. Patients with PH and RV dysfunction had lower cardiac index and pulmonary artery oxygen saturation, larger and more hypertrophic RV, larger right atrial dimensions, smaller RV and LV ejection fractions, and more frequent LGE.

RV strain analysis

Global RV strain and strain rate values for the whole sample and Groups A-C are presented in Table 2. All strain and strain rates were reduced in patients with PH and impaired RVEF in comparison with those without PH and those with preserved RVEF. In addition, GCSR was significantly reduced in the group with PH and preserved RVEF group compared to the control group (Table 2 and Fig. 3).

Table 2 Global right ventricular strain and strain rate
Fig. 3
figure 3

Right ventricular global circumferential strain rate (GCSR) values in the three groups. Among patients with preserved right ventricular ejection fraction, there was a significantly lower value for those with pulmonary hypertension (PH) (Group B) compared to those without PH (Group A)

CMR-FT strain measurements demonstrated low inter- and intra-observer variability as indicated by the intraclass correlation coefficient (0.96/0.99; 0.94/0.98; 0.96/0.97; and 0.96/0.98 for GLS, GCS, GLSR, and GCSR, respectively; Additional file 1: Table S1 and Additional file 2: Figure S1).

Correlations between RV deformation and disease severity

Correlations between RV strain/strain rates and hemodynamic and CMR measurements are shown in Additional file 3: Table S2. For these analyses, the absolute values (without the negative sign) of deformation indices were used. Overall, there were strong and moderate positive correlations of global strains and strain rates, respectively, with RVEF, and moderate negative corrrelations with RV mass index. There were also moderate negative correlations with RV volumes and mean pulmonary artery pressures, and weaker correlations with right atrial size and remaining hemodynamic parameters (except between GCS/GCSR and cardiac index). In addition, all deformation indices were reduced in patients with LGE (Additional file 4: Table S3).

Global RV strain parameters and clinical outcome

Median duration of follow-up was 730 days (range 190 to 1585 days). There were 78 events: 38 deaths (34.5 %), 2 lung transplantations (1.8 %), and 38 cases of NYHA class worsening (34.5 %). Patients with PH and RVEF <50 % had a greater incidence of events (p < 0.001) compared to those with either PH and normal RVEF or no PH. There was no statistically significant difference regarding the composite endpoint between Groups A and B (Table 1).

The associations of all invasive and non-invasive variables with the combined endpoint are reported in Table 3. In brief, the medication use and functional class differed between those who did or did not have events. Among RHC parameters, systemic and pulmonary pressures, cardiac index, pulmonary vascular resistance index, and pulmonary oxygen saturation were associated with the primary endpoint. Similarly, among CMR measurements, reduced biventricular ejection fractions and LV end-diastolic volume, enlarged RV and right atrium, RV hypertrophy and presence of LGE each were more common in patients experiencing events. All strain and strain rates were also decreased in the event-positive group.

Table 3 Univariate analysis of the associations with the combined endpoint

From non–strain-related variables associated with outcomes (Table 3) the following were selected for inclusion in the Cox analysis based on their univariate associations and known clinical value: NYHA class ≥2, cardiac index, pulmonary artery oxygen saturation, mean pulmonary artery pressure, RV mass index, and RVEF. We did not include RV end-diastolic volume index, systolic pulmonary pressure, or pulmonary vascular resistance for multivariate analysis since they showed collinearity with previous variables, and to avoid model overfitting. Table 4 shows the final predictive models with the corresponding adjusted hazard ratios obtained when each global strain parameter was introduced into the model. When the strain parameters were analyzed as continuous variables, GLS (together with cardiac index, RV mass index and mean pulmonary artery pressure), GLSR (together with cardiac index and RV mass index), and GCSR (together with cardiac index, RV mass index and mean pulmonary artery pressure) were each independently associated with the combined endpoint, whereas GCS was not. GLS, GLSR, GCR and GCSR cut-off points, obtained from optimal values in the ROC curves (GLS -17 %, GLSR -1.1 s-1, GCS -12 %, and GCSR -0.8 s-1) were able to differentiate patients that presented the combined end-point during the follow-up, as shown in Kaplan-Meyers curves (Fig. 4). However, only GCSR cut-off point remained as independent outcome predictor after adjusting in the multivariate model (Fig. 5).

Table 4 Cox proportional adjusted hazard ratio for final multivariate models of each RV strain parameter
Fig. 4
figure 4

Kaplan Meyer survival curves for survival free of the combined endpoint according to CMR global RV strain parameters

Fig. 5
figure 5

Cox survival curve for survival free of the combined endpoint according to CMR global RV strain parameter cutoff point. GCSR cut off point > -0.8 s-1 was able to predict events in the multivariate model

An additional multivariate analysis was also performed including guideline-recommended variables for prognostication in PH [21]: NYHA class ≥2, cardiac index, pulmonary artery oxygen saturation, right atrial pressure and RVEF. GLS (together with cardiac index), GLSR (together with RVEF), and GCSR (together with cardiac index) were each independently associated with the combined endpoint, and again GCS was not (Additional file 5: Table S4). Finally, when all strain parameters were forced together into the first model, the final variables retained were GCSR, cardiac index, mean pulmonary artery pressure, and RV mass index (Additional file 6: Table S5).


The major findings of our study are as follows: 1) Systematic evaluation of RV strain in PH with CMR-FT is feasible and reproducible, 2) measures of RV strain correlate with disease severity, 3) reduced RV strain indices are associated with subsequent clinical deterioration, and 4) GCSR specifically is associated with outcomes even after adjustment for other prognostic predictors in PH.

RV function is the main predictor of survival in PH patients [13]. While echocardiography remains the most commonly used modality to assess this factor in clinical practice, CMR is today considered the gold standard for the quantification of RV volumes and ejection fraction [5]; however, other measures of myocardial deformation may be of interest and provide additive information regarding RV performance. Myocardial tissue tagging with CMR is considered the noninvasive reference method to assess LV regional deformation [22, 23]. Although feasible [24], its application to the right chamber is limited mostly because of the relatively thin chamber wall, and the method requires laborious offline post-processing analysis. In contrast to conventional tagging, strain-encoded CMR is an alternative technique that provides direct myocardial strain imaging, and its capability of quantifying RV strain in PH has been reported [25, 26]. However, the sequence may not be widely available, and additional scanning is required. Therefore, a method like FT that can be applied in routinely acquired cine images can be potentially more practical and time saving.

CMR-FT of the RV has been successfully used to demonstrate RV dysfunction in congenital heart disease [13, 27] and arrhythmogenic RV dysplasia [28] with good inter-study and inter-observer reproducibility for RV global longitudinal strain measures. More recently, Ohyama et al. [19] quantified RV longitudinal strain from cine images in 26 patients with PH using a different pixel-based multimodality tissue tracking method. They validated it against strain encoded imaging and also reported high reproducibility, although they could not show differences with normal volunteers. In our much larger series of 116 patients, strain could not be evaluated only in 6 patients (5 %). In addition, we were able to demonstrate significant reductions in GLS, GCS, GLSR, and GCSR compared with individuals lacking PH. Similar to Ohyama and colleagues, reductions in strain parameters correlated with disease severity whether measured as impairment in hemodynamics or abnormalities in right heart chambers. Importantly, decreased myocardial deformation was already noted in the presence of PH alongside preserved RVEF, although this reduction reached statistical significance in comparison with individuals without PH only for GCSR, perhaps due to sample size limitations. This finding indicates a potential role for this parameter for noninvasive detection of PH or of early RV decompensation.

To our knowledge, this is the first study to show the relationship between CMR-FT-derived strain analysis and outcomes in patients with PH. After adjustment for 6 covariates with known clinical and prognostic value (namely functional class, cardiac index, pulmonary artery oxygen saturation, mean pulmonary artery pressure, RV mass index, and RVEF), GLS, GLSR, and GCSR remained independently associated with the combined endpoint of death, lung transplantation, or decrease in NYHA class. Our findings are consistent with those reported in the echo literature using speckle-tracking based analysis of RV deformation. Haeck et al evaluated 150 patients with PH and demonstrated that RV GLS is a significant tripled determinant of all-cause mortality risk compared to patients with higher GLS (threshold of < −19 %) [7]. Sachdev et al also found that patients with PH and lower RV strain and strain rate values had poorer survival [29]. In a prospective study with a large cohort of 575 patients, RV strain predicted outcome independent of other clinical and echocardiographic variables [8]. Interestingly, Hardegree et al performed a study to assess whether serial quantitative assessment of RV strain by speckle-tracking was affected by PH-specific treatments and concluded that strain imaging independently predicts clinical deterioration and mortality in PH patients after the institution of medical therapy [30]. Although GLS and GLSR were also predictors of events in our study, probably reflecting the larger dependence of the RV on longitudinal shortening [31], survival Cox analysis demonstrated GCSR to be superior in identifying patients at heightened risk of events. Strain rate describes the degree of change in myocardial deformation with respect to time and has been found to reflect myocardial contractility better than strain parameters, which are more pre-load and after-load dependent and may change with ventricular dimensions [32, 33]. Normal RV function is highly dependent on longitudinal shortening [34]; however, it has been suggested that transverse wall motion may be a better reflection of systolic function in PH than longitudinal motion, and it is possible that as PH advances and the RV becomes more hypertrophic, circumferential deformation becomes relatively more important resembling LV contractility patterns [35]. In addition, we have previously demonstrated heterogeneity of regional RV function in PH even before RVEF decreases [36]. Thus, evaluating RV myocardial deformation may identify subclinical RV dysfunction before the development of abnormalities on conventional measures of RV performance. Further evidence that strain reflects RV performance beyond RVEF (similar to findings in the LV) is the independent value of strain when adjusting for RVEF shown in this study.


This was a single-center, retrospective study performed at a large tertiary hospital, so the inherent limitations of this study design cannot be avoided, specifically the possibilities of referral bias and residual confounding. We tried to limit sample heterogeneity by including only patients in PH Group 1 and Group 5 but etiologies of PH and underlying conditions still varied, thus we cannot exclude the possibility that other diseases influenced RV strain. However, this reflects the clinical reality of a relatively rare disorder. Although the number of events was rather large at 78, we were unable to include every single variable associated with outcomes in the multivariate models. Nonetheless, we could find significant associations between RV strain parameters and events risk after adjusting for 6 clinically and prognostically relevant covariates. The control group was not free of disease, a fact that could interfere with the strain values. Nevertheless, biventricular volumes and ejection fractions were normal, and no LGE was noted in these patients except one focal case at the RV insertion point in a patient with scleroderma. CMR and RHC were not performed simultaneously, which may reduce the associations between hemodynamic parameters and RV deformation; however, the interval between the 2 exams was not more than 1 month. The temporal resolution of the CMR cine images (33–45 ms) was lower than commonly recommended for strain rate evaluation; therefore, true strain rate values are likely underestimated. Finally, the CMR-FT strain analysis software used in this study was designed primarily for LV strain analysis and adapted for use on the RV; strain analysis software specific to RV has yet to be developed.


CMR-FT assessment of RV strain can be successfully incorporated into a comprehensive CMR exam protocol for PH and may represent a valuable noninvasive method of evaluation that could be widely adopted at PH centers. Routine RV strain measurement using CMR-FT is feasible in PH patients, correlates with disease severity, and is associated with clinically relevant outcomes.


CMR, cardiovascular magnetic resonance; FT, feature-tracking; GCS, global circumferential strain; GCSR, global circumferential strain rate; GLS, global longitudinal strain; GLSR, global longitudinal strain rate; LGE, late gadolinium enhancement; LV, left ventricular; NYHA, New York Heart Association; PH, pulmonary hypertension; RHC, right heart catheterization; RV, right ventricular; RVEF, right ventricular ejection fraction; SSFP, steady-state free precession


  1. McLaughlin VV, Presberg KW, Doyle RL, Abman SH, McCrory DC, Fortin T, Ahearn G, American College of Chest P. Prognosis of pulmonary arterial hypertension: Accp evidence-based clinical practice guidelines. Chest. 2004;126:78S–92S.

    Article  PubMed  Google Scholar 

  2. Simonneau G, Gatzoulis MA, Adatia I, Celermajer D, Denton C, Ghofrani A, Gomez Sanchez MA, Krishna Kumar R, Landzberg M, Machado RF, Olschewski H, Robbins IM, Souza R. Updated clinical classification of pulmonary hypertension. Journal of the American College of Cardiology. 2013;62:D34–41.

    Article  PubMed  Google Scholar 

  3. Farber HW, Loscalzo J. Pulmonary arterial hypertension. The New England journal of medicine. 2004;351:1655–65.

    CAS  Article  PubMed  Google Scholar 

  4. Haddad F, Hunt SA, Rosenthal DN, Murphy DJ. Right ventricular function in cardiovascular disease, part i: Anatomy, physiology, aging, and functional assessment of the right ventricle. Circulation. 2008;117:1436–48.

    Article  PubMed  Google Scholar 

  5. Vonk-Noordegraaf A, Souza R. Cardiac magnetic resonance imaging: What can it add to our knowledge of the right ventricle in pulmonary arterial hypertension? The American journal of cardiology. 2012;110:25S–31S.

    Article  PubMed  Google Scholar 

  6. Smith BC, Dobson G, Dawson D, Charalampopoulos A, Grapsa J, Nihoyannopoulos P. Three-dimensional speckle tracking of the right ventricle: Toward optimal quantification of right ventricular dysfunction in pulmonary hypertension. Journal of the American College of Cardiology. 2014;64:41–51.

    Article  PubMed  Google Scholar 

  7. Haeck ML, Scherptong RW, Marsan NA, Holman ER, Schalij MJ, Bax JJ, Vliegen HW, Delgado V. Prognostic value of right ventricular longitudinal peak systolic strain in patients with pulmonary hypertension. Circulation Cardiovascular imaging. 2012;5:628–36.

    Article  PubMed  Google Scholar 

  8. Fine NM, Chen L, Bastiansen PM, Frantz RP, Pellikka PA, Oh JK, Kane GC. Outcome prediction by quantitative right ventricular function assessment in 575 subjects evaluated for pulmonary hypertension. Circulation Cardiovascular imaging. 2013;6:711–21.

    Article  PubMed  Google Scholar 

  9. Hor KN, Baumann R Fau - Pedrizzetti G, Pedrizzetti G Fau - Tonti G, Tonti G Fau - Gottliebson WM, Gottliebson Wm Fau - Taylor M, Taylor M Fau - Benson W, Benson W Fau - Mazur W, Mazur W. Magnetic resonance derived myocardial strain assessment using feature tracking. Lid - 10.3791/2356 [doi] lid - 2356 [pii]. J Vis Exp. 2011;12.

  10. Maret E, Todt T, Brudin L, Nylander E, Swahn E, Ohlsson JL, Engvall JE. Functional measurements based on feature tracking of cine magnetic resonance images identify left ventricular segments with myocardial scar. Cardiovascular ultrasound. 2009;7:53.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Hor KN, Gottliebson WM, Carson C, Wash E, Cnota J, Fleck R, Wansapura J, Klimeczek P, Al-Khalidi HR, Chung ES, Benson DW, Mazur W. Comparison of magnetic resonance feature tracking for strain calculation with harmonic phase imaging analysis. JACC Cardiovascular imaging. 2010;3:144–51.

    Article  PubMed  Google Scholar 

  12. Wu L, Germans T, Guclu A, Heymans MW, Allaart CP, van Rossum AC. Feature tracking compared with tissue tagging measurements of segmental strain by cardiovascular magnetic resonance. Journal of cardiovascular magnetic resonance : official journal of the Society for Cardiovascular Magnetic Resonance. 2014;16:10.

    Article  Google Scholar 

  13. Kempny A, Diller GP, Orwat S, Kaleschke G, Kerckhoff G, Bunck A, Maintz D, Baumgartner H. Right ventricular-left ventricular interaction in adults with tetralogy of fallot: A combined cardiac magnetic resonance and echocardiographic speckle tracking study. International journal of cardiology. 2012;154:259–64.

    Article  PubMed  Google Scholar 

  14. Galie N, Hoeper MM, Humbert M, Torbicki A, Vachiery JL, Barbera JA, Beghetti M, Corris P, Gaine S, Gibbs JS, Gomez-Sanchez MA, Jondeau G, Klepetko W, Opitz C, Peacock A, Rubin L, Zellweger M, Simonneau G, Guidelines ESCCfP. Guidelines for the diagnosis and treatment of pulmonary hypertension: The 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 the international society of heart and lung transplantation (ishlt). European heart journal. 2009;30:2493–537.

    Article  PubMed  Google Scholar 

  15. Tandri H, Daya SK, Nasir K, Bomma C, Lima JA, Calkins H, Bluemke DA. Normal reference values for the adult right ventricle by magnetic resonance imaging. The American journal of cardiology. 2006;98:1660–4.

    Article  PubMed  Google Scholar 

  16. Maceira AM, Cosin-Sales J, Roughton M, Prasad SK, Pennell DJ. Reference right atrial dimensions and volume estimation by steady state free precession cardiovascular magnetic resonance. Journal of cardiovascular magnetic resonance : official journal of the Society for Cardiovascular Magnetic Resonance. 2013;15:29.

    Article  Google Scholar 

  17. Sanz J, Dellegrottaglie S, Kariisa M, Sulica R, Poon M, O’Donnell TP, Mehta D, Fuster V, Rajagopalan S. Prevalence and correlates of septal delayed contrast enhancement in patients with pulmonary hypertension. The American journal of cardiology. 2007;100:731–5.

    Article  PubMed  Google Scholar 

  18. Rudski LG, Lai WW, Afilalo J, Hua L, Handschumacher MD, Chandrasekaran K, Solomon SD, Louie EK, Schiller NB. Guidelines for the echocardiographic assessment of the right heart in adults: A report from the american society of echocardiography endorsed by the european association of echocardiography, a registered branch of the european society of cardiology, and the canadian society of echocardiography. Journal of the American Society of Echocardiography : official publication of the American Society of Echocardiography. 2010;23:685–713. quiz 786-688.

    Article  Google Scholar 

  19. Ohyama Y, Ambale-Venkatesh B, Chamera E, Shehata ML, Corona-Villalobos CP, Zimmerman SL, Hassoun PM, Bluemke DA, Lima JA. Comparison of strain measurement from multimodality tissue tracking with strain-encoding mri and harmonic phase mri in pulmonary hypertension. International journal of cardiology. 2015;182C:342–8.

    Article  Google Scholar 

  20. Kempny A, Fernandez-Jimenez R, Orwat S, Schuler P, Bunck AC, Maintz D, Baumgartner H, Diller GP. Quantification of biventricular myocardial function using cardiac magnetic resonance feature tracking, endocardial border delineation and echocardiographic speckle tracking in patients with repaired tetralogy of fallot and healthy controls. Journal of cardiovascular magnetic resonance : official journal of the Society for Cardiovascular Magnetic Resonance. 2012;14:32.

    Article  Google Scholar 

  21. Galie N, Humbert M, Vachiery JL, Gibbs S, Lang I, Torbicki A, Simonneau G, Peacock A, Vonk Noordegraaf A, Beghetti M, Ghofrani A, Gomez Sanchez MA, Hansmann G, Klepetko W, Lancellotti P, Matucci M, McDonagh T, Pierard LA, Trindade PT, Zompatori M, Hoeper M. 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). The European respiratory journal. 2015;46:903–75.

    CAS  Article  PubMed  Google Scholar 

  22. Lima JA, Jeremy R, Guier W, Bouton S, Zerhouni EA, McVeigh E, Buchalter MB, Weisfeldt ML, Shapiro EP, Weiss JL. Accurate systolic wall thickening by nuclear magnetic resonance imaging with tissue tagging: Correlation with sonomicrometers in normal and ischemic myocardium. Journal of the American College of Cardiology. 1993;21:1741–51.

    CAS  Article  PubMed  Google Scholar 

  23. Yeon SB, Reichek N, Tallant BA, Lima JA, Calhoun LP, Clark NR, Hoffman EA, Ho KK, Axel L. Validation of in vivo myocardial strain measurement by magnetic resonance tagging with sonomicrometry. Journal of the American College of Cardiology. 2001;38:555–61.

    CAS  Article  PubMed  Google Scholar 

  24. Fayad ZA, Ferrari VA, Kraitchman DL, Young AA, Palevsky HI, Bloomgarden DC, Axel L. Right ventricular regional function using mr tagging: Normals versus chronic pulmonary hypertension. Magn Reson Med. 1998;39:116–23.

    CAS  Article  PubMed  Google Scholar 

  25. Oyama-Manabe N, Sato T, Tsujino I, Kudo K, Manabe O, Kato F, Osman NF, Terae S. The strain-encoded (senc) mr imaging for detection of global right ventricular dysfunction in pulmonary hypertension. The international journal of cardiovascular imaging. 2013;29:371–8.

    Article  PubMed  Google Scholar 

  26. Shehata ML, Harouni AA, Skrok J, Basha TA, Boyce D, Lechtzin N, Mathai SC, Girgis R, Osman NF, Lima JA, Bluemke DA, Hassoun PM, Vogel-Claussen J. Regional and global biventricular function in pulmonary arterial hypertension: A cardiac mr imaging study. Radiology. 2013;266:114–22.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Schmidt R, Orwat S, Kempny A, Schuler P, Radke R, Kahr PC, Hellige A, Baumgartner H, Diller GP. Value of speckle-tracking echocardiography and mri-based feature tracking analysis in adult patients after fontan-type palliation. Congenital heart disease. 2014;9:397–406.

    Article  PubMed  Google Scholar 

  28. Heermann P, Hedderich DM, Paul M, Schulke C, Kroeger J, Baessler B, Wichter T, Maintz D, Waltenberger J, Heindel W, Bunck AC. Biventricular myocardial strain analysis in patients with arrhythmogenic right ventricular cardiomyopathy (arvc) using cardiovascular magnetic resonance feature tracking. Journal of cardiovascular magnetic resonance : official journal of the Society for Cardiovascular Magnetic Resonance. 2014;16:75.

    Article  Google Scholar 

  29. Sachdev A, Villarraga HR, Frantz RP, McGoon MD, Hsiao JF, Maalouf JF, Ammash NM, McCully RB, Miller FA, Pellikka PA, Oh JK, Kane GC. Right ventricular strain for prediction of survival in patients with pulmonary arterial hypertension. Chest. 2011;139:1299–309.

    Article  PubMed  Google Scholar 

  30. Hardegree EL, Sachdev A, Villarraga HR, Frantz RP, McGoon MD, Kushwaha SS, Hsiao JF, McCully RB, Oh JK, Pellikka PA, Kane GC. Role of serial quantitative assessment of right ventricular function by strain in pulmonary arterial hypertension. The American journal of cardiology. 2013;111:143–8.

    Article  PubMed  Google Scholar 

  31. Kukulski T, Hubbert L, Arnold M, Wranne B, Hatle L, Sutherland GR. Normal regional right ventricular function and its change with age: A doppler myocardial imaging study. Journal of the American Society of Echocardiography : official publication of the American Society of Echocardiography. 2000;13:194–204.

    CAS  Article  Google Scholar 

  32. Greenberg NL, Firstenberg MS, Castro PL, Main M, Travaglini A, Odabashian JA, Drinko JK, Rodriguez LL, Thomas JD, Garcia MJ. Doppler-derived myocardial systolic strain rate is a strong index of left ventricular contractility. Circulation. 2002;105:99–105.

    Article  PubMed  Google Scholar 

  33. Sutherland GR, Di Salvo G, Claus P, D’Hooge J, Bijnens B. Strain and strain rate imaging: A new clinical approach to quantifying regional myocardial function. Journal of the American Society of Echocardiography : official publication of the American Society of Echocardiography. 2004;17:788–802.

    Article  Google Scholar 

  34. Carlsson M, Ugander M, Heiberg E, Arheden H. The quantitative relationship between longitudinal and radial function in left, right, and total heart pumping in humans. American journal of physiology. Heart and circulatory physiology. 2007;293:H636–644.

    CAS  Article  PubMed  Google Scholar 

  35. Kind T, Mauritz GJ, Marcus JT, van de Veerdonk M, Westerhof N, Vonk-Noordegraaf A. Right ventricular ejection fraction is better reflected by transverse rather than longitudinal wall motion in pulmonary hypertension. Journal of cardiovascular magnetic resonance : official journal of the Society for Cardiovascular Magnetic Resonance. 2010;12:35.

    Article  Google Scholar 

  36. Fernandez-Friera L, Garcia-Alvarez A, Guzman G, Bagheriannejad-Esfahani F, Malick W, Nair A,Fuster V, Garcia MJ, Sanz J. Apical right ventricular dysfunction in patients with pulmonary hypertension demonstrated with magnetic resonance. Heart. 2011;97:1250–6.

    Article  PubMed  Google Scholar 

Download references

Authors’ contributions

MEMS was involved in acquiring and interpreting the data. She was primarily responsible for drafting the manuscript. JS conceived of the idea for the project and helped critically revise the manuscript for intellectual content. PS, AC, JN, VF contributed valuable comments and formulations. EPO performed the statistical analysis and revised the manuscript. KM was responsible for the interobserver variability. SG, VF, CBC helped in the data analysis. All authors gave final approval of this version to be published.

Competing interests

Dr Sengupta is an advisor for Saffron Technology and Heart Test Laboratories and a consultant for Edwards Lifesciences. The other authors report no conflicts.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Javier Sanz.

Additional files

Additional file 1: Table S1.

Inter- and intra-observer agreement. (DOCX 48 kb)

Additional file 2: Figure S1.

Bland-Altman plots for inter-observer agreement. (PPTX 509 kb)

Additional file 3: Table S2.

Correlation coefficients of right ventricular strain indices (absolute values) with CMR and RHC-derived parameters. (DOCX 82 kb)

Additional file 4: Table S3.

Strain parameters in patients with and without late gadolinium enhancement. (DOCX 47 kb)

Additional file 5: Table S4.

Cox proportional adjusted hazard ratio for final multivariate alternative model of each RV strain parameter. (DOCX 45 kb)

Additional file 6: Table S5.

Cox proportional adjusted hazard ratio for final multivariate model including all RV strain parameters. (DOCX 39 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

de Siqueira, M.E.M., Pozo, E., Fernandes, V.R. et al. Characterization and clinical significance of right ventricular mechanics in pulmonary hypertension evaluated with cardiovascular magnetic resonance feature tracking. J Cardiovasc Magn Reson 18, 39 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Pulmonary hypertension
  • Strain
  • Right ventricle