Multiple 3D inversion recovery imaging for volume T1 mapping of the heart

The future of cardiovascular magnetic resonance: All-in-one vs. real-time (Part 1)

Cardiovascular magnetic resonance (CMR) protocols can be lengthy and complex, which has driven the research community to develop new technologies to make these protocols more efficient and patient-friendly. Two different approaches to improving CMR have been proposed, specifically "all-in-one" CMR, where several contrasts and/or motion states are acquired simultaneously, and "real-time" CMR, in which the examination is accelerated to avoid the need for breathholding and/or cardiac gating. The goal of this two-part manuscript is to describe these two different types of emerging rapid CMR. To this end, the vision of each is described, along with techniques which have been devised and tested along the pathway of clinical implementation. The pros and cons of the different methods are presented, and the remaining open needs of each are detailed. Part 1 will tackle the "all-in-one" approaches, and Part 2 the "real-time" approaches along with an overall summary of these emerging methods.

Recent advances in PET-MRI for cardiac sarcoidosis

The diagnosis of cardiac sarcoidosis (CS) remains challenging. While only a small fraction of patients with systemic sarcoidosis present with clinically symptomatic CS, cardiac involvement has been associated with adverse outcomes, such as ventricular arrhythmia, heart block, heart failure and sudden cardiac death. Despite the clinical relevance of having an early and accurate diagnosis of CS, there is no gold-standard technique available for the assessment of CS. Non-invasive PET and MR imaging have shown promise in the detection of different histopathological features of CS. More recently, the introduction of hybrid PET-MR scanners has enabled the acquisition of these hallmarks in a single scan, demonstrating higher sensitivity and specificity for CS detection and risk stratification than with either imaging modality alone. This article describes recent developments in hybrid PET-MR imaging for improving the diagnosis of CS and discusses areas of future development that could make cardiac PET-MRI the preferred diagnostic tool for the comprehensive assessment of CS.

SAturation-recovery and Variable-flip-Angle-based three-dimensional free-breathing cardiovascular magnetic resonance T1 mapping at 3 T

The purpose of the current study was to develop and validate a three-dimensional (3D) free-breathing cardiac T₁ -mapping sequence using SAturation-recovery and Variable-flip-Angle (SAVA). SAVA sequentially acquires multiple electrocardiogram-triggered volumes using a multishot spoiled gradient-echo sequence. The first volume samples the equilibrium signal of the longitudinal magnetization, where a flip angle of 2° is used to reduce the time for the magnetization to return to equilibrium. The succeeding three volumes are saturation prepared with variable delays, and are acquired using a 15° flip angle to maintain the signal-to-noise ratio. A diaphragmatic navigator is used to compensate the respiratory motion. T₁ is calculated using a saturation-recovery model that accounts for the flip angle. We validated SAVA by simulations, phantom, and human subject experiments at 3 T. SAVA was compared with modified Look-Locker inversion recovery (MOLLI) and saturation-recovery single-shot acquisition (SASHA) in vivo. In phantoms, T₁ by SAVA had good agreement with the reference (R²  = 0.99). In vivo 3D T₁ mapping by SAVA could achieve an imaging resolution of 1.25 × 1.25 × 8 mm³ . Both global and septal T₁ values by SAVA (1347 ± 37 and 1332 ± 42 ms) were in between those by SASHA (1612 ± 63 and 1618 ± 51 ms) and MOLLI (1143 ± 59 and 1188 ± 65 ms). According to the standard deviation (SD) and coefficient of variation (CV), T₁ precision measured by SAVA (SD: 99 ± 14 and 60 ± 8 ms; CV: 7.4% ± 0.9% and 4.5% ± 0.6%) was comparable with MOLLI (SD: 99 ± 25 and 46 ± 12 ms; CV: 8.8% ± 2.5% and 3.9% ± 1.1%) and superior to SASHA (SD: 222 ± 89 and 132 ± 33 ms; CV: 13.8% ± 5.5% and 8.1% ± 2.0%). It was concluded that the proposed free-breathing SAVA sequence enables more efficient 3D whole-heart T₁ estimation with good accuracy and precision.

Free-breathing simultaneous myocardial T1 and T2 mapping with whole left ventricle coverage

PURPOSE

To develop a free-breathing sequence, that is, Multislice Joint T₁ -T₂ , for simultaneous measurement of myocardial T₁ and T₂ for multiple slices to achieve whole left-ventricular coverage.

METHODS

