Temporally resolved parametric assessment of Z-magnetization recovery (TOPAZ): Dynamic myocardial T1 mapping using a cine steady-state look-locker approach

Deep image prior cine MR fingerprinting with B1 + spin history correction

PURPOSE

To develop a deep image prior (DIP) reconstruction for B1 + -corrected 2D cine MR fingerprinting (MRF).

METHODS

The proposed method combines low-rank (LR) modeling with a DIP to generate cardiac phase-resolved parameter maps without motion correction, employing self-supervised training to enforce consistency with undersampled spiral k-space data. Two implementations were tested: one approach (DIP) for cine T1 , T2 , and M0 mapping, and a second approach (DIP with effective B1 + estimation [DIP-B1]) that also generated an effective B1 + map to correct for errors due to RF transmit inhomogeneities, through-plane motion, and blood flow. Cine MRF data were acquired in 14 healthy subjects and four reconstructions were compared: LR, low-rank motion-corrected (LRMC), DIP, and DIP-B1. Results were compared to diastolic ECG-triggered MRF, MOLLI, and T2 -prep bSSFP. Additionally, bright-blood and dark-blood images calculated from cine MRF maps were used to quantify ventricular function and compared to reference cine measurements.

RESULTS

DIP and DIP-B1 outperformed other cine MRF reconstructions with improved noise suppression and delineation of high-resolution details. Within-segment variability in the myocardium (reported as the coefficient of variation for T1 /T2 ) was lowest for DIP-B1 (2.3/8.3%) followed by DIP (2.7/8.7%), LRMC (3.5/10.5%), and LR (15.3/39.6%). Spatial homogeneity improved with DIP-B1 having the lowest intersegment variability (2.6/4.1%). The mean bias in ejection fraction was -1.1% compared to reference cine scans.

CONCLUSION

A DIP reconstruction for 2D cine MRF enabled cardiac phase-resolved mapping of T1 , T2 , M0 , and the effective B1 + with improved noise suppression and precision compared to LR and LRMC.

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.

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.

Free-running 3D whole-heart T1 and T2 mapping and cine MRI using low-rank reconstruction with non-rigid cardiac motion correction

PURPOSE

To introduce non-rigid cardiac motion correction into a novel free-running framework for the simultaneous acquisition of 3D whole-heart myocardial <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> maps and cine images, enabling a <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>∼</mml:mo></mml:mrow> <mml:annotation>$$ \sim $$</mml:annotation></mml:semantics> </mml:math> 3-min scan.

METHODS

Data were acquired using a free-running 3D golden-angle radial readout interleaved with inversion recovery and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> -preparation pulses. After correction for translational respiratory motion, non-rigid cardiac-motion-corrected reconstruction with dictionary-based low-rank compression and patch-based regularization enabled 3D whole-heart <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> mapping at any given cardiac phase as well as whole-heart cardiac cine imaging. The framework was validated and compared with established methods in 11 healthy subjects.

RESULTS

Good quality 3D <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> maps and cine images were reconstructed for all subjects. Septal <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> values using the proposed approach ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mn>1200</mml:mn> <mml:mo>±</mml:mo> <mml:mn>50</mml:mn></mml:mrow> <mml:annotation>$$ 1200\pm 50 $$</mml:annotation></mml:semantics> </mml:math> ms) were higher than those from a 2D MOLLI sequence ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mn>1063</mml:mn> <mml:mo>±</mml:mo> <mml:mn>33</mml:mn></mml:mrow> <mml:annotation>$$ 1063\pm 33 $$</mml:annotation></mml:semantics> </mml:math> ms), which is known to underestimate <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> , while <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> values from the proposed approach ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mn>51</mml:mn> <mml:mo>±</mml:mo> <mml:mn>4</mml:mn></mml:mrow> <mml:annotation>$$ 51\pm 4 $$</mml:annotation></mml:semantics> </mml:math> ms) were in good agreement with those from a 2D GraSE sequence ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mn>51</mml:mn> <mml:mo>±</mml:mo> <mml:mn>2</mml:mn></mml:mrow> <mml:annotation>$$ 51\pm 2 $$</mml:annotation></mml:semantics> </mml:math> ms).

