0901

3D Magnetic Resonance Fingerprinting with Quadratic Phase in Mouse Brain on Preclinical 7T System
Charlie Wang1, Yuning Gu2, Rasim Boyacioglu2, Charlie Androjna3, Mark Alan Griswold2, and Xin Yu2
1Metrohealth Hospital, Cleveland, OH, United States, 2Case Western Reserve University, Cleveland, OH, United States, 3Cleveland Clinic, Cleveland, OH, United States

Synopsis

Magnetic Resonance Fingerprinting with Quadratic Phase (qRF-MRF) was previously validated in 2D clinical imaging for simultaneous off-resonance, T1, T2, and T2* mapping. Translation of qRF-MRF to high-field preclinical systems for small animal imaging is challenging due to the higher field inhomogeneity and the higher spatial resolution required. Here, a 3D qRF-MRF method was explored to address these challenges. High-resolution simultaneous mapping of off-resonance, T1, T2, and T2* on in vivo mouse brain at 7T was demonstrated. Computational limitations for large dictionary parameter space and reconstruction times were addressed using randomized SVD time compression and quadratic fitting methods.

Introduction

MRI has been widely used in preclinical investigations employing animal models to study human diseases. For example, iron nanoparticle labeled transgenic mouse models has been used to delineate the progression of multiple sclerosis1. Multi-parametric quantitative MR methods enables comprehensive characterization of these models. Magnetic resonance fingerprinting (MRF) has been developed to map multiple tissue properties and system properties rapidly. Previously, a 2D MRF with quadratic RF phase (qRF-MRF) method was developed and validated on clinical systems to map off-resonance (ΔB0), T1, T2, and T2* simultaneously2. Translating this method to high-field preclinical imaging is nontrivial due to additional challenges such as higher imaging resolution, lower SNR, larger field inhomogeneity, and physiological noise. In particular, the large field inhomogeneity at high-field can mask the subtle differences in intrinsic tissue T2*. To overcome these challenges, we aimed to develop a 3D approach for preclinical imaging of mouse brain at 7T. The computational challenge of extending the qRF-MRF method to 3D was explored previously3. In this work, time compression by randomized singular value decomposition (rSVD)4 and tissue property compression by quadratic fit5 were used to reduce computation time. Further, the qRF-MRF sequence parameters were tailored for improve sensitivity to tissue parameters at high-field. Quantitative T1, T2, ΔB0 and T2* maps were obtained from in vivo mouse brain and compared to literature values.

Methods

The 3D preclinical qRF-MRF pulse sequence and sequence parameters are shown in Figure 1. 1758 imaging frames were acquired. Two inversion pulses were implemented at the beginning and middle of the sequence for enhanced T1 sensitivity. A relatively low flip angle (FA) was used in each imaging frame to reduce SAR and improve excitation profile homogeneity (Fig. 1a). RF phase was continuously varied in a quadratic approach (Fig. 1b), resulting in the resonance band response shown in Figure 1c. A constant TR of 5 ms was used for all frames. 3D spatial encoding was performed using a stack-of-spirals sampling pattern6. In-plane spatial encoding used variable density spiral trajectories that required 48 and 12 interleaves to fully sample the peripheral and center of k-space, respectively (Fig. 1c). 192 repetitions were acquired, each with a different combination of rotated in-plane spiral trajectories and 16 phase-encoding steps, to cover an FOV of 20x20x8 mm3. This acquisition scheme achieved a nominal spatial resolution of 156x156x500 μm. Each repetition required 8.8 s, with a 2 s delay allowed for partial magnetization recovery. Total acquisition time was 35 min. Prior to reconstruction, both data and target dictionary were compressed to rank 200 from 1758 time points. As was shown previously3, the truncation matrix was calculated with rSVD of a coarse dictionary that spanned the expected range of physiologic ΔB0, T1, T2, and intravoxel field dispersion (Γ) values. Following tissue property matching with this coarse dictionary, high resolution dictionary tissue property (in the tissue property dimension) matches were recovered with interpolation with the previously proposed quadratic fit method. Here, the inner product values of the closest dictionary entries in the tissue property space with the set of correct tissue properties are known to form a quadratic curve. The neighboring entries in each of the four tissue property dimensions (81 combinations in total) were fit in a quadratic function and the maximum along the fitted curve was taken as the new match. The 3D qRF-MRF was tested on mouse brain, and ROIs were compared to literature values of mouse brain at 7T.

Results

