Recent advances in MRI technique such as the multiband sequence have made fast sampling available up to 100-200 ms TR for whole brain acquisition. However, its usefulness is limited because there is not found any significant functional changes needed to be sampled with that short TR. In this study, we challenged to detect task-induced fast signals with the ultra-high temporal resolution imaging with 100 Hz sampling rates.
(Pulse sequence & image reconstruction) A typical GRE sequence with multiple echoes was modified to obtain radial trajectories. Each session consisted of consecutive three scans: initial dummy scan for steady-state magnetization (20 sec.), reference scan, and projection scans in order. The image reconstruction was made by inverse imaging method [4].
(fMRI protocol) To derive the activation in BOLD signal in the primary visual area, a checkerboard pattern flickering at 8 Hz was used. Each session consisted of initial rest period of 20 s, and 4 repetitions of the stimulus period of 20 s and the resting period of 20 s as shown in Figure 1(a), which made a single session 180 s. Four fMRI sessions were made for each subject.
(Data acquisition) Six healthy volunteers participated in this study. All MRI experiments were performed using a 3T Skyra-fit system (Siemens Healthcare GmbH, Erlangen, Germany) with a standard, 16-channel head matrix coil. For imaging, multi-echo projection sequence was used with three TEs of 4.0, 5.3, 6.7 ms, TR of 10 ms, FA of 8°, FOV of 220 mm, in-plane resolution of 3.4 mm, slice thickness of 5.0 mm, and BW of 2kHz/pixel. A single slice was taken to include the calcarine sulcus of the brain. The orientation angles of the read-out gradient were 0°, 120°, and 60°.
(Analysis) The BOLD activation map was estimated by general linear model (GML) analysis. A voxel with the highest beta value and adjacent eight voxels in the primary visual area were selected as a region of interest (10.2 x 10.2 x 5.0 mm3, ROI). To evaluate task-induced fast signal changes, we examined the power spectrum density (PSD) in six frequency bands corresponding to typical EEG signals such as α and β waves that were defined in Table 1. The signals were evaluated in the nine periods as depicted in Figure 1(b). To normalize PSDs, the criterion was taken as the PSD value obtained from the initial period of each session. For the statistical analysis, paired t-test was performed between PSDs of stimulus and resting periods.
[1] Steen Moeller et al., “Multiband multislice GE-EPI at 7 tesla, with 16-fold acceleration using partial parallel imaging with application to high spatial and temporal whole-brain fMRI,” Magn. Reson. Med., vol. 63, no. 5, pp. 1144–1153, 2010.
[2] F.-H. Lin, et al., “Dynamic magnetic resonance inverse imaging of human brain function,” Magn. Reson. Med., vol. 56, no. 4, pp. 787–802, 2006.
[3] J. Hennig, et al., “MR-Encephalography?: Fast multi-channel monitoring of brain physiology with magnetic resonance,” vol. 34, pp. 212–219, 2007.
[4] Kang et al., “2D inverse imaging with multiple-echo scheme for ultra-high sampling in fMRI,” submitted in ISMRM 2017.