Multislice Joint T₁ -T₂ adopts slice-interleaved acquisition to collect 10 single-shot electrocardiogram-triggered images for each slice prepared by saturation and T₂ preparation to simultaneously estimate myocardial T₁ and T₂ and achieve whole left-ventricular coverage. Prospective slice-tracking using a respiratory navigator and retrospective image registration are used to reduce through-plane and in-plane motion, respectively. Multislice Joint T₁ -T₂ was validated through numerical simulations and phantom and in vivo experiments, and compared with saturation-recovery single-shot acquisition and T₂ -prepared balanced Steady-State Free Precession (T₂ -prep SSFP) sequences.

RESULTS

Phantom T₁ and T₂ from Multislice Joint T₁ -T₂ had good accuracy and precision, and were insensitive to heart rate. Multislice Joint T₁ -T₂ yielded T₁ and T₂ maps of nine left-ventricular slices in 1.4 minutes. The mean left-ventricular T₁ difference between saturation-recovery single-shot acquisition and Multislice Joint T₁ -T₂ across healthy subjects and patients was 191 ms (1564 ± 60 ms versus 1373 ± 50 ms; P < .05) and 111 ms (1535 ± 49 ms vs 1423 ± 49 ms; P < .05), respectively. The mean difference in left-ventricular T₂ between T₂ -prep SSFP and Multislice Joint T₁ -T₂ across healthy subjects and patients was -6.3 ms (42.4 ± 1.4 ms vs 48.7 ± 2.5; P < .05) and -5.7 ms (41.6 ± 2.5 ms vs 47.3 ± 2.7; P < .05), respectively.

CONCLUSION

Multislice Joint T₁ -T₂ enables quantification of whole left-ventricular T₁ and T₂ during free breathing within a clinically feasible scan time of less than 2 minutes.

Free-breathing whole-heart multi-slice myocardial T1 mapping in 2 minutes

PURPOSE

To develop and validate a saturation-delay-inversion recovery preparation, slice tracking and multi-slice based sequence for measuring whole-heart native T₁ .

METHOD

The proposed free-breathing sequence performs T₁ mapping of multiple left-ventricular slices by slice-interleaved acquisition to collect 10 electrocardiogram-triggered single-shot slice-selective images for each slice. A saturation-delay-inversion recovery pulse is used for T₁ preparation. Prospective slice tracking by the diaphragm navigator and retrospective registration are used to reduce through-plane and in-plane motion, respectively. The proposed sequence was validated in both phantom and human subjects (12 healthy subjects and 15 patients who were referred for a clinical cardiac MR exam) and compared with saturation recovery single-shot acquisition (SASHA) and modified Look-Locker inversion recovery (MOLLI).

RESULTS

Phantom T₁ measured by the proposed sequence had excellent agreement (R²  = 0.99) with the ground-truth T₁ and was insensitive to heart rate. In both healthy subjects and patients, the proposed sequence yielded nine left-ventricular T₁ maps per volume in less than 2 minutes (healthy volunteers: 1.8 ± 0.4 minutes; patients: 1.9 ± 0.2 minutes). The average T₁ of whole left ventricle for all healthy subjects and patients were 1560 ± 61 and 1535 ± 49 ms by SASHA, 1208 ± 42 and 1233 ± 56 ms by MOLLI5(3)3, and 1397 ± 34 and 1433 ± 56 ms by the proposed sequence, respectively. The corresponding coefficient of variation of T₁ were 6.2 ± 1.4% and 5.8 ± 1.6%, 5.3 ± 1.1% and 5.1 ± 0.8%, and 4.9 ± 0.8% and 4.5 ± 0.8%, respectively.

CONCLUSION

The proposed sequence enables quantification of whole heart T₁ with good accuracy and precision in less than 2 minutes during free breathing.

Inflow artifact reduction using an adaptive flip-angle navigator restore pulse for late gadolinium enhancement of the left atrium

PURPOSE

Late gadolinium enhancement (LGE) of the left atrium is susceptible to artifacts arising from the right pulmonary veins, caused by inflowing blood tagged by the navigator restore pulse. The purpose of this study was to evaluate a new method to reduce the inflow artifact using an adaptive flip-angle restore pulse.

METHODS

A low-restore angle reduces the inflow artifact but may lead to a poor navigator SNR. The proposed approach aims to determine the patient-specific restore angle, which optimizes the trade-off between inflow artifacts and navigator SNR. Three-dimensional LGE with adaptive navigator restore (3D LGE_A ) was implemented by incrementing the flip angle of the restore pulse from a starting value of 0°, based on the navigator normalized cross-correlation. Magnetic resonance imaging experiments were performed on a 1.5T scanner. The value of 3D LGE_A was compared with 3D LGE with a constant 180° restore pulse (3D LGE₁₈₀ ) in 22 patients with heart diseases. The values of 3D LGE_A and 3D LGE₁₈₀ were compared in terms of pulmonary vein blood signal relative to reference blood in the descending aorta (PV_rel ) and visual scoring to determine level of motion artifacts using a 4-point scale (1 = severe artifacts; 4 = no artifacts).

