Mapping right ventricular myocardial mechanics using 3D cine DENSE cardiovascular magnetic resonance
© Auger et al; licensee BioMed Central Ltd. 2012
Received: 13 September 2011
Accepted: 11 January 2012
Published: 11 January 2012
The mechanics of the right ventricle (RV) are not well understood as studies of the RV have been limited. This is, in part, due to the RV's thin wall, asymmetric geometry and irregular motion. However, the RV plays an important role in cardiovascular function. This study aims to describe the complex mechanics of the healthy RV using three dimensional (3D) cine displacement encoding with stimulated echoes (DENSE) cardiovascular magnetic resonance (CMR).
Whole heart 3D cine DENSE data were acquired from five healthy volunteers. Tailored post-processing algorithms for RV mid-wall tissue tracking and strain estimation are presented. A method for sub-dividing the RV into four regions according to anatomical land marks is proposed, and the temporal evolution of strain was assessed in these regions.
The 3D cine DENSE tissue tracking methods successfully capture the motion and deformation of the RV at a high spatial resolution in all volunteers. The regional Lagrangian peak surface strain and time to peak values correspond with previous studies using myocardial tagging, DENSE and strain encoded CMR. The inflow region consistently displays lower peak strains than the apical and outflow regions, and the time to peak strains suggest RV mechanical activation in the following order: inflow, outflow, mid, then apex.
Model-free techniques have been developed to study the myocardial mechanics of the RV at a high spatial resolution using 3D cine DENSE CMR. The consistency of the regional RV strain patterns across healthy subjects is encouraging and the techniques may have clinical utility in assessing disrupted RV mechanics in the diseased heart.
Keywordsdisplacement encoding DENSE right ventricle cardiac function strain
Right ventricular (RV) function may be impaired in a number of heart conditions, including myocardial infarction, congenital heart disease and cardiomyopathy . The function of the RV may also be affected in diseases of the left ventricle (LV) where it is difficult to ignore the complex nature of ventricular interaction [2, 3]. In the past, the importance of the LV in cardiac research has overshadowed the study of the RV. This neglect is, in part, because the RV is difficult to image. The wall of the RV myocardium is thin (2-5 mm) when compared to that of the LV (7-11 mm) . Furthermore, the RV has a complex geometry, eccentric motion [1, 4] and it is heavily trabeculated, thus it does not offer the clearly defined endocardial margins typically seen in the LV.
Various cardiovascular magnetic resonance (CMR) studies have assessed properties such as strain, motion and volumes of the RV. The CMR techniques used include turbo gradient echo and steady state free precession (SSFP) , myocardial tagging [4, 6, 7], phase contrast velocity encoding , strain encoded (SENC) CMR  and 2D multi-slice displacement encoding with stimulated echoes (DENSE) .
The standard cardiac imaging planes are based on the relatively uncomplicated geometry of the LV. However, no standardized localization exists for the RV. Two dimensional strain estimates of the LV are reliable in short axis views because the through plane motion is relatively uniform . This does not apply to the RV because of its complex shape and motion, so 2D imaging may be insufficient to accurately assess strain in the RV . To date, no detailed studies of the RV have been presented using DENSE CMR.
DENSE is a quantitative CMR technique used for measuring myocardial displacement and strain. DENSE encodes tissue displacement directly into the image phase (typically with reference to end diastole), thus allowing for the extraction of motion data at a pixel resolution [12, 13]. A recently developed free breathing 3D cine DENSE sequence  is well suited for quantifying the complex behavior of the RV. A high signal to noise ratio (SNR), which is necessary for imaging the thin RV walls, is achieved using a spiral k-space trajectory and three point phase cycling. Furthermore, DENSE is inherently a black-blood technique, which together with fat suppression, provides a more accurate delineation of the RV walls.
The purpose of this study was to develop tailored processing techniques for assessing detailed 3D RV motion and surface strain using 3D cine DENSE, and to quantify these parameters for the healthy human heart.
Whole heart 3D cine DENSE data were acquired from five healthy male volunteers (age range 21 - 45) on a 1.5T MRI scanner (Siemens MAGNETOM Avanto, Erlangen, Germany) using a two-channel anterior body array coil and an eight-channel spine array coil.
The entire heart was imaged with the imaging volume aligned along the heart's long axis and at a 2.8 × 2.8 × 5 mm3 spatial resolution and 32 ms temporal resolution. Displacement was encoded in three orthogonal directions and a spiral k-space trajectory was used with 6 interleaves per 3D partition. Spiral aliasing artifacts existed on the edge of images, but not on the heart, which was placed at the center of the field of view. Images were acquired during a 20.5 ± 5.7 min scan time with prospective ECG gating, diaphragmatic navigator respiratory gating, and no multichannel parallel imaging. Other imaging parameters include: field of view (FOV) = 360 × 360 × 140 mm3 (in-plane dimensions defined by the diameter of the spiral), displacement encoding frequency ke = 0.06 cycles/mm, ramped flip angle up to 20 degrees, TR = 16 ms, TE = 1.3 ms, cardiac phases = 22 and number of acquired 3D partitions = 14 (zero-padded to 28 during image reconstruction). After the Fourier Transform in the partition direction, three partitions at each end of the volume were discarded to avoid aliasing. The resulting 3D image matrix size was 128 × 128 × 22. The number of partitions chosen per subject in order to include the RV was approximately 14, which corresponds to an overall stack thickness of 70 mm. All imaging was conducted after informed consent and in accordance with protocols approved by the University of Virginia Institutional Review Board.
A number of tailored post processing steps were implemented to analyze the RV myocardial DENSE data. All software development was performed using MATLAB (The Mathworks, Natick, MA).
Contouring, phase unwrapping and tissue tracking
Both the RV and LV myocardium was demarcated from surrounding structures by manually drawing each epicardial and endocardial contour on the reconstructed short-axis DENSE magnitude images.
There are typically too few transmural pixels spanning the RV to acquire a full 3D strain tensor, so 2D strain at the RV mid-wall surface was calculated. The mid-line between epicardial and endocardial contours was used to produce tissue tracking seed points, which were spaced at pixel-distance intervals on the first cardiac phase. The mid-line contours were spatially smoothed within and across partitions using 4th order polynomial fitting to ensure a continuous RV surface. An example of the RV mid-line points is shown in green in Figure 1(a).
A 2D tissue tracking algorithm involving interpolation of the displacement fields  was directly extended to 3D for this application. The position of each mid-line point along its motion trajectory throughout the cardiac cycle was estimated using 3D distance weighted linear interpolation. The full 3D mid-line motion trajectories were calculated by subtracting the interpolated vectors of successive frames from one another. Prospective gating only allowed about 90 percent of the cardiac cycle to be imaged, therefore, temporal fitting was done for each ordinate direction of each of the trajectories using a 10th order polynomial, as 5th order Fourier basis functions have been shown to sufficiently describe cardiac motion .
Lagrangian surface strain
where T denotes the transpose operation, and I represent the identity matrix. Note that this construction of a 3D tensor using a unit normal vector assumes that the cardiac tissue is non-compressible and has been implemented here to simplify the strain calculations. The 3D strain tensor was decomposed into its corresponding eigenvalues (strain) and eigenvectors (direction of deformation). The eigenvalue corresponding to the orthogonal unit vector (E1) is aligned with the direction normal to the surface (radial) and is zero because of the aforementioned unit normal vector assumption. The remaining two eigenvalues (E2 and E3) provide an estimate of the 2D tangential surface strain. The eigenvalues E2 and E3 are of similar magnitude, and there are regional differences in the myofiber (and eigenvector) orientation across the RV surface, so we were unable to discern between E2 and E3. The results of the surface strain analysis are thus presented as a mean Lagrangian 2D principal strain, which has been averaged for the various co-planar vector configurations.
It has been shown that tangential strain in the circumferential and longitudinal directions are good indicators of LV dyssynchrony . In order to infer details about the strain directions in the RV, one dimensional (1D) strains were calculated in the longitudinal and circumferential directions in the same manner as usually done for the LV, and relative to the orientation of the LV. Each 1D strain was calculated from the motion trajectories by only considering points above and below (longitudinal) or to the left and right (circumferential) of each mid-line point. This is also illustrated in Figure 2.
Anatomical sub-divisions for strain-time analyses
Unlike the LV, there is currently no standard method for dividing the RV into anatomical sub-regions. Previous studies of the RV using myocardial tagging have adopted different conventions for sub-dividing the RV. Klein et al. divides each RV basal, mid and apical short axis slice into three regions: superior, mid and inferior . Haber et al. , divided the RV into four regions according to anatomical landmarks. The parietal and septal bands were used to demarcate the outflow region, while the free wall was identified based on the normalized height of the septum in a long axis view . Fayad et al. divides the RV using the supraventricular crest (SC), moderator band, papillary muscles and the tricuspid valve as landmarks, creating four distinct regions: the inflow, outflow, mid and apical regions .
Coordinate 1: The SC is an accentuation of the muscular wall demarcating the inflow and outflow tracts within the RV. Navigate from apex to base along the short axis magnitude images until the SC appears in the image. The first coordinate is marked on the RV free wall at the level of and adjacent to the SC.
Coordinates 2 and 3: The moderator band is a muscular band connecting the interventricular septum to the anterior papillary muscle of the tricuspid valve. Navigate from base to apex along the short axis magnitude images and define the slice where the moderator band becomes evident as a segment of myocardium crossing the RV cavity. The second and third coordinates are defined as the anterior and inferior RV-LV insertion points at this level, respectively.
The apical RV region is defined as the myocardium lying apical to the slice identified by the moderator band. The inflow, mid and outflow RV regions are defined using two planes: (i) a plane defined by coordinate 1, coordinate 2 and a point displaced normal to the RV mid-line surface in this region, and (ii) a plane defined by coordinate 1, coordinate 3 and a point displaced normal to the RV mid-line surface in this region. The strain values (2D and 1D Lagrangian) were averaged in each of the 4 segments and assessed through the duration of the cardiac cycle for all 5 volunteers.
Right ventricular regional Lagrangian strain and time to peak strain
Strain quantification method and anatomical region
Time to peak (ms)
-0.10 ± 0.04
-0.15 ± 0.05
-0.17 ± 0.06
-0.16 ± 0.06
-0.16 ± 0.06
-0.22 ± 0.08
-0.20 ± 0.07
-0.18 ± 0.06
Mean principal strain
-0.12 ± 0.05
-0.18 ± 0.07
-0.17 ± 0.06
-0.16 ± 0.05
All strain values are consistently negative, indicating muscle shortening or contraction, and the longitudinal shortening is greater for all regions when compared to circumferential shortening.
The inflow region has the lowest peak strain value of all the regions, while the outflow region generally reaches the highest peak strain. However, Figure 6 illustrates that during tangential shortening in the circumferential direction, the apical region tends to have a slightly greater strain than the outflow region. Good consistency between the patterns of strain evolution is shown for the various strain estimates in Figures 5 and 6. The four regions each arrive at a maximum strain at slightly different time points in the cardiac cycle. The time sequence of regional myocardial shortening, as confirmed by all strain estimates is: inflow, outflow, mid then apical.
Additional file 1: Three dimensional mid-line motion trajectories in a plane view of the LV and RV seeded from the mid-line points, for each partition/slice. Motion trajectories begin at end diastole. The dots represent the motion trajectory positions, while the color represents mean principal 2D Lagrangian surface strain. (MPEG 824 KB)
Additional file 2: Three dimensional mid-line motion trajectories in an oblique view of the LV and RV seeded from the mid-line points, for each partition/slice. Motion trajectories begin at end diastole. The dots represent the motion trajectory positions, while the color represents mean principal 2D Lagrangian surface strain. (MPEG 819 KB)
This paper introduces post processing algorithms developed for 3D tissue tracking and strain analysis of the RV from 3D cine DENSE data. Unlike previous studies of 3D RV motion using tagging [4, 7], no models were used to describe the underlying data. Furthermore, no spatial smoothing was applied to the data, and the only form of filtering was the temporal polynomial fitting applied to the mid-line motion trajectories. This temporal fitting provides a more realistic trajectory behavior and allows for motion to be estimated at any temporal resolution.
Both the peak strain and the time to peak strain were found to vary across the different RV regions. The inflow region consistently demonstrated lower strains and early strain peak times, while the apex consistently showed larger strain values and later peak times. Some variation does exist in the outflow region when comparing the different strain calculation types.
As there is no current standard for dividing the RV, and multiple approaches for calculating strain have been adopted, it is difficult to make a direct comparison with results from previous studies. In general, however, the magnitude trends and timing parameters of the strain measurements from this study compare favorably with previous RV studies using myocardial tagging, SENC and DENSE. From our study the peak mean principal strains values in the apical, mid and outflow regions are -0.17, -0.16 and -0.18, respectively. These values are similar to the corresponding values of -0.17, -0.19 and -0.22 reported in a tagging study by Haber et al. .
In the SENC study by Hamdan et al. , they show longitudinal strains to be largest at the base, with a lower strain at the apex, and the lowest strain in the mid region. For circumferential strain they show the largest strain at the apex, with a lower strain in the mid and basal regions. These regional strain comparisons show excellent agreement with those presented in Table 1. In our results, during longitudinal contraction, the largest strain is found in the outflow region (basal), followed by the apical and mid regions, while circumferential strain in the apical region is larger, followed by the mid and outflow (basal) regions.
The 1D peak longitudinal strains in Table 1 are all larger than the corresponding 1D peak circumferential and mean principal strains, a result confirmed in tagging studies by both Haber et al  and Fayad et al. . Fayad et al., who used a very similar method to ours for sub-dividing the RV, found regional peak circumferential strains to range from largest to smallest in the following order: apical, mid, inflow and outflow. For peak longitudinal strains the corresponding order was: outflow, inflow, apical and mid regions. The peak strain values in Table 1 also correlate well with the segmental shortening presented by Fayad et al., except for the inflow region where our values are consistently lower.
Wen et al.  present circumferential strain results from 2D DENSE data as an average of the entire RV free wall, in three slices: basal, mid and apical. Here, the strain is derived from the lengthening of the RV contour. They also found decreasing peak circumferential strain from apex to base, but the results plotted in  have a consistently larger magnitude than those presented here.
Limitations of our study include a lengthy scan time, the need for manual contouring, and the difficulty to discern between the minimum and intermediate principal strains. A reduction in the scan time could be accomplished by incorporating parallel imaging and/or outer volume suppression. The time taken to contour the images may be considerably reduced with the extension of previous work on motion guided segmentation for 2D cine DENSE images . Even with the use of the eigenvectors, discerning between the two principal strains was not possible. A geometric or finite element model of the RV could solve this problem, whereby the equivalent of an RLC coordinate system for the LV is defined according to the RV geometry.
The 2D surface strain estimates are expected to have a lower accuracy in the through-plane direction because the through-plane image resolution is lower than the in-plane resolution. This should lead to a less accurate representation of strain in the longitudinal direction compared to the circumferential direction, but this is not evident in the standard deviation curves in Figure 6.
An inherent limitation of DENSE is the reduction in SNR associated with the stimulated echo. The SNR was typically sufficient to clearly identify the RV because of the short echo times associated with the spiral readout. However, signal loss due to through-plane dephasing is particularly pronounced in the RV because of its large through-plane strain . This reduction in SNR was evident in a few instances, and these data had to be excluded from the analyses. For future work on the RV, this through-plane dephasing could be minimized by lowering the displacement encoding frequency, but at the expense of displacement sensitivity.
Three dimensional cine DENSE CMR is a promising technique for assessing RV mechanics. The methods presented here enable the motion and deformation of the entire RV to be captured at a pixel resolution. The RV mechanics measured by 3D cine DENSE is comparable to that of previous studies using tagging, SENC and 2D DENSE.
We further describe a method for dividing RV into regions according to anatomical land marks. The differing temporal strain evolution in the proposed RV anatomical sub-divisions indicates that these regions may also be functionally distinct. Future work will involve applying these techniques to studying variations in regional RV function in cardiac disease.
The authors would like to thank Aaron Hess, Ph.D. for the helpful discussions about the surface strain algorithm, and Leon Axel M.D., Ph.D. for his valuable and insightful recommendations on the RV regional anatomical divisions. The authors further appreciate the clinical input from Bongani Mayosi M.D., Ph.D.
- Oldershaw P: Assessment of right ventricular function and its role in clinical practice. Br Heart J. 1992, 68: 12-15. 10.1136/hrt.68.7.12.PubMed CentralView ArticlePubMedGoogle Scholar
- 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-1448. 10.1161/CIRCULATIONAHA.107.653576.View ArticlePubMedGoogle Scholar
- McKenna WJ, Thiene G, Nava A, Fontaliran F, Blomstrom-Lundqvist C, Fontaine G, Camerini F: Diagnosis of arrhythmogenic right ventricular dysplasia/cardiomyopathy. Br Heart J. 1994, 71: 215-218. 10.1136/hrt.71.3.215.PubMed CentralView ArticlePubMedGoogle Scholar
- Haber I, Metaxas DN, Geva T, Axel L: Three-dimensional systolic kinematics of the right ventricle. Am J Physiol Heart Circ Physiol. 2005, 289: 1826-1833. 10.1152/ajpheart.00442.2005.View ArticleGoogle Scholar
- Alfakih K, Plein S, Thiele H, Jones T, Ridgway JP, Sivananthan MU: Normal human left and right ventricular dimensions for MRI as assessed by turbo gradient echo and steady-state free precession imaging sequences. J Magn Reson Imaging. 2003, 17: 323-329. 10.1002/jmri.10262.View ArticlePubMedGoogle Scholar
- Haber I, Metaxas DN, Axel L: Three-dimensional motion reconstruction and analysis of the right ventricle using tagged MRI. Med Imaging Analysis. 2004, 4: 335-355.View ArticleGoogle Scholar
- Young AA, Fayad ZA, Axel L: Right ventricular midwall surface motion and deformation using magnetic resonance tagging. Am J Physiol. 1996, 271 (6 Pt 2): H2677-H2688.PubMedGoogle Scholar
- Kayser HWM, van der Geest RJ, van der Wall EE, Duchateau C, de Roos A: Right Ventricular Function in Patients After Acute Myocardial Infarction Assessed With Phase Contrast MR Velocity Mapping Encoded in Three Directions. J Magn Reson Imaging. 2000, 11: 471-475. 10.1002/(SICI)1522-2586(200005)11:5<471::AID-JMRI1>3.0.CO;2-D.View ArticlePubMedGoogle Scholar
- Paetsch I, Gebker R, Wellnhofer E, Schnackenbur B, Fahmy AS, Osman NF, Fleck E, Hamdan A, Thouet T, Sebastian K: Regional right ventricular function and timing of contraction in healthy volunteers evaluated by strain-encoded MRI. J Magn Reson Imaging. 2008, 28: 1379-1385. 10.1002/jmri.21526.View ArticlePubMedGoogle Scholar
- Wen H, Marsolo KA, Bennett EE, Kutten KS, Lewis RP, Lipps DB, Epstein ND, Plehn JF, Croisille P: Adaptive Postprocessing Techniques for Myocardial Tissue Tracking with Displacement-encoded MR Imaging. Radiology. 2008, 246 (1): 229-240.PubMed CentralView ArticlePubMedGoogle Scholar
- Tecelão SRR, Zwanenburg JJM, Kuijer JPA, de Cock CC, Germans T, van Rossum AC, Marcus JT: Quantitative comparison of 2D and 3D circumferential strain using MRI tagging in normal and LBBB hearts. Magn Reson Med. 2007, 57: 485-493. 10.1002/mrm.21142.View ArticlePubMedGoogle Scholar
- Aletras AH, Ding S, Balaban RS, Wen H: DENSE: Displacement Encoding with Stimulated Echoes in Cardiac Functional MRI. J Magn Reson. 1999, 137: 247-252. 10.1006/jmre.1998.1676.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim D, Gilson WD, Kramer CM, Epstein FH: Myocardial tissue tracking with two-dimensional cine displacement-encoded MR imaging: development and initial evaluation. Radiology. 2004, 230: 862-871. 10.1148/radiol.2303021213.View ArticlePubMedGoogle Scholar
- Zhong X, Spottiswoode BS, Meyer CH, Kramer CM, Epstein FH: Imaging Three-Dimensional Myocardial Mechanics Using Navigator-Gated Volumetric Spiral Cine DENSE MRI. Magn Reson Med. 2010, 64: 1089-1097. 10.1002/mrm.22503.PubMed CentralView ArticlePubMedGoogle Scholar
- Spottiswoode BS, Zhong X, Hess AT, Kramer CM, Meintjes EM, Mayosi BM, Epstein FH: Tracking myocardial motion from cine DENSE images using spatiotemporal phase unwrapping and temporal fitting. IEEE Trans Med Imaging. 2007, 26: 15-30.View ArticlePubMedGoogle Scholar
- Young AA: Epicardial deformations from coronary cineangiograms. Theory of the Heart: Biomechanics, Biophysics, and Nonlinear Dynamics of Cardiac Function. Edited by: Glass LL, Hunter P, McCulloch A. 1991, New York: Springer-Verlag, 175-207.View ArticleGoogle Scholar
- Helm RH, Leclercq C, Faris OP, Ozturk C, McVeigh E, Lardo AC, Kass DA: Cardiac dyssynchrony analysis using circumferential versus longitudinal strain: implications for assessing cardiac resynchronization. Circulation. 2005, 111: 2760-2767. 10.1161/CIRCULATIONAHA.104.508457.PubMed CentralView ArticlePubMedGoogle Scholar
- Klein SS, Graham TP, Lorenz CH: Noninvasive delineation of normal right ventricular contractile motion with magnetic resonance imaging myocardial tagging. Annals Biomed Eng. 1998, 26: 756-63.View ArticleGoogle Scholar
- 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. J Magn Reson Med. 1998, 39: 116-123. 10.1002/mrm.1910390118.View ArticleGoogle Scholar
- Spottiswoode BS, Zhong X, Lorenz CH, Mayosi BM, Meintjes EM, Epstein FH: Motion-guided segmentation for cine DENSE MRI. Med Image Anal. 2009, 13: 105-115. 10.1016/j.media.2008.06.016.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.