Impact of abnormal longitudinal rotation on the assessment of right ventricular systolic function in patients with severe pulmonary hypertension
Introduction
Right ventricular (RV) systolic function is an important determinant of prognosis in patients with pulmonary hypertension (PH) (1,2). Indeed, current guidelines recommend that echocardiographic assessment of RV systolic function be carried out as an important component of the follow-up care of patients with PH (1). Given the well-documented complexity involved in echocardiographic assessment of RV systolic function (2), it is not surprising that a number of echocardiographic methods for evaluating RV systolic function have been used in clinical practice. Frequently used, recommended parameters include fractional area change (FAC), tricuspid annular plane systolic excursion (TAPSE), and pulsed wave tissue Doppler tricuspid annular systolic longitudinal excursion velocity (RVS’) (3). Both TAPSE and RVS’ assess RV longitudinal function (with the inference that the greater the descent or speed of descent of the base in systole, the better the RV systolic function), and when these parameters decreased, it has been shown to reflect RV dysfunction, and to be of prognostic value in PH (4-6). RV strain (the percentage change in myocardial deformation) has emerged as a novel marker of RV systolic function, serving as a marker of both regional and global RV contractility. It has the theoretical advantage of being able to distinguish between myocardial deformation and myocardial displacement (3).
Recently, in our echocardiographic laboratory, we observed the occurrence of an abnormal RV rocking motion in patients with PH that bore some similarities to our prior findings of abnormal left ventricular (LV) longitudinal rotation in the setting of LV dilatation and dysfunction (7-9). The aims of this study were to: (I) assess longitudinal rotation of the RV (RVLR) in patients with PH; (II) determine whether such RVLR might have an impact on the conventional echocardiographic parameters used in the assessment of overall RV systolic function (i.e., TAPSE and RVS’); and (III) determine whether RVLR has an impact on prognosis.
Methods
Study sample
We studied patients in whom an echocardiogram was performed in an outpatient setting, as a part of the initial visit or subsequent follow-up evaluation, at the request of one of the PH specialists of the Cleveland Clinic Respiratory Institute. All consecutive patients seen during the period between April 1st 2011 and April 1st 2012, who fulfilled the inclusion criteria were included. The inclusion criteria were: an estimated systolic pulmonary arterial pressure of >60 mmHg on a formal echocardiogram report, and performance of the echocardiogram study using either Vivid 7 or Vivid 9 ultrasound systems (GE Medical, Milwaukee, USA). This was necessary, because only images acquired using these machines were amenable to the relevant analysis using the proprietary automated software (EchoPAC, Version 110.0.0. Advanced Analysis Technologies; GE Medical, Milwaukee, USA). Exclusion criteria were: PH due to left heart disease [i.e., Group 2 according to World Health Organization classification of PH (10)]; incomplete study (limited acquisition of RV and LV images); and suboptimal echocardiographic image quality precluding accurate measurement of RV and LV parameters (defined as inadequate visualization of more than two ventricular segments, or significant artifact).
The study protocol was approved by the Institutional Review Board (IRB) of the Cleveland Clinic Foundation (IRB number: 12-347). The medical records of all included patients were further reviewed to gather clinical data.
Data analysis
Transthoracic echocardiograms were reviewed and analyzed by two independent, experienced echocardiography readers. Standard echocardiographic measurements of the LV and RV were made in accordance with current guidelines (3,11). Quantitative estimation of RV systolic function was carried out using FAC. Maximal systolic excursion of the tricuspid annulus in the longitudinal plane was measured using 2-dimensional (2D) images. RVS’ was measured by tissue Doppler imaging from an apical 4-chamber window, with the pulsed Doppler sample volume placed at the tricuspid annulus. Measurement of right ventricular global longitudinal strain (RVGLS) was performed by speckle-tracking using EchoPAC (Version 110.0.0. Advanced Analysis Technologies; GE Medical, Chicago, Illinois, USA) as previously described (12). Estimated peak RV systolic pressures (RVSP) were calculated by summing the tricuspid gradient (estimated by the modified Bernoulli equation from the tricuspid regurgitant velocity) and the right atrial pressure (estimated by inferior vena cava size and changes with respiration). RVLR was calculated by off-line analysis using EchoPAC (Version 110.0.0. Advanced Analysis Technologies; GE Medical, Chicago, Illinois, USA). In brief, following selection of the optimal quality 2D digital images in a dedicated RV long-axis view, the endocardium was traced in an optimal frame, from which a speckle tracking region of interest was automatically selected to approximate the myocardium between the endocardium and epicardium. The region of interest width was adjusted as needed to fit the wall thickness. Measurements were deemed unacceptable, if two or more segments did not track adequately after two manual adjustments of the endocardial borders. The software then automatically selected suitable, stable objects for tracking and searched for them in the next frame using the sum of absolute differences algorithm before computing an RVLR profile defining the ventricular centroid from the mid-myocardial line on each frame. Peak systolic RVLR was defined as the maximal rotation value of the RV during systole. In accordance with engineering notation, the negative sign indicated clockwise rotation, and a positive sign indicated counterclockwise rotation.
In vitro studies
Because speckle tracking rotation quantitation software was designed to measure angular motion of rotationally symmetrical object occurring relative to the point of reference that coincides with its center of gravity, we performed an in vitro experiment.
Commercially obtained porcine hearts weighing between 200 and 250 grams were transected in order to obtain a slice (~15 mm in thickness) corresponding to the apical four-chamber view. The hearts were then mounted on a circular plate that rotated around the axle passing either through the point corresponding to the base of the left ventricle, or to the point corresponding to the base of the right ventricle, and placed into a water tank. The plate was connected via a system of opposing springs to the rotating spindle that caused the motion of the plate in a see-saw manner, thus imitating the “rocking” motion of the heart. The diameter of the circle of rotating spindle was set to zero, maximum, and midway between zero and maximum. A total of 6 different setups were run, each of them imaged with three different frame rates (at 30 frames per second, 46 frames per second, and 57 frames per second). At each setup, two dimensional, M mode, and pulsed wave tissue Doppler data were collected using a Vivid 9 echocardiographic machine (GE Medical, Milwaukee, USA) with a standard transthoracic probe held by a holder in a water tank.
For the gold standard of angular displacement, we analyzed the fiducial points on images obtained by the video camera (30 frames per second) placed over the water tank. Data were digitized using the Engauge Digitizer v.15 (Mark Mitchell, 2012; https://markummitchell.github.io/engauge-digitizer) with the change of plate angle used for reference. Correlations were assessed by Pearson correlation coefficient.
The in vitro studies confirmed the accuracy of RV rotation measurement with good agreement between the gold standard (in vitro experiment) and echocardiography-derived rotation, with an overall correlation coefficient r of 0.98 (P<0.01).
Statistical analysis
Data were presented as either mean ± standard deviations (SDs) or median [interquartile range (IQR)] for normally and non-normally distributed continuous variables, respectively, whereas frequencies and percentages (in parentheses) summarized categorical variables. Multiple linear regression analyses were performed to assess whether the dependent variables, TAPSE and RVS’ could be predicted from a linear combination of the independent variables, RVLR and RVGLS. To assess inter-observer variability, RVLR was measured by two independent blinded investigators in a subgroup of 15 consecutive study patients using the Bland-Altman analysis. Intra-observer variability was measured by repeated assessment of RVLR in the same group of patients by one investigator at a separate time using Bland-Altman analysis. A Cox regression proportional hazards model was used to assess survival related to baseline covariates of TAPSE, RVGLS and RVLR. To control for the known impact of clinical parameters on prognosis, we first calculated the REVEAL score (13) of individual patients. Then, we performed non-parsimonious modelling of proportional hazards by first applying the REVEAL score and PH type (with patients dichotomized according to the presence or absence of type 1 PH) as predictors. In the second step, we performed forward stepwise regression using TAPSE, RVGLS, RVS’ and RVLR as additional prognostic variables. Because cross-correlation between TAPSE, RVGLS, RVS’ and RVLR was expected, we attempted separate stepwise regression for each of these four parameters into the model.
Statistical analysis was performed using the statistical software packages, Sigma Plot (Version 11.0, Systat Software, Inc., San Jose, USA) and SPSS Statistics for Windows (Version 21.0., IBM Corp, Armonk, NY, USA). A P value of <0.05 was considered significant.
Results
We identified 117 patients who fulfilled the clinical inclusion criteria. Out of this group, 11 patients were excluded due to inadequate image quality and a further 6 patients were excluded due to incomplete studies. Therefore, a total of 100 consecutively identified patients formed the final study population. The mean patient age was: 58±15 years. Baseline demographics and medication use are summarized in Table 1. According to the classification of PH based on current guidelines, the etiologies of PH in the study population were as follows: Group 1—pulmonary arterial hypertension (80%), Group 2—PH due to left heart disease (0%), Group 3—PH due to lung disease and/or hypoxia (5%), Group 4—chronic thromboembolic PH (10%), and Group 5—PH with unclear and/or multifactorial mechanisms (5%) (1). For 88 of the 100 study patients (88%), invasive hemodynamic data were available, with a mean right atrial pressure of 12±6 mmHg, a mean pulmonary artery pressure of 55±18 mmHg, a mean wedge pressure of 12±7 mmHg, and a mean cardiac index (Fick method) of 2.4±0.9 L/min/m2.
Full table
The tracking of rotational motion of the right ventricle in a set of systolic images obtained in apical 4-chamber views obtained in a patient with a prominent clockwise RVLR is demonstrated in Figure 1, while the video clip obtained in the same patient is available as Figures S1 and S2. All echocardiographic variables are summarized in Table 2. In keeping with inclusion criteria of severe PH, right heart catheterization data revealed a mean pulmonary artery pressure of 55±18 mmHg, while median (IQR) estimated RVSP measured 87 [74–99] mmHg. Most patients had RV enlargement with median RV end-diastolic area of 32 [25–43] cm2 [compared to lower reference value (LRV) ±95% confidence intervals (CIs) of 25% (24–27%) as per guidelines] (3). This study population also had documented RV systolic dysfunction with decreased mean FAC values of 27.3%±8.9% [LRV ±95% CIs of 35% (32–38%) as per guidelines] (3). Consistent with this finding, RV deformation was also abnormally decreased, with mean RVGLS values of −13.9%±4.6% [LRV ±95% CIs of −20% (−15% to −24%)], which was approximately half of published normal values (3). In contrast, TAPSE and RVS’ were relatively preserved, with mean TAPSE values of 17.6±4.4 mm {compared to LRV ±95% CIs of 16 [15–18] mm} and mean RVS’ values of 11.2±3.1 cm/s [compared to LRV ±95% CIs of 6 mm (5–7 mm)] (3). While 37 patients had identifiable pericardial effusions, the majority of these pericardial effusions were either trivial (n=16) or small (n=17). Peak RVLR measured −4.2°±3.7°. Absolute inter-observer variability was −0.6°±1.2° (95% CIs of mean of differences: −1.3° to 0.1°). The error was random with no systematic trend or proportional bias observed (P=0.36), while the inter-observer coefficient of repeatability was 2.4°. Absolute intra-observer variability was −0.1°±1.0° (95% CIs of mean of differences: −0.4° to 0.6°). Again, the error was random with no systematic trend or proportional bias observed (P=0.19), while the intra-observer coefficient of repeatability was 2.0°.
Full table
Multiple linear regression demonstrated that TAPSE can be predicted from a linear combination of peak RVLR and peak RVGLS as follows: TAPSE =10.3 − (0.3 × RVLR) − (0.4 × RVGLS); R=0.56 (P<0.001) (Figure 2A,B). Similarly, RVS’ can be predicted from a linear combination of peak RVLR and peak RVGLS as follows: RVS’ =7.1 − (0.3 × RVLR) − (0.2 × RVGLS); R=0.52 (P<0.001) (Figure 3A,B). While no association was found between RVLR and RV size, estimated RVSP or the presence of a pericardial effusion, a mild correlation was noted between longitudinal RV rotation and QRS duration (R=0.25, P=0.01).
Over a median follow-up period of 996 [533–1,098] days, 41 patients died. In univariate analysis, RVGLS, TAPSE and RVLR correlated with patient survival, while RVS’ showed trends towards significance. In a multivariate analysis, only RVGLS had significant incremental prognostic value (P=0.048) over and above the REVEAL score, and the etiology of PH (Table 3).
Full table
Discussion
This study investigated a previously undescribed movement of the RV in the setting of severe PH, namely clockwise longitudinal rotation. We demonstrated that such rotatory (or rocking) motion influenced both TAPSE and RVS’, both of which are important conventional measures of RV systolic function, and markers of prognosis (4-6). In particular, greater degrees of clockwise rotation were associated with greater amounts of both tricuspid annular displacements (as measured by TAPSE) and systolic velocity (RVS’). However, the values of these parameters may be falsely reassuring despite the presence of significant RV dysfunction.
Prior studies that investigated the different parameters of RV systolic function have shown that FAC, TAPSE and RVS’ were generally concordant (14). In contrast, in the present cohort of consecutive patients with severe PH, we found FAC and TAPSE to be somewhat discordant, i.e., FAC was markedly reduced, while TAPSE and RVS’ levels were found to be above the published LRV (even if still below the means) (3). Consideration of these findings with reference to the three discrete elements known to incorporate RV systolic function may however offer further insight. Contraction of spiral longitudinal muscles (which produces systolic basal-apical shortening), inter-ventricular septal motion (affected by biventricular interdependence and conduction (and motion of the RV free wall are important determinants of FAC, whereas longitudinal function alone is the predominant determinant of TAPSE and RVS’.
Peak RVLR and RVGLS were found to be predictive of both TAPSE and RVS’, although regression coefficients (R=0.56 and 0.52, respectively) were relatively modest. However, RV remodeling in PH is clearly a complex phenomenon involving a structure with an already complex geometry. Indeed, a multitude of modifying factors may exist, not the least of which includes differences in afterload, ventricular hypertrophy, fibrosis, ischemia, wall stress and the degree of flattening of the inter-ventricular septum. One such factor that may amplify the effects of RV rotation on tricuspid annular motion is the degree of RV dilation. This is known to occur relatively early, in response to increased afterload on the right side of the heart compared to the left, due to the relatively increased compliance and lower mass/volume ratio of the RV (15-17). However, we did not find a relationship between peak RVLR and RV size, which may relate to the fact that most patients already had significant RV dilation. Although pericardial effusions are commonly found in association with PH, our data suggested that the absence or presence of a pericardial effusion appears to have little or no role in augmenting rotatory motion, and this may relate to our previous findings that only a relatively small amount of energy is actually required to induce longitudinal rotation (7). Unique to the RV is the premier contribution played by longitudinal fibers to overall RV systolic function (15). It is not surprising therefore that, among non-volumetric echocardiographic indices, speckle tracking longitudinal RV strain was shown to exhibit the closest association with RV ejection fraction as determined by cardiac magnetic resonance imaging (the reference standard of RV systolic function) (18). Thus, a noteworthy finding in the current study is that the peak RVLR was of similar importance to speckle tracking longitudinal RV strain, in terms of the impact on TAPSE and RVS’.
The etiology of abnormal rotatory motion of the RV in the setting of PH without intrinsic LV dysfunction can only be speculated. One potential mechanism for abnormal rotation may relate to an imbalance of ventricular strains. We have previously shown that in patients with LV dysfunction, clockwise LV longitudinal rotation occurred in patients with high strains of lateral LV wall (7,8). Furthermore, we have previously shown that in PH, longitudinal LV lateral wall strain was typically preserved, while RV free wall and inter-ventricular septal longitudinal strains were typically reduced (proportional to the degree of pressure overload as assessed by either estimated RVSP or by the LV eccentricity index—a measure of septal displacement or flattening) (19). Thus, the presence of a ventricular strain gradient between the LV lateral wall (with preserved strain) and RV free wall may contribute to systolic clockwise rotation of the dilated, pressure-overloaded RV.
Although overall RV systolic function should never be assessed by annular motion alone, from a clinical practice point of view, TAPSE is a parameter that is more commonly obtained than other markers of RV systolic function, a point that is highlighted in the current echocardiographic guidelines (3). A likely explanation is that the RV free wall cannot be completely visualized, yet annular motion may still be clearly evident. In this scenario, however, increased clockwise rotation can still be qualitatively identified by assessing the relative angulation between the annular plane and the inter-ventricular septum. When RV contraction is preserved, this angle normally decreases. Our findings suggest that when this angle remains fixed, it could be speculated that “relative preservation” of annular motion may result from RVLR, rather than preserved RV contraction, raising the suspicion for possible significant underlying RV dysfunction (Figures 4,5). Awareness of the potential confounding contribution of RVLR to the “relative preservation” of annular motion may therefore limit the possibility of missing RV dysfunction, particularly where such a finding may have prognostic importance.
Limitations
Although this study population represented a consecutive series of patients, subjects were retrospectively identified and recruited into this study. Many of these patients with PH already had advanced RV systolic dysfunction. However, we would argue that the demonstration of relative preservation of TAPSE and RVS’ in this context further emphasized our conclusions regarding the potential confounding contribution of RVLR. Additionally, the software used was not designed to specifically assess change in RV orientation in a horizontal plane. Possibly, better software, perhaps by separating lateral and septal walls, would have better accuracy to quantitate rotational motion of RV silhouette.
Further, the adjustment for clinical factors in the survival analysis was performed using the REVEAL score (13), which was developed specifically for the initial assessment of survival in type 1 PH. On the other hand, almost all of the parameters incorporated into the REVEAL score (e.g., functional class, age, kidney function, and blood pressure) were prognostic in other types of PH. A larger dataset is needed to adequately assess the hazard associated with excessive RVLR. Nevertheless, our survival analysis was hypothesis-generating, and may serve as the rationale for a future evaluation of RVLR as a potential survival predictor for patients with PH.
Conclusions
Clockwise angular displacement of the right ventricle (RVLR) was observed in a population of patients with PH and RV systolic dysfunction. RVLR influenced the measurements of the TAPSE and RVS’ in these patients. In contrast, RVLR did not affect RV global longitudinal strain, which may indicate its potential usefulness as a prognostic marker in PH.
Acknowledgements
None.
Footnote
Conflicts of Interest: The authors have no conflicts of interest to declare.
Ethical Statement: The study protocol was approved by the Institutional Review Board (IRB) of the Cleveland Clinic Foundation (IRB number: 12-347).
References
- Galie N, Hoeper MM, Humbert M, et al. Guidelines for the diagnosis and treatment of pulmonary hypertension. Eur Respir J 2009;34:1219-63. [Crossref] [PubMed]
- Jiang L, Levine RA, Weyman AE. Echocardiographic Assessment of Right Ventricular Volume and Function. Echocardiography 1997;14:189-206. [Crossref] [PubMed]
- Rudski LG, Lai WW, Afilalo J, et al. 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. J Am Soc Echocardiogr 2010;23:685-713. [Crossref] [PubMed]
- Forfia PR, Fisher MR, Mathai SC, et al. Tricuspid annular displacement predicts survival in pulmonary hypertension. Am J Respir Crit Care Med 2006;174:1034-41. [Crossref] [PubMed]
- Lindqvist P, Waldenstrom A, Henein M, et al. Regional and global right ventricular function in healthy individuals aged 20-90 years: a pulsed Doppler tissue imaging study: Umea General Population Heart Study. Echocardiography 2005;22:305-14. [Crossref] [PubMed]
- Mathai SC, Sibley CT, Forfia PR, et al. Tricuspid annular plane systolic excursion is a robust outcome measure in systemic sclerosis-associated pulmonary arterial hypertension. J Rheumatol 2011;38:2410-8. [Crossref] [PubMed]
- Popovic ZB, Grimm RA, Ahmad A, et al. Longitudinal rotation: an unrecognised motion pattern in patients with dilated cardiomyopathy. Heart 2008;94. [Crossref] [PubMed]
- Phillips KP, Popovic ZB, Lim P, et al. Opposing wall mechanics are significantly influenced by longitudinal cardiac rotation in the assessment of ventricular dyssynchrony. JACC Cardiovasc Imaging 2009;2:379-86. [Crossref] [PubMed]
- Stankovic I, Aarones M, Smith HJ, et al. Dynamic relationship of left-ventricular dyssynchrony and contractile reserve in patients undergoing cardiac resynchronization therapy. Eur Heart J 2014;35:48-55. [Crossref] [PubMed]
- Simonneau G, Gatzoulis MA, Adatia I, et al. Updated clinical classification of pulmonary hypertension. J Am Coll Cardiol 2013;62:D34-41. [Crossref] [PubMed]
- Lang RM, Bierig M, Devereux RB, et al. Recommendations for chamber quantification: a report from the American Society of Echocardiography's Guidelines and Standards Committee and the Chamber Quantification Writing Group, developed in conjunction with the European Association of Echocardiography, a branch of the European Society of Cardiology. J Am Soc Echocardiogr 2005;18:1440-63. [Crossref] [PubMed]
- Chow PC, Liang XC, Cheung EW, et al. New two-dimensional global longitudinal strain and strain rate imaging for assessment of systemic right ventricular function. Heart 2008;94:855-9. [Crossref] [PubMed]
- Benza RL, Gomberg-Maitland M, Miller DP, et al. The REVEAL Registry risk score calculator in patients newly diagnosed with pulmonary arterial hypertension. Chest 2012;141:354-62. [Crossref] [PubMed]
- Tamborini G, Pepi M, Galli CA, et al. Feasibility and accuracy of a routine echocardiographic assessment of right ventricular function. Int J Cardiol 2007;115:86-9. [Crossref] [PubMed]
- Haddad F, Hunt SA, Rosenthal DN, et al. Right ventricular function in cardiovascular disease, part I: Anatomy, physiology, aging, and functional assessment of the right ventricle. Circulation 2008;117:1436-48. [Crossref] [PubMed]
- Chin KM, Rubin LJ. Pulmonary arterial hypertension. J Am Coll Cardiol 2008;51:1527-38. [Crossref] [PubMed]
- Haddad F, Doyle R, Murphy DJ, et al. Right ventricular function in cardiovascular disease, part II: pathophysiology, clinical importance, and management of right ventricular failure. Circulation 2008;117:1717-31. [Crossref] [PubMed]
- Leong DP, Grover S, Molaee P, et al. Nonvolumetric echocardiographic indices of right ventricular systolic function: validation with cardiovascular magnetic resonance and relationship with functional capacity. Echocardiography 2012;29:455-63. [Crossref] [PubMed]
- Puwanant S, Park M, Popovic ZB, et al. Ventricular geometry, strain, and rotational mechanics in pulmonary hypertension. Circulation 2010;121:259-66. [Crossref] [PubMed]
- Collier P, Xu B, Kusunose K, et al. Grayscale video of patient shown in —displaying clockwise right ventricular rotation. Asvide 2018;5:706. Available online: http://www.asvide.com/article/view/26622
- Collier P, Xu B, Kusunose K, et al. Video of the same patient with color-coded tracking of right ventricular rotation. Asvide 2018;5:707. Available online: http://www.asvide.com/article/view/26623