RESULTS

The value of PV_rel was significantly lower for 3D LGE_A than for 3D LGE₁₈₀ (1.16 ± 0.23 vs. 1.59 ± 0.29, P < .001). Furthermore, visual scoring of the motion artifacts yielded no difference (P = .78).

CONCLUSION

Adaptively adjusting the navigator restore flip angle based on the navigator normalized cross-correlation reduces the 3D LGE inflow artifact without affecting image quality or the scan time.

FASt single-breathhold 2D multislice myocardial T1 mapping (FAST1) at 1.5T for full left ventricular coverage in three breathholds

Background

Conventional myocardial T₁ mapping techniques such as modified Look–Locker inversion recovery (MOLLI) generate one T₁ map per breathhold. T₁ mapping with full left ventricular coverage may be desirable when spatial T₁ variations are expected. This would require multiple breathholds, increasing patient discomfort and prolonging scan time.

Purpose

To develop and characterize a novel FASt single‐breathhold 2D multislice myocardial T₁ mapping (FAST1) technique for full left ventricular coverage.

Study Type

Prospective.

Population/Phantom

Numerical simulation, agarose/NiCl₂ phantom, 9 healthy volunteers, and 17 patients.

Field Strength/Sequence

1.5T/FAST1.

Assessment

Two FAST1 approaches, FAST1‐BS and FAST1‐IR, were characterized and compared with standard 5‐(3)‐3 MOLLI in terms of accuracy, precision/spatial variability, and repeatability.

Statistical Tests

Kruskal‐Wallis, Wilcoxon signed rank tests, intraclass correlation coefficient analysis, analysis of variance, Student's t‐tests, Pearson correlation analysis, and Bland–Altman analysis.

Results

In simulation/phantom, FAST1‐BS, FAST1‐IR, and MOLLI had an accuracy (expressed as T₁ error) of 0.2%/4%, 6%/9%, and 4%/7%, respectively, while FAST1‐BS and FAST1‐IR had a precision penalty of 1.7/1.5 and 1.5/1.4 in comparison with MOLLI, respectively. In healthy volunteers, FAST1‐BS/FAST1‐IR/MOLLI led to different native myocardial T₁ times (1016 ± 27 msec/952 ±22 msec/987 ± 23 msec, P < 0.0001) and spatial variability (66 ± 10 msec/57 ± 8 msec/46 ± 7 msec, P < 0.001). There were no statistically significant differences between all techniques for T₁ repeatability (P = 0.18). In vivo native and postcontrast myocardial T₁ times in both healthy volunteers and patients using FAST1‐BS/FAST1‐IR were highly correlated with MOLLI (Pearson correlation coefficient ≥0.93).

Data Conclusion

FAST1 enables myocardial T₁ mapping with full left ventricular coverage in three separated breathholds. In comparison with MOLLI, FAST1 yield a 5‐fold increase of spatial coverage, limited penalty of T₁ precision/spatial variability, no significant difference of T₁ repeatability, and highly correlated T₁ times. FAST1‐IR provides improved T₁ precision/spatial variability but reduced accuracy when compared with FAST1‐BS.

Level of Evidence: 1

Technical Efficacy: Stage 3

J. Magn. Reson. Imaging 2020;51:492–504.

Accelerated free-breathing 3D T1ρ cardiovascular magnetic resonance using multicoil compressed sensing

BACKGROUND

Endogenous contrast T1ρ cardiovascular magnetic resonance (CMR) can detect scar or infiltrative fibrosis in patients with ischemic or non-ischemic cardiomyopathy. Existing 2D T1ρ techniques have limited spatial coverage or require multiple breath-holds. The purpose of this project was to develop an accelerated, free-breathing 3D T1ρ mapping sequence with whole left ventricle coverage using a multicoil, compressed sensing (CS) reconstruction technique for rapid reconstruction of undersampled k-space data.

METHODS

We developed a cardiac- and respiratory-gated, free-breathing 3D T1ρ sequence and acquired data using a variable-density k-space sampling pattern (A = 3). The effect of the transient magnetization trajectory, incomplete recovery of magnetization between T1ρ-preparations (heart rate dependence), and k-space sampling pattern on T1ρ relaxation time error and edge blurring was analyzed using Bloch simulations for normal and chronically infarcted myocardium. Sequence accuracy and repeatability was evaluated using MnCl2 phantoms with different T1ρ relaxation times and compared to 2D measurements. We further assessed accuracy and repeatability in healthy subjects and compared these results to 2D breath-held measurements.