CONCLUSION

The proposed technique provides 3D <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>1</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_1 $$</mml:annotation></mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow> <mml:mrow><mml:mn>2</mml:mn></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$$ {T}_2 $$</mml:annotation></mml:semantics> </mml:math> maps and cine images with isotropic spatial resolution in a single <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>∼</mml:mo></mml:mrow> <mml:annotation>$$ \sim $$</mml:annotation></mml:semantics> </mml:math> 3.3-min scan.

Simultaneous multi-parametric acquisition and reconstruction techniques in cardiac magnetic resonance imaging: Basic concepts and status of clinical development

Simultaneous multi-parametric acquisition and reconstruction techniques (SMART) are gaining attention for their potential to overcome some of cardiovascular magnetic resonance imaging's (CMR) clinical limitations. The major advantages of SMART lie within their ability to simultaneously capture multiple "features" such as cardiac motion, respiratory motion, T1/T2 relaxation. This review aims to summarize the overarching theory of SMART, describing key concepts that many of these techniques share to produce co-registered, high quality CMR images in less time and with less requirements for specialized personnel. Further, this review provides an overview of the recent developments in the field of SMART by describing how they work, the parameters they can acquire, their status of clinical testing and validation, and by providing examples for how their use can improve the current state of clinical CMR workflows. Many of the SMART are in early phases of development and testing, thus larger scale, controlled trials are needed to evaluate their use in clinical setting and with different cardiac pathologies.

Cardiac phase-resolved late gadolinium enhancement imaging

Late gadolinium enhancement (LGE) with cardiac magnetic resonance (CMR) imaging is the clinical reference for assessment of myocardial scar and focal fibrosis. However, current LGE techniques are confined to imaging of a single cardiac phase, which hampers assessment of scar motility and does not allow cross-comparison between multiple phases. In this work, we investigate a three step approach to obtain cardiac phase-resolved LGE images: (1) Acquisition of cardiac phase-resolved imaging data with varying T₁ weighting. (2) Generation of semi-quantitative T1* maps for each cardiac phase. (3) Synthetization of LGE contrast to obtain functional LGE images. The proposed method is evaluated in phantom imaging, six healthy subjects at 3T and 20 patients at 1.5T. Phantom imaging at 3T demonstrates consistent contrast throughout the cardiac cycle with a coefficient of variation of 2.55 ± 0.42%. In-vivo results show reliable LGE contrast with thorough suppression of the myocardial tissue is healthy subjects. The contrast between blood and myocardium showed moderate variation throughout the cardiac cycle in healthy subjects (coefficient of variation 18.2 ± 3.51%). Images were acquired at 40–60 ms and 80 ms temporal resolution, at 3T and 1.5, respectively. Functional LGE images acquired in patients with myocardial scar visualized scar tissue throughout the cardiac cycle, albeit at noticeably lower imaging resolution and noise resilience than the reference technique. The proposed technique bears the promise of integrating the advantages of phase-resolved CMR with LGE imaging, but further improvements in the acquisition quality are warranted for clinical use.

Quantification of changes in myocardial T1 * values with exercise cardiac MRI using a free-breathing non-electrocardiograph radial imaging

PURPOSE

To develop and evaluate a free breathing non-electrocardiograph (ECG) myocardial T₁ * mapping sequence using radial imaging to quantify the changes in myocardial T₁ * between rest and exercise (T₁ *_reactivity ) in exercise cardiac MRI (Ex-CMR).

METHODS

A free-running T₁ * sequence was developed using a saturation pulse followed by three Look-Locker inversion-recovery experiments. Each Look-Locker continuously acquired data as radial trajectory using a low flip-angle spoiled gradient-echo readout. Self-navigation was performed with a temporal resolution of ∼100 ms for retrospectively extracting respiratory motion. The mid-diastole phase for every cardiac cycle was retrospectively detected on the recorded electrocardiogram signal using an empirical model. Multiple measurements were performed to obtain mean value to reduce effects from the free-breathing acquisition. Finally, data acquired at both mid-diastole and end-expiration are picked and reconstructed by a low-rank plus sparsity constraint algorithm. The performance of this sequence was evaluated by simulations, phantoms, and in vivo studies at rest and after physiological exercise.

