3396

The effect of black blood and fat suppression prepulses on signal modelling in T1-mapping and dynamic contrast enhanced MRI
Wilhelm Stehling1, Myrte Wennen1,2, Eric Schrauben1, Pim van Ooij1, Kak Khee Yeung3, Aart Nederveen1, and Oliver Gurney-Champion1,4
1Radiology and Nuclear Medicine, Amsterdam University Medical Centers, Amsterdam, Netherlands, 2Department of Intensive Care, Erasmus Medical Centre, Rotterdam, Netherlands, 3Vascular Surgery, Amsterdam University Medical Centers, Amsterdam, Netherlands, 4Cancer Center Amsterdam, Imaging and Biomarkers, Amsterdam, Netherlands

Synopsis

Keywords: System Imperfections: Measurement & Correction, Simulations

In this work, we derived an equation describing the signal evolution during an SGRE sequence with preparation pulses. The signal changes especially during the earlier RF pulses and differs from the values for later RF pulses, where it approaches the signal from an SPGR sequence without prepulses. Using the derived equation might enable more accurate DCE examinations.

Introduction

Abdominal aortic aneurysm (AAA) is a degenerative inflammatory disease of the aortic wall which when getting ruptured becomes life-threatening. AAA progression may be assessed using dynamic contrast-enhancement (DCE) MRI1. Imaging the vessel walls of AAA patients requires suppression of blood and fat signal to accurately image the contrast uptake. However, the usage of blood and fat suppression disrupts the steady state of the spoiled gradient echo (SGRE), which is generally used to acquire DCE MRI. Therefore, conventional approaches to convert signal into contrast concentration may be inadequate. Hence, to investigate the effect of black blood and fat suppression prepulses on the signal, we modelled the signal evolution for an SGRE combined with preparation prepulses and validated it in a phantom with known T1 properties.

Methods

Theory
A GRE with fat suppression and black blood typically consists of a block of N (Turbo Field Factor) regular GRE pulses, interleaved with a block with a SAR and fat suppression block (to apply fat suppression and limit SAR), a black blood block2 and an additional spoiling block (Fig.1). During the fat suppression and spoiling blocks, the spins of interest undergo T1 relaxation, whereas during the black blood block, the spins of interest undergo T2 relaxation as they are in the x-y plane. To accommodate for this additional T1 and T2 relaxation, we extended the formula of a conventional SGRE to include these pre-pulses.
The longitudinal magnetization available before each of the N RF-pulses within a shot is given by:

$$ M_n^-=M_0 (1-E_1 )*\frac{1-[E_1*cos⁡(α) ]^n}{1-E_1*cos⁡(α) }+M_0^-*cos^n⁡(α)*E_1^n $$

Where n is the number of the RF pulses starting from 0 until N-1, N is the total number of RF pulses within one shot and is the longitudinal magnetization at the end of the preparation just before the first RF pulse ($$$\alpha$$$ is the flip angle and $$$M_0$$$ the net magnetization in z-direction; See Fig.1). Instead of having a steady-state during each RF pulse, the steady-state will repeat itself during each block of RF pulses3,4, which leads to $$$M_0^-$$$:

$$ M_0^-=((M_N^-*E_{Fat-sat}+M_0 (1-E_{Fat-sat}))*E_{Black-blood} ) E_{Spoil}+M_0 (1-E_{Spoil}) $$

During the first break, in which a fat-saturation pulse can be performed, the magnetization of interest is parallel to the B0 field and only T1 relaxation takes place ($$$E_{Fat-sat} = exp(-T_{Fat-saturation}/T_1)$$$). The magnetization available before the black-blood pulse only undergoes T2 decay as it is in the x-y plane (due to the 90°-180°-180°-90° pulses). During the last gap, in which the black-blood spoiling gradient is played out, only T1 relaxation occurs ($$$E_{Spoil} = exp(-T_{Spoil}/T_1)$$$). Inserting eq.2 into eq.1 provides the signal evolution depending on T1, T2 and the imaging parameters:

$$$ {M_n}^- = M_{0} (1-E_{1})\frac{1-[E_{1}cos(\alpha)]^{n}}{1-E_{1}cos(\alpha)}+[E_{1}cos(\alpha)]^{n}*M_{0}\frac{(1-E_{Spoil})+(1-E_{1})E_{Fat-sat}E_{Black-blood}E_{Spoil}\frac{1-[E_{1}cos(\alpha)]^{N}}{1-E_{1}cos(\alpha)}+(1-E_{Fat-sat})E_{Black-blood}E_{Spoil}}{1-cos(\alpha)^{N}{E_{1}}^{N}E_{Fat-sat}E_{Black-blood}E_{Spoil}} $$$

