3D myocardial T1 mapping using saturation recovery

Free-breathing 3D whole-heart joint T1/T2 mapping and water/fat imaging at 0.55 T

PURPOSE

To develop and validate a highly efficient motion compensated free-breathing isotropic resolution 3D whole-heart joint T1/T2 mapping sequence with anatomical water/fat imaging at 0.55 T.

METHODS

The proposed sequence takes advantage of shorter T1 at 0.55 T to acquire three interleaved water/fat volumes with inversion-recovery preparation, no preparation, and T2 preparation, respectively. Image navigators were used to facilitate nonrigid motion-compensated image reconstruction. T1 and T2 maps were jointly calculated by a dictionary matching method. Validations were performed with simulation, phantom, and in vivo experiments on 10 healthy volunteers and 1 patient. The performance of the proposed sequence was compared with conventional 2D mapping sequences including modified Look-Locker inversion recovery and T2-prepared balanced steady-SSFP sequence.

RESULTS

The proposed sequence has a good T1 and T2 encoding sensitivity in simulation, and excellent agreement with spin-echo reference T1 and T2 values was observed in a standardized T1/T2 phantom (R2 = 0.99). In vivo experiments provided good-quality co-registered 3D whole-heart T1 and T2 maps with 2-mm isotropic resolution in a short scan time of about 7 min. For healthy volunteers, left-ventricle T1 mean and SD measured by the proposed sequence were both comparable with those of modified Look-Locker inversion recovery (640 ± 35 vs. 630 ± 25 ms [p = 0.44] and 49.9 ± 9.3 vs. 54.4 ± 20.5 ms [p = 0.42]), whereas left-ventricle T2 mean and SD measured by the proposed sequence were both slightly lower than those of T2-prepared balanced SSFP (53.8 ± 5.5 vs. 58.6 ± 3.3 ms [p < 0.01] and 5.2 ± 0.9 vs. 6.1 ± 0.8 ms [p = 0.03]). Myocardial T1 and T2 in the patient measured by the proposed sequence were in good agreement with conventional 2D sequences and late gadolinium enhancement.

CONCLUSION

The proposed sequence simultaneously acquires 3D whole-heart T1 and T2 mapping with anatomical water/fat imaging at 0.55 T in a fast and efficient 7-min scan. Further investigation in patients with cardiovascular disease is now warranted.

Native myocardial T1 mapping using inversion recovery T1-weighted turbo field echo sequence

This study proposes the use of the inversion recovery T1-weighted turbo field echo (IR-T1TFE) sequence for myocardial T1 mapping and compares the results obtained with those of the modified Look-Locker inversion recovery (MOLLI) method for accuracy, precision, and reproducibility. A phantom containing seven vials with different T1 values was imaged, thereby comparing the T1 measurements between the inversion recovery spin-echo (IR-SE) technique, MOLLI, and the IR-T1TFE. The accuracy, precision, and reproducibility of the T1-mapping sequences were analyzed in a phantom study. Fifteen healthy subjects were recruited for the in vivo comparison of native myocardial T1 mapping using MOLLI and IR-T1TFE sequences. After myocardium segmentation, the T1 value of the entire myocardium was calculated. In the phantom study, excellent accuracy was achieved using IR-T1TFE for all T1 ranges. MOLLI displayed lower accuracy than IR-T1TFE (p =0.016), substantially underestimating T1 at large T1 values (> 1000 ms). In the in vivo study, the first mean myocardial T1 values ± SD using MOLLI and IR-T1TFE were 1306 ± 70 ms and 1484 ± 28 ms, respectively, and the second were 1297 ± 68 ms and 1474 ± 43 ms, respectively. The native myocardial T1 obtained with MOLLI was lower than that of IR-T1TFE (p < 0.001). The reproducibility of native myocardial T1 mapping within the same sequence was not statistically significant (p = 0.11). This study demonstrates the utility and validity of myocardial T1 mapping using IR-T1TFE, which is a common sequence. This method was found to have high accuracy and reproducibility.

Utility of 1.5 Tesla MRI Scanner in the Management of Small Sample Sizes Driven from 3D Breast Cell Culture