RESULTS

The error in T1ρ due to incomplete recovery of magnetization between T1ρ-preparations was T1ρhealthy = 6.1% and T1ρinfarct = 10.8% at 60 bpm and T1ρhealthy = 13.2% and T1ρinfarct = 19.6% at 90 bpm. At a heart rate of 60 bpm, error from the combined effects of readout-dependent magnetization transients, k-space undersampling and reordering was T1ρhealthy = 12.6% and T1ρinfarct = 5.8%. CS reconstructions had improved edge sharpness (blur metric = 0.15) compared to inverse Fourier transform reconstructions (blur metric = 0.48). There was strong agreement between the mean T1ρ estimated from the 2D and accelerated 3D data (R2 = 0.99; P < 0.05) acquired on the MnCl2 phantoms. The mean R1ρ estimated from the accelerated 3D sequence was highly correlated with MnCl2 concentration (R2 = 0.99; P < 0.05). 3D T1ρ acquisitions were successful in all human subjects. There was no significant bias between undersampled 3D T1ρ and breath-held 2D T1ρ (mean bias = 0.87) and the measurements had good repeatability (COV2D = 6.4% and COV3D = 7.1%).

CONCLUSIONS

This is the first report of an accelerated, free-breathing 3D T1ρ mapping of the left ventricle. This technique may improve non-contrast myocardial tissue characterization in patients with heart disease in a scan time appropriate for patients.

Three-dimensional free breathing whole heart cardiovascular magnetic resonance T1 mapping at 3 T

BACKGROUND

This study demonstrates a three-dimensional (3D) free-breathing native myocardial T1 mapping sequence at 3 T.

METHODS

The proposed sequence acquires three differently T1-weighted volumes. The first two volumes receive a saturation pre-pulse with different recovery time. The third volume is acquired without magnetization preparation and after a significant recovery time. Respiratory navigator gating and volume-interleaved acquisition are adopted to mitigate misregistration. The proposed sequence was validated through simulation, phantom experiments and in vivo experiments in 12 healthy adult subjects.

RESULTS

In phantoms, good agreement on T1 measurement was achieved between the proposed sequence and the reference inversion recovery spin echo sequence (R2 = 0.99). Homogeneous 3D T1 maps were obtained from healthy adult subjects, with a T1 value of 1476 ± 53 ms and a coefficient of variation (CV) of 6.1 ± 1.4% over the whole left-ventricular myocardium. The averaged septal T1 was 1512 ± 60 ms with a CV of 2.1 ± 0.5%.

CONCLUSION

Free-breathing 3D native T1 mapping at 3 T is feasible and may be applicable in myocardial assessment. The proposed 3D T1 mapping sequence is suitable for applications in which larger coverage is desired beyond that available with single-shot parametric mapping, or breath-holding is unfeasible.

Single breath-hold 3D cardiac T1 mapping using through-time spiral GRAPPA

The quantification of cardiac T1 relaxation time holds great potential for the detection of various cardiac diseases. However, as a result of both cardiac and respiratory motion, only one two-dimensional T1 map can be acquired in one breath-hold with most current techniques, which limits its application for whole heart evaluation in routine clinical practice. In this study, an electrocardiogram (ECG)-triggered three-dimensional Look-Locker method was developed for cardiac T1 measurement. Fast three-dimensional data acquisition was achieved with a spoiled gradient-echo sequence in combination with a stack-of-spirals trajectory and through-time non-Cartesian generalized autocalibrating partially parallel acquisition (GRAPPA) acceleration. The effects of different magnetic resonance parameters on T1 quantification with the proposed technique were first examined by simulating data acquisition and T1 map reconstruction using Bloch equation simulations. Accuracy was evaluated in studies with both phantoms and healthy subjects. These results showed that there was close agreement between the proposed technique and the reference method for a large range of T1 values in phantom experiments. In vivo studies further demonstrated that rapid cardiac T1 mapping for 12 three-dimensional partitions (spatial resolution, 2 × 2 × 8 mm3 ) could be achieved in a single breath-hold of ~12 s. The mean T1 values of myocardial tissue and blood obtained from normal volunteers at 3 T were 1311 ± 66 and 1890 ± 159 ms, respectively. In conclusion, a three-dimensional T1 mapping technique was developed using a non-Cartesian parallel imaging method, which enables fast and accurate T1 mapping of cardiac tissues in a single short breath-hold.

