1151

Data acquisition strategies to mitigate cardiac-induced noise in quantitative R2* maps of the brain
Quentin Raynaud1, Thomas Dardano1,2, Christopher Roy3, Jérôme Yerly3,4, Tobias Kober3,5,6, Ruud B. van Heeswijk3, and Antoine Lutti1
1Department for Clinical Neuroscience, Laboratory for Research in Neuroimaging, Lausanne University Hospital and University of Lausanne, Lausanne, Switzerland, 2Physics section, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland, 3Department of Diagnostic and Interventional Radiology, Lausanne University Hospital and University of Lausanne, Lausanne, Switzerland, 4Center for Biomedical Imaging (CIBM), Lausanne, Switzerland, 5Advanced Clinical Imaging Technology, Siemens Healthineers International AG, Lausanne, Switzerland, 6LTS5, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland

Synopsis

Keywords: Artifacts, Relaxometry

Cardiac pulsation enhances the noise level in MR images of the brain and reduces the sensitivity of the data in studies of brain disease. We propose two data acquisition strategies that mitigate cardiac-induced noise in quantitative brain maps of the MRI parameter R2*. The first strategy sets the number of samples at each k-space location according to the local level of cardiac-induced noise. The second strategy adjusts data acquisition in real-time to acquire the data most sensitive to cardiac-induced noise during the diastolic period of the cardiac cycle.

Introduction

The transverse relaxation rate R2* correlates with e.g. iron and myelin concentration in brain tissue and enables the monitoring of microscopic disease-related changes in patients1-3. However cardiac pulsation enhances the noise level in R2* maps of the brain and reduces the sensitivity of the data to brain tissue changes4. No data acquisition strategy has been shown to mitigate cardiac-induced noise in R2* maps of the brain. Raynaud et al.5 quantified the effect of cardiac-related noise on R2* estimates and identified the affected k-space regions. From these characteristics, we propose data acquisition strategies optimized to mitigate cardiac-induced noise in R2* maps of the brain. The first mitigation strategy sets the number of samples at each k-space location according to the local level of cardiac-induced noise6. The alternative strategy synchronizes prospectively in real-time the acquisition of the most cardiac-sensitive data around the diastolic period of the cardiac cycle7-10.

Methods

Multi-echo data was acquired on 5 participants using a 3T scanner (MAGNETOM Prisma, Siemens Healthcare, Erlangen, Germany) and a customized 3D Cartesian FLASH sequence (repetition time TR =40ms; echo times =2.34ms to 35.10ms with 2.34ms spacing). The pixel size was 2mm and 4mm along the readout and phase-encode directions. Data was acquired continuously for one hour using a sampling kernel optimized to mitigate spurious non-cardiac noise. The data were retrospectively distributed across 12 bins according to the phase of the cardiac cycle at which they were acquired, resulting in a fully sampled five-dimensional datasets (three spatial dimensions, echo time, phase of the cardiac cycle) that enabled the analysis of cardiac pulsation effects on brain R2* maps5.
We conducted numerical simulations of cardiac pulsation effects in multi-echo 3D (i.e. 4D) data normally acquired to compute brain R2* maps. This data was obtained by picking data from the 5D model of cardiac-induced noise according to the k-space sampling trajectory (i.e. k-space coordinates of the data acquired consecutively in time) and a log of cardiac pulsation obtained experimentally using a pulse-oximeter. The ability of the sampling trajectories to mitigate cardiac-induced noise was assessed from the standard deviation of the R2* maps and their goodness of fit (RMSE) across repetitions with the same acquisition but different recordings of cardiac pulsation.
Figure 1 shows characteristics of cardiac-induced noise in k-space data. The amplitude of cardiac-induced noise is radially distributed and is highest around the origin5 (Figure 1A). Within a given radius, cardiac-induced noise follows a Gaussian distribution (Figure 1B). Figure 1C shows the variance of this distribution, estimated within all circles centred on the origin. From this result, we defined a mitigation strategy from the number of samples averaged at each k-space location (Nsamples) that produces the most efficient reduction in effective noise level. This is illustrated in figure 2A for a voxel with cardiac-induced noise variance 10 times that from thermal noise (decrease in effective noise level ~1/Nsamples). At Nsamples =4 an additional sample decreases the effective noise level by 0.5 (‘efficiency’), identical to a voxel dominated by thermal noise where only 1 sample is acquired. The most efficient number of samples in this case is therefore Nsamples =4. Figure 2B shows the distribution of the number of samples for efficient reduction in the effective noise level, which can be best reproduced with a Cartesian pseudo-spiral trajectory with 68 spiral arms with 63 multi-echoes readouts each, and a radial sampling density matched to the optimal number of samples (Figure 2C).
With a standard 2D linear trajectory the k-space centre, most sensitive to cardiac-induced noise, is acquired across all phases of the cardiac cycle (Figure 3A). While the cardiac phase varies smoothly along the fast encoding direction where data points are acquired consecutively every TR =40ms, abrupt changes in the cardiac phase take place along the slow encoding direction, which leads to aliasing in the reconstructed images11,12. However, cardiac-induced noise is strongest around the systolic period of the cardiac cycle ($$$\phi_c=[\pi/2,3\pi/2]$$$) and lowest around the diastolic period ($$$\phi_c=[3\pi/2,\pi/2]$$$, Figure 1D). From this, we defined a real-time mitigation strategy in which the k-space centre data is acquired during the diastole (Figure 3B), minimizing both the level of cardiac-induced noise and aliasing.

Results