The aim of this work was to use and optimize a 1.5 Tesla magnetic resonance imaging (MRI) system for three-dimensional (3D) images of small samples obtained from breast cell cultures in vitro. The basis of this study was to design MRI equipment to enable imaging of MCF-7 breast cancer cell cultures (about 1 million cells) in 1.5 and 2 mL glass tubes and/or bioreactors with an external diameter of less than 20 mm. Additionally, the development of software to calculate longitudinal and transverse relaxation times is described. Imaging tests were performed using a clinical MRI scanner OPTIMA 360 manufactured by GEMS. Due to the size of the tested objects, it was necessary to design additional receiving circuits allowing for the study of MCF-7 cell cultures placed in glass bioreactors. The examined sample's volume did not exceed 2.0 mL nor did the number of cells exceed 1 million. This work also included a modification of the sequence to allow for the analysis of T1 and T2 relaxation times. The analysis was performed using the MATLAB package (produced by MathWorks). The created application is based on medical MR images saved in the DICOM3.0 standard which ensures that the data analyzed are reliable and unchangeable in an unintentional manner that could affect the measurement results. The possibility of using 1.5 T MRI systems for cell culture research providing quantitative information from in vitro studies was realized. The scanning resolution for FOV = 5 cm and the matrix was achieved at a level of resolution of less than 0.1 mm/pixel. Receiving elements were built allowing for the acquisition of data for MRI image reconstruction confirmed by images of a phantom with a known structure and geometry. Magnetic resonance sequences were modified for the saturation recovery (SR) method, the purpose of which was to determine relaxation times. An application in MATLAB was developed that allows for the analysis of T1 and T2 relaxation times. The relaxation times of cell cultures were determined over a 6-week period. In the first week, the T1 time value was 1100 ± 40 ms, which decreased to 673 ± 59 ms by the sixth week. For T2, the results were 171 ± 10 ms and 128 ± 12 ms, respectively.

Motion-compensated T1 mapping in cardiovascular magnetic resonance imaging: a technical review

T 1 mapping is becoming a staple magnetic resonance imaging method for diagnosing myocardial diseases such as ischemic cardiomyopathy, hypertrophic cardiomyopathy, myocarditis, and more. Clinically, most T 1 mapping sequences acquire a single slice at a single cardiac phase across a 10 to 15-heartbeat breath-hold, with one to three slices acquired in total. This leaves opportunities for improving patient comfort and information density by acquiring data across multiple cardiac phases in free-running acquisitions and across multiple respiratory phases in free-breathing acquisitions. Scanning in the presence of cardiac and respiratory motion requires more complex motion characterization and compensation. Most clinical mapping sequences use 2D single-slice acquisitions; however newer techniques allow for motion-compensated reconstructions in three dimensions and beyond. To further address confounding factors and improve measurement accuracy, T 1 maps can be acquired jointly with other quantitative parameters such as T 2 , T 2 ∗ , fat fraction, and more. These multiparametric acquisitions allow for constrained reconstruction approaches that isolate contributions to T 1 from other motion and relaxation mechanisms. In this review, we examine the state of the literature in motion-corrected and motion-resolved T 1 mapping, with potential future directions for further technical development and clinical translation.

Pilot tone-based prospective correction of respiratory motion for free-breathing myocardial T1 mapping

OBJECTIVE

To provide respiratory motion correction for free-breathing myocardial T1 mapping using a pilot tone (PT) and a continuous golden-angle radial acquisition.

MATERIALS AND METHODS

During a 45 s prescan the PT is acquired together with a dynamic sagittal image covering multiple respiratory cycles. From these images, the respiratory heart motion in head-feet and anterior-posterior direction is estimated and two linear models are derived between the PT and heart motion. In the following scan through-plane motion is corrected prospectively with slice tracking based on the PT. In-plane motion is corrected for retrospectively. Our method was evaluated on a motion phantom and 11 healthy subjects.

RESULTS

Non-motion corrected measurements using a moving phantom showed T1 errors of 14 ± 4% (p < 0.05) compared to a reference measurement. The proposed motion correction approach reduced this error to 3 ± 4% (p < 0.05). In vivo the respiratory motion led to an overestimation of T1 values by 26 ± 31% compared to breathhold T1 maps, which was successfully corrected to an average difference of 3 ± 2% (p < 0.05) between our free-breathing approach and breathhold data.

DISCUSSION

Our proposed PT-based motion correction approach allows for T1 mapping during free-breathing with the same accuracy as a corresponding breathhold T1 mapping scan.

3D model-based super-resolution motion-corrected cardiac T1 mapping

Objective. To provide 3D high-resolution cardiac T1 maps using model-based super-resolution reconstruction (SRR).Approach. Due to signal-to-noise ratio limitations and the motion of the heart during imaging, often 2D T1 maps with only low through-plane resolution (i.e. slice thickness of 6-8 mm) can be obtained. Here, a model-based SRR approach is presented, which combines multiple stacks of 2D acquisitions with 6-8 mm slice thickness and generates 3D high-resolution T1 maps with a slice thickness of 1.5-2 mm. Every stack was acquired in a different breath hold (BH) and any misalignment between BH was corrected retrospectively. The novelty of the proposed approach is the BH correction and the application of model-based SRR on cardiac T1 Mapping. The proposed approach was evaluated in numerical simulations and phantom experiments and demonstrated in four healthy subjects.Main results. Alignment of BH states was essential for SRR even in healthy volunteers. In simulations, respiratory motion could be estimated with an RMS error of 0.18 ± 0.28 mm. SRR improved the visualization of small structures. High accuracy and precision (average standard deviation of 69.62 ms) of the T1 values was ensured by SRR while the detectability of small structures increased by 40%.Significance. The proposed SRR approach provided T1 maps with high in-plane and high through-plane resolution (1.3 × 1.3 × 1.5-2 mm³). The approach led to improvements in the visualization of small structures and precise T1 values.

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 3D cardiac T1 mapping with transmit B1 correction at 3T