T1-refBlochi: high resolution 3D post-contrast T1 myocardial mapping based on a single 3D late gadolinium enhancement volume, Bloch equations, and a reference T1

BACKGROUND

High resolution 3D T1 mapping is important for assessment of diffuse myocardial fibrosis in left atrium or other thin-walled structures. In this work, we investigated a fast single-TI 3D high resolution T1 mapping method that directly transforms a 3D late gadolinium enhancement (LGE) volume to a 3D T1 map.

METHODS

The proposed method, T1-refBlochi, is based on Bloch equation modeling of the LGE signal, a single-point calibration, and assumptions that proton density and T2* are relatively uniform in the heart. Several sources of error of this method were analyzed mathematically and with simulations. Imaging was performed in phantoms, eight swine and five patients, comparing T1-refBlochi to a standard spin-echo T1 mapping, 3D multi-TI T1 mapping, and 2D ShMOLLI, respectively.

RESULTS

The method has a good accuracy and adequate precision, even considering various sources of error. In phantoms, over a range of protocols, heart-rates and T1 s, the bias ±1SD was -3 ms ± 9 ms. The porcine studies showed excellent agreement between T1-refBlochi and the multi-TI method (bias ±1SD = -6 ± 22 ms). The proton density and T2* weightings yielded ratios for scar/blood of 0.94 ± 0.01 and for myocardium/blood of 1.03 ± 0.02 in the eight swine, confirming that sufficient uniformity of proton density and T2* weightings exists among heterogeneous tissues of the heart. In the patients, the mean T1 bias ±1SD in myocardium and blood between T1-refBlochi and ShMOLLI was -9 ms ± 21 ms.

CONCLUSION

T1-refBlochi provides a fast single-TI high resolution 3D T1 map of the heart with good accuracy and adequate precision.

Simultaneous multislice imaging for native myocardial T1 mapping: Improved spatial coverage in a single breath-hold

PURPOSE

To develop a saturation recovery myocardial T1 mapping method for the simultaneous multislice acquisition of three slices.

METHODS

Saturation pulse-prepared heart rate independent inversion recovery (SAPPHIRE) T1 mapping was implemented with simultaneous multislice imaging using FLASH readouts for faster coverage of the myocardium. Controlled aliasing in parallel imaging (CAIPI) was used to achieve minimal noise amplification in three slices. Multiband reconstruction was performed using three linear reconstruction methods: Slice- and in-plane GRAPPA, CG-SENSE, and Tikhonov-regularized CG-SENSE. Accuracy, spatial variability, and interslice leakage were compared with single-band T1 mapping in a phantom and in six healthy subjects.

RESULTS

Multiband phantom T1 times showed good agreement with single-band T1 mapping for all three reconstruction methods (normalized root mean square error <1.0%). The increase in spatial variability compared with single-band imaging was lowest for GRAPPA (1.29-fold), with higher penalties for Tikhonov-regularized CG-SENSE (1.47-fold) and CG-SENSE (1.52-fold). In vivo multiband T1 times showed no significant difference compared with single-band (T1 time ± intersegmental variability: single-band, 1580 ± 119 ms; GRAPPA, 1572 ± 145 ms; CG-SENSE, 1579 ± 159 ms; Tikhonov, 1586 ± 150 ms [analysis of variance; P = 0.86]). Interslice leakage was smallest for GRAPPA (5.4%) and higher for CG-SENSE (6.2%) and Tikhonov-regularized CG-SENSE (7.9%).

CONCLUSION

Multiband accelerated myocardial T1 mapping demonstrated the potential for single-breath-hold T1 quantification in 16 American Heart Association segments over three slices. A 1.2- to 1.4-fold higher in vivo spatial variability was observed, where GRAPPA-based reconstruction showed the highest homogeneity and the least interslice leakage. Magn Reson Med 78:462-471, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Three-dimensional Cardiac MR Imaging: Related Techniques and Clinical Applications

Three-dimensional (3D) cardiac magnetic resonance (MR) imaging has several advantages, including the easy coverage of the entire heart without misregistration, reduction of breath-holding times, and availability for postprocessing reconstruction. These advantages are associated with some techniques such as breath-hold or navigator gating and parallel imaging. However, the image quality of 3D cardiac MR images is compromised by the use of a shorter repetition time and parallel imaging. Thus, a steady-state free precession sequence, contrast agent administration, and presaturation pulses are used to maintain the image quality. In this review, we introduce the MR imaging techniques used in 3D cardiac MR imaging and demonstrate the typical 3D cardiac MR images, followed by discussion about their advantages and disadvantages.