RESULTS

Numerical simulation demonstrated that changes in T₁ * are related to the changes in T₁ ; however, other factors such as breathing motion could influence T₁ * measurements. Phantom T₁ * values measured using free-running T₁ * mapping sequence had good correlation with spin-echo T₁ values and was insensitive to heart rate. In the Ex-CMR study, the measured T₁ * reactivity was 10% immediately after exercise and declined over time.

CONCLUSION

The free-running T₁ * mapping sequence allows free-breathing non-ECG quantification of changes in myocardial T₁ * with physiological exercise. Although, absolute myocardial T₁ * value is sensitive to various confounders such as B₁ and B₀ inhomogeneity, quantification of its change may be useful in revealing myocardial tissue properties with exercise.

Dual flip-angle IR-FLASH with spin history mapping for B1+ corrected T1 mapping: Application to T1 cardiovascular magnetic resonance multitasking

PURPOSE

To develop a single-scan method for B 1 + -corrected T1 mapping and apply it for free-breathing (FB) cardiac MR multitasking without electrocardiogram (ECG) triggering.

METHODS

One dual flip-angle (2FA) inversion recovery (IR)-FLASH scan provides two observations of T 1 ∗ (apparent T1 ) corresponding to two distinct combinations of the nominal FA α and B 1 + . Spatiotemporally coregistered T1 and B 1 + spin history maps are obtained by fitting the 2FA signal model. T1 estimate accuracy and repeatability for single flip-angle (1FA) and 2FA IR-FLASH sequence MR multitasking were evaluated at 3T. A T1 phantom was first imaged on the scanner table, then on two human subjects' thoraxes in both breath-hold (BH) and FB conditions. IR-turbo spin echo (IR-TSE) static phantom T1 measurements served as reference. In 10 healthy subjects, myocardial T1 was evaluated with ECG-free, FB multitasking sequences alongside ECG-triggered BH MOLLI.

RESULTS

For phantom-on-table T1 estimates, 2FA agreed better with IR-TSE (intraclass correlation coefficient [ICC] = 0.996, mean error ± SD = -1.6% ± 1.9%) than did 1FA (ICC = 0.922; mean error ± SD = -4.3% ± 12%). For phantom-on-thorax, 2FA was more repeatable and robust to respiration than 1FA (coefficient of variation [CoV] = 1.2% 2FA, = 11.3% 1FA). In vivo, in intrasession T1 repeatability, 2FA (septal CoV = 2.4%, six-segment CoV = 4.4%) outperformed 1FA (septal CoV = 3.1%, six-segment CoV = 5.5%). In six-segment T1 homogeneity, 2FA (CoV = 7.9%) also outperformed 1FA (CoV = 11.1%).

CONCLUSION

The 2FA IR-FLASH improves T1 estimate accuracy and repeatability over 1FA IR-FLASH, and enables single-scan B 1 + -corrected T1 mapping without BHs or ECG when used with MR multitasking.

Development, validation, qualification, and dissemination of quantitative MR methods: Overview and recommendations by the ISMRM quantitative MR study group

On behalf of the International Society for Magnetic Resonance in Medicine (ISMRM) Quantitative MR Study Group, this article provides an overview of considerations for the development, validation, qualification, and dissemination of quantitative MR (qMR) methods. This process is framed in terms of two central technical performance properties, i.e., bias and precision. Although qMR is confounded by undesired effects, methods with low bias and high precision can be iteratively developed and validated. For illustration, two distinct qMR methods are discussed throughout the manuscript: quantification of liver proton-density fat fraction, and cardiac T₁ . These examples demonstrate the expansion of qMR methods from research centers toward widespread clinical dissemination. The overall goal of this article is to provide trainees, researchers, and clinicians with essential guidelines for the development and validation of qMR methods, as well as an understanding of necessary steps and potential pitfalls for the dissemination of quantitative MR in research and in the clinic.