PURPOSE

To develop a cardiac T₁ mapping method for free-breathing 3D T₁ mapping of the whole heart at 3 T with transmit B₁ ( B 1 + ) correction.

METHODS

A free-breathing, electrocardiogram-gated inversion-recovery sequence with spoiled gradient-echo readout was developed and optimized for cardiac T₁ mapping at 3 T. High-frame-rate dynamic images were reconstructed from sparse (k,t)-space data acquired along a stack-of-stars trajectory using a subspace-based method for accelerated imaging. Joint T₁ and flip-angle estimation was performed in T₁ mapping to improve its robustness to B 1 + inhomogeneity. Subject-specific timing of data acquisition was used in the estimation to account for natural heart-rate variations during the imaging experiment.

RESULTS

Simulations showed that accuracy and precision of T₁ mapping can be improved with joint T₁ and flip-angle estimation and optimized electrocardiogram-gated spoiled gradient echo-based inversion-recovery acquisition scheme. The phantom study showed good agreement between the T₁ maps from the proposed method and the reference method. Three-dimensional cardiac T₁ maps (40 slices) were obtained at a 1.9-mm in-plane and 4.5-mm through-plane spatial resolution from healthy subjects (n = 6) with an average imaging time of 14.2 ± 1.6 minutes (heartbeat rate: 64.2 ± 7.1 bpm), showing myocardial T₁ values comparable to those obtained from modified Look-Locker inversion recovery. The proposed method generated B 1 + maps with spatially smooth variation showing 21%-32% and 11%-15% variations across the septal-lateral and inferior-anterior regions of the myocardium in the left ventricle.

CONCLUSION

The proposed method allows free-breathing 3D T₁ mapping of the whole heart with transmit B₁ correction in a practical imaging time.

Recent advances in multimodality imaging of the tricuspid valve

INTRODUCTION

The tricuspid valve (TV) and the right heart chambers are complex three-dimensional structures that are difficult to assess using tomographic imaging techniques. The progressive aging of the general population and the advancements in treating left-sided heart diseases by transcatheter procedures have contributed to the tricuspid regurgitation (TR) becoming a major public health problem associated with progression to refractory heart failure and poor outcome. Recent advances in multimodality cardiac imaging allow a better understanding of the pathophysiology of TR that may translate in better management of patients.

AREAS COVERED

Three-dimensional echocardiography, cardiac magnetic resonance, and computed tomography provide complementary information to i. assess the TV complex; ii. identify the etiology and the mechanisms of TR; iii. evaluate its severity and hemodynamic consequences; iv. explore the remodeling of the right heart chambers; and v. properly plan, guide, and monitor the transcatheter interventions aimed to reduce the severity of TR.

EXPERT OPINION

We need thorough understanding of both the TV and the right heart chamber geometry and function to understand the pathophysiology of TR. The integrated use of multimodality cardiac imaging is pivotal to assess patients with TR and to identify tailored and timely treatment of TR in properly selected patients.

Generalized low-rank nonrigid motion-corrected reconstruction for MR fingerprinting

PURPOSE

Develop a novel low-rank motion-corrected (LRMC) reconstruction for nonrigid motion-corrected MR fingerprinting (MRF).

METHODS

Generalized motion-corrected (MC) reconstructions have been developed for steady-state imaging. Here we extend this framework to enable nonrigid MC for transient imaging applications with varying contrast, such as MRF. This is achieved by integrating low-rank dictionary-based compression into the generalized MC model to reconstruct MC singular images, reducing motion artifacts in the resulting parametric maps. The proposed LRMC reconstruction was applied for cardiac motion correction in 2D myocardial MRF (T₁ and T₂ ) with extended cardiac acquisition window (~450 ms) and for respiratory MC in free-breathing 3D myocardial and 3D liver MRF. Experiments were performed in phantom and 22 healthy subjects. The proposed approach was compared with reference spin echo (phantom) and with 2D electrocardiogram-triggered/breath-hold MOLLI and T₂ gradient-and-spin echo conventional maps (in vivo 2D and 3D myocardial MRF).

RESULTS

