Radiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Tseng, W.-Y. I.
Right arrow Articles by Wedeen, V. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Tseng, W.-Y. I.
Right arrow Articles by Wedeen, V. J.
(Radiology. 2000;216:128-139.)
© RSNA, 2000


Cardiac Imaging

Myocardial Fiber Shortening in Humans: Initial Results of MR Imaging1

Wen-Yih I. Tseng, MD, PhD, Timothy G. Reese, PhD, Robert M. Weisskoff, PhD, Thomas J. Brady, MD and Van J. Wedeen, MD

1 From the Department of Radiology, National Taiwan University Medical College, Laser Medical Research Center, Taipei, Taiwan (W.Y.I.T.); and the Department of Radiology, NMR Center, Massachusetts General Hospital-East, Bldg 149, 13th St, Charlestown, MA 02129 (T.G.R., R.M.W., T.J.B., V.J.W). Received July 16, 1999; revision requested August 17; revision received October 7; accepted October 25. Supported by National Institutes of Health 1RO1-HL56737 grant (V.J.W., W.Y.I.T.), the American Heart Association Established Investigator Grant 9740208N (V.J.W.), and grants from the New York Cardiac Center and the Sol Goldman Charitable Trust (V.J.W.). Address correspondence to V.J.W. (e-mail: van@nmr.mgh.harvard.edu).


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 REFERENCES
 
PURPOSE: To use diffusion-sensitive magnetic resonance (MR) imaging to obtain images of fiber orientation in vivo and to map fiber shortening in humans by means of integrating such data with strain images.

MATERIALS AND METHODS: Images of fiber shortening for midventricular short-axis sections were acquired in eight healthy subjects. Fiber orientation maps obtained by means of diffusion-sensitive MR imaging were coregistered with systolic strain maps obtained by means of velocity-sensitive MR imaging. Fiber shortening was quantified by use of the component of systolic strain in the fiber direction.

RESULTS: The results were reproducible among subjects and were consistent with published values. MR imaging of myocardial fibers showed axisymmetric progression of fiber angles from -90° epicardially to +90° endocardially, with maxima near 0°. Fiber shortening (mean, 0.12 ± 0.01 [SD]) was more uniform than radial, circumferential, longitudinal, or cross-fiber strain or any principal strain. Fiber orientation coincided with the direction of maximum contraction epicardially, with that of minimum contraction endocardially, and varied between these extremes linearly with wall depth (r = 0.6).

CONCLUSION: Registered diffusion and strain MR imaging can be used quantitatively to map fiber orientation and its relations to myocardial deformation in humans.

Index terms: Heart, function, 51.91, 51.92 • Heart, MR, 51.121411, 51.121416, 51.12144 • Magnetic resonance (MR), diffusion study, 51.12144 • Magnetic resonance (MR), phase imaging, 51.121411, 51.121416 • Myocardium, MR, 511.121411, 511.121416, 511.12144


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 REFERENCES
 
Myocardial function at the local scale emerges from a relationship between fiber orientation and local wall deformation (1,2). We recently described a methodology based on diffusion-sensitive magnetic resonance (MR) imaging to map myocardial fiber architecture in vivo (3,4). In the present work, this method is combined with MR imaging of local wall deformation to provide images of myocardial fiber shortening and of its relation to myocardial strain in healthy human subjects. This MR imaging methodology, by providing data that are (a) three-dimensional (3D), (b) spatially resolved, (c) able to be used to define fiber orientation and wall deformation with identical cardiac configuration and in identical physiologic conditions, and (d) obtainable in humans noninvasively affords an opportunity to refine and extend analyses of myocardial fiber shortening traditionally assessed by means of invasive methods.

The present study had three aims: first, to define an MR imaging methodology to obtain registered images of myocardial structure and function and derive from them quantitative maps of myocardial fiber shortening and related parameters; second, to validate this method in a series of healthy subjects by means of showing reproducibility and agreement with published values obtained by using conventional methods; and third, to show that MR imaging affords improved recognition of general characteristics of cardiac function.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 REFERENCES
 
Overview
The present method is based on joint analysis of registered MR images of myocardial fiber orientation, represented by the myocardial diffusion tensor field D, and of local wall deformation, represented by the myocardial systolic strain tensor field Ssys. Structure-function relations between these two tensor fields may be converted into quantitative images. First among these is the myocardial fiber shortening SF. Fiber orientation f is equated with the direction of maximum diffusion, which is the leading eigenvector of the diffusion tensor (5,6). Fiber shortening is then the component of the systolic strain tensor along the fiber:

The immediate goal of the present MR imaging acquisition is therefore to obtain accurate, spatiotemporally registered images of the myocardial diffusion tensor D and systolic strain tensor Ssys.

The image acquisition and reconstruction procedure is as follows: (a) Acquire an MR movie of myocardial two-dimensional (2D) strain rates. (b) Compute from the strain-rate movie the optimum cardiac phase delay {Psi} for diffusion MR imaging. (c) At t = {Psi}, acquire an image of the cardiac diffusion tensor D and a registered image of the myocardial 3D strain-rate tensor S'({Psi}). (d) By using the image of 3D strain rate S'({Psi}) and the strainrate movie, compute the tensor image of net systolic 3D strain Ssys. (e) Compute from the tensor images D and Ssys myocardial fiber shortening and associated quantities.

The details concerning computation of optimum cardiac phase delay in step b and of the net systolic strain in step d are given in the Appendix.

MR Imaging Acquisitions
Eight healthy volunteers (six men, two women; age range, 26–40 years) without previous history of heart disease were recruited. Each subject provided written informed consent, which was approved by the hospital human study committee. All experiments were performed with a 1.5-T MR imager (Signa; GE Medical Systems, Milwaukee, Wis) with Instascan echo-planar imaging. After electrocardiographic leads were applied, a 15 x 28-cm rectangular receive-only radio-frequency surface coil was positioned precordially, and the subject was put in the magnet lying in a 45° left anterior oblique orientation.

To suppress the effects of respiratory motion, all of the images were acquired at the end of expiration by using synchronized breathing. With imaging electrocardiographically triggered to every fifth heartbeat, subjects were requested to take a breath after the sound of each imaging pulse, then passively exhale, await the next image acquisition, and repeat. The resultant respiratory rate was within normal limits and could be sustained comfortably for many minutes. Both diffusion and strain-rate images were spatially encoded by using single-shot echo-planar imaging with a 20 x 40-cm field of view and a 64 x 128 image matrix, with a resultant in-plane resolution of 3 x 3 mm.

After a sagittal T1-weighted study was performed to establish left ventricular location and axis, a midventricular short-axis section was defined for the studies of myocardial strain rate (a 2D strain-rate movie to determine {Psi}, followed by 3D strain-rate tensor imaging at t = {Psi}) and the registered study of myocardial diffusion. Registration of the diffusion image and the strain-rate image required temporal alignment among three corresponding components of the respective images: (a) section-selection for section registration, (b) diffusion and velocity encoding pulses for registration of architecture and function, and (c) readout period for in-plane registration. The registered pulse sequences of diffusion and strain rate are illustrated in Figure 1. To ensure coregistration, the same phase-encoding direction and readout direction were used in diffusion and strain-rate imaging sequences.