3D myocardial T1 mapping using saturation recovery

Purpose

To propose a 3D quantitative high‐resolution T₁ mapping technique, called 3D SASHA (saturation‐recovery single‐shot acquisition), which combines a saturation recovery pulse with 1D‐navigator‐based‐respiratory motion compensation to acquire the whole volume of the heart in free breathing. The sequence was tested and validated both in a T₁ phantom and in healthy subjects.

Materials and Methods

The 3D SASHA method was implemented on a 1.5T scanner. A diaphragmatic navigator was used to allow free‐breathing acquisition and the images were acquired with a resolution of 1.4 × 1.4 × 8 mm³. For assessment of accuracy and precision the sequence was compared with the reference gold‐standard inversion‐recovery spin echo (IRSE) pulse sequence in a T₁ phantom, while for the in vivo studies (10 healthy volunteers) 3D SASHA was compared with the clinically used 2D MOLLI (3‐3‐5) and 2D SASHA protocols.

Results

There was good agreement between the T₁ values measured in a T₁ phantom with 3D SASHA and the reference IRSE pulse sequences (1111.6 ± 31 msec vs. 1123.6 ± 8 msec, P = 0.9947). Mean and standard deviation of the myocardial T₁ values in healthy subjects measured with 2D MOLLI, 2D SASHA, and 3D SASHA sequences were 881 ± 40 msec, 1181.3 ± 32 msec, and 1153.6 ± 28 msec respectively.

Conclusion

The proposed 3D SASHA sequence allows for high‐resolution free‐breathing whole‐heart T₁‐mapping with T₁ values in good agreement with the 2D SASHA and improved precision.

Level of Evidence: 2

Technical Efficacy: Stage 1

J. MAGN. RESON. IMAGING 2017;46:218–227

Fast and robust 3D T1 mapping using spiral encoding and steady RF excitation at 7 T: application to cardiac manganese enhanced MRI (MEMRI) in mice

Mapping longitudinal relaxation times in 3D is a promising quantitative and non-invasive imaging tool to assess cardiac remodeling. Few methods are proposed in the literature allowing us to perform 3D T1 mapping. These methods often require long scan times and use a low number of 3D images to calculate T1 . In this project, a fast 3D T1 mapping method using a stack-of-spirals sampling scheme and regular RF pulse excitation at 7 T is presented. This sequence, combined with a newly developed fitting procedure, allowed us to quantify T1 of the whole mouse heart with a high spatial resolution of 208 × 208 × 315 µm(3) in 10-12 min acquisition time. The sensitivity of this method for measuring T1 variations was demonstrated on mouse hearts after several injections of manganese chloride (doses from 25 to 150 µmol kg(-1) ). T1 values were measured in vivo in both pre- and post-contrast experiments. This protocol was also validated on ischemic mice to demonstrate its efficiency to visualize tissue damage induced by a myocardial infarction. This study showed that combining spiral gradient shape and steady RF excitation enabled fast and robust 3D T1 mapping of the entire heart with a high spatial resolution.

Simultaneous three-dimensional myocardial T1 and T2 mapping in one breath hold with 3D-QALAS

BACKGROUND

Quantification of the longitudinal- and transverse relaxation time in the myocardium has shown to provide important information in cardiac diagnostics. Methods for cardiac relaxation time mapping generally demand a long breath hold to measure either T1 or T2 in a single 2D slice. In this paper we present and evaluate a novel method for 3D interleaved T1 and T2 mapping of the whole left ventricular myocardium within a single breath hold of 15 heartbeats.

METHODS

The 3D-QALAS (3D-quantification using an interleaved Look-Locker acquisition sequence with T2 preparation pulse) is based on a 3D spoiled Turbo Field Echo sequence using inversion recovery with interleaved T2 preparation. Quantification of both T1 and T2 in a volume of 13 slices with a resolution of 2.0x2.0x6.0 mm is obtained from five measurements by using simulations of the longitudinal magnetizations Mz. This acquisition scheme is repeated three times to sample k-space. The method was evaluated both in-vitro (validated against Inversion Recovery and Multi Echo) and in-vivo (validated against MOLLI and Dual Echo).

RESULTS

