Toward high-resolution myocardial tagging

Myocardial tagging by cardiovascular magnetic resonance: evolution of techniques--pulse sequences, analysis algorithms, and applications

Cardiovascular magnetic resonance (CMR) tagging has been established as an essential technique for measuring regional myocardial function. It allows quantification of local intramyocardial motion measures, e.g. strain and strain rate. The invention of CMR tagging came in the late eighties, where the technique allowed for the first time for visualizing transmural myocardial movement without having to implant physical markers. This new idea opened the door for a series of developments and improvements that continue up to the present time. Different tagging techniques are currently available that are more extensive, improved, and sophisticated than they were twenty years ago. Each of these techniques has different versions for improved resolution, signal-to-noise ratio (SNR), scan time, anatomical coverage, three-dimensional capability, and image quality. The tagging techniques covered in this article can be broadly divided into two main categories: 1) Basic techniques, which include magnetization saturation, spatial modulation of magnetization (SPAMM), delay alternating with nutations for tailored excitation (DANTE), and complementary SPAMM (CSPAMM); and 2) Advanced techniques, which include harmonic phase (HARP), displacement encoding with stimulated echoes (DENSE), and strain encoding (SENC). Although most of these techniques were developed by separate groups and evolved from different backgrounds, they are in fact closely related to each other, and they can be interpreted from more than one perspective. Some of these techniques even followed parallel paths of developments, as illustrated in the article. As each technique has its own advantages, some efforts have been made to combine different techniques together for improved image quality or composite information acquisition. In this review, different developments in pulse sequences and related image processing techniques are described along with the necessities that led to their invention, which makes this article easy to read and the covered techniques easy to follow. Major studies that applied CMR tagging for studying myocardial mechanics are also summarized. Finally, the current article includes a plethora of ideas and techniques with over 300 references that motivate the reader to think about the future of CMR tagging.

Total removal of unwanted harmonic peaks (TruHARP) MRI for single breath-hold high-resolution myocardial motion and strain quantification

Current MRI methods for myocardial motion and strain quantification have limited resolution because of Fourier space spectral peak interference. Methods have been proposed to remove this interference in order to improve resolution; however, these methods are clinically impractical due to the prolonged imaging times. In this paper, we propose total removal of unwanted harmonic peaks (TruHARP); a myocardial motion and strain quantification methodology that uses a novel single breath-hold MR image acquisition protocol. In post-processing, TruHARP separates the spectral peaks in the acquired images, enabling high-resolution motion and strain quantification. The impact of high resolution on calculated circumferential and radial strains is studied using realistic Monte Carlo simulations, and the improvement in strain maps is demonstrated in six human subjects.

Transmural myocardial strain in mouse: quantification of high-resolution MR tagging using harmonic phase (HARP) analysis

MR tagging allows noninvasive examination of regional myocardial function with high accuracy and reproducibility. The current tagging method is limited by low tagging resolution for accurate transmural strain quantification. Previously, a spatial modulation of magnetization (SPAMM)-based method was proposed to increase the tagging resolution by combining two or more tagged images with different tagging grid positions. However, there has been limited application due to the challenge in image processing of multiple data sets. In the current study, we propose a harmonic phase (HARP)-based method for automated and fast analysis of high tagging resolution images. First-order harmonic peaks from low tagging resolution images were combined to generate the composite second-order harmonic peak for strain computation. The combined images reached a tagging resolution of 0.3 mm. The proposed method was applied to the quantification of transmural myocardial wall strain in seven normal C57BL/6 mice. Principal strains, as well as radial and circumferential strains, were quantified using the current method.

TRUHARP: SINGLE BREATH-HOLD MRI FOR HIGH RESOLUTION CARDIAC MOTION AND STRAIN QUANTIFICATION

MRI techniques for tissue motion and strain quantifications have limited resolution because of interference from the conjugate echo or spectral peak in Fourier space. Methods have been proposed to remove this interference in order to improve resolution; however, these methods are clinically impractical due to long image acquisition time. In this paper, we propose TruHARP, an MRI motion and strain quantification methodology that involves a novel single breath-hold imaging protocol. In post-processing, TruHARP separates the spectral peaks in the acquired datasets, enabling high resolution motion and strain quantification. The impact of high resolution on circumferential and radial strain is studied using a realistic simulation and the improvement in strain maps is demonstrated in an in-vivo human study.