Phantom results were in general agreement with reference spin-echo measurements, presenting relative errors of approximately 5.4% and 5.5% for T₁ and short T₂ (<100 ms), respectively. The proposed LRMC MRF reduced residual blurring artifacts with respect to no MC for cardiac or respiratory motion in all cases (2D and 3D myocardial, 3D abdominal). In 2D myocardial MRF, left-ventricle T₁ values were 1150 ± 41 ms for LRMC MRF and 1010 ± 56 ms for MOLLI; T₂ values were 43.8 ± 2.3 ms for LRMC MRF and 49.5 ± 4.5 ms for T₂ gradient and spin echo. Corresponding measurements for 3D myocardial MRF were 1085 ± 30 ms and 1062 ± 29 ms for T₁ , and 43.5 ± 1.9 ms and 51.7 ± 1.7 ms for T₂ . For 3D liver, LRMC MRF measured liver T₁ at 565 ± 44 ms and liver T₂ at 35.4 ± 2.4 ms.

CONCLUSION

The proposed LRMC reconstruction enabled generalized (nonrigid) MC for 2D and 3D MRF, both for cardiac and respiratory motion. The proposed approach reduced motion artifacts in the MRF maps with respect to no motion compensation and achieved good agreement with reference measurements.

Accelerating whole-heart 3D T2 mapping: Impact of undersampling strategies and reconstruction techniques

PURPOSE

We aim to determine an advantageous approach for the acceleration of high spatial resolution 3D cardiac T2 relaxometry data by comparing the performance of different undersampling patterns and reconstruction methods over a range of acceleration rates.

METHODS

Multi-volume 3D high-resolution cardiac images were acquired fully and undersampled retrospectively using 1) optimal CAIPIRINHA and 2) a variable density random (VDR) sampling. Data were reconstructed using 1) multi-volume sensitivity encoding (SENSE), 2) joint-sparsity SENSE and 3) model-based SENSE. Four metrics were calculated on 3 naïve swine and 8 normal human subjects over a whole left-ventricular region of interest: root-mean-square error (RMSE) of image signal intensity, RMSE of T2, the bias of mean T2, and standard deviation (SD) of T2. Fully sampled data and volume-by-volume SENSE with standard equally spaced undersampling were used as references. The Jaccard index calculated from one swine with acute myocardial infarction (MI) was used to demonstrate preservation of segmentation of edematous tissues with elevated T2.

RESULTS

In naïve swine and normal human subjects, all methods had similar performance when the net reduction factor (Rnet) <2.5. VDR sampling with model-based SENSE showed the lowest RMSEs (10.5%-14.2%) and SDs (+1.7-2.4 ms) of T2 when Rnet>2.5, while VDR sampling with the joint-sparsity SENSE had the lowest bias of mean T2 (0.0-1.1ms) when Rnet>3. The RMSEs of parametric T2 values (9.2%-24.6%) were larger than for image signal intensities (5.2%-18.4%). In the swine with MI, VDR sampling with either joint-sparsity or model-based SENSE showed consistently higher Jaccard index for all Rnet (0.71-0.50) than volume-by-volume SENSE (0.68-0.30).

CONCLUSIONS

Retrospective exploration of undersampling and reconstruction in 3D whole-heart T2 parametric mapping revealed that maps were more sensitive to undersampling than images, presenting a more stringent limiting factor on Rnet. The combination of VDR sampling patterns with model-based or joint-sparsity SENSE reconstructions were more robust for Rnet>3.

Extracellular volume-guided late gadolinium enhancement analysis for non-ischemic cardiomyopathy: The Women's Interagency HIV Study

BACKGROUND

Quantification of non-ischemic myocardial scar remains a challenge due to the patchy diffuse nature of fibrosis. Extracellular volume (ECV) to guide late gadolinium enhancement (LGE) analysis may achieve a robust scar assessment.

METHODS

Three cohorts of 80 non-ischemic-training, 20 non-ischemic-validation, and 10 ischemic-validation were prospectively enrolled and underwent 3.0 Tesla cardiac MRI. An ECV cutoff to differentiate LGE scar from non-scar was identified in the training cohort from the receiver-operating characteristic curve analysis, by comparing the ECV value against the visually-determined presence/absence of the LGE scar at the highest signal intensity (SI) area of the mid-left ventricle (LV) LGE. Based on the ECV cutoff, an LGE semi-automatic threshold of n-times of standard-deviation (n-SD) above the remote-myocardium SI was optimized in the individual cases ensuring correspondence between LGE and ECV images. The inter-method agreement of scar amount in comparison with manual (for non-ischemic) or full-width half-maximum (FWHM, for ischemic) was assessed. Intra- and inter-observer reproducibility were investigated in a randomly chosen subset of 40 non-ischemic and 10 ischemic cases.

RESULTS