In-vitro, a strong relation was found between 3D-QALAS and Inversion Recovery (R = 0.998; N = 10; p < 0.01) and between 3D-QALAS and Multi Echo (R = 0.996; N = 10; p < 0.01). The 3D-QALAS method showed no dependence on e.g. heart rate in the interval of 40-120 bpm. In healthy myocardium, the mean T1 value was 1083 ± 43 ms (mean ± SD) for 3D-QALAS and 1089 ± 54 ms for MOLLI, while the mean T2 value was 50.4 ± 3.6 ms 3D-QALAS and 50.3 ± 3.5 ms for Dual Echo. No significant difference in in-vivo relaxation times was found between 3D-QALAS and MOLLI (N = 10; p = 0.65) respectively 3D-QALAS and Dual Echo (N = 10; p = 0.925) for the ten healthy volunteers.

CONCLUSIONS

The 3D-QALAS method has demonstrated good accuracy and intra-scan variability both in-vitro and in-vivo. It allows rapid acquisition and provides quantitative information of both T1 and T2 relaxation times in the same scan with full coverage of the left ventricle, enabling clinical application in a broader spectrum of cardiac disorders.

Cardiac MRI: a central prognostic tool in myocardial fibrosis

Fibrotic remodelling of the extracellular matrix is a healing mechanism necessary immediately after myocardial injury. However, prolonged increase in myocardial fibrotic activity results in stiffening of the myocardium and heralds adverse outcomes related to systolic and diastolic dysfunction, as well as arrhythmogenesis. Cardiac MRI provides a noninvasive phenotyping tool for accurate and easy detection and quantification of myocardial fibrosis by probing the retention of gadolinium-contrast agent in myocardial tissue. Late-gadolinium enhancement (LGE) cardiac MRI has been used extensively in a large number of studies for measurement of myocardial scarring. T1 mapping, a fairly new technique that can be used to identify the exact T1 value of the tissue, provides a direct measurement of the extracellular volume fraction of the myocardium. In contrast to LGE, T1 mapping can be used to measure diffuse myocardial fibrosis and differentiate between disease processes. In this Review, we describe the basic principles of imaging myocardial fibrosis using contrast-enhanced MRI and summarize its use for prognostic purposes.

Free-breathing combined three-dimensional phase sensitive late gadolinium enhancement and T1 mapping for myocardial tissue characterization

PURPOSE

To develop a novel MR sequence for combined three-dimensional (3D) phase-sensitive (PS) late gadolinium enhancement (LGE) and T1 mapping to allow for simultaneous assessment of focal and diffuse myocardial fibrosis.

METHODS

In the proposed sequence, four 3D imaging volumes are acquired with different T1 weightings using a combined saturation and inversion preparation, after administration of a gadolinium contrast agent. One image is acquired fully sampled with the inversion time selected to null the healthy myocardial signal (the LGE image). The other three images are three-fold under-sampled and reconstructed using compressed sensing. An acquisition scheme with two interleaved imaging cycles and joint navigator-gating of those cycles ensures spatial registration of the imaging volumes. T1 maps are generated using all four imaging volumes. The signal-polarity in the LGE image is restored using supplementary information from the T1 fit to generate PS-LGE images. The accuracy of the proposed method was assessed with respect to a inversion-recovery spin-echo sequence. In vivo T1 maps and LGE images were acquired with the proposed sequence and quantitatively compared with 2D multislice Modified Look-Locker inversion recovery (MOLLI) T1 maps. Exemplary images in a patient with focal scar were compared with conventional LGE imaging.

RESULTS

The deviation of the proposed method and the spin-echo reference was < 11 ms in phantom for T1 times between 250 and 600 ms, regardless of the inversion time selected in the LGE image. There was no significant difference in the in vivo T1 times of the proposed sequence and the 2D MOLLI technique (myocardium: 292 ± 75 ms versus 310 ± 49 ms, blood-pools: 191 ± 75 ms versus 182.0 ± 33). The LGE images showed proper nulling of the healthy myocardium in all subjects and clear depiction of scar in the patient.

CONCLUSION

The proposed sequence enables simultaneous acquisition of 3D PS-LGE images and spatially registered 3D T1 maps in a single scan.

Adaptive registration of varying contrast-weighted images for improved tissue characterization (ARCTIC): application to T1 mapping

PURPOSE

To propose and evaluate a novel nonrigid image registration approach for improved myocardial T1 mapping.

METHODS

Myocardial motion is estimated as global affine motion refined by a novel local nonrigid motion estimation algorithm. A variational framework is proposed, which simultaneously estimates motion field and intensity variations, and uses an additional regularization term to constrain the deformation field using automatic feature tracking. The method was evaluated in 29 patients by measuring the DICE similarity coefficient and the myocardial boundary error in short axis and four chamber data. Each image series was visually assessed as "no motion" or "with motion." Overall T1 map quality and motion artifacts were assessed in the 85 T1 maps acquired in short axis view using a 4-point scale (1-nondiagnostic/severe motion artifact, 4-excellent/no motion artifact).

