Examining signal changes at the level of cortical layers using T2*-weighted fMRI suffers from a systematic bias caused by extravascular effects of draining veins. T2-weighting promises higher specificity to signal changes near the origin of neural activity. In this work, we use a T2-prepared, multi-echo FLASH sequence at 7 T in a flickering checkerboard experiment targeting the visual cortex to investigate the feasibility of this sequence for laminar fMRI. The layer profiles indicate the desired T2-weighted response in grey matter, but also a peak at the pial surface likely caused by intravascular effects in the pial veins.
Methods
The echo time of the T2-prep module was 50 ms for optimal sensitivity in GM (T2=50ms @ 7 T). This module attenuates the longitudinal magnetisation that is available for the FLASH readout following it. As venous blood has a very short T22 (non-stimulated condition: oxygenation level (Y)=60%, T2=8ms; active condition: Y=70%, T2=13.5ms), the blood signal is expected to be close to fully attenuated by the T2-prep module (exp(-50/13.5) = 2%). As the FLASH multi-echo readout has TEs ranging from 2.6ms to 8.8ms, the extravascular contribution is expected to be attenuated as well, to approximately 25%-60% of its maximum response at TE=T2*GM=30ms. Taken together, the expectation is that the T2-weighting maximises sensitivity to capillary changes in GM, whilst suppressing/attenuating unwanted venous effects. Experiments were performed on a 7T Siemens scanner with a 32ch head coil. A 256x256x32 matrix was acquired per echo, TEs were: 2.88ms, 4.84ms, 6.8ms, and 8.76ms. TR was 11.2ms. We used an undersampling factor of 8 in the phase-encode dimension, but shifted the sampling patterns of the echoes w.r.t. each other to aid the joint reconstruction process3. This method uses all echoes to inform the reconstruction of a single one resulting in lower g-factors. Given the additional phase-encoding between echoes, one could argue this is 32-fold segmented EPI, but as the echoes are individually reconstructed, the FLASH designation seems more appropriate. The volume TR was 11.4 seconds. As T1 relaxation would certainly nullify the T2 preparation over such a long time, 8-fold segmentation was applied, resulting in 1.42 seconds of acquisition time per segment. To futher reduce T1-relaxation artefacts, the encoding was performed centre-out. Activity was induced using a flashing checkerboard in a 40-minute run in a single subject. Structural data (mp-rage) were processed in FreeSurfer to be able to perform cortical depth sampling in the functional data. Standard fMRI processing steps (registration, temporal filtering, etc.) were carried out using SPM12.1. J. Hua et al.: Whole-Brain Three-Dimensional T2-Weighted BOLD Functional Magnetic Resonance Imaging at 7 Tesla. MRM. 72: 1530 - 40 (2014).
2. K. Uludag et al.: An integrative model for neuronal activity-induced signal changes for gradient and spin echo functional imaging. NeuroImage 48: 150 - 65 (2009)
3. B. Bilgic et al.: Improving Parallel Imaging by Jointly Reconstructing Multi-Contrast Data. Proc. MRM. 80: 619 - 32 (2018).