Dual-excitation flip-angle simultaneous cine and T1 mapping using spiral acquisition with respiratory and cardiac self-gating

PURPOSE

To develop a free-breathing cardiac self-gated technique that provides cine images and B1+ slice profile-corrected T₁ maps from a single acquisition.

METHODS

Without breath-holding or electrocardiogram gating, data were acquired continuously on a 3T scanner using a golden-angle gradient-echo spiral pulse sequence, with an inversion RF pulse applied every 4 seconds. Flip angles of 3° and 15° were used for readouts after the first four and second four inversions. Self-gating cardiac triggers were extracted from heart image navigators, and respiratory motion was corrected by rigid registration on each heartbeat. Cine images were reconstructed from the steady-state portion of 15° readouts using a low-rank plus sparse reconstruction. The T₁ maps were fit using a projection onto convex sets approach from images reconstructed using slice profile-corrected dictionary learning. This strategy was evaluated in a phantom and 14 human subjects.

RESULTS

The self-gated signal demonstrated close agreement with the acquired electrocardiogram signal. The image quality for the proposed cine images and standard clinical balanced SSFP images were 4.31 (±0.50) and 4.65 (±0.30), respectively. The slice profile-corrected T₁ values were similar to those of the inversion-recovery spin-echo technique in phantom, and had a higher global T₁ value than that of MOLLI in human subjects.

CONCLUSIONS

Cine and T₁ mapping using spiral acquisition with respiratory and cardiac self-gating successfully acquired T₁ maps and cine images in a single acquisition without the need for electrocardiogram gating or breath-holding. This dual-excitation flip-angle approach provides a novel approach for measuring T₁ while accounting for B1+ and slice profile effect on the apparent T1∗ .

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.

Cardiac cine magnetic resonance fingerprinting for combined ejection fraction, T1 and T2 quantification

This study introduces a technique called cine magnetic resonance fingerprinting (cine-MRF) for simultaneous T₁ , T₂ and ejection fraction (EF) quantification. Data acquired with a free-running MRF sequence are retrospectively sorted into different cardiac phases using an external electrocardiogram (ECG) signal. A low-rank reconstruction with a finite difference sparsity constraint along the cardiac motion dimension yields images resolved by cardiac phase. To improve SNR and precision in the parameter maps, these images are nonrigidly registered to the same phase and matched to a dictionary to generate T₁ and T₂ maps. Cine images for computing left ventricular volumes and EF are also derived from the same data. Cine-MRF was tested in simulations using a numerical relaxation phantom. Phantom and in vivo scans of 19 subjects were performed at 3 T during a 10.9 seconds breath-hold with an in-plane resolution of 1.6 x 1.6 mm² and 24 cardiac phases. Left ventricular EF values obtained with cine-MRF agreed with the conventional cine images (mean bias -1.0%). Average myocardial T₁ times in diastole/systole were 1398/1391 ms with cine-MRF, 1394/1378 ms with ECG-triggered cardiac MRF (cMRF) and 1234/1212 ms with MOLLI; and T₂ values were 30.7/30.3 ms with cine-MRF, 32.6/32.9 ms with ECG-triggered cMRF and 37.6/41.0 ms with T₂ -prepared FLASH. Cine-MRF and ECG-triggered cMRF relaxation times were in good agreement. Cine-MRF T₁ values were significantly longer than MOLLI, and cine-MRF T₂ values were significantly shorter than T₂ -prepared FLASH. In summary, cine-MRF can potentially streamline cardiac MRI exams by combining left ventricle functional assessment and T₁ -T₂ mapping into one time-efficient acquisition.

Low-Rank Tensor Models for Improved Multi-Dimensional MRI: Application to Dynamic Cardiac T 1 Mapping

