- Open Access
Left ventricular torsion shear angle volume analysis in patients with hypertension: a global approach for LV diastolic function
Journal of Cardiovascular Magnetic Resonance volume 16, Article number: 70 (2014)
Torsion shear angle φ is an important measure of left ventricular (LV) systolic and diastolic functions. Here we provide a novel index utilizing LV normalized torsion shear angle volume loop to assess LV diastolic functional properties. We defined the area within loop as torsion hysteresis area, and hypothesized that it may be an important global parameter of diastolic function. We evaluated the changes to increased during early diastole as a potential measure of LV suction.
Sixty resistant hypertension patients (HTN), forty control volunteers were studied using cardiovascular magnetic resonance with tissue tagging. Volumetric and torsional parameters were evaluated.
HTN demonstrated concentric remodeling with preserved ejection fraction. HTN had significantly decreased normalized early filling rate, early diastolic mitral annulus velocity and E/A (1.33 ± 1.13 vs. 2.19 ± 1.07, P < 0.0001) vs. control. Torsion hysteresis area was greater (0.11 ± 0.07 vs. 0.079 ± 0.045, P < 0.001) and peak at early diastole was higher (10.46 ± 8.51 vs. 6.29 ± 3.85, P = 0.002) than control. Torsion hysteresis area was significantly correlated with E/A (r = −0.23, P = 0.025). Thirteen HTN patients had both E/A ratio < 1.12 (Control mean E/A-1SD) and torsion hysteresis area > 0.12 (Control mean torsion hysteresis area + 1SD).
Torsion hysteresis area and peak early diastolic were significantly increased in hypertensive concentric remodeling. The loop takes into account the active and passive recoil processes of LV diastolic and systolic phases, therefore provides a new global description of LV diastolic function.
Left ventricular (LV) diastolic dysfunction is characterized by abnormal myocardial mechanical properties that include impaired diastolic distensibility, impaired LV filling and slow or delayed myocardial relaxation . Up to 50% of patients with heart failure have predominant diastolic dysfunction in the presence of preserved LV ejection fraction (EF) .
Invasive hemodynamic measurements, such as the time constant of relaxation(1), τ, and LV end diastolic pressure volume relationship (EDPVR)  are considered useful to assess diastolic function. However, they are not frequently performed in routine clinical practice. In contrast, echocardiographic tissue Doppler measurement of early diastolic mitral annular velocity with other appropriate parameters is frequently used as a non-invasive diagnostic tool of diastolic dysfunction . However these methodologies may have significant limitations .
LV twist, measured as the myocardial rotation gradient from the base to apex along a longitudinal axis , is an important mechanical property of the myocardium that results from the helical fiber arrangement of LV. Torsion shear angle is twist normalized to long-axis length and LV radius . Contraction of myocardial bundles and their interaction with extracellular matrix during systole result in storage of torsional potential energy. Torsional recoil during isovolumic relaxation and early diastole releases the potential energy stored in the deformed matrix during systole ,. Previous work with cardiovascular magnetic resonance (CMR) in an animal model  and recently in humans , showed that LV early untwist rate correlated closely with τ.
In this paper, a global approach utilizing normalized LV torsion shear angle volume (normalized φ is indicated as and normalized V is indicated as ) loop is proposed to assess LV diastolic function. Torsion hysteresis area (THA) is quantified as the area within the loop. This ‘hysteresis’ concept was inspired by the history-dependent stress strain relationships in viscoelastic materials undergoing cyclic loading, representing energy loss or work. As it is difficult to measure instantaneous LV wall stress non-invasively, temporal changes in LV volume throughout the cardiac cycle provide an indirect measure of changes in ventricular stress. Moreover, LV wall tensile stress is directly related to wall geometry and wall thickness. Torsion shear angle is also affected by these factors, representing a global measure of ventricular strain. We therefore evaluated the relationship of changes in torsion shear angle to LV volumetric changes over the cardiac cycle and develop THA as a non-invasive measurement which quantifies the relationship between myocardial strain and volume change. Normalizing torsion shear angle and volume allow for comparison among groups with various geometry remodeling. Furthermore, since untwist is an important mechanism of LV suction, we explored the relationship of peak changes to increase in during early diastolic phase as a potential measure of LV suction.
We hypothesized that loop represented a new global approach to assess diastolic function. THA and peak at the early diastole derived from the loop are higher in patients with resistant hypertension (HTN).
The study population consisted of 40 normal control volunteers and 60 HTN, defined as uncontrolled hypertension (ambulatory BP > 140/90 mmHg at two or more clinic visits in spite of the use of three or more antihypertensive medications at optimal doses). Detailed methodology for enrollment of HTN patients in this prospective study design was previously described . The normal control volunteers had no history of cardiovascular disease and were not taking any cardiovascular medicines. Cine and tagged CMR was performed on all HTN patients and control volunteers for comparison purposes. Early morning ambulatory brain natriuretic peptide (BNP) was measured in HTN. The study protocol was approved by the Institutional Review Boards of University of Alabama at Birmingham and Auburn University and informed consent was obtained from all participants.
Cardiovascular Magnetic Resonance (CMR)
Cine CMR was performed on a 1.5-T MRI scanner (Signa, GE, Milwaukee, Wisconsin) optimized for cardiac application. Electrocardiographically gated breath-hold steady-state free precision technique was used to obtain standard (2-, 3-, and 4-chamber long axis and short-axis) views using the following general parameters: slice thickness of the imaging planes 8 mm, field of view 40 × 40 cm, scan matrix 256 × 128, flip angle 45°, repetition/echo times 3.8/1.6 ms, views per segment 8–10, number of reconstructed cardiac phases 20.
Tagged CMR was done on exact slice prescriptions as above by applying grid tagging to the short axis views and stripe tagging to long axis views using spatial modulation of magnetization encoding gradients method (FGR-SPAMM) as previously described  with following general parameters: prospective ECG triggering, trigger time 10 ms from R wave, slice thickness 8 mm, zero interslice gap, field of view 40 × 40 cm, scan matrix 256 × 128, flip angle 10°, repetition/echo times 8.0/4.2 ms, views per segment 8–10, tag spacing 7 mm and number of reconstructed cardiac phases 20. Because the tag lines faded with time due to T1 relaxation, tagged image derived parameters were only valid throughout systole and the first 67% of diastole.
LV geometric parameters were measured from endocardial and epicardial contours manually traced on cine images acquired near end-diastole (ED) and end-systole (ES). These contours were propagated throughout the cardiac cycle using in-house software . LV volume, volume time curve and its derived peak ejection rates, peak early and late filling rates were calculated as described . LVED mass was measured excluding the papillary muscle. Peak early and late diastolic mitral annular velocities were calculated using non-rigid registration to track a manually-selected point on the mitral annulus through the cardiac cycle . The 3D radius of curvature to wall thickness ratio (R/T) was computed by the reciprocal of the product of the endocardial circumferential curvature and 3D wall thickness as previously described .
Two-dimensional (2D) strain at each timeframe and rates were measured using harmonic phase (HARP) analysis . 2D basal and apical rotations at each timeframe were measured by tracking a circular mesh of points in the basal and apical slice of that timeframe. The mesh was identified in the first time based on user-defined contours and tracked through the remaining imaged phases using improved HARP tracking . 2D twist was computed as the apical rotation minus the basal rotation at the same timeframe. Twist time curve was constructed and differentiated to obtain the twist rate time curve. Peak early diastolic untwist rate was defined as the maximum twist rate at the early diastolic phase. Torsion shear angle φ at timeframe t was computed as 
where ρ(t) is the epicardial radius at time t and L is the distance between the basal and distal slices at ED timeframe. 2D φt curve was therefore constructed for each subject.
Torsion shear angle volume loop
For each subject, the φt curve was normalized by its maximum to generate curve and LV Vt curve was normalized by its maximum V to generate curve. Then by forcing the maximum normalized torsion shear angel and minimum normalized volume of each subject to be at the same time point, the systolic and diastolic phases of the t curve and t curve were interpolated with 10 and 18 time points respectively to create the loop as shown in Figure 1. This process was called ‘extrema matching’ and it was implemented to account for the differences in temporal resolution and changes in heart rate in cine vs tag MR images. As MR tag lines fade during mid to late diastole and torsion shear angle cannot be reliably measured, for consistency, we chose the first 12 interpolated data points (67% of diastole) in diastole for torsion hysteresis area computation. Area under the systolic and diastolic arms of the curve were computed numerically using the trapezoidal rule (Figure 1). Systolic and diastolic areas were computed over the volume interval between ES and 67% diastole. THA was computed by subtracting the diastolic area from the systolic area. Peak at the early diastole phase was calculated as the negative peak slope of the diastolic arm of the curve. The peak slope of the diastolic arm of the curve was defined as the slope of a fitted linear regression model to the first four points of the diastolic arm. Both THA and Peak were computed for each subject and averaged over each group.
Student’s two-sample t test (for continuous variables) and Fisher’s exact test (for categorical variables) were conducted to compare control and HTN in terms of demographic, geometric, and functional variables. Pearson’s correlation analysis was THA and other CMR-derived indices of LV diastolic function (normalized peak early filling rate, E/A ratio, normalized peak early diastolic mitral annulus velocity) as well as LV mass/volume ratio and mean atrial pressure (MAP). We also created multivariate model for torsional and strain rate parameters to adjust for age using linear regression models.
Data are in mean ± standard deviation (SD). A P < 0.05 was considered statistically significant. All statistical analysis was performed using SAS version 9.3.
Control’s mean age was significantly lower than that of HTN (Table 1). Thus, the comparisons between control and HTN in torsional and strain parameters were adjusted for age. HTN had higher systolic and diastolic blood pressures vs. controls at the time of CMR. BNP level in HTN was measured as 35.25 ± 34.32 pg/mL, ranging from 2 to 169.8 pg/mL.
CMR-derived LV geometric and systolic functional parameters
Controls and HTN did not differ in volumes and stroke volume indices (normalized to BSA, as shown in Table 2). LV mass and LV mass/volume ratio were significantly increased while LVED R/T ratio was significantly decreased in HTN vs. controls. Systolic function measured as LVEF and normalized peak ejection rate were higher in HTN.
CMR-derived LV diastolic functional parameters
HTN demonstrated decreased CMR-derived volumetric peak early filling rate and greater peak late filling rate expressed as EDV/s, and lower E/A ratio (1.33 ± 1.13 vs. 2.19 ± 1.07, P < 0.0001) vs. controls (Table 3). Moreover, HTN depicted lower peak early diastolic mitral annular velocity and higher peak late diastolic mitral annular velocity vs. controls. Peak early diastolic circumferential and longitudinal strain rates were significantly lower in HTN vs. controls. Adjusted for age, these differences were still significant. Peak untwist rate in the early diastole was significantly increased in HTN with and without adjusted for age. Time to peak early untwist rate did not differ in HTN vs. controls.
LV torsion shear angle volume loop
The volume vs. time curve (Figure 2A) indicates an increase in peak ejection rate (Table 2) and a decrease in early filling rate (Table 3); while the torsion shear angle vs. time curve (Figure 2B) indicates an increase of peak torsion shear angle (9.18 ± 2.33 vs. 6.72 ± 1.66°, P < 0.0001) in the HTN group vs. controls. HTN curve (Figure 2C) slightly shifted to the left of the control loop with a steeper diastolic arm. THA and peak at the early diastole phase were significantly increased in HTN vs. controls (Table 4). After adjusted for age, these significant differences still persisted.
THA demonstrated significant correlation with E/A ratio (r = −0.23, P = 0.025). Twenty two percent HTN patients (13 out of 60) had both E/A ratio < 1.12 (Control mean E/A-1SD) and THA > 0.12 (Control mean THA + 1SD) as shown in Figure 3. Furthermore, THA was not significantly correlated with LV mass/volume ratio nor MAP.
This study provides a novel approach utilizing normalized torsion shear angle volume loop, which takes into account the global cardiac cycle, to assess diastolic function. We demonstrate, 1) the phenomenon of hysteresis, based on loop where the systolic component is not identical to the diastolic component of the loop; 2) that the loop derived THA is increased in HTN vs. controls; 3) HTN has increased peak early diastolic vs. control; 4) THA is significantly correlated with E/A ratio, a conventional parameter for diastolic function. These findings suggest stiffer ventricles with impaired suction and LV early filling in HTN. Therefore, loop may be a useful approach of non-invasive global assessment of diastolic function that takes into account not only active and passive recoil processes of the LV diastolic phase but also the systolic phase.
HTN is one of the most important risk factors for developing diastolic dysfunction and heart failure. This is attributed to maladaptive concentric LV remodeling, disorganization of sarcomeres, alteration in calcium handling, and increase in extracellular matrix . Our HTN demonstrates concentric LV remodeling by decreased LVED R/T ratio and increased LV mass/volume ratio. Due to increased LV wall thickness with relatively greater epicardial radius compared to endocardial radius and other factors associated with the remodeling, there is likely a greater gradient of torque from epicardium to the endocardium in HTN vs. controls, resulting in a higher peak torsion shear angle in HTN that will therefore be associated with greater stored potential energy in the myocardium at ES. During diastole, some of the stored potential energy is used as work against the viscoelastic properties of myocardium and the rest likely contributes to LV filling and suction, myocardial relaxation and enhanced early diastolic filling ,,.
In the current study, we find that peak early diastolic is significantly increased in HTN, indicating reduced LV filling to the same changes of untwist in HTN vs. controls. This therefore indicates impaired LV suction and early filling in HTN. Moreover, we find significantly higher THA in HTN. Therefore the abnormal torsion shear angle volume relationship obtained by changes in torsion shear angle to LV volumetric changes over the cardiac cycle in HTN group likely indicates inefficiency in ventricular relaxation due to stiffer/relatively non-compliant ventricle. This is further supported by reduced early diastolic strain rates in HTN vs. controls.
A number of echocardiographic parameters have been proposed for evaluating diastolic function , such as normalized peak early filling rate, E/A ratio, normalized peak early diastolic mitral annulus velocity. In our study, THA is significantly correlated with CMR-derived E/A ratio. Twenty-six percent HTN patients (13 out of 60) had both E/A ratio < 1.12 (control’s mean E/A-1SD) and THA > 0.12 (control’s mean THA + 1SD) as shown in Figure 3. HTN patients with small E/A ratio are not highly overlapped with patients with large THA, indicating that other factors may play an important role in determining THA, which is not expressed in E/A ratio. It is now increasingly recognized that patients with preserved ejection fraction and diastolic dysfunction may demonstrate impaired systolic strains at rest or blunted systolic response with exertion. Early diastolic E/A ratio is purely a diastolic filling parameter. The loop framework we propose is unique in that it takes into account both the systolic and diastolic phases and provides insight into ventricular stiffness and suction which are important determinants of diastolic dysfunction.
Previous study by Takeuchi M, et al. reported decreased early diastolic LV untwisting rate during isovolumic relaxation period in hypertensive patients utilizing 2D speckle tracking imaging . Due to the relative lower temporal resolution of tagged CMR, peak untwist rate during isovolumic relaxation period is not available in our study. Instead, peak untwist rate during the early diastole phase is computed. However, in our study, we used widely available tagged CMR technique to define a new approach that wasn’t likely affected by the somewhat lower temporal resolution compared to echocardiography. MR tag lines faded during mid to late diastole and torsion shear angle cannot be reliably measured. In our study, we were able to evaluate 2/3 of the diastolic phase for all subjects. Since major portion of diastolic untwist happened quickly during early-diastolic phase and minimal untwisting occurred during late diastole, substantial change in current conclusion due to the limitation of our technique was unlikely. Regardless, our data demonstrates consistent differences amongst the two groups in the conventional non-invasive parameters frequently used for evaluating diastolic function.
In conclusion, loop may be an important global approach to assess diastolic function. To our knowledge, this is the first study to propose the concept of THA and early diastole as measures of diastolic functional property and their application may provide greater insight into heart failure. Future study will perform invasive hemodynamic measurements and simultaneous echocardiography for assessing routine measures of diastolic function and compare them to the loop derived measures to comprehensively assess all the factors responsible for this phenomenon.
Cardiovascular magnetic resonance
Torsion hysteresis area
LV end diastolic pressure volume relationship
Cardiovascular magnetic resonance
Brain natriuretic peptide
Zile MR, Baicu CF, Gaasch WH: Diastolic heart failure—abnormalities in active relaxation and passive stiffness of the left ventricle. N Engl J Med. 2004, 350: 1953-9. 10.1056/NEJMoa032566.
Redfield MM, Jacobsen SJ, Borlaug BA, Rodeheffer RJ, Kass DA: Age- and gender-related ventricular-vascular stiffening: a community-based study. Circulation. 2005, 112: 2254-62. 10.1161/CIRCULATIONAHA.105.541078.
Maurer MS, Spevack D, Burkhoff D, Kronzon I: Diastolic dysfunction: can it be diagnosed by doppler echocardiography?. J Am Coll Cardiol. 2004, 44: 1543-9. 10.1016/j.jacc.2004.07.034.
Ommen SR, Nishimura RA, Appleton CP, Miller FA, Oh JK, Redfield MM, Tajik AJ: Clinical utility of doppler echocardiography and tissue doppler imaging in the estimation of left ventricular filling pressures: a comparative simultaneous doppler-catheterization study. Circulation. 2000, 102: 1788-94. 10.1161/01.CIR.102.15.1788.
Tschope C, Paulus WJ: Is echocardiographic evaluation of diastolic function useful in determining clinical care? Doppler echocardiography yields dubious estimates of left ventricular diastolic pressures. Circulation. 2009, 120: 810-20. 10.1161/CIRCULATIONAHA.109.869628.
Young A, Cowan B: Evaluation of left ventricular torsion by cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2012, 14: 49-10.1186/1532-429X-14-49.
Bell SP, Nyland L, Tischler MD, McNabb M, Granzier H, LeWinter MM: Alterations in the determinants of diastolic suction during pacing tachycardia. Circ Res. 2000, 87: 235-40. 10.1161/01.RES.87.3.235.
Rademakers FE, Buchalter MB, Rogers WJ, Zerhouni EA, Weisfeldt ML, Weiss JL, Shapiro EP: Dissociation between left ventricular untwisting and filling. Accentuation by catecholamines. Circulation. 1992, 85: 1572-81. 10.1161/01.CIR.85.4.1572.
Dong S-J, Hees PS, Siu CO, Weiss JL, Shapiro EP: MRI assessment of lv relaxation by untwisting rate: a new isovolumic phase measure of tau. Am J Physiol Heart Circ Physiol. 2001, 281: H2002-9.
Burns AT, La Gerche A, Prior DL, MacIsaac AI: Left ventricular untwisting is an important determinant of early diastolic function. J Am Coll Cardiol Img. 2009, 2: 709-16. 10.1016/j.jcmg.2009.01.015.
Gaddam K, Corros C, Pimenta E, Ahmed M, Denney T, Aban I, Inusah S, Gupta H, Lloyd SG, Oparil S, Husain A, Dell'Italia LJ, Calhoun DA: Rapid reversal of left ventricular hypertrophy and intracardiac volume overload in patients with resistant hypertension and hyperaldosteronism a prospective clinical study. Hypertension. 2010, 55: 1137-42. 10.1161/HYPERTENSIONAHA.109.141531.
Castillo E, Osman NF, Rosen BD, El-Shehaby I, Pan L, Jerosch-Herold M, Lai S, Bluemke DA, Lima JA: Quantitative assessment of regional myocardial function with mr-tagging in a multi-center study: interobserver and intraobserver agreement of fast strain analysis with harmonic phase (harp) mri. J Cardiovasc Magn Reson. 2005, 7: 783-91. 10.1080/10976640500295417.
Feng W, Nagaraj H, Gupta H, Lloyd SG, Aban I, Perry GJ, Calhoun DA, Dell'Italia LJ, Denney TS: A dual propagation contours technique for semi-automated assessment of systolic and diastolic cardiac function by cmr. J Cardiovasc Magn Reson. 2009, 11: 30-10.1186/1532-429X-11-30.
Schiros CG, Dell'Italia LJ, Gladden JD, Clark D, Aban I, Gupta H, Lloyd SG, McGiffin DC, Perry G, Denney TS, Ahmed MI: Magnetic resonance imaging with 3-dimensional analysis of left ventricular remodeling in isolated mitral regurgitation: Implications beyond dimensions. Circulation. 2012, 125: 2334-42. 10.1161/CIRCULATIONAHA.111.073239.
Osman NF, Kerwin WS, McVeigh ER, Prince JL: Cardiac motion tracking using cine harmonic phase (harp) magnetic resonance imaging. Magn Reson Med. 1999, 42: 1048-60. 10.1002/(SICI)1522-2594(199912)42:6<1048::AID-MRM9>3.0.CO;2-M.
Khalifa A, Youssef AB, Osman N: Improved harmonic phase (harp) method for motion tracking a tagged cardiac mr images. Conf Proc IEEE Eng Med Biol Soc. 2005, 4: 4298-301.
Hill JA, Olson EN: Cardiac plasticity. N Engl J Med. 2008, 358: 1370-80. 10.1056/NEJMra072139.
Yun KL, Niczyporuk MA, Daughters GT, Ingels NB, Stinson EB, Alderman EL, Hansen DE, Miller DC: Alterations in left ventricular diastolic twist mechanics during acute human cardiac allograft rejection. Circulation. 1991, 83: 962-73. 10.1161/01.CIR.83.3.962.
Oh JK, Hatle L, Tajik AJ, Little WC: Diastolic heart failure can be diagnosed by comprehensive two-dimensional and doppler echocardiography. J Am Coll Cardiol. 2006, 47: 500-6. 10.1016/j.jacc.2005.09.032.
Takeuchi M, Borden WB, Nakai H, Nishikage T, Kokumai M, Nagakura T, Otani S, Lang RM: Reduced and delayed untwisting of the left ventricle in patients with hypertension and left ventricular hypertrophy: a study using two-dimensional speckle tracking imaging. Eur Heart J. 2007, 28: 2756-62. 10.1093/eurheartj/ehm440.
NIH NHLBI R01-HL104018
The authors declare that they have no competing interests.
CGS and HG were responsible for conception and design of the paper, and manuscript drafting and finalizing. CGS, BAV and TSD were responsible for tagged image analysis, algorithm development and data generation. DAC, LJD’I, SGL, and HG were responsible for patient population design and recruiting. CGS, LJD’I, TSD, SGL, RVD, KKG and SA were responsible for data generation and interpretation. CGS, TSD, SGL, HG were responsible for revising paper critically for important intellectual content. Final approval of the manuscript was done by HG. All authors read and approved the final manuscript.
About this article
Cite this article
Schiros, C.G., Desai, R.V., Venkatesh, B.A. et al. Left ventricular torsion shear angle volume analysis in patients with hypertension: a global approach for LV diastolic function. J Cardiovasc Magn Reson 16, 70 (2014). https://doi.org/10.1186/s12968-014-0070-4
- Torsion hysteresis area
- Diastolic function
- Cardiovascular magnetic resonance