Compared to standard linear sampling, the spiral and cardiac-triggered trajectories reduce the variability of R2* by 31% and 52% and of RMSE by 30% and 50% respectively, on average over the whole brain (Figure 4). The increase in scan time for the two mitigation strategies was 72% and 29%. Both strategies lead to a strong reduction in aliasing of cardiac-induced noise from blood vessels (e.g., circle of Willis) into the brain (Figure 4, blue arrows).

Discussion & conclusion

From the characteristics of cardiac-induced noise in k-space, we propose two sampling strategies to mitigate cardiac pulsation effects on brain R2* maps. The first strategy adjusts the number of samples at each k-space location to reduce the effective noise level efficiently. The alternative strategy restricts the acquisition of the sensitive k-space regions during the diastolic period of the cardiac cycle. Our results suggest that the cardiac-triggered strategy reduces cardiac-induced noise in brain R2* maps most efficiently.

Acknowledgements

No acknowledgement found.

References

1. Barbosa JHO, Santos AC, Tumas V, et al. Quantifying brain iron deposition in patients with Parkinson’s disease using quantitative susceptibility mapping, R2 and R2*. Magn. Reson. Imaging 2015;33:559–565 doi: 10.1016/j.mri.2015.02.021.

2. Ulla M, Bonny JM, Ouchchane L, Rieu I, Claise B, Durif F. Is R2* a New MRI Biomarker for the Progression of Parkinson’s Disease? A Longitudinal Follow-Up. PLoS One 2013;8:1–8 doi: 10.1371/journal.pone.0057904.

3. Khalil M, Enzinger C, Langkammer C, et al. Quantitative assessment of brain iron by R2* relaxometry in patients with clinically isolated syndrome and relapsing-remitting multiple sclerosis. Mult. Scler. 2009;15:1048–1054 doi: 10.1177/1352458509106609.

4. Hermes D, Wu H, Kerr AB, Wandell B. Measuring brain beats: cardiac-aligned fast fMRI signals. bioRxiv 2022:2022.02.18.480957.

5. Raynaud Q, Yerly J, Van Heeswijk RB, Lutti A. Characterization of cardiac noise in brain quantitative relaxometry MRI data. ISMRM 2022.

6. Tijssen RHN, Jenkinson M, Brooks JCW, Jezzard P, Miller KL. Optimizing RetroICor and RetroKCor corrections for multi-shot 3D FMRI acquisitions. Neuroimage 2014;84:394–405 doi: 10.1016/j.neuroimage.2013.08.062.

7. Skare S, Andersson JLR. On the effects of gating in diffusion imaging of the brain using single shot EPI. Magn. Reson. Imaging 2001;19:1125–1128 doi: 10.1016/S0730-725X(01)00415-5.

8. Nunes RG, Jezzard P, Clare S. Investigations on the efficiency of cardiac-gated methods for the acquisition of diffusion-weighted images. J. Magn. Reson. 2005;177:102–110 doi: 10.1016/j.jmr.2005.07.005.

9. Habib J, Auer DP, Morgan PS. A quantitative analysis of the benefits of cardiac gating In practical diffusion tensor imaging of the brain. Magn. Reson. Med. 2010;63:1098–1103 doi: 10.1002/mrm.22232.

10. Bopp MHA, Yang J, Nimsky C, Carl B. The effect of pulsatile motion and cardiac-gating on reconstruction and diffusion tensor properties of the corticospinal tract. Sci. Rep. 2018;8:1–12 doi: 10.1038/s41598-018-29525-0.

11. Tijssen RHN, Okell TW, Miller KL. Real-time cardiac synchronization with fixed volume frame rate for reducing physiological instabilities in 3D FMRI. Neuroimage 2011;57:1364–1375 doi: 10.1016/j.neuroimage.2011.05.070.

12. Kristoffersen A, Goa PE. Cardiac-induced physiological noise in 3D gradient echo brain imaging: Effect of k-space sampling scheme. J. Magn. Reson. 2011;212:74–85 doi: 10.1016/j.jmr.2011.06.012.

Figures

Figure 1: Cardiac-induced noise characteristics. A) Amplitude of the modeled cardiac noise (reproduced from Raynaud et al.5). B) Distribution of the modeled cardiac noise within a given radius from the origin of k-space. C) Variance of this distribution estimated within all circles centred on the k-space origin. D) Mean cardiac noise amplitude in k-space as a function of the cardiac phase, averaged across all participants.

Figure 2: Cardiac noise mitigation by efficient averaging. A) Effective noise level and averaging efficiency in a voxel with cardiac-induced noise variance that is 10 times larger than the thermal noise. The most efficient number of samples in this case is Nsamples =4. B) Distribution of the number of samples that allows an efficient reduction in the effective noise level. This distribution can best be reproduced with a Cartesian pseudo-spiral trajectory (C), where all missing k-space lines are in the periphery and acquired once in a linear sampling strategy after the spiral.

Figure 3: Cardiac noise mitigation by synchronizing data acquisition with the cardiac cycle. A) With a 2D linear trajectory the phase of the cardiac cycle at the time of acquisition of the data shows abrupt changes along the slow encoding direction, which leads to aliasing in the reconstructed images. B) With the cardiac-triggering strategy, the k-space centre is acquired during the diastolic period of the cardiac cycle and the changes in the cardiac phase at the time of acquisition of the data are smooth. The red line indicates the k-space data at which cardiac triggering is performed.

Figure 4: Variability of R2* and RMSE maps across 100 repetitions, for the linear, spiral, and cardiac-triggered samplings. The blue arrows point at a region of high R2* variability caused by the aliasing of cardiac-induced noise from blood vessels present in standard 2D linear sampling.

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