The non-ischemic groups were all female with the HIV positive rate of 73.8% (training) and 80% (validation). The ischemic group was all male with reduced LV function. An ECV cutoff of 31.5% achieved optimum performance (sensitivity: 90%, specificity: 86.7% in training; sensitivity: 100%, specificity: 81.8% in validation dataset). The identified n-SD threshold varied widely (range 3 SD-18 SD), and was independent of scar amount (β = -0.01, p = 0.92). In the non-ischemic cohorts, results suggested that the manual LGE assessment overestimated scar (%) in comparison to ECV-guided analysis [training: 4.5 (3.2-6.4) vs. 0.92 (0.1-2.1); validation: 2.5 (1.2-3.7) vs. 0.2 (0-1.6); P < 0.01 for both]. Intra- and inter-observer analyses of global scar (%) showed higher reproducibility in ECV-guided than manual analysis with CCC = 0.94 and 0.78 versus CCC = 0.86 and 0.73, respectively (P < 0.01 for all). In ischemic validation, the ECV-guided LGE analysis showed a comparable scar amount and reproducibility with the FWHM.

CONCLUSIONS

ECV-guided LGE analysis is a robust scar quantification method for a non-ischemic cohort. Trial registration ClinicalTrials.gov; NCT00000797, retrospectively-registered 2 November 1999; NCT02501811, registered 15 July 2015.

Single breath-hold saturation recovery 3D cardiac T1 mapping via compressed SENSE at 3T

Objectives

To propose and validate a novel imaging sequence that uses a single breath-hold whole-heart 3D T1 saturation recovery compressed SENSE rapid acquisition (SACORA) at 3T.

Methods

The proposed sequence combines flexible saturation time sampling, compressed SENSE, and sharing of saturation pulses between two readouts acquired at different RR intervals. The sequence was compared with a 3D saturation recovery single-shot acquisition (SASHA) implementation with phantom and in vivo experiments (pre and post contrast; 7 pigs) and was validated against the reference inversion recovery spin echo (IR-SE) sequence in phantom experiments.

Results

Phantom experiments showed that the T1 maps acquired by 3D SACORA and 3D SASHA agree well with IR-SE. In vivo experiments showed that the pre-contrast and post-contrast T1 maps acquired by 3D SACORA are comparable to the corresponding 3D SASHA maps, despite the shorter acquisition time (15s vs. 188s, for a heart rate of 60 bpm). Mean septal pre-contrast T1 was 1453 ± 44 ms with 3D SACORA and 1460 ± 60 ms with 3D SASHA. Mean septal post-contrast T1 was 824 ± 66 ms and 824 ± 60 ms.

Conclusion

3D SACORA acquires 3D T1 maps in 15 heart beats (heart rate, 60 bpm) at 3T. In addition to its short acquisition time, the sequence achieves good T1 estimation precision and accuracy.

Electronic supplementary material

The online version of this article (10.1007/s10334-020-00848-2) contains supplementary material, which is available to authorized users.

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.

3D whole-heart isotropic-resolution motion-compensated joint T1 /T2 mapping and water/fat imaging

PURPOSE

To develop a free-breathing isotropic-resolution whole-heart joint T₁ and T₂ mapping sequence with Dixon-encoding that provides coregistered 3D T₁ and T₂ maps and complementary 3D anatomical water and fat images in a single ~9 min scan.

METHODS

Four interleaved dual-echo Dixon gradient echo volumes are acquired with a variable density Cartesian trajectory and different preparation pulses: 1) inversion recovery-preparation, 2) and 3) no preparations, and 4) T₂ preparation. Image navigators are acquired to correct each echo for 2D translational respiratory motion; the 8 echoes are jointly reconstructed with a low-rank patch-based reconstruction. A water/fat separation algorithm is used to obtain water and fat images for each acquired volume. T₁ and T₂ maps are generated by matching the signal evolution of the water images to a simulated dictionary. Complementary bright-blood and fat volumes for anatomical visualization are obtained from the T₂ -prepared dataset. The proposed sequence was tested in phantom experiments and 10 healthy subjects and compared to standard 2D MOLLI T₁ mapping, 2D balance steady-state free precession T₂ mapping, and 3D T₂ -prepared Dixon coronary MR angiography.

RESULTS

High linear correlation was found between T₁ and T₂ quantification with the proposed approach and phantom spin echo measurements (y = 1.1 × -11.68, R² = 0.98; and y = 0.85 × +5.7, R² = 0.99). Mean myocardial values of T₁ /T₂ = 1116 ± 30.5 ms/45.1 ± 2.38 ms were measured in vivo. Biases of T₁ /T₂ = 101.8 ms/-0.77 ms were obtained compared to standard 2D techniques.

CONCLUSION

The proposed joint T₁ /T₂ sequence permitted the acquisition of motion-compensated isotropic-resolution 3D T₁ and T₂ maps and complementary coronary MR angiography and fat volumes, showing promising results in terms of T₁ and T₂ quantification and visualization of cardiac anatomy and pericardial fat.