View larger version (16K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 1. Diagram shows pulse sequences for acquisition of (top) registered diffusion MR images and (bottom) strain-rate MR images. The diffusion image is obtained with use of a stimulated echo with two electrocardiographic (ECG) triggers to locate the diffusion-encoding gradient pulses gD at identical phase delays {Psi} in consecutive cardiac cycles. The strain-rate image is obtained with use of a velocity-sensitive phase-contrast spin echo. Spatial and temporal registration of diffusion and strain images requires equal delays from the R wave of the section-select pulses (last 90° pulse of the diffusion image and 90° pulse of the strain-rate image), of the diffusion- and velocity-encoding gradient pulses, and of readout echo centers. EPI = echo-planar imaging.

 
Strain-Rate MR Imaging
Strain-rate MR imaging reconstructs from velocity-sensitive phase-contrast MR imaging an image in which the value at each location is the myocardial 2D or 3D strain-rate tensor, given by a symmetric 2 x 2 or 3 x 3 matrix, respectively, at each pixel (7). Images of myocardial strain-rate tensor fields are computed from the velocity-sensitive data by using

where vi and Sij' are the cartesian components of the velocity field v and strain-rate tensor field S', i and j are coordinate indices, and {delta}xi is the differential between pixels adjacent in the xi direction. The strain-rate image obtained with use of Equation (2) has spatial resolution identical to that of the original velocity-sensitive phase-contrast images. Images of 2D strain rates were computed by using velocity-sensitive data for a single section, whereas 3D strain rates required velocity-sensitive data for two contiguous sections offset in the through-plane, or z direction, to define z components of strain.

Strain-rate MR imaging data were acquired by using a single-shot spin-echo echo-planar imaging sequence augmented by velocity-encoding bipolar gradient pulses. Velocity-encoding gradient pulses were applied with amplitude |gv| = 10 mT · m-1 and pulse duration {delta}v = 6 msec for a velocity sensitivity |kv{approx} 130 radian · m-1, with tetrahedral orientations corresponding to the nonopposed corners of a cube (8). Velocity sensitivity was selected to yield myocardial phase shifts that had differentials that did not exceed {pi}/2 radians per 3-mm voxel at peak ejection rate. A movie of 2D strain rates was acquired with a 3.0-mm section thickness, progressive cardiac delay of 30 msec, four velocity encoding directions, one signal acquired, one section, a repetition time of five R-R intervals with synchronized breathing, and an echo time of 48 msec, for an imaging time of 7–9 minutes.

After determination of the optimum time delay t = {Psi} (Appendix), an image of the 3D strain-rate tensor field S'({Psi}) was acquired at this same section by using two contiguous 3.0-mm sections with opposite 1.5-mm offsets from the section center. For this acquisition, we used four velocity-encoding directions, four signals acquired, two sections, and a repetition time of five R-R intervals for an imaging time of about 2 minutes, which typically yielded a signal-to-noise ratio of at least 40:1 for strain rates. This corresponded to a root-mean-square error of approximately 11° in the principal thickening direction and of approximately 13° in the principal shortening direction.

Diffusion MR Imaging
Diffusion tensor MR imaging reconstructs the symmetric diffusion tensor at each location in an MR image. It is based on the diffusion-dependent reduction in the MR signal intensity that results from the application of a reversible spatial modulation of magnetization by a pair of temporally separated magnetic gradient pulses of opposite signs (9). The measured signal intensities are related to the diffusion tensor D by

where Io and Ik represent the signal intensities of the nonattenuated and attenuated images, respectively; {Delta} is the diffusion time; and kD is the gradient-induced spatial modulation of magnetization produced in a time {delta}D << {Delta}. By measuring this attenuation for spatial modulations in six directions and an image of null gradient (kD {approx} 0), we computed the diffusion tensor at each pixel by using standard algebraic inversion of Equation (3) (10).

In the present study, a double-gated stimulated-echo pulse sequence was used to acquire diffusion tensor MR images at the required time delay t = {Psi} with 6.0-mm section thickness. Diffusion-sensitizing gradient pulses were applied in six nonopposed edge-centers of a cube, |1,±1,0|, |0,1,±1|, and |±1,0,1|, and were of amplitude |gD| = 10 mT · m-1 duration {delta}D = 8.6 msec, which corresponded to a spatial modulation

{gamma} being the proton gyromagnetic ratio, and diffusion sensitivity b = {int}{Delta} |kD|2 d{tau} {approx} |kD|2 {Delta} {approx} 420 sec · mm-2, given a diffusion time {Delta} {approx} 800 msec (one R-R interval). With 16 signals acquired and a repetition time of five R-R intervals, the total acquisition time was approximately 7.5 minutes. With 75% k-space readout and an echo time of 48 msec, the net signal-to-noise ratio was 45:1 for the nonattenuated diffusion image (b {approx} 0 sec/mm2) and 25:1 for the attenuated (b {approx} 420 sec/mm2) images, which corresponded to a root-mean-square error in the direction of principal diffusion of approximately 8°.

Strains in Local Fiber and Cardiac Coordinates
Strain components of interest were computed by transforming the strain tensor Ssys from the lab frame |x, y, z| into local cardiac coordinates or local fiber coordinates by using

where R is a rotation matrix defined at each pixel for each coordinate system.

For each short-axis section, local fiber orientation f was defined as the first eigenvector of the diffusion tensor d1 with sign assigned to produce a counterclockwise orientation relative to the ventricular centroid

where {theta} is the gradient of the polar angle about the left ventricular centroid. In each section, a ventricular axis L was defined as the direction about which the fiber orientations have maximum symmetry, specifically, as the vector L that maximized the net axial contribution

where the sum extends over the set M of all left ventricular myocardial pixels, ri (L) is the radial direction normal to the axis L at the ith pixel, and fi is the fiber orientation at that pixel. Local cardiac coordinates were given by the unit vectors of cylindric coordinates coaxial with the axis L, at each point an orthonormal frame consisting of circumferential, longitudinal, and radial unit vectors: |c, l, r| = RTcard, respectively. Local fiber coordinates were an orthonormal frame that included fiber orientation, cross-fiber, and radial directions, |fxr^| = RTfiber,  where

and

The unit vector r^ is the best approximation of r that is orthogonal to local fiber orientation f, and local cross-fiber direction x is perpendicular to f and r^(Fig 2).



View larger version (26K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 2. Diagram shows geometric relations between local cardiac coordinates and local fiber coordinates. At each position in the middle of the ventricle, axes of local cardiac coordinates (open arrows) and local fiber coordinates (solid arrows) can be defined as follows. The longitudinal axis {ell} parallels the long axis of the left ventricle L. The radial axis r points away from the ventricular centroid and perpendicular to {ell}. The circumferential axis c is determined by the vector product {ell} x r. This diagram shows a voxel in the inner layer where fibers run obliquely from upper right to lower left. The direction of the fiber at this voxel is given so that the fiber axis f always orients counterclockwise about L. With f defined, r^ is the projection of r on the plane perpendicular to f. Since fibers lie almost tangential to the wall surface, the difference between r^ and r is negligible. Cross-fiber axis x is determined by the vector product r^ x f. The fiber helix angle h(f) is the angle between c and f with the sign of the angle consistent with the sign of the longitudinal component of f (Eq [9]). According to this convention, h(f) at this voxel is positive.

 
The helix angle h(v) of any vector v was defined as the angle to the circumferential direction of the projection of v to the circumferential-longitudinal plane, which was computed by using the complex form

with vc and vl the local circumferential and longitudinal components of v, respectively, and i = (-1)1/2. In particular, fiber helix angles were defined by h(f) (Fig 2).

Statistics
For each strain component in each subject, the mean and SD for the left ventricle in the imaged section were computed. The mean for each subject and the mean and SD for all subjects together were then computed for every averaged strain component to evaluate the intersubject variability. To study spatial uniformity of fiber shortening, we computed 95% variance ratio CIs of normal strains with respect to fiber shortening. In addition, linear regression of fiber shortening and of other normal strains against wall depth was performed. The resultant slopes and residuals indicating transmural variability and local variability, respectively, were compared among strain components.

Geometric relations between fibers and principal shortening strains were studied by computing the regression slopes, intercepts, and correlation coefficients between fiber helix angles and helix angles of principal shortening directions. The angles between fibers and directions of principal thickening also were analyzed in terms of the mean and SD.

In this study, MR data acquisition was performed by T.G.R. and W.Y.I.T. together, and data analysis was performed by both W.Y.I.T. and V.J.W.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 REFERENCES
 
MR Imaging of Myocardial Fiber Architecture and 3D Strain
Results of diffusion MR imaging of myocardial fiber architecture in healthy subjects were mutually consistent and in good agreement with known myocardial fiber architecture. Figure 3a is a graphic rendering of diffusion MR imaging of myofiber orientation in one subject. Fiber orientations all appeared perpendicular to the radial direction and showed smooth circular symmetry about the ventricular axis. Fiber helix angles showed a uniform transmural progression from approximately -90° epicardially through 0° in the middle of the wall to approximately +90° endocardially. The graphs of helix angles for all eight subjects (Fig 3b) had maxima at the circumferential orientation (0°) and decreased from this maximum smoothly and symmetrically (11).



View larger version (140K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 3a. (a) Map of myocardial fiber orientation of a midventricular short-axis section in a healthy subject; the anterior wall is at the top, and the interventricular septum is at the left. Fiber orientation at each pixel is represented by a small cylinder whose axis parallels the leading eigenvector of the diffusion tensor. Color indicates fiber helix angle according to the code inset at lower right. Fiber orientation shows a transmural progression from right-handed helicity in the subepicardial layers (turquoise) through circumferential orientation in the middle of the wall (blue) to left-handed helicity in the subendocardial layers (pink). (b) Graphs of fiber helix angles, one for each subject, consistently show predominant circumferential fiber population and a balanced distribution among populations of positive and negative helix angles.

 


View larger version (14K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 3b. (a) Map of myocardial fiber orientation of a midventricular short-axis section in a healthy subject; the anterior wall is at the top, and the interventricular septum is at the left. Fiber orientation at each pixel is represented by a small cylinder whose axis parallels the leading eigenvector of the diffusion tensor. Color indicates fiber helix angle according to the code inset at lower right. Fiber orientation shows a transmural progression from right-handed helicity in the subepicardial layers (turquoise) through circumferential orientation in the middle of the wall (blue) to left-handed helicity in the subendocardial layers (pink). (b) Graphs of fiber helix angles, one for each subject, consistently show predominant circumferential fiber population and a balanced distribution among populations of positive and negative helix angles.

 
Images of the myocardial 3D strain tensor fields showed accurate registration with the corresponding images of myocardial diffusion. They were mutually consistent and were in good agreement with results from previous studies (12). Figure 4a displays a 3 x 3 matrix of images of the components of Ssys in local cardiac coordinates |c, l, r|. The maps of normal components of strains showed shortening in circumferential and longitudinal directions and thickening in the radial direction—that is, a 2D shortening machine (13). Closer inspection showed that circumferential shortening had a larger transmural gradient, increasing toward the endocardium, than longitudinal shortening did. Radial strains were mostly positive and had a moderate transmural gradient, but they were negative in some subendocardial regions (adjacent to the anterior and posterior papillary muscle groups), which was consistent with systolic compression of the trabecular myocardium. In this study, we used a black-blood phase-contrast technique, so the phase signal came mainly from the myocardium. A partial volume effect from blood signal in the trabecular space is negligible.



View larger version (56K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 4a. Components of myocardial 3D strain in local cardiac coordinates registered with the fiber image in Figure 3a. SC, SL, and SR denote normal strains along circumferential, longitudinal, and radial directions, respectively; Scl, Scr, and Slr denote shear strains in circumferential-longitudinal, circumferential-radial, and longitudinal-radial planes, respectively. (a) Schematic shows each strain component displayed with a common gray scale in which dark gray indicates negative values, medium gray indicates zero, and light gray to white indicates positive values (see left side of image). The maximum root-mean-square error of strain measurement at each pixel estimated from MR noise is about 0.01 and is indicated by the width of a small rectangle in the gray-level bar. Strain component maps are arranged in a 3 x 3 matrix form that is consistent with the matrix form of a strain tensor. See text for a detailed description of individual strain patterns. (b) Graph shows strain components (mean ± SD) in eight healthy subjects. Different symbols indicate the means of individual subjects. Vertical bars indicate SD about the means, and • is the mean of the means among the subjects. The means of absolute shears were computed because of sign inconsistency in the same section. Reproducibility of this method is evident in the small variability of the mean strain components among subjects.

 


View larger version (13K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 4b. Components of myocardial 3D strain in local cardiac coordinates registered with the fiber image in Figure 3a. SC, SL, and SR denote normal strains along circumferential, longitudinal, and radial directions, respectively; Scl, Scr, and Slr denote shear strains in circumferential-longitudinal, circumferential-radial, and longitudinal-radial planes, respectively. (a) Schematic shows each strain component displayed with a common gray scale in which dark gray indicates negative values, medium gray indicates zero, and light gray to white indicates positive values (see left side of image). The maximum root-mean-square error of strain measurement at each pixel estimated from MR noise is about 0.01 and is indicated by the width of a small rectangle in the gray-level bar. Strain component maps are arranged in a 3 x 3 matrix form that is consistent with the matrix form of a strain tensor. See text for a detailed description of individual strain patterns. (b) Graph shows strain components (mean ± SD) in eight healthy subjects. Different symbols indicate the means of individual subjects. Vertical bars indicate SD about the means, and • is the mean of the means among the subjects. The means of absolute shears were computed because of sign inconsistency in the same section. Reproducibility of this method is evident in the small variability of the mean strain components among subjects.

 
Like normal strains, shear strains were qualitatively and quantitatively consistent among subjects. Circumferential-longitudinal shear, called "epicardial shear," was positive and uniform throughout the ventricle of the same section and corresponded to myocardial twist, which is the rotation of the ventricular apex about the ventricular axis relative to the base (14). The other two shears, called "transverse shears," were larger than epicardial shear and had more complicated spatial patterns. Circumferential-radial shear presented variable sign changes across the wall, whereas the longitudinal-radial shear showed a reversal of sign between the septum and the free wall, as reported in animal models (15,16).

The means and SDs of cardiac strain components within and among subjects are plotted in Figure 4b. Despite the large variance within each strain component in each subject, the intersubject variance of the component means was small.

In Figure 5, the strain tensor field Ssys is graphically rendered with boxes so that 3D patterns of the strain tensor field orientation are appreciated easily. We see that the direction of principal thickening was predominantly radial. The lengths of two principal shortenings were unequal, and the skew of the greater shortening corresponded to myocardial twist. Furthermore, the shortenings were seen to rotate about the thickening direction from epi- to endocardium, with greatest shortening longitudinal more epicardially and circumferential more endocardially.



View larger version (141K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 5. Graphic rendering of myocardial 3D strain tensor data of Figure 4. The strain tensor at each pixel is represented by a box that has axes parallel to the directions of the principal strains and that has lengths proportional to the stretch in each direction. The longest axis (perpendicular to the pink surface of the box) of each box indicates the orientation of the principal thickening, the shortest axis (perpendicular to the green surface of the box) indicates the maximum shortening direction, and the third axis (perpendicular to the blue surface of the box) indicates the minimum shortening direction. See text for characteristic orientations of the two principal strains.

 
MR Imaging of Myocardial Fiber Shortening and Structure-Function Relations
Images of fiber shortening and related strain components were computed by using Equation (4) to transform the systolic strain tensor Ssys into local fiber coordinates |f, x, r^|. The matrix of images of the systolic strain tensor components in local fiber coordinates is illustrated in Figure 6a. The image of fiber shortening SF is at the upper left. This image qualitatively shows fiber shortening was uniform throughout the middle level of the left ventricle and shows no evidence of regions of exceptional or anomalous shortening. Cross-fiber shortening, the middle image in the middle column in Figure 6a, showed a steep and uniform transmural gradient that was small at the epicardium but gradually exceeded fiber shortening in the inner half of the ventricular wall to attain a high maximum in the area beneath the endocardium. Myocardial fiber shears, represented by the off-diagonal images located at the upper right and lower left triangles of the 3 x 3 matrix, showed a preponderance of shear transverse to the fiber in the cross-fiber–radial plane; this component also showed a sign reversal between the free wall and the septum. Shears in fiber–cross-fiber and fiber-radial planes were much smaller.



View larger version (58K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 6a. Myocardial fiber shortening and other related strain components in fiber coordinates determined from registered images of fibers (Fig 3) and strains (Fig 4). The format of the illustration is similar to that of Figure 4. SF, SR^, and SX, denote normal strains along fiber, radial, and cross-fiber directions, respectively; Sfr{vee}, Sfx, and Sxr^ denote shear strains in fiber-radial, fiber-cross-fiber, and cross-fiber-radial planes, respectively. (a) Schematic shows that the maximum root-mean-square error of strain measurement at each pixel estimated from MR noise is about 0.02 and is indicated by the width of a small rectangle in the gray-level bar. See text for a detailed description of individual strain patterns. (b) Graph shows summary plot of strain components in the fiber coordinates with the same symbols used in Figure 4b. The means among all subjects are tightly clustered. Note that in both intersubject and intrasubject variabilities, fiber shortening is the smallest among all strains.

 


View larger version (13K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 6b. Myocardial fiber shortening and other related strain components in fiber coordinates determined from registered images of fibers (Fig 3) and strains (Fig 4). The format of the illustration is similar to that of Figure 4. SF, SR^, and SX, denote normal strains along fiber, radial, and cross-fiber directions, respectively; Sfr{vee}, Sfx, and Sxr^ denote shear strains in fiber-radial, fiber-cross-fiber, and cross-fiber-radial planes, respectively. (a) Schematic shows that the maximum root-mean-square error of strain measurement at each pixel estimated from MR noise is about 0.02 and is indicated by the width of a small rectangle in the gray-level bar. See text for a detailed description of individual strain patterns. (b) Graph shows summary plot of strain components in the fiber coordinates with the same symbols used in Figure 4b. The means among all subjects are tightly clustered. Note that in both intersubject and intrasubject variabilities, fiber shortening is the smallest among all strains.

 
Mean fiber strains were quantitatively highly consistent among subjects (Fig 6b). Among all subjects, mean systolic fiber shortening was 0.12 ± 0.01, and mean cross-fiber shortening was 0.11 ± 0.03. Fiber shortening showed a substantially narrower distribution of values than any other normal strain, as shown in Figure 6b. The widths of strain components may be represented by the variances of the mean-normalized distributions. As shown in Table 1, the variances of other normal strains were always larger than that of fiber shortening, with the ratios ranging from 1.6 to 4.3.


View this table:
[in this window]
[in a new window]

 
TABLE 1. Variance Ratio CI of Normal Strains with Respect to Fiber Shortening
 
Variation of strain components in the transmural direction also was assessed by means of linear regression. An example of these regressions for one subject is shown in Figure 7. In all subjects, the transmural slope of fiber shortening was consistently the smallest of the measured components (Table 2). Even though normalized transmural slopes of SR and SR^ were comparable to that of SF, the regression residual of SF was smaller than that of SR and SR^ (Table 2).



View larger version (47K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 7. Graphs show linear regressions of normal strain components versus wall depth, with 95% CIs. The points displayed here are sequentially subsampled from the total sorted pixels. Negative strains, SC, SF, SL, SX, S1, and S2, are made positive and shown with the common positive gray scale indicated at the lower right. Despite additional noise from the diffusion measurement, fiber shortening still shows a small transmural gradient and small regression residual. For quantitative comparison with other normal strains, see Table 2.

 

View this table:
[in this window]
[in a new window]

 
TABLE 2. Slopes and Residuals of Linear Regression of Normal Strains versus Wall Depth
 
The geometric relation of fiber orientation to 3D strain was studied by measuring fiber orientations relative to the three principal strain directions, namely maximum shortening s1, minimum shortening s2, and principal thickening s3. Fibers were approximately perpendicular to the direction of principal thickening, as illustrated in the image at the top of Figure 8. Among all subjects, the mean fiber angle to principal thickening was 93° ± 3, with SD in each subject averaging 11° ± 2 (Table 3).



View larger version (61K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure 8. Top: Images show the angle of fibers to principal thickening color-coded according to the scale at the right. Fibers and principal thickening are approximately orthogonal except those in the trabecular endocardium. Bottom: Color-coded images show the helix angles of myocardial fiber orientation h(f), maximum principal shortening h(s1), and minimum principal shortening h(s2). These images show smooth variation of fiber angles from the direction of maximum shortening subepicardially to that of minimum shortening subendocardially.

 

View this table:
[in this window]
[in a new window]

 
TABLE 3. Relations between Fiber and Principal Strains
 
The bottom of Figure 8 shows images of fiber helix angles and the helix angles of maximum shortening and minimum shortening. We see that fiber orientation rotated by 90° relative to the shortening strains from epi- to endocardium: Fibers were aligned with the direction of maximum shortening s1 in the epicardium and with that of minimum shortening s2 at the endocardium. Figure 8 also shows that this transmural variation of fiber-to-shortening angle is characteristic of the entire ventricular myocardium, except the papillary muscle insertions. By computing the linear regressions of fiber helix angles h(f) versus helix angles of principal shortening h(s1), a mean correlation coefficient r = 0.58 ± 0.11 was found over all subjects (Table 3).


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 REFERENCES
 
Fiber Orientation
Results of the present MR imaging study of myocardial fiber architecture in vivo are in good agreement with those of previous studies in which conventional histologic examination was performed (11). Classic features of fiber orientation demonstrated by means of MR imaging include circular symmetry, the orthogonality of the fibers to the radial direction, the smooth transmural progression of fiber helix angles from approximately -90° epicardially to 0° (circumferential) at the middle of the wall to approximately +90° endocardially, with a distribution of fiber angles that is symmetric about its maximum, which is near 0°.

Advantages of MR imaging are illustrated in the trabecular myocardium. Histologic studies of myocardial fiber architecture were initially equivocal about the trabecular zones but ultimately showed trabecular regions to contain fiber angles in a range of +70° to +90° and also showed that this contribution substantially improved the overall symmetry of the ventricular fiber angle graph (11). MR imaging results, conversely, indicated that trabecular and compact myocardium were equally accessible. MR imaging immediately obtains a full and symmetric fiber angle graph. Furthermore, images of fiber orientation show a uniform progression of fiber angles across the interface between trabecular and compact myocardium, which is direct evidence of the structural unity of compact and noncompact myocardium.

Fiber Shortening
In the present study, mean fiber shortening in healthy human subjects was 0.12 ± 0.01. This is smaller than the shortening of 0.15 reported by MacGowan et al (17). In that study, fiber shortening was computed by using strain measurements obtained by means of MR imaging tagging in vivo in conjunction with fiber orientations determined in unrelated autopsy specimens. Although the variability of fiber angles among autopsy specimens was small, it is still unclear to what extent these angles faithfully represent fiber geometry in situ. Their endocardial fiber angles were approximately +70°, which is substantially less than the +90° of our in vivo results. As shown in Figure 8, in vivo imaging shows endocardial fiber angles up to +90°, at which point they are orthogonal to the direction of maximum shortening. Therefore, calculation of fiber shortening by using fiber angles less than +90° would lead to overestimation of shortening by an amount cos(20°)-2 {approx} 1.13, or about half of the noted disparity.

The hypothesis that myocardial fiber shortening is transmurally uniform was proposed by Arts et al (18) and subsequently was derived as an optimum solution to a general model of myocardial fiber orientation and function (19). This hypothesis has been supported by experimental observations in the dog, including those of Waldman et al (2), who found fiber shortenings in the inner and outer walls of 0.06 ± 0.06 and 0.09 ± 0.04, respectively, by using radiopaque marker strain measurement and postmortem histologic examination for fiber architecture. Rademakers et al (20), using MR imaging tagging and postmortem histologic examination, found shortening in the endocardium and epicardium of 0.09 ± 0.10 and 0.06 ± 0.11, respectively. Bloomgarden et al (21), using MR imaging tagging and published histologic examination results, found fiber shortenings of 0.12 and 0.13, respectively, in the inner and outer halves of the ventricular wall.

To our knowledge, the present quantitative analysis of the uniformity of fiber shortening and its relation to other strain components has not been reported previously. Its results show that variance ratios of strain components are all substantially greater than 1, as compared with those of fiber shortening. Results of analysis of transmural profiles further showed that fiber shortening had the narrowest distribution and smallest transmural gradient among all strain components (Tables 1, 2), which verifies that fiber shortening is the most uniform of all cardiac strain components.

Cross-Fiber Shortening
Cross-fiber shortening in normal human hearts was 0.10 ± 0.03 and had the largest variance ratio with respect to fiber shortening (Table 1). The reason for the large variance in cross-fiber shortening was mainly because it had a steep transmural gradient, with 0.05 ± 0.02 at the epicardium and 0.20 ± 0.02 at the endocardium (Fig 7). This prominent transmural gradient was found previously in canine and human studies. In canine left ventricles, cross-fiber shortening in the outer and the inner walls, respectively, were reported to be 0.04 ± 0.04 and 0.17 ± 0.03 by Waldman et al (2), 0.006 ± 0.082 and 0.251 ± 0.103 by Rademakers et al (20), and 0.09 and 0.22 by Bloomgarden et al (21). In healthy human subjects, MacGowan et al (17) found a similar transmural trend: 0.08 ± 0.01 and 0.26 ± 0.01 in the areas beneath the epicardium and beneath the endocardium, respectively.

Advantages of Measuring Myocardial Strains in Local Fiber Coordinates
Having characterized transmural patterns of fiber and cross-fiber shortening, advantages of transforming the strain tensor from local cardiac coordinates to local fiber coordinates becomes evident. Whereas circumferential and longitudinal strains both showed transmural variations, fiber and cross-fiber strains took the opposite modes of shortening patterns across the wall: Fiber shortening was virtually constant, whereas cross-fiber shortening showed a steep gradient (Table 2). With these two drastically different transmural patterns as norms, it is easy to detect any deviation from these two standard patterns in diseased hearts.

Another advantage of transforming the strain tensor into local fiber coordinates is that two different myocardial functions can be separated. Fiber shortening indicates myocardial contractility, whereas cross-fiber shortening is believed to reflect a dynamic rearrangement of the myocardial sheet structure to facilitate systolic wall thickening (15). Different cardiac disease states may have different effects on these two functions (22); these changes would not be shown clearly in local cardiac coordinates.

The present MR imaging method demonstrated its capability to provide images of fiber orientation and myocardial strains with identical cardiac configurations in identical physiologic conditions; it directly mapped the fiber shortening and its relation to local wall deformation with anatomic details not previously available. This imaging capability should be of particular value in the study of disease in which the characterization of regional abnormalities and their change over time in the living individual is essential (23). Furthermore, for disease detection and characterization, conventional cardiac coordinates and corresponding strain components are likely to become equivocal, since anatomic distortions affect these geometry-based coordinates. Fiber shortening and related functional parameters, however, should be less affected by these distortions, and they thus better reflect disease location and severity.

Fiber and Principal Shortening
Since clarification by Streeter (11) of the basic helical geometry of myocardial fiber orientations, several explanations of its functional significance have been offered (18,20). Arts et al (1) emphasized the capacity of helical fiber orientations to make fiber stress and fiber strain transmurally more uniform and constructed a model in which observed patterns of structure and function emerged purely through local feedback. Teleologically, this uniformity may facilitate physiologic tuning by means of concomitant uniformity of cellular workload, of local tissue perfusion, and of myocyte metabolism and morphology (2426). Waldman et al (2) and Rademakers et al (20) suggested that the helical fiber geometry serves to transmit fiber contraction from the outer to the inner wall through the connective tissue network and produces the remarkable cross-fiber shortening in the inner wall that in turn generates substantial radial thickening. Present MR imaging findings in humans support both of these views, as well as their mutual synergy.

As was noted earlier, MR imaging demonstrated the change with transmural depth of fiber orientation from the direction of maximum to that of minimum shortening. This arrangement allowed principal shortening to increase with wall depth while the fibers, by rotating away from the direction of greatest shortening, continued to shorten by constant amounts. MR imaging showed that this rotation was linear with transmural depth (Fig 8, bottom) and that it exhausted the meaningful angular ranges of fiber orientation and of its relation to shortening, with fiber orientation rotating anatomically 180° relative to the heart wall even as it rotated 90° relative to the direction of maximum shortening. We see therefore that the observed fiber geometry in healthy humans accomplishes an equipartition of fibers across the entire range of relevant states. This condition of pushing to possible limits lends new support to the hypothesis that fiber geometry and its relation to myocardial 3D strain are objects of close physiologic optimization (27).

Technical Considerations
Multisection ventricular coverage from apex to base is constrained currently by the inherently low signal-to-noise ratio of present cardiac diffusion MR imaging methodology and consequent long imaging times. This study presents myocardial structure-function relations only at the midventricular level, and information concerning more apical or basal regions is lacking. In the future, imaging speed may be increased by replacing synchronized breathing with one of the automated MR imaging methods now being developed for respiratory artifact suppression (28). Acquisition economies also may be obtainable by use of intersecting sections in the cardiac long axis.

Present methodology for identification of optimum cardiac phase delay to acquire accurate diffusion data incurs certain systematic errors, as described previously (3). Among these are potential inaccuracies in the assumed linearity between myocardial strain and the observed diffusion, and sweet-spot localization errors arising from R-R variability or errors of myocardial tracking in the present sampling strategy. Within the range of normal cardiac function, we estimated that these systematic errors led to diffusion error of 15% or less, or uncertainty in fiber direction of 5° or less.

In the present study, acquisition of accurate diffusion images with minimum strain effects depended on cardiac synchrony. Application of this method to disease would be expected to be less accurate owing to mechanical dyssynergy. Partial relief from this problem may result from the fact that dyssynergic walls are most often hypokinetic and so would have reduced strains and correspondingly reduced measurement errors. Furthermore, new MR imaging pulse sequences that provide cardiac diffusion MR imaging completely free of myocardial strain effects have been proposed (29). Methods such as this would provide accurate maps of myocardial fiber architecture in the normal or the diseased heart and at any point in the cardiac cycle.

MR imaging of diffusion and strain makes possible quantitative noninvasive imaging of myocardial fiber shortening and its relation to local myocardial deformation in healthy human subjects. This method addresses the problems of in vivo access and of metric distortion that hindered prior attempts to relate myocardial structure and function. Present data confirm key quantitative features of fiber shortening previously established invasively in animal models and indirectly in humans; but as image data, MR imaging resolves these measurements with anatomic details not previously available. Recognition of robust regularities in myocardial fiber architecture and function with the present method may provide new insight into myocardial design and furnish new normative parameters to facilitate evaluation of myocardial dysfunction.


    APPENDIX
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 REFERENCES
 
Determination of Optimum Cardiac Phase Delay for Diffusion MR Imaging
We have shown previously that cardiac diffusion MR imaging can be performed by using diffusion encoding spanning a cardiac cycle (30), that the resultant data are affected by myocardial strains (4,31), and that this strain effect disappears in the normal heart if the diffusion encoding is performed at particular phase delays in the cardiac cycle (3). Because the effect of strain on diffusion is approximately linear in strain, these phase delays are the points in the cardiac cycle when the myocardial strain state equals its time average during the cardiac cycle; because of the synchrony of the normal heart, such times coincide for all strain components and all locations and are called "sweet spots" (32). Observed diffusion Dobs is related to true diffusion D measured without motion: Dobs = <U-1 D U-1> {approx} D - <S>D - D<S>, where angle brackets denote the time-average during the diffusion encoding time of one cardiac cycle, U is the time-dependent material stretch tensor, and S is the corresponding strain tensor.

To have a more intuitive picture of the sweet spot, consider a measurement of diffusion in a rubber band that we deform cyclically. Diffusion encoding and decoding gradient pulses are applied at the same phase of the deformation cycle to eliminate signal dropout owing to bulk motion. If these pulses are applied when the band is shortest, then the subsequent stretching of the band also will stretch the phase modulation inscribed on it and lead to a reduction of the phase modulation. As this phase modulation is reduced during the diffusion time, the diffusion-induced signal attenuation will be less than it would have been had there been no stretching, and measured diffusion would appear smaller than expected. Conversely, if diffusion encoding is applied when the band is maximally stretched, the apparent diffusion will be larger than expected. In general, the deviation of the apparent diffusion from the true diffusion is proportional to the time integral of the strain curve during a full cycle (3). This deviation varies periodically from positive to negative, and the sweet spots are the time points when the deviation crosses zero. Diffusion measurement at the sweet spots is therefore free of strain effects.

In each subject, sweet-spot location was determined by analyzing the mean radial strain, a robust and easily measured scalar index of strain state, as described (3). From the strain-rate movie, we measured the time-course of the mean radial strain rate S'rr({tau}), where the radial was defined by the ventricular centroid and the underscore denotes the mean over ventricular pixels. This was integrated numerically to construct a scalar function F(t) that varies linearly with the mean radial strain

The requisite times {Psi} when the myocardium occupied its mean strain state were then the solutions of

where the angle brackets denote the time-average in {tau} during one cardiac cycle. A typical curve F(t) and its relation to {Psi} is shown in Figure A1. The normal cardiac cycle contains two time points that satisfy Equation (A2), one in midejection and the other in rapid filling. In this work, we always used the {Psi} that occurs in midsystole to match the cardiac configuration of strain images acquired in systole.



View larger version (19K):
[in this window]
[in a new window]
[Download PPT slide]
 
Figure A1. Graph shows time course of myocardial mean radial strain F(t) computed from an MR imaging movie of myocardial 2D strain rate. Two sweet-spot time points {Psi} are defined by the condition that F({Psi}) equals the time-average of F(t) during the cardiac cycle. In this subject, these points are located at 206 and 547 msec after the R wave. The reference line is set to let F({Psi}) = 0 so that the integral of F(t) during a cardiac cycle equals zero.

 
Computation of the Myocardial Net Systolic Strain
The image of the myocardial systolic strain tensor was constructed by means of an indirect integration of strain rates by using synchrony and the function F(t) of Equation (A1). An image of the myocardial systolic strain tensor S{Psi} is equal to the integral of the strain rates

where {eta} refers to the time of end systole. Synchrony of strain rates implies S'(t1) L(t2) = S'(t2) L(t1), for L(t) any linear functional of strain rate and any times t1 and t2. Setting L(t) = S'rr({tau}) and t1 = {Psi}, one obtains

which gives, according to Equations (A1) and (A3),

where the image S{Psi} expresses net systolic strains relative to cardiac configuration at t = {Psi}.

According to standard practice, we computed net systolic strain relative to cardiac configuration at end-diastole Ssys by change of coordinates Ssys = U S{Psi} U for

where U is the stretch from t = {Psi} to end-diastole in image coordinates.


    FOOTNOTES
 
Abbreviations: 2D = two-dimensional, 3D = three-dimensional

Author contributions: Guarantors of integrity of entire study, W.Y.I.T., V.J.W., R.M.W.; study concepts and design, W.Y.I.T., V.J.W., T.J.B.; definition of intellectual content, W.Y.I.T., V.J.W.; literature research, W.Y.I.T.; clinical studies, W.Y.I.T., T.G.R., V.J.W.; experimental studies, W.Y.I.T., T.G.R., V.J.W.; data acquisition, T.G.R., W.Y.I.T.; data analysis, W.Y.I.T., V.J.W.; statistical analysis, W.Y.I.T.; manuscript preparation, W.Y.I.T.; manuscript editing, V.J.W., R.M.W., T.J.B.; manuscript review, V.J.W., T.G.R., R.M.W., T.J.B.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 REFERENCES
 

  1. Arts T, Prinzen FW, Reneman RS. Mechanics of the wall of the left ventricle. In: Strackee J, Westerhof N, eds. The physics of heart and circulation. London, England: Institute of Physics, 1993; 153-173.
  2. Waldman LK, Nosan D, Villarreal F, Covell JK. Relation between transmural deformation and local fiber direction in canine left ventricle. Circ Res 1988; 63:550-562.[Abstract/Free Full Text]
  3. Tseng WI, Reese TG, Weisskoff RM, Wedeen VJ. Cardiac diffusion tensor MRI in vivo without strain correction. Magn Reson Med 1999; 42:393-403.[Medline]
  4. Reese TG, Weisskoff RM, Smith RN, Rosen BR, Dinsmore RE, Wedeen VJ. Imaging myocardial fiber architecture in vivo with magnetic resonance. Magn Reson Med 1995; 34:786-791.[Medline]
  5. Hsu EW, Muzikant AL, Matulevicius SA, Penland RC, Henriquez CS. Magnetic resonance myocardial fiber-orientation mapping with direct histological correlation. Am J Physiol 1998; 274:H1627-1634.
  6. Scollan DF, Holmes A, Winslow R, Forder J. Histological validation of myocardial microstructure obtained from diffusion tensor magnetic resonance imaging. Am J Physiol 1998; 275:H2308-2318.
  7. Wedeen V. Magnetic resonance imaging of myocardial kinematics: technique to detect, localize and quantify the strain rates of the active human myocardium. Magn Reson Med 1992; 27:52-67.[Medline]
  8. Pelc NJ, Bernstein MA, Shimakawa A, Glover GH. Encoding strategies for three-direction phase-contrast MR imaging of flow. J Magn Reson Imaging 1991; 1:405-413.[Medline]
  9. Tanner JE. Use of the stimulated echo in NMR diffusion studies. Chem Phys 1969; 52:2523-2526.
  10. Basser PJ, Pierpaoli C. A simplified method to measure the diffusion tensor from seven MR images. Magn Reson Med 1998; 39:928-934.[Medline]
  11. Streeter DD. Gross morphology and fiber geometry of the heart. In: Berne RM, eds. Handbook of physiology. Bethesda, Md: American Physiological Society, 1979; 61-112.
  12. Young AA, Kramer CM, Ferrari VA, Axel L, Reichek N. Three-dimensional left ventricular deformation in hypertrophic cardiomyopathy. Circulation 1994; 90:854-867.[Abstract/Free Full Text]
  13. Waldman LK. Multidimensional measurement of regional strains in the intact heart. In: Glass L, Hunter P, McCulloch A, eds. Theory of heart. New York, NY: Springer-Verlag, 1991; 145-174.
  14. Buchalter MB, Weiss JL, Rogers WJ, Zerhouni EA. Noninvasive quantification of left ventricular rotational deformation in normal humans using magnetic resonance imaging myocardial tagging. Circulation 1990; 81:1236-1244.[Abstract/Free Full Text]
  15. LeGrice IJ, Takayama Y, Covell JW. Transverse shear along myocardial cleavage planes provides a mechanism for normal systolic wall thickening. Circ Res 1995; 77:182-193.[Abstract/Free Full Text]
  16. Costa KD. The structural basis of three-dimensional ventricular mechanics San Diego, Calif: University of California, 1996; Dissertation..
  17. MacGowan GA, Shapiro EP, Azhari H, et al. Noninvasive measurement of shortening in the fiber and cross-fiber directions in the normal human left ventricle and in idiopathic dilated cardiomyopathy. Circulation 1997; 96:535-541.[Abstract/Free Full Text]
  18. Arts T, Reneman RS, Veenstra PC. A model of the mechanics of the left ventricle. Ann Biomed Eng 1979; 7:299-318.[Medline]
  19. Rijcken J, Arts T, Bovendeerd P, Schoofs B, van Campen D. Optimization of left ventricular fibre orientation of the normal heart for homogeneous sarcomere length during ejection. Eur J Morph 1996; 34:39-46.
  20. Rademakers FE, Rogers WJ, Guier WH, et al. Relation of regional cross-fiber shortening to wall thickening in the intact heart: three-dimensional strain analysis by NMR tagging. Circulation 1994; 89:1174-1182.[Abstract/Free Full Text]
  21. Bloomgarden DC, Young AA, Kraitchman DL, et al. Noninvasive assessment of regional myocardial fiber and cross-fiber strain using finite element modeling and magnetic resonance tagging (abstr) In: Proceedings of the Fifth Meeting of the International Society for Magnetic Resonance in Medicine. Berkeley, Calif: International Society for Magnetic Resonance in Medicine, 1997.
  22. Wedeen VJ, Tseng WY, Reese TG, Kantor HA, Weisskoff RM. MRI finds normal myocardial fiber shortening and reduced cross-fiber shortening in hypertrophic cardiomyopathy; Presented at the 71st Scientific Session of the American Heart Association, Dallas, Texas, November 21–25, 1998..
  23. Bogaert J, Maes A, Van de Werf F, et al. Functional recovery of subepicardial myocardial tissue in transmural myocardial infarction after successful reperfusion: an important contribution to the improvement of regional and global left ventricular function. Circulation 1999; 99:36-43.[Abstract/Free Full Text]
  24. Feigl EO. Coronary physiology. Physiol Rev 1983; 63:1-205.[Abstract/Free Full Text]
  25. Prinzen FW. Gradients in fiber shortening and metabolism across the ischemic left ventricular wall. Am J Physiol 1986; 250:H255-264.
  26. Meggs LG, Coupet J, Huang H, et al. Regulation of angiotensin II receptors on ventricular myocytes after myocardial infarction in rats. Circ Res 1993; 72:1149-1162.[Abstract/Free Full Text]
  27. Arts T, Prinzen FW, Snoeckx LH, Rijcken JM, Reneman RS. Adaptation of cardiac structure by mechanical feedback in the environment of the cell: a model study. Biophys J 1994; 66:953-961.[Medline]
  28. Reese TG, Poncelet BP, Weisskoff RM, Brady TJ. Improved accuracy in EPI measurement of myocardial perfusion using motion-tracking images (abstr) In: Proceedings of the Sixth Meeting of the International Society for Magnetic Resonance in Medicine. Berkeley, Calif: International Society for Magnetic Resonance in Medicine, 1998.
  29. Wedeen VJ, Reese TG, Tseng WI, Weisskoff RM. Strain-free cardiac diffusion MRI (abstr) In: Proceedings of the Seventh Meeting of the International Society for Magnetic Resonance in Medicine. Berkeley, Calif: International Society for Magnetic Resonance in Medicine, 1999.
  30. Edelman RR, Gaa J, Wedeen VJ, et al. In-vivo measurement of water diffusion in the human heart. Magn Reson Med 1994; 32:423-428.[Medline]
  31. Reese TG, Wedeen VJ, Weisskoff RM. Measuring diffusion in the presence of material strain. J Magn Reson B 1996; 112:253-258.[Medline]
  32. Wedeen VJ, Weisskoff RM, Reese TG, et al. Motionless movies of myocardial strain rate using stimulated echoes. Magn Reson Med 1995; 33:401-408.[Medline]



This article has been cited by other articles:


Home page
Exp PhysiolHome page
S. G. Campbell, E. Howard, J. Aguado-Sierra, B. A. Coppola, J. H. Omens, L. J. Mulligan, A. D. McCulloch, and R. C. P. Kerckhoffs
Effect of transmurally heterogeneous myocyte excitation-contraction coupling on canine left ventricular electromechanics
Exp Physiol, May 1, 2009; 94(5): 541 - 552.
[Abstract] [Full Text] [PDF]


Home page
Circ Cardiovasc ImagingHome page
M.-T. Wu, M.-Y. M. Su, Y.-L. Huang, K.-R. Chiou, P. Yang, H.-B. Pan, T. G. Reese, V. J. Wedeen, and W.-Y. I. Tseng
Sequential Changes of Myocardial Microstructure in Patients Postmyocardial Infarction by Diffusion-Tensor Cardiac MR: Correlation With Left Ventricular Structure and Function
Circ Cardiovasc Imaging, January 1, 2009; 2(1): 32 - 40.
[Abstract] [Full Text] [PDF]


Home page
Phil Trans R Soc AHome page
S. G Campbell, S. N Flaim, C. H Leem, and A. D McCulloch
Mechanisms of transmurally varying myocyte electromechanics in an integrated computational model
Phil Trans R Soc A, September 28, 2008; 366(1879): 3361 - 3380.
[Abstract] [Full Text] [PDF]


Home page
RadiologyHome page
H. Wen, K. A. Marsolo, E. E. Bennett, K. S. Kutten, R. P. Lewis, D. B. Lipps, N. D. Epstein, J. F. Plehn, and P. Croisille
Adaptive Postprocessing Techniques for Myocardial Tissue Tracking with Displacement-encoded MR Imaging
Radiology, January 1, 2008; 246(1): 229 - 240.
[Abstract] [Full Text] [PDF]


Home page
Eur. J. Cardiothorac. Surg.Home page
S. H. Gilbert, A. P. Benson, P. Li, and A. V. Holden
Regional localisation of left ventricular sheet structure: integration with current models of cardiac fibre, sheet and band structure
Eur. J. Cardiothorac. Surg., August 1, 2007; 32(2): 231 - 249.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
S. M. Felton, T. A. Gaige, T. G. Reese, V. J. Wedeen, and R. J. Gilbert
Mechanical basis for lingual deformation during the propulsive phase of swallowing as determined by phase-contrast magnetic resonance imaging
J Appl Physiol, July 1, 2007; 103(1): 255 - 265.
[Abstract] [Full Text] [PDF]


Home page
CirculationHome page
M.-T. Wu, W.-Y. I. Tseng, M.-Y. M. Su, C.-P. Liu, K.-R. Chiou, V. J. Wedeen, T. G. Reese, and C.-F. Yang
Diffusion Tensor Magnetic Resonance Imaging Mapping the Fiber Architecture Remodeling in Human Myocardium After Infarction: Correlation With Viability and Wall Motion
Circulation, September 5, 2006; 114(10): 1036 - 1045.
[Abstract] [Full Text] [PDF]


Home page
Eur. J. Cardiothorac. Surg.Home page
F. Carreras, M. Ballester, S. Pujadas, R. Leta, and G. Pons-Llado
Morphological and functional evidences of the helical heart from non-invasive cardiac imaging
Eur. J. Cardiothorac. Surg., April 1, 2006; 29(Suppl_1): S50 - S55.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
P. P. Sengupta, B. K. Khandheria, J. Korinek, J. Wang, and M. Belohlavek
Biphasic tissue Doppler waveforms during isovolumic phases are associated with asynchronous deformation of subendocardial and subepicardial layers
J Appl Physiol, September 1, 2005; 99(3): 1104 - 1111.
[Abstract] [Full Text] [PDF]


Home page
HeartHome page
D J Sahn and G W Vick III
Review of new techniques in echocardiography and magnetic resonance imaging as applied to patients with congenital heart disease
Heart, December 1, 2001; 86(90002): ii41 - 53.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Tseng, W.-Y. I.
Right arrow Articles by Wedeen, V. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Tseng, W.-Y. I.
Right arrow Articles by Wedeen, V. J.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
RADIOLOGY RADIOGRAPHICS RSNA JOURNALS ONLINE