Quadra-FSE: A Multi-Platform Pulse Sequence for Multispectral qMRI (PD, T1, T2)
Hernan Jara1, Arnaud Guidon2, Jorge A Soto1, and Osamu Sakai1

1Radiology, Boston University, Boston, MA, United States, 2Global MR Applications and Workflow, GE Healthcare, Boston, MA, United States

Synopsis

Purpose: To describe the quadra fast spin-echo (quadra-FSE) pulse sequence, which is the concatenation of two dual-echo FSE acquisitions differing only in TR and to describe the matching qMRI algorithms for mapping T1, T2, and PD. Methods: quadra-FSE was tested at 3T with a multi-compartment agarose phantom and relaxometry was compared to gold standard relaxometry scans using qMRI algorithms developed in house. Results: PD, T1, and T2 maps generated with the quadra-FSE scans are accurate and of excellent image quality. Conclusion: Concatenation of two DE-FSE scans with different TRs can be used for combined and accurate PD, T1, and T2 mapping.

Introduction

Recent interest in the development of multispectral (MS-) qMRI pulse sequences (1-10) has been motivated by the increasing awareness of the clinical benefits afforded by: a) continuous contrast navigation via on-demand MR image synthesis, and b) the improved diagnostic accuracy possible with tissue integrity and hydration measurements; e.g. by the relaxation times (T1, T2) and the water-normalized proton density (PD). Therefore, there is great need for platform-independent MS-qMRI pulse sequences that are scantime efficient and produce superior image quality. Here we report the development of the quadra fast spin-echo (quadra-FSE), which is the concatenation of two dual-echo FSE acquisitions differing solely in the repetition time (Figure 1) and describe the theoretical underpinnings of the accompanying qMRI algorithms for mapping T1, T2, and PD. In addition, we validate the quantitative accuracy of this dual-saturation and dual-echo FSE pulse sequence and the corresponding qMRI algorithms, with a multi-compartment qMRI phantom (Figure 2) independently calibrated with the reference “gold standard” conventional pulse sequences for T1 (multi-TI inversion recovery single-shot-FSE) and multi-TE conventional spin echo for T2.

Methods and Materials

All qMRI algorithms were programmed in Mathcad (PTC, Needham, MA). The T1 mapping algorithm is based on the T1-implicit equation:

see Eq. 1 in Figure 4

which leads to a root-finding problem solvable the Ridder and Brent numerical methods (standard Mathcad functions); is the pixel value ratio between the first echoes of each DE-FSE acquisition and the function F is given by (11, 12):

see Eq. 2 in Figure 4

The T2 mapping algorithm combines the exact dual-echo formula applicable for short and intermediate T2 values and a maximum-likelihood algorithm valid for long T2s. Finally, the PD mapping algorithm is based on the principle of pulse sequence relaxation time unweighting (T1 and T2), followed by low spatial frequency demodulation to correct for receiver coil profile inhomogeneities, and pixel value calibration relative to water. Phantom experimentation was performed at room temperature (20±0.5ºC) using a 3T clinical scanner (Discovery MR750w, GE Healthcare, Waukesha, WI).

Results

The logical architecture of the MS-qMRI algorithms is shown in the flowchart of Figure 3. The top row shows the four directly acquired images per slice generated by the quadra-FSE pulse sequence; these are labelled as Sat1_E1&2 and Sat2_E1&2. In the current implementation, Sat1_E1 and Sat2_E1 are used for T1 mapping via Eq. 1. Furthermore, Sat2_E1 and Sat2_E2 are used for T2 mapping. In the final steps, the residual T1 and T2 weightings of the Sat1_E1 image are reversed, the low spatial frequencies are demodulated, and the pixel values are normalized to that of water. Typical resulting maps are displayed in the bottom row. Additionally, the quadra-FSE T1 and T2 values are in very good agreement (within 10%) with the reference values afforded the gold standard conventional multi-TI and multi-TE conventional pulse sequences.