Multi-dimensional, multi-contrast magnetic resonance imaging (MRI) has become increasingly available for comprehensive and time-efficient evaluation of various pathologies, providing large amounts of data and offering new opportunities for improved image reconstructions. Recently, a cardiac phase-resolved myocardial T 1 mapping method has been introduced to provide dynamic information on tissue viability. Improved spatio-temporal resolution in clinically acceptable scan times is highly desirable but requires high acceleration factors. Tensors are well-suited to describe inter-dimensional hidden structures in such multi-dimensional datasets. In this study, we sought to utilize and compare different tensor decomposition methods, without the use of auxiliary navigator data. We explored multiple processing approaches in order to enable high-resolution cardiac phase-resolved myocardial T 1 mapping. Eight different low-rank tensor approximation and processing approaches were evaluated using quantitative analysis of accuracy and precision in T 1 maps acquired in six healthy volunteers. All methods provided comparable T 1 values. However, the precision was significantly improved using local processing, as well as a direct tensor rank approximation. Low-rank tensor approximation approaches are well-suited to enable dynamic T 1 mapping at high spatio-temporal resolutions.

Robust Single-Shot T2 Mapping via Multiple Overlapping-Echo Acquisition and Deep Neural Network

Quantitative magnetic resonance imaging (MRI) is of great value to both clinical diagnosis and scientific research. However, most MRI experiments remain qualitative, especially dynamic MRI, because repeated sampling with variable weighting parameter makes quantitative imaging time-consuming and sensitive to motion artifacts. A single-shot quantitative T₂ mapping method based on multiple overlapping-echo acquisition (dubbed MOLED-4) was proposed to obtain reliable T₂ mapping in milliseconds. Different from traditional MRI acceleration methods, such as compressed sensing and parallel imaging, MOLED-4 accelerates quantitative T₂ mapping via synchronized multisampling and then deep learning to map the complex nonlinear relationship that is difficult to solve by traditional optimization-based methods. The results of simulation, phantom, and in vivo human brain experiments show the great performance of the proposed method. The principle of MOLED-4 may be extended to other ultrafast quantitative parameter mappings and potentially lead to new dynamic MRI with high efficiency to catch quantitative variation of tissue properties.

Functional LGE Imaging: Cardiac Phase-Resolved Assessment of Focal Fibrosis

Cardiac Magnetic Resonance Imaging (CMR) is a central tool for diagnosis of various ischemic and non-ischemic cardiomyopathies. CMR protocols commonly comprise assessment of functional properties using cardiac phase-resolved CINE MRI and characterization of myocardial viability using late gadolinium enhancement (LGE) imaging. Conventional LGE imaging requires inversion recovery preparation with a specific inversion time to null the healthy myocardium, which restricts the acquisition to a single cardiac phase. In turn, this necessitates separate scans for cardiac function and viability. In this work, we develop a new method for functional LGE imaging in a single breath-hold using a three-step approach: 1) ECG-triggered multi-contrast data is acquired for each cardiac phase, 2) semi-quantitative relaxation maps are generated, 3) LGE imaging contrast is synthesized based on the semi-quantitative maps. The proposed functional LGE method is evaluated in four healthy subject and 20 patients at 1.5T and 3T. Thorough suppression of the healthy myocardium, as well as 40-80ms temporal resolution are achieved, with no visually apparent temporal blurring at tissue interfaces. Functional LGE in patients with focal scar demonstrates robust hyperenhancement in the scar area throughout all cardiac phases, allowing for visual assessment of scar motility. The proposed technique bears the potential to simplify and speedup common cardiac imaging protocols, while enabling improved data fusion of functional and viability information for improved evaluation of CMR.

Free-breathing, non-ECG, continuous myocardial T1 mapping with cardiovascular magnetic resonance multitasking

PURPOSE

To evaluate the accuracy and repeatability of a free-breathing, non-electrocardiogram (ECG), continuous myocardial T₁ and extracellular volume (ECV) mapping technique adapted from the Multitasking framework.

METHODS