Preliminary validation of angle-independent myocardial elastography using MR tagging in a clinical setting

Myocardial elastography (ME), a radio-frequency (RF) based speckle tracking technique, was employed in order to image the entire two-dimensional (2D) transmural deformation field in full echocardiographic views and was validated against tagged magnetic resonance imaging (tMRI) in normal as well as reperfused (i.e., treated myocardial infarction [MI]) human left ventricles. RF ultrasound and tMRI frames were acquired at the papillary muscle level in 2D short-axis (SA) views at the frame rates of 136 (fps; real-time) and 33 fps (electrocardiogram [ECG]-gated), respectively. In ME, in-plane, 2D (lateral and axial) incremental displacements were iteratively estimated using one-dimensional (1D) cross-correlation and recorrelation techniques in a 2D search with a 1D matching kernel. In tMRI, cardiac motion was estimated by a template-matching algorithm on a 2D grid-shaped mesh. In both ME and tMRI, cumulative 2D displacements were obtained and then used to estimate 2D Lagrangian finite systolic strains, from which polar (i.e., radial and circumferential) strains, namely angle-independent measures, were further obtained through coordinate transformation. Principal strains, which are angle-independent and less centroid-dependent than polar strains, were also computed and imaged based on the 2D finite strains using methodology previously established. Both qualitatively and quantitatively, angle-independent ME is shown to be capable of (1) estimating myocardial deformation in good agreement with tMRI estimates in a clinical setting and of (2) differentiating abnormal from normal myocardium in a full left-ventricular view. The principal strains were concluded to be a potential diagnostic measure for detection of cardiac disease with reduced centroid dependence.

3D myocardial tissue tracking with slice followed cine DENSE MRI

PURPOSE

To track three-dimensional (3D) myocardial tissue motion using slice followed cine displacement encoded imaging with stimulated echoes (DENSE).

MATERIALS AND METHODS

Slice following (SF) has previously been developed for 2D myocardial tagging to compensate for the effect of through-plane motion on 2D tissue tracking. By incorporating SF into a cine DENSE sequence, and applying displacement encoding in three orthogonal directions, we demonstrate the ability to track discrete elements of a slice of myocardium in 3D as the heart moves through the cardiac cycle. The SF cine DENSE tracking algorithm was validated on a moving phantom, and the effects of through-plane motion on 2D cardiac strain were investigated in six healthy subjects.

RESULTS

A through-plane tracking accuracy of 0.46 +/- 0.32 mm was measured for a typical range of myocardial motion using a rotating phantom. In vivo 3D measurements of cardiac motion were consistent with prior myocardial tagging results. Through-plane rotation in a mid-ventricularshort-axis view was shown to decrease the magnitude of the 2D end-systolic circumferential strain by 3.91 +/- 0.43% and increase the corresponding radial strain by 6.01 +/- 1.07%.

CONCLUSION

Slice followed cine DENSE provides an accurate method for 3D tissue tracking.

Myocardial tagging with steady state free precession techniques and semi-automatic postprocessing--impact on diagnostic value

Our aim was to determine the diagnostic value of myocardial tagging sequences with regard to the evaluable share of the cardiac cycle. Thirty-three patients were examined at 1.5 T using tagging sequences with gradient-echo (GRE) readout, 18 patients at 1.5 T with steady-state free precession (SSFP), and 11 patients at 3 T using GRE. Two observers graded image quality and determined the share of the cardiac cycle for which postprocessing could be performed (1, optimal; 2, little interaction; 3, whole cycle assessable; 4, diastole non-assessable; 5, systole incomplete; 6, non-diagnostic). With GRE at 1.5 T, median image quality was 4.0 (95% CI 4.0-5.0), while it was significantly better with 2.0 (2.0-3.0) using the SSFP technique and similar at 3 T with 2.9 (1.7-3.5). With GRE at 1.5 T, systole could be assessed in 69% of patients, and an evaluation of the whole cardiac cycle was not possible. With the SSFP sequence at 1.5 T and GRE at 3 T, an evaluation of the whole cardiac cycle was possible in 71% and 70% of the patients, respectively, and systole was assessable in all patients. Tagging sequences with SSFP readout at 1.5 T make a semi-automatic evaluation of the whole cardiac cycle feasible in a large share of patients.