Conclusion

A new MS-qMRI pulse sequence and accompanying mapping T1, T2, and PD algorithms have been developed. This pulse sequence termed quadra-FSE is simple, does not use inversion pulses which are prone to flip angle errors, is scan-time efficient, and nearly platform independent. This work could have implications for synthetic MRI, clinical and research applications of qMRI, and for correlation time diffusion MRI.

Acknowledgements

This work was supported in part by a research grant from GE Heathcare

References

1. Schmitt P, Griswold MA, Jakob PM, et al. Inversion recovery TrueFISP: quantification of T(1), T(2), and spin density. Magnetic Resonance in Medicine. 2004;51(4):661-7.

2. Farraher SW, Jara H, Chang KJ, Hou A, Soto JA. Liver and spleen volumetry with quantitative MR imaging and dual-space clustering segmentation. Radiology. 2005;237(1):322-8.

3. Suzuki S, Sakai O, Jara H. Combined volumetric T1, T2 and secular-T2 quantitative MRI of the brain: age-related global changes (preliminary results). Magn Reson Imaging. 2006;24(7):877-87.

4. Warntjes JB, Dahlqvist O, Lundberg P. Novel method for rapid, simultaneous T1, T*2, and proton density quantification. Magnetic Resonance in Medicine. 2007;57(3):528-37.

5. Deoni SC, Rutt BK, Arun T, Pierpaoli C, Jones DK. Gleaning multicomponent T1 and T2 information from steady-state imaging data. Magn Reson Med. 2008;60(6):1372-87.

6. Saito N, Sakai O, Ozonoff A, Jara H. Relaxo-volumetric multispectral quantitative magnetic resonance imaging of the brain over the human lifespan: global and regional aging patterns. Magn Reson Imaging. 2009;27(7):895-906.

7. Liu X, Feng Y, Lu ZR, Morrell G, Jeong EK. Rapid simultaneous acquisition of T1 and T2 mapping images using multishot double spin echo EPI and automated variations of TR and TE (ms DSEPI T12). NMR in Biomedicine. 2010;23(1):97-104.

8. Newbould RD, Skare ST, Alley MT, Gold GE, Bammer R. Three-dimensional T1, T2 and proton density mapping with inversion recovery balanced SSFP. Magnetic resonance imaging. 2010;28(9):1374-82.

9. Ehses P, Seiberlich N, Ma D, et al. IR TrueFISP with a golden-ratio-based radial readout: Fast quantification of T1, T2, and proton density. Magnetic Resonance in Medicine. 2012:n/a-n/a.

10. Weiskopf N, Suckling J, Williams G, et al. Quantitative multi-parameter mapping of R1, PD(*), MT, and R2(*) at 3T: a multi-center validation. Frontiers in Neuroscience. 2013;7:95.

11. Lee JN, Riederer SJ. A modified saturation-recovery approximation for multiple spin-echo pulse sequences. Magnetic Resonance in Medicine. 1986;3(1):132-4.

12. Rydberg JN, Riederer SJ, Rydberg CH, Jack CR. Contrast optimization of fluid-attenuated inversion recovery (flair) imaging. Magnetic Resonance in Medicine. 1995;34(6):868-77.

Figures

Figure 1: The herein termed quadra-FSE pulse sequence in its simplest form is the concatenation of two DE-FSE acquisitions with different TR values that are run with the same pre-scan parameters. This pulse sequence uses two different levels of T1 saturation and two levels of T2 weighting thus enables the generation of T1, T2, and PD maps. It is also easy to implement in any platform, produces excellent image quality, and does not use inversion pulses.

Figure 2: qMRI phantom containing multiple compartments with various agarose gels in distilled water (2%, 3%, and 4%), distilled water, and olive oil to represent fat.

Figure 3: Logical flowchart of qMRI mapping algorithms (see text)

Figure 4: Main equations (see text)



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
1817