Faster 3D saturation-recovery based myocardial T1 mapping using a reduced number of saturation points and denoising

Purpose

To accelerate the acquisition of free-breathing 3D saturation-recovery-based (SASHA) myocardial T1 mapping by acquiring fewer saturation points in combination with a post-processing 3D denoising technique to maintain high accuracy and precision.

Methods

3D SASHA T1 mapping acquires nine T1-weighted images along the saturation recovery curve, resulting in long acquisition times. In this work, we propose to accelerate conventional cardiac T1 mapping by reducing the number of saturation points. High T1 accuracy and low standard deviation (as a surrogate for precision) is maintained by applying a 3D denoising technique to the T1-weighted images prior to pixel-wise T1 fitting. The proposed approach was evaluated on a T1 phantom and 20 healthy subjects, by varying the number of T1-weighted images acquired between three and nine, both prospectively and retrospectively. Following the results from the healthy subjects, three patients with suspected cardiovascular disease were acquired using five T1-weighted images. T1 accuracy and precision was determined for all the acquisitions before and after denoising.

Results

In the T1 phantom, no statistical difference was found in terms of accuracy and precision for the different number of T1-weighted images before or after denoising (P = 0.99 and P = 0.99 for accuracy, P = 0.64 and P = 0.42 for precision, respectively). In vivo, both prospectively and retrospectively, the precision improved considerably with the number of T1-weighted images employed before denoising (P<0.05) but was independent on the number of T1-weighted images after denoising.

Conclusion

We demonstrate the feasibility of accelerating 3D SASHA T1 mapping by reducing the number of acquired T1-weighted images in combination with an efficient 3D denoising, without affecting accuracy and precision of T1 values.

Free-Breathing Three-Dimensional T1 Mapping of the Heart Using Subspace-Based Data Acquisition and Image Reconstruction

Mapping the longitudinal relaxation time constant (T1) of the myocardium using Magnetic Resonance Imaging (MRI) is an emerging technique for quantitative assessment of the morphology and viability of the myocardium. However, three-dimensional (3D) T1 mapping of the heart is challenging due to the high dimensionality of the signal and the presence of cardiac and respiratory motions. We propose a subspace-based method for free-breathing 3D T1 mapping of the heart without respiratory gating. The image function is represented as a high-order partially separable (PS) function to explore the inherent spatiotemporal correlations of the underlying signal. A special data acquisition scheme enabled by the high-order PS model is used for sparse sampling of the (k,t)-space, where complementary sparse datasets are acquired, each covering only a small portion of the (k,t)-space to characterize a single subspace (spatial or temporal). High-resolution dynamic MR images are reconstructed from the highly undersampled (k,t)-space using low-rank tensor and sparsity constraints. We demonstrate the feasibility of our proposed method using in vivo data obtained from healthy subjects on a 3T MR scanner. The proposed method can enable new clinical applications of T1 mapping in cardiac MR.

Multimodality Imaging of the Tricuspid Valve and Right Heart Anatomy

The characterization of tricuspid valve and right-heart anatomy has been gaining significant interest in the setting of new percutaneous transcatheter interventions for tricuspid regurgitation. Multimodality cardiac imaging provides a wealth of information about the anatomy and function of the tricuspid valve apparatus, right ventricle, and right atrium, which is pivotal for diagnosis and prognosis and for planning of percutaneous interventions. The present review describes the role of echocardiography, cardiac magnetic resonance, and multidetector row cardiac computed tomography for right heart and tricuspid valve assessment.

Molecular and Nonmolecular Magnetic Resonance Coronary and Carotid Imaging

Atherosclerosis is the leading cause of cardiovascular morbidity and mortality. Over the past 2 decades, increasing research attention is converging on the early detection and monitoring of atherosclerotic plaque. Among several invasive and noninvasive imaging modalities, magnetic resonance imaging (MRI) is emerging as a promising option. Advantages include its versatility, excellent soft tissue contrast for plaque characterization and lack of ionizing radiation. In this review, we will explore the recent advances in multicontrast and multiparametric imaging sequences that are bringing the aspiration of simultaneous arterial lumen, vessel wall, and plaque characterization closer to clinical feasibility. We also discuss the latest advances in molecular magnetic resonance and multimodal atherosclerosis imaging.

Myocardial T1 and ECV Measurement: Underlying Concepts and Technical Considerations