Automatic 3D tracking of cardiac material markers using slice-following and harmonic-phase MRI

A method to track a grid of cardiac material points in three dimensions using slice-following (SF) tagged magnetic resonance imaging (MRI) and harmonic-phase MRI is presented. A three-dimensional grid of material points on the lines of intersections of short-axis (SA) and long-axis (LA) planes is automatically tracked by combining two-dimensional pathlines that are computed on both SA and LA image planes. This process yields the true three-dimensional motion of points originating on the image plane intersections. Experimental data from normal volunteers, each obtained in four short breath-holds using the SF harmonic phase MRI pulse sequence, is presented. A validation of two-dimensional in-plane tracks using this pulse sequence on a moving phantom is also presented.

Tagged magnetic resonance imaging of the heart: a survey

Magnetic resonance imaging (MRI) of the heart with magnetization tagging provides a potentially useful new way to assess cardiac mechanical function, through revealing the local motion of otherwise indistinguishable portions of the heart wall. While still an evolving area, tagged cardiac MRI is already able to provide novel quantitative information on cardiac function. Exploiting this potential requires developing tailored methods for both imaging and image analysis. In this paper, we review some of the progress that has been made in developing such methods for tagged cardiac MRI, as well as some of the ways these methods have been applied to the study of cardiac function.

Spiral MR myocardial tagging

In the present study, complementary spatial modulation of magnetization (CSPAMM) myocardial tagging was extended with an interleaved spiral imaging sequence. The use of a spiral sequence enables the acquisition of grid-tagged images with a tagline distance as low as 4 mm in a single breath-hold. Alternatively, a high temporal resolution of 77 frames per second was obtained with 8-mm grid spacing. Ten healthy adult subjects were studied. With this new approach, high-quality images can be obtained and the tags persist throughout the entire cardiac cycle.

Steady-state free precession MR imaging: improved myocardial tag persistence and signal-to-noise ratio for analysis of myocardial motion

Tagging with balanced steady-state free-precession (SSFP) magnetic resonance (MR) imaging by using a steady-state storage scheme for myocardial motion analysis was evaluated. Signal-to-noise ratio (SNR), blood-tissue contrast, and tag persistence in volunteers and phantoms showed improved performance of SSFP imaging with tagging compared with that of radiofrequency spoiled gradient-echo (SPGR) MR imaging with tagging. Choice of flip angle with SSFP imaging involved a trade-off among SNR, blood-tissue contrast, and tag persistence. Increased SNR and tag persistence can be achieved simultaneously with SSFP imaging compared with SPGR tagging methods. As a result, the proposed technique may be useful for analysis of diastolic ventricular function.

Importance of k-space trajectory in echo-planar myocardial tagging at rest and during dobutamine stress

Hybrid fast gradient echo/echo-planar imaging (FGRE-EPI) can be used to increase temporal resolution, enhance tag contrast, and/or decrease scan time for breathhold myocardial tagging. However, off-resonance effects and motion can lead to local phase discontinuities in FGRE-EPI raw data when a conventional interleaved bottom-up k-space trajectory is used. These discontinuities can be particularly problematic for myocardial tagging, where the image energy is not only concentrated near the k-space origin, but is also concentrated in multiple spectral peaks centered throughout k-space. In this study, tag distortion artifacts in FGRE-EPI tagging due to off-resonance and velocity-induced phase discontinuities were characterized at rest and dobutamine stress, and the flyback and gradient moment smoothing (GMS) methods were shown to reduce these artifacts. For the specific parameters used in this study, flyback and GMS resulted in improved image quality at rest and stress, increased myocardium-tag contrast-to-noise ratio (11.4 +/- 2.1 vs. 10.0 +/- 2.9, P < 0.01 at rest; 11.1 +/- 1.8 vs. 8.1 +/- 2.4, P < 0.01 at stress), and reduced full width at half maximum of the tag profile (3.6 vs. 3.8 pixels at rest; 4.0 vs. 5.1 pixels at stress) compared to the conventional trajectory. A limitation of the improved trajectory is a parameter-dependent decrease in data acquisition efficiency. For the specific imaging protocol used, the repetition time of the improved trajectory increased by 36% compared to the conventional trajectory.