Validation:
To validate the equation we scanned the NIST phantom (CaliberMRI, T1 reference values 10ms-1879ms) on a 3T system (Ingenia, Philips, Best, Netherlands) using a 16-channel head coil. Imaging was performed with a SGRE sequence with 24 RF pulses and four different flip angles (2°, 7°, 15°, 25°) (TR=10ms, TBlack-blood=13ms, TSpoil=6.3ms, TFat-saturation=82ms). K-space was fully sampled using the PROUD patch5 and each k-line was sampled 24 times before going to the next line. Separate images with data only from one RF pulse during each shot were then reconstructed using ReconFrame (Gyrotools GmbH, Switzerland). For the RF pulse for which the model differs most from SGRE (first pulse in the shot), T1 maps are acquired and a Bland-Altman plot is generated showing the difference between the estimated and ground truth T1 values for conventional SGRE fitting and our new approach. Additionally, the carotid artery of a healthy volunteer was imaged using the same acquisition technique to visually assess the strength of black-blood and fat suppression during the shot.

Results

Spheres with shorter T1 and T2 values show a major signal change during the earlier RF pulses (Fig.2). For later RF pulses, the signal reaches a plateau which is similar to the signal from an SGRE. For the spheres with lower T1 values, the fit using the derived equation is closer to the reference values than using a normal SGRE signal equation (Fig.3). For higher T1 values, both fits deviate from the reference values. Although the measured data is described well by both models (Fig.4) for the flip angles measured, the T1 values estimated by SGRE deviate substantially.

Discussion

In this work, we derived an equation describing the signal evolution during an SGRE sequence with preparation pulses. The signal changes especially during the earlier RF pulses and differs from the values for later RF pulses, where it approaches signal from an SPGR sequence without prepulses. Depending on the acquisition settings, in particular for protocols that fill low k-lines first, it might be useful to take our new equation into account. One large disadvantage of our equation is the additional T2-dependence of the signal. Skipping early RF pulses during each shot returns the signal to (almost) SGRE and hence can overcome this dependency on T2, rendering T1 mapping less challenging. Important to consider then is the decreasing effect of preparation pulses as shown in the animation in Fig.5, where fat signal returns, while the blood suppression is visible throughout the entire shot.

Acknowledgements

No acknowledgement found.

References

1. Nguyen VL, Backes WH, Kooi ME, Wishaupt MC, Hellenthal FA, Bosboom EM, et al. Quantification of abdominal aortic aneurysm wall enhancement with dynamic contrast-enhanced MRI: feasibility, reproducibility, and initial experience. J Magn Reson Imaging. 2014;39(6):1449-56.

2. Wang J, Yarnykh VL, Yuan C. Enhanced image quality in black-blood MRI using the improved motion-sensitized driven-equilibrium (iMSDE) sequence. J Magn Reson Imaging. 2010;31(5):1256-63.

3. Myrte Wennen MM, Tim Marcus, Joost Kuijer, Leo Heunks, Christina Lavini, Gustav Strijkers, Aart Nederveen, and Oliver Gurney-Champion. A model for fat-suppressed variable flip angle T1-mapping and dynamic contrast enhanced MRI. ISMRM 20222022.

4. Qi H, Huang F, Zhou Z, Koken P, Balu N, Zhang B, et al. Large coverage black-bright blood interleaved imaging sequence (LaBBI) for 3D dynamic contrast-enhanced MRI of vessel wall. Magn Reson Med. 2018;79(3):1334-44.

5. Gottwald LM, Peper ES, Zhang Q, Coolen BF, Strijkers GJ, Nederveen AJ, et al. Pseudo-spiral sampling and compressed sensing reconstruction provides flexibility of temporal resolution in accelerated aortic 4D flow MRI: A comparison with k-t principal component analysis. NMR Biomed. 2020;33(4):e4255.

Figures

Figure 1: Sequence design: After N RF-pulses with flip angle α and separated by TR a first temporal gap $$$T_{Fat-saturation}$$$ takes place in which a fat saturation could be played out (or just a temporal gap to stick to SAR limits) and T1 decay occurs. After that, an iMSDE black-blood suppression2 is applied (modeled by T2 decay), followed by a second temporal gap $$$T_{Spoil}$$$ in which the iMSDE spoiling gradient takes place.

Figure 2: Simulated signal evolution for the 14 spheres of the T1 layer in the NIST phantom for flip anlge 2° and 15°. Imaging acq. parameters are: TR=10ms, 24 RF pulses, $$$T_{Black-blood}$$$=13 ms, $$$T_{Spoil}$$$=6.3ms, $$$T_{Fat-saturation}$$$=82 ms.

Figure 3: T1 fitting results using the presented equation (red) and a normal SGRE signal equation without any preparation (blue).

Figure 4: For three different T1 reference values (1184ms, 690ms, 86ms) four different measurements were performed (2°,7°,15°,25°). The derived and normal SGRE equations without any preparation were both fit to these sampled points. Fitted T1 values are shown in the legend.

Figure 5: Using the descriped sampling scheme separate images with k-space data acquired only after one specific RF pulse were reconstructed. The preparation consists of a SPIR fat suppression prepulse and an iMSDE prepulse with 13ms duration. As can be seen, the fat signal returns during the shot while the blood remains dark. Shown is an axial slice of the neck region of a healthy volunteer.

Proc. Intl. Soc. Mag. Reson. Med. 31 (2023)
3396
DOI: https://doi.org/10.58530/2023/3396