Myocardial native T1 and extracellular volume fraction (ECV) mapping have emerged as cardiac magnetic resonance biomarkers providing unique insight into cardiac pathophysiology. Single breath-hold acquisition techniques, available on clinical scanners across multiple vendor platforms, have made clinical T1 and ECV mapping a reality. Although the relationship between changes in native T1 and alterations in cardiac microstructure is complex, an understanding of how edema, blood volume, myocyte and interstitial expansion, lipids, and paramagnetic substances affect T1 and ECV can provide insight into how and why these parameters change in various cardiac pathologies. The goals of this state-of-the-art review will be to review factors influencing native T1 and ECV, to describe how native T1 and ECV are measured, to discuss potential challenges and pitfalls in clinical practice, and to describe new T1 mapping techniques on the horizon.

3D Whole-heart free-breathing qBOOST-T2 mapping

PURPOSE

To develop an accelerated motion corrected 3D whole-heart imaging approach (qBOOST-T2) for simultaneous high-resolution bright- and black-blood cardiac MR imaging and quantitative myocardial T2 characterization.

METHODS

Three undersampled interleaved balanced steady-state free precession cardiac MR volumes were acquired with a variable density Cartesian trajectory and different magnetization preparations: (1) T2-prepared inversion recovery (T2prep-IR), (2) T2-preparation, and (3) no preparation. Image navigators were acquired prior the acquisition to correct for 2D translational respiratory motion. Each 3D volume was reconstructed with a low-rank patch-based reconstruction. The T2prep-IR volume provides bright-blood anatomy visualization, the black-blood volume is obtained by means of phase sensitive reconstruction between first and third datasets, and T2 maps are generated by matching the signal evolution to a simulated dictionary. The proposed sequence has been evaluated in simulations, phantom experiments, 11 healthy subjects and compared with 3D bright-blood cardiac MR and standard 2D breath-hold balanced steady-state free precession T2 mapping. The feasibility of the proposed approach was tested on 4 patients with suspected cardiovascular disease.

RESULTS

High linear correlation (y = 1.09 × -0.83, R² = 0.99) was found between the proposed qBOOST-T2 and T2 spin echo measurements in phantom experiment. Good image quality was observed in vivo with the proposed 4x undersampled qBOOST-T2. Mean T2 values of 53.1 ± 2.1 ms and 55.8 ± 2.7 ms were measured in vivo for 2D balanced steady-state free precession T2 mapping and qBOOST-T2, respectively, with linear correlation of y = 1.02x+1.46 (R² = 0.61) and T2 bias = 2.7 ms.

CONCLUSION

The proposed qBOOST-T2 sequence allows the acquisition of 3D high-resolution co-registered bright- and black-blood volumes and T2 maps in a single scan of ~11 min, showing promising results in terms of T2 quantification.

Free-running 3D whole heart myocardial T1 mapping with isotropic spatial resolution

PURPOSE

To develop a free-running (free-breathing, retrospective cardiac gating) 3D myocardial T1 mapping with isotropic spatial resolution.

METHODS

The free-running sequence is inversion recovery (IR)-prepared followed by continuous 3D golden angle radial data acquisition. 1D respiratory motion signal is extracted from the k-space center of all spokes and used to bin the k-space data into different respiratory states, enabling estimation and correction of 3D translational respiratory motion, whereas cardiac motion is recorded using electrocardiography and synchronized with data acquisition. 3D translational respiratory motion compensated T1 maps at diastole and systole were generated with 1.5 mm isotropic spatial resolution with low-rank inversion and high-dimensionality patch-based undersampled reconstruction. The technique was validated against conventional methods in phantom and 9 healthy subjects.

RESULTS

Phantom results demonstrated good agreement (R2 = 0.99) of T1 estimation with reference method. Homogeneous systolic and diastolic 3D T1 maps were reconstructed from the proposed technique. Diastolic septal T1 estimated with the proposed method (1140 ± 36 ms) was comparable to the saturation recovery single-shot acquisition (SASHA) sequence (1153 ± 49 ms), but was higher than the modified Look-Locker inversion recovery (MOLLI) sequence (1037 ± 33 ms). Precision of the proposed method (42 ± 8 ms) was comparable to MOLLI (41 ± 7 ms) and improved with respect to SASHA (87 ± 19 ms).

CONCLUSIONS

The proposed free-running whole heart T1 mapping method allows for reconstruction of isotropic resolution 3D T1 maps at different cardiac phases, serving as a promising tool for whole heart myocardial tissue characterization.

Whole-heart T1 mapping using a 2D fat image navigator for respiratory motion compensation

Purpose

To combine a 3D saturation‐recovery‐based myocardial T₁ mapping (3D SASHA) sequence with a 2D image navigator with fat excitation (fat‐iNAV) to allow 3D T₁ maps with 100% respiratory scan efficiency and predictable scan time.

Methods