The Multitasking framework is adapted to quantify both myocardial native T₁ and ECV with a free-breathing, non-ECG, continuous acquisition T₁ mapping method. We acquire interleaved high-spatial resolution image data and high-temporal resolution auxiliary data following inversion-recovery pulses at set intervals and perform low-rank tensor imaging to reconstruct images at 344 inversion times, 20 cardiac phases, and 6 respiratory phases. The accuracy and repeatability of Multitasking T₁ mapping in generating native T₁ and ECV maps are compared with conventional techniques in a phantom, a simulation, 12 healthy subjects, and 10 acute myocardial infarction patients.

RESULTS

In phantoms, Multitasking T₁ mapping correlated strongly with the gold-standard spin-echo inversion recovery (R² = 0.99). A simulation study demonstrated that Multitasking T₁ mapping has similar myocardial sharpness to the fully sampled ground truth. In vivo native T₁ and ECV values from Multitasking T₁ mapping agree well with conventional MOLLI values and show good repeatability for native T₁ and ECV mapping for 60 seconds, 30 seconds, or 15 seconds of data. Multitasking native T₁ and ECV in myocardial infarction patients correlate positively with values from MOLLI.

CONCLUSION

Multitasking T₁ mapping can quantify native T₁ and ECV in the myocardium with free-breathing, non-ECG, continuous scans with good image quality and good repeatability in vivo in healthy subjects, and correlation with MOLLI T₁ and ECV in acute myocardial infarction patients.

Saturation-Recovery Myocardial T1-Mapping during Systole: Accurate and Robust Quantification in the Presence of Arrhythmia

Myocardial T₁-mapping, a cardiac magnetic resonance imaging technique, facilitates a quantitative measure of fibrosis which is linked to numerous cardiovascular symptoms. To overcome the problems of common techniques, including lack of accuracy and robustness against partial-voluming and heart-rate variability, we introduce a systolic saturation-recovery T₁-mapping method. The Saturation-Pulse Prepared Heart-rate independent Inversion-Recovery (SAPPHIRE) T₁-mapping method was modified to enable imaging during systole. Phantom measurements were used to evaluate the insensitivity of systolic T₁-mapping towards heart-rate variability. In-vivo feasibility and accuracy were demonstrated in ten healthy volunteers with native and post-contrast T₁-mappping during systole and diastole. To show benefits in the presence of RR-variability, six arrhythmic patients underwent native T₁-mapping. Resulting systolic SAPPHIRE T₁-values showed no dependence on arrhythmia in phantom (CoV < 1%). In-vivo, significantly lower T₁ (1563 ± 56 ms, precision: 84.8 ms) and ECV-values (0.20 ± 0.03) than during diastole (T₁ = 1580 ± 62 ms, p = 0.0124; precision: 60.2 ms, p = 0.03; ECV = 0.21 ± 0.03, p = 0.0098) were measured, with a strong correlation of systolic and diastolic T₁ (r = 0.89). In patients, mis-triggering-induced motion caused significant imaging artifacts in diastolic T₁-maps, whereas systolic T₁-maps displayed resilience to arrythmia. In conclusion, the proposed method enables saturation-recovery T₁-mapping during systole, providing increased robustness against partial-voluming compared to diastolic imaging, for the benefit of T₁-measurements in arrhythmic patients.

Locally Low-Rank Tensor Regularization for High-Resolution Quantitative Dynamic MRI

Quantitative dynamic MRI acquisitions have the potential to diagnose diffuse diseases in conjunction with functional abnormalities. However, their resolutions are limited due to the long acquisition time. Such datasets are multi-dimensional, exhibiting interactions between ≥ 4 dimensions, which cannot be easily identified using sparsity or low-rank matrix methods. Hence, low-rank tensors are a natural fit to model such data. But in the presence of multitude of different tissue types in the field-of-view, it is difficult to find an appropriate value of tensor rank, which avoids under- or over-regularization. In this work, we propose a locally low-rank tensor regularization approach to enable high-resolution quantitative dynamic MRI. We show this approach successfully enables dynamic T 1 mapping at high spatio-temporal resolutions.