Estimating Motion From MRI Data

INVITED PAPER: Magnetic resonance imaging (MRI) is an ideal imaging modality to measure blood flow and tissue motion. It provides excellent contrast between soft tissues, and images can be acquired at positions and orientations freely defined by the user. From a temporal sequence of MR images, boundaries and edges of tissues can be tracked by image processing techniques. Additionally, MRI permits the source of the image signal to be manipulated. For example, temporary magnetic tags displaying a pattern of variable brightness may be placed in the object using MR saturation techniques, giving the user a known pattern to detect for motion tracking. The MRI signal is a modulated complex quantity, being derived from a rotating magnetic field in the form of an induced current. Well-defined patterns can also be introduced into the phase of the magnetization, and could be thought of as generalized tags. If the phase of each pixel is preserved during image reconstruction, relative phase shifts can be used to directly encode displacement, velocity and acceleration. New methods for modeling motion fields from MRI have now found application in cardiovascular and other soft tissue imaging. In this review, we shall describe the methods used for encoding, imaging, and modeling motion fields with MRI.

Magnitude image CSPAMM reconstruction (MICSR)

Image reconstruction of tagged cardiac MR images using complementary spatial modulation of magnetization (CSPAMM) requires the subtraction of two complex datasets to remove the untagged signal. Although the resultant images typically have sharper and more persistent tags than images formed without complementary tagging pulses, handling the complex data is problematic and tag contrast still degrades significantly during diastole. This article presents a magnitude image CSPAMM reconstruction (MICSR) method that is simple to implement and produces images with improved contrast and tag persistence. The MICSR method uses only magnitude images - i.e., no complex data - but yields tags with zero mean, sinusoidal profiles. A trinary display of MICSR images emphasizes their long tag persistence and demonstrates a novel way to visualize myocardial deformation. MICSR contrast and contrast-to-noise ratios (CNR) were evaluated using simulations, a phantom, and two normal volunteers. Tag contrast 1000 msec after the R wave trigger was 3.0 times better with MICSR than with traditional CSPAMM reconstruction techniques, while CNRs were 2.0 times better.

Real-time imaging of two-dimensional cardiac strain using a harmonic phase magnetic resonance imaging (HARP-MRI) pulse sequence

The harmonic phase (HARP) method provides automatic and rapid analysis of tagged magnetic resonance (MR) images for quantification and visualization of myocardial strain. In this article, the development and implementation of a pulse sequence that acquires HARP images in real time are described. In this pulse sequence, a CINE sequence of images with 1-1 spatial modulation of magnetization (SPAMM) tags are acquired during each cardiac cycle, alternating between vertical and horizontal tags in successive heartbeats. An incrementing train of imaging RF flip angles is used to compensate for the decay of the harmonic peaks due to both T(1) relaxation and the applied imaging pulses. The magnitude images displaying coarse anatomy are automatically reconstructed and displayed in real time after each heartbeat. HARP strain images are generated offline at a rate of four images per second; real-time processing should be possible with faster algorithms or computers. A comparison of myocardial contractility in non-breath-hold and breath-hold experiments in normal humans is presented.

Myocardial tagging with SSFP

This work presents the first implementation of myocardial tagging with refocused steady-state free precession (SSFP) and magnetization preparation. The combination of myocardial tagging (a noninvasive method for quantitative measurement of regional and global cardiac function) with the high tissue signal-to-noise ratio (SNR) obtained with SSFP is shown to yield improvements in terms of the myocardium-tag contrast-to-noise ratio (CNR) and tag persistence when compared to the current standard fast gradient-echo (FGRE) tagging protocol. Myocardium-tag CNR and tag persistence were studied using numerical simulations as well as phantom and human experiments. Both quantities were found to decrease with increasing imaging flip angle (alpha) due to an increased tag decay rate and a decrease in myocardial steady-state signal. However, higher alpha yielded better blood-myocardium contrast, indicating that optimal alpha is dependent on the application: higher alpha for better blood-myocardium boundary visualization, and lower alpha for better tag persistence. SSFP tagging provided the same myocardium-tag CNR as FGRE tagging when acquired at four times the bandwidth and better tag- and blood-myocardium CNRs than FGRE tagging when acquired at equal or twice the receiver bandwidth (RBW). The increased acquisition efficiency of SSFP allowed decreases in breath-hold duration, or increases in temporal resolution, as compared to FGRE.