Data from T₁ phantom and 10 subjects were acquired at 1.5T. For respiratory motion compensation, a 2D fat‐iNAV was acquired before each 3D SASHA k‐space segment to correct for 2D translational motion in a beat‐to‐beat fashion. The effect of the fat‐iNAV on the 3D SASHA T1 estimation was evaluated on the T₁ phantom. For 3 representative subjects, the proposed free‐breathing 3D SASHA with fat‐iNAV was compared to the original implementation with the diaphragmatic navigator. The 3D SASHA with fat‐iNAV was compared to the breath‐hold 2D SASHA sequence in terms of accuracy and precision.

Results

In the phantom study, the Bland‐Altman plot shows that the 2D fat‐iNAVs does not affect the T₁ quantification of the 3D SASHA acquisition (0 ± 12.5 ms). For the in vivo study, the 2D fat‐iNAV permits to estimate the respiratory motion of the heart, while allowing for 100% scan efficiency, improving the precision of the T₁ measurement compared to non‐motion‐corrected 3D SASHA. However, the image quality achieved with the proposed 3D SASHA with fat‐iNAV is lower compared to the original implementation, with reduced delineation of the myocardial borders and papillary muscles.

Conclusions

We demonstrate the feasibility to combine the 3D SASHA T₁ mapping imaging sequence with a 2D fat‐iNAV for respiratory motion compensation, allowing 100% respiratory scan efficiency and predictable 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.

Cardiac T1 mapping: Techniques and applications

A key advantage of cardiac magnetic resonance (CMR) imaging over other cardiac imaging modalities is the ability to perform detailed tissue characterization. CMR techniques continue to evolve, with advanced imaging sequences being developed to provide a reproducible, quantitative method of tissue interrogation. The T₁ mapping technique, a pixel-by-pixel method of quantifying T₁ relaxation time of soft tissues, has been shown to be promising for characterization of diseased myocardium in a wide variety of cardiomyopathies. In this review, we describe the basic principles and common techniques for T₁ mapping and its use for native T₁ , postcontrast T₁ , and extracellular volume mapping. We will review a wide range of clinical applications of the technique that can be used for identification and quantification of myocardial edema, fibrosis, and infiltrative diseases with illustrative clinical examples. In addition, we will explore the current limitations of the technique and describe some areas of ongoing development. Level of Evidence: 5 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2020;51:1336-1356.

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.

Cardiac magnetic resonance T1 mapping. Part 1: Aspects of acquisition and evaluation

While an enormous number of studies have documented pathological alterations of the myocardial native longitudinal relaxation time (T1) and the fraction of the extracellular myocardial volume (ECV), it has also become clear that continuously evolving T1 mapping sequence, acquisition and evaluation techniques have a substantial impact on quantitative results, making the translation of reported findings into routine clinical use particularly challenging. To provide a basis for the discussion of pathological myocardial T1 and ECV alterations, the present review aims to summarize the methodological aspects of myocardial T1 mapping along with technical and physiological factors influencing results and normal ranges of myocardial native T1 and ECV reported across studies.

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.

3D SASHA myocardial T1 mapping with high accuracy and improved precision

Purpose

To improve the precision of a free-breathing 3D saturation-recovery-based myocardial T1 mapping sequence using a post-processing 3D denoising technique.

Methods

A T1 phantom and 15 healthy subjects were scanned on a 1.5 T MRI scanner using 3D saturation-recovery single-shot acquisition (SASHA) for myocardial T1 mapping. A 3D denoising technique was applied to the native T1-weighted images before pixel-wise T1 fitting. The denoising technique imposes edge-preserving regularity and exploits the co-occurrence of 3D spatial gradients in the native T1-weighted images by incorporating a multi-contrast Beltrami regularization. Additionally, 2D modified Look-Locker inversion recovery (MOLLI) acquisitions were performed for comparison purposes. Accuracy and precision were measured in the myocardial septum of 2D MOLLI and 3D SASHA T1 maps and then compared. Furthermore, the accuracy and precision of the proposed approach were evaluated in a standardized phantom in comparison to an inversion-recovery spin-echo sequence (IRSE).

Results

For the phantom study, Bland–Altman plots showed good agreement in terms of accuracy between IRSE and 3D SASHA, both on non-denoised and denoised T1 maps (mean difference −1.4 ± 18.9 ms and −4.4 ± 21.2 ms, respectively), while 2D MOLLI generally underestimated the T1 values (69.4 ± 48.4 ms). For the in vivo study, there was a statistical difference between the precision measured on 2D MOLLI and on non-denoised 3D SASHA T1 maps (P = 0.005), while there was no statistical difference after denoising (P = 0.95).

Conclusion

The precision of 3D SASHA myocardial T1 mapping was substantially improved using a 3D Beltrami regularization based denoising technique and was similar to that of 2D MOLLI T1 mapping, while preserving the higher accuracy and whole-heart coverage of 3D SASHA.

Electronic supplementary material

The online version of this article (10.1007/s10334-018-0703-y) contains supplementary material, which is available to authorized users.

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.