Coarse dictionary match generated in vivo ΔB0, T1, T2, Γ maps from 3D qRF-MRF are shown in Figure 2a-d. T2* calculated from matched T2 and Γ are shown in Figure 2e. Additionally, the relative M0 map from the coarse dictionary match is included (Fig. 2f). Corresponding matches with high parameter resolution after quadratic fit are shown in Fig. 3. Consistent with prior qRF-MRF implementations, parameter mapping in mouse brain using 3D qRF-MRF showed robust ΔB0 mapping. The increased field inhomogeneity at 7T is also evident, particularly near the auditory canals and sinuses due to the susceptibility effects at air-tissue boundaries. Average matching values from manually drawn ROIs are shown compared with literature values in Table 1. While T1 and T2 showed good agreement with the literature, T2* was significantly shorter. This discrepancy might be due to the 3-fold difference in voxel size between our current study and the study published in literature. Experimental conditions may have also contributed to this disagreement as the T2* values reported in literature were significantly longer than T2 values reported in other studies performed at the same field strength.

Discussion and Conclusion

An initial implementation of a 3D qRF-MRF method for simultaneous mapping of T1, T2, off-resonance and T2* on preclinical high-field system was explored. The method can be further improved by increasing the spatial resolution to further reduce the effect of high background field inhomogeneity masking the underlying intrinsic T2*. Additionally, the current method requires long acquisition time as the substantial undersampling potential seen from prior 3D MRF methods (such as 144 fold6) have not yet been fully explored. High resolution 3D qRF-MRF mapping would be a valuable tool to study the vast array of preclinical animal models of human disease and their treatments.

Acknowledgements

This work was supported in part by a grant from NIH (R01-EB023704).

References

1. Nathoo N, Yong V, Dunn J. Understanding disease processes in multiple sclerosis through magnetic resonance imaging studies in animal models. NeuroImage Clin. 2014 4:743-756.

2. Wang C, Coppo S, Mehta B, et al. Magnetic resonance fingerprinting with quadratic RF phase for measurement of T2* simultaneously with δf, T1, and T2. Magn. Reson. Med. 2019 81: 1849– 1862.

3. Boyacioglu R, Wang C, Ma D, et al. 3D Magnetic Resonance Fingerprinting with Quadratic RF Phase. 27th annual ISMRM Proceedings, Montreal (2019).

4. Yang M, Ma D, Jiang Y, et al. Low rank approximation methods for MR fingerprinting with large scale dictionaries. Magn Reson Med. 2018 79(4):2392-2400.

5. McGivney D, Boyacioglu R, Jiang Y, et al. Towards Continuous Tissue Property Resolution in MR Fingerprinting using a Quadratic Inner Product Model. 27th annual ISMRM Proceedings, Montreal (2019).

6. Ma D, Jiang Y, Chen Y, et al. Fast 3D magnetic resonance fingerprinting for a whole-brain coverage. Magn Reson Med. 2018 79(4):2190-2197.

7. Guilfoyle D, Dyakin V, O'Shea J, et al. Quantitative measurements of proton spin‐lattice (T1) and spin–spin (T2) relaxation times in the mouse brain at 7.0 T. Magn. Reson. Med. 2003 49: 576-580.

8. Cong L, Muir E, Chen C, et al. Multimodal MRI Evaluation of the MitoPark Mouse Model of Parkinson’s Disease. PLoS One. 2016 11(3):e0151884.

Figures

Figure 1. Pulse sequence parameters of the preclinical 3D qRF-MRF method. For each repetition, 1758 frames were acquired with varied FA (a), and RF phase (b). The resulting excitation resonance bands (c) are shown. Each frame acquired a rotated variable density spiral trajectory arm (d) which required 48 and 12 trajectories to full sample the outside and inside of k-space, respectively.

Figure 2. Maps from in vivo mouse brain obtained through dictionary matching with coarse dictionary. Dictionary match results from a) off-resonance (ΔB0), b) T1, c) T2, and d) intravoxel frequency dispersion (Γ) are shown. e) Resulting T2* calculated from T2 and Γ. f) Resulting relative M0 calculated from respective matching dictionary entries.

Figure 3. Corresponding high parameter resolution maps from in vivo mouse brain obtained after quadratic fit of coarse dictionary results. a) off-resonance (ΔB0), b) T1, c) T2, and d) intravoxel frequency dispersion (Γ) are shown. e) Resulting T2* calculated from T2 and Γ maps.

Table 1. Average tissue property values from high resolution tissue property maps from ROIs. Literature values7,8 for corresponding anatomy at 7T mouse brain are also shown.

Proc. Intl. Soc. Mag. Reson. Med. 28 (2020)
0901