RESULTS

Increased DICE similarity coefficient (0.78 ± 0.14 to 0.87 ± 0.03, P < 0.001), reduced myocardial boundary error (1.29 ± 0.72 mm to 0.84 ± 0.20 mm, P < 0.001), improved overall T1 map quality (2.86 ± 1.04 to 3.49 ± 0.77, P < 0.001), and reduced T1 map motion artifacts (2.51 ± 0.84 to 3.61 ± 0.64, P < 0.001) were obtained after motion correction of "with motion" data (∼56% of data).

CONCLUSIONS

The proposed nonrigid registration approach reduces the respiratory-induced motion that occurs during breath-hold T1 mapping, and significantly improves T1 map quality.

Accuracy, precision, and reproducibility of four T1 mapping sequences: a head-to-head comparison of MOLLI, ShMOLLI, SASHA, and SAPPHIRE

PURPOSE

To compare accuracy, precision, and reproducibility of four commonly used myocardial T1 mapping sequences: modified Look-Locker inversion recovery (MOLLI), shortened MOLLI (ShMOLLI), saturation recovery single-shot acquisition (SASHA), and saturation pulse prepared heart rate independent inversion recovery (SAPPHIRE).

MATERIALS AND METHODS

This HIPAA-compliant study was approved by the institutional review board. All subjects provided written informed consent. Accuracy, precision, and reproducibility of the four T1 mapping sequences were first compared in phantom experiments. In vivo analysis was performed in seven healthy subjects (mean age ± standard deviation, 38 years ± 19; four men, three women) who were imaged twice on two separate days. In vivo reproducibility of native T1 mapping and extracellular volume (ECV) were measured. Differences between the sequences were assessed by using Kruskal-Wallis and Wilcoxon rank sum tests (phantom data) and mixed-effect models (in vivo data).

RESULTS

T1 mapping accuracy in phantoms was lower with ShMOLLI (62 msec) and MOLLI (44 msec) than with SASHA (13 msec; P < .05) and SAPPHIRE (12 msec; P < .05). MOLLI had similar precision to ShMOLLI (4.0 msec vs 5.6 msec; P = .07) but higher precision than SAPPHIRE (6.8 msec; P = .002) and SASHA (8.7 msec; P < .001). All sequences had similar reproducibility in phantoms (P = .1). The four sequences had similar in vivo reproducibility for native T1 mapping (∼25-50 msec; P > .05) and ECV quantification (∼0.01-0.02; P > .05).

CONCLUSION

SASHA and SAPPHIRE yield higher accuracy, lower precision, and similar reproducibility compared with MOLLI and ShMOLLI for T1 measurement. Different sequences yield different ECV values; however, all sequences have similar reproducibility for ECV quantification.

Free-breathing post-contrast three-dimensional T1 mapping: Volumetric assessment of myocardial T1 values

PURPOSE

To develop a three-dimensional (3D) free-breathing myocardial T1 mapping sequence for assessment of left ventricle diffuse fibrosis after contrast administration.

METHODS

In the proposed sequence, multiple 3D inversion recovery images are acquired in an interleaved manner. A mixed prospective/retrospective navigator scheme is used to obtain the 3D Cartesian k-space data with fully sampled center and randomly undersampled outer k-space. The resulting undersampled 3D k-space data are then reconstructed using compressed sensing. Subsequently, T1 maps are generated by voxel-wise curve fitting of the individual interleaved images. In a phantom study, the accuracy of the 3D sequence was evaluated against two-dimensional (2D) modified Look-Locker inversion recovery (MOLLI) and spin-echo sequences. In vivo T1 times of the proposed method were compared with 2D multislice MOLLI T1 mapping. Subsequently, the feasibility of high-resolution 3D T1 mapping with spatial resolution of 1.7 × 1.7 × 4 mm(3) was demonstrated.

RESULTS

The proposed method shows good agreement with 2D MOLLI and the spin-echo reference in phantom. No significant difference was found in the in vivo T1 times estimated using the proposed sequence and the 2D MOLLI technique (myocardium, 330 ± 66 ms versus 319 ± 93 ms; blood pools, 211 ± 68 ms versus 210 ± 98 ms). However, improved homogeneity, as measured using standard deviation of the T1 signal, was observed in the 3D T1 maps.

CONCLUSION

The proposed sequence enables high-resolution 3D T1 mapping after contrast injection during free-breathing with volumetric left ventricle coverage.