MRI cardiac tagging using a sinc-modulated RF pulse train

A method of spatial modulation of magnetization using the RF pulse train is described for improved cardiac tagging. The method is based on a sinc-modulated RF pulse train in the presence of constant gradient. Such an RF pulse train yields a rectangular profile of saturated magnetization, and the bandwidth of the sinc envelope controls the width of the saturation tag. Compared to DANTE and SPAMM tagging techniques, this method offers sharper-edged tags and the ability to change the ratio of tag width to tag separation. The approach was implemented on a 9.4 T vertical bore NMR system and demonstrated with both phantom and in vivo mouse heart studies.

Phase contrast MRI of myocardial 3D strain by encoding contiguous slices in a single shot

Quantitative measurements of inherently three-dimensional (3D) cardiac strain and strain rate require 3D data; MRI provides uniquely high sensitivity to material strain by combining phase contrast with single-shot acquisition methods, such as echo-planar imaging (EPI). Previous MRI methods applied to 3D strain used multiple two-dimensional (2D) acquisitions and suffered loss of sensitivity due to magnification within the strain calculation of physiologic noise related to cardiac beat-to-beat variability. In the present work, each single-shot acquisition generates 3D image data by acquiring two contiguous 2D Fourier transform (FT) images in a single echo train of an EPI readout. Although strain encoding divides across multiple EPI shots, each strain component is computed only within single-shot data, avoiding noise magnification. Strain tensor maps are displayed using iconic 3D graphics or a simple color code of tensor shape. In a deforming gel phantom, gradient-recalled echo (GRE) MRI movies of 3D strain rates match expected strain fields. In normal human subjects, 3D strain rate tensor movies of heart and brain comprising seven slices in each of seven cardiac phases were completed in 56 heartbeats. Stimulated echo (STE) MRI of net systolic 3D strain was also demonstrated. Two-slices-in-one-shot spatial encoding permits a complete quantitative survey of ventricular 3D strain in under a minute, with routine patient supervision and turnkey image processing.

Myocardial wall tagging with undersampled projection reconstruction

Azimuthally undersampled projection reconstruction (PR) acquisition is investigated for use in myocardial wall tagging with MR using grid tags to provide increased temporal and spatial resolution. PR can provide the high-resolution images required for tagging with very few projections, at the expense of artifact. Insight is provided into the PR undersampling artifact, in the context of measuring myocardial motion with tags. For Fourier transform imaging, at least 112 phase-encodings must be collected to image tagging grids spaced 7 pixels apart. PR requires about 80 projections, a 1.4-fold reduction in scan time. Magn Reson Med 45:562-567, 2001. Published 2001 Wiley-Liss, Inc.

Variance components of two-dimensional strain parameters in the left-ventricular heart wall obtained by magnetic resonance tagging

This study quantifies variance components of two-dimensional strains in the left-ventricular heart wall assessed by magnetic resonance (MR) tagging in 18 healthy xxvolunteers. For a 7-mm tagging grid and homogeneous strain analysis, the intersubject variability and measurement error were estimated, as well as the intra- and interobserver variability. The variance components were calculated for the mean strain of a circumferential sector. The results show that the measurement error was almost equal to the intra-observer variability. With four circumferential sectors of 90 degrees each, approximately 65% of the total variance in epsilonr and epsilonc was due to intersubject variability, the remaining 35% was due to measurement error. With 12 sectors of 30 degrees each, the intersubject variability and measurement error both contributed 50% to the total variance. With 18 sectors of 20 degrees each, only 40% of the total variance was due to intersubject variability. The total variability increased with the number of sectors and therefore the number of sectors used in a study will be a trade-off between segment size (defining spatial resolution) and variability.

Ultrafast MR grid-tagging sequence for assessment of local mechanical properties of the lungs

While MR imaging with tagged magnetization has shown great utility in the study of muscle mechanics, the evaluation of pulmonary mechanics has long been hindered by the technical difficulties in MR imaging of lung parenchyma. In this study, a fast MR grid-tagging technique is described for dynamic assessment of regional pulmonary deformation. The method is based on a fast FLASH sequence with short TR and short TE. Tagging was achieved by using double DANTE pulse train or inversion pulses. Our results show that this technique is able to detect changes of the tagging grid caused by physiological deformation of the lung. Quantitative analysis of the data shows that this method is capable of assessing local pulmonary mechanics. The application of this technique could improve our understanding of ventilatory control, and thus provide a unique metric for assessing pulmonary disorders. Magn Reson Med 45:24-28, 2001.

Ultrafast pulse sequence techniques for cardiac magnetic resonance imaging

Cardiac magnetic resonance imaging is a rapidly emerging field that has seen tremendous advances in the past decade. Central to the development of effective imaging strategies has been the advent of high-performance gradient hardware and the exploitation of their speed characteristics through specialized pulse sequences well suited for cardiac imaging. These advances have facilitated unprecedented acquisition times that now approach echocardiographic frame rates, while maintaining excellent image quality. This article provides a detailed overview of advanced pulse sequence technology and approaches currently taken to maximize speed performance and image quality. In particular, segmented K-space techniques that include single-echo and multiecho spoiled gradient-echo imaging as well as steady-state free precession imaging are discussed. Finally, spiral and fast spin-echo techniques are explored. Examples of common applications of these pulse sequences are presented.

Detailed motion analysis of the left ventricular myocardium using an MR tagging method with a dense grid

Detailed analysis of myocardial deformation through a whole cardiac cycle was accomplished using a tagging method with a high-density grid. Four sets of tagged images with a 4-mm-spacing grid were measured by generating four tagging pulses arranged at regular intervals in the cardiac cycle. Through each set of images, tag intersections were tracked semi-automatically. The estimated motions of tag intersections were concatenated so that sequential positions of myocardium were connected through a whole cardiac cycle. In vitro evaluation of the precision of this technique showed that the mean error of tracked 4-mm tag intersections was less than 0.47 +/- 0.17 mm, even on the quite low-contrast images, and the concatenation error caused by double concatenation was comparable to the interpolation error in the subendocardial area obtained with 8-mm tag intersection motion. The small difference between the two mean distance curves of the in vivo evaluation indicated that the method is useful for analyzing heart wall abnormalities. Magn Reson Med 44:73-82, 2000.

Optimization of fast cardiac imaging using an echo-train readout

Fast gradient-echo sequences that use an echo-train readout are becoming more widely used, particularly for imaging the heart. An important issue for these sequences involves determining the optimal duration for the echo-train readout. In normal volunteer scans and theoretically the echo-train readout duration was varied from 2.4 to 32.8 msec. Myocardial signal-to-noise ratio (SNR), myocardium-tag signal difference-to-noise ratio (SDNR), flow artifact-to-noise ratio (FNR), and geometric distortion were measured and/or calculated. Our results showed that to obtain high SNR, SDNR, and data acquisition efficiency while minimizing FNR and geometric distortion, the readout duration should be 10-15 msec at 1.5 T.

MRI myocardial tagging

MRI myocardial tagging is now a well-developed method for evaluation of regional myocardial contraction. A series of progressively more refined imaging strategies, combined with advances in analytic strategies have provided a strong armamentarium of methods. Important insights into normal human physiology of left ventricular systolic and diastolic function have been developed using one-dimensional, two-dimensional and three-dimensional analyses of myocardial deformation. In disease states, improved understanding and detection of early alterations in myocardial function in hypertensive heart disease has been possible. In addition, improved understanding of effects of ischemia and infarction on regional function has been possible. Further, after acute myocardial infarction, clearer definition of the natural history of contractile dysfunction in the infarct region and the zone adjacent to the infarct have been possible. Similarly, effects on regional function of a number of important pharmacologic agents used for treatment, such as angiotensin converting enzyme inhibitors, beta blockers and angiotensin receptor blockers have been characterized. In the cardiomyopathies, myocardial tagging has permitted more reliable assessment of heterogeneity of segmental function, especially in hypertrophic cardiomyopathy. Finally, initial applications of myocardial tagging to assessment of right ventricular regional function in hypertrophied hearts with and without major congenital abnormalities have generated advances in understanding of effects of hypertrophy on right ventricular function.J. Magn. Reson. Imaging 1999;10:609-616.