Whole brain measurements of the hemodynamic response function variability during a finger tapping task at 7T show regional differences in hrf profiles.
Yohan Boillat1, Rolf Gruetter1,2,3, and Wietske van der Zwaag1,4

1Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland, 2Universite de Lausanne, Lausanne, Switzerland, 3Universite de Geneve, Geneva, Switzerland, 4Spinoza Centre for Neuroimaging, Amsterdam, Switzerland

Synopsis

The temporal properties of the BOLD response, including inter-trial variability, were studied using a 3D-EPI-CAIPI acquisition at 7T with TR=400ms. The HRFs of six different brain regions in the motor network were characterized, showing a reduced post-stimulus undershoot in the cerebellar regions of interest, as well as differences in peak height, with higher response amplitude in M1 than any other region and onset time, with the cerebellar lobule VIII response starting later than the other ROIs. Trial-to-trial variability was highest in CVIII and lowest in SI.

Target Audience

Physicists and neuroscientists interested in fast sampling sequences for fMRI and/or hrf properties.

Purpose

Temporal resolution in fMRI is typically limited to 2-3s per volume. However, recent developments in acquisition schemes, such as multiband1 and 3D-EPI-CAIPI2 , are overcoming this constraint. The 3D-EPI-CAIPI sequence can provide whole brain coverage including the cerebellum, in 400 milliseconds. This allows a precise hemodynamic response function (HRF) characterization without having to jitter the stimulus onsets.

The purpose of this study was to investigate the temporal properties of the BOLD response, including inter-trial variability, using a high-temporal resolution acquisition at 7T. Several different HRF parameters, such as time of onset, time to peak and undershoot amplitude, as well as inter-trial variability, were measured and compared between six brain regions involved in motor tasks.

Methods

Seven participants (3 females; 19-24 yrs) were scanned at 7T (Siemens, Germany). Participants performed a bilateral finger tapping task following a visual cue. Functional ROIs were defined based on a block task localizer (15s ON, 30s OFF alternated). HRF measurements were then taken from a 7-minute event-related run in which the subjects moved their fingers once every 30s.

For the localizer, a 2D EPI (2*2*2mm voxels, 30 coronal-oblique slices, TR/TE: 2500/26ms) was used. For the event-related task, a 3D-EPI-CAIPI (2*2*2mm voxels, 60 coronal-oblique slices, TR/TE:400/27ms, GRAPPAz:6, ΔCAIPIRINHA:2, partial Fourier:6/8) was used. Respiratory and cardiac traces were recorded for physiological noise removal using RETROICOR.

Localizer data were GLM analyzed (spm12, p<0.05) to define six bilateral ROIs: primary motor cortex (M1), primary sensory cortex (S1), supplementary motor area (SMA), secondary sensory cortex (S2), cerebellar lobule V (CV) and cerebellar lobule VIII (CVIII). The S1/M1 boundary was drawn manually. The timecourses of the within-ROI voxels were extracted from the realigned 3D-EPI-CAIPI data (CAIPI). Timecourses were also downsampled for comparison (dsCAIPI; five point average; TR=2sec) and the following was performed on both: Voxel timecourses were filtered, physiological noise was regressed out and subsequently, timecourses were averaged within ROIs, then normalized to baseline (two seconds before stimulus onset). Three inverse logit functions3,4 were fitted to the mean timecourses using the Nelder-Mead Simplex method (Figure 1A,B) to obtain positive peak height, onset time, time-to-peak, full-width at half maximum (FWHM) and undershoot amplitude (Figure 1B). Repeated measures statistics, using a linear mixed model with ROIs as repeated factors, were performed independently on CAIPI and dsCAIPI results, followed by post-hoc analysis (significance level p<.05; FWE ).

Results

The six different ROIs all clearly demonstrate an HRF (Figure 2). However, the variability of the HRF illustrated by the shaded band (i.e. mean standard deviation) is important especially for the cerebellar regions, and is also notably higher in the post-peak period than prior to stimulus onset. For the CAIPI data, significant differences (p<.05) between ROIs were found for peak height, onset time, FWHM and the undershoot amplitude (Figure 3). Cerebral regions showed higher peak amplitude compared to the cerebellar ones, with M1 highest overall. The undershoot amplitude in CV/CVIII and S2 was significantly smaller than in M1, and also somewhat smaller than in S1 and SMA. The absence of a post-stimulus undershoot in the cerebellum is also visible in the time courses (Figure 2). In terms of timing, the HRF started later for the cerebellum, especially CVIII, than most of the cerebral ROIs. Trial-to-trial variability (Figure 2) was highest in CVIII and lowest in SI. Concerning the dsCAIPI, a main effect was observed only for the height and the FWHM (p<.05) and no significant post-hoc differences were found.

Discussion

Several significant differences in terms of amplitudes and timings were observed between different brain regions involved in motor control when using a TR=400ms 3D-EPI-CAIPI acquisition. This temporal resolution means trial-to-trial fluctuations can also be studied, showing different behavior even between S1 and M1. In contrast, the TR=2 of dsCAIPI did not provide enough information to obtain measurements allowing a good characterization of the HRF.

In the 400ms CAIPI data, the post-stimulus undershoot in cerebellar regions was consistently smaller and even absent in several subjects. The underlying cause could be a vascular phenomenon, such as different vascular compliance5 in the cerebellum, or a different post-stimulus neuronal modulation6 from the cerebrum. The lack of HRF studies targeting the cerebellum and its very different anatomy from the forebrain still leave both hypotheses open.

Conclusion

We conclude that the HRFs of six different brain regions exhibit between-region differences, such as a smaller cerebellar post-stimulus undershoot and high trial-to-trial stability in the primary sensory area.

Acknowledgements

This work was supported by the Swiss National Science Fondation

References

1. Chen, L. et al. Evaluation of highly accelerated simultaneous multi-slice EPI for fMRI. Neuroimage 104, 452–459 (2015).

2. Mayur Narsude, Daniel Gallichan, Wietske van der Zwaag, Rolf Gruetter and Jose P. Marques. Three-Dimensional Echo Planar Imaging with Controlled Aliasing: A Sequence for High Temporal Resolution Functional MRI. Magn. Reson. Med. (2015).

3. Lindquist, M. a., Meng Loh, J., Atlas, L. Y. & Wager, T. D. Modeling the hemodynamic response function in fMRI: efficiency, bias and mis-modeling. Neuroimage 45, S187–S198 (2009).

4. Lindquist, M. a. & Wager, T. D. Validity and power in hemodynamic response modeling: A comparison study and a new approach. Hum. Brain Mapp. 28, 764–784 (2007).

5. Van Zijl, P. C. M., Hua, J. & Lu, H. The BOLD post-stimulus undershoot, one of the most debated issues in fMRI. Neuroimage 62, 1092–1102 (2012).

6. Mullinger, K. J., Mayhew, S. D., Bagshaw, A. P., Bowtell, R. & Francis, S. T. Poststimulus undershoots in cerebral blood flow and BOLD fMRI responses are modulated by poststimulus neuronal activity. Proc. Natl. Acad. Sci. U. S. A. 110, 13636–41 (2013).

Figures

Figure 1. A) The solid lines represent the three inverse logit function fitted on the HRF time course (dotted line). B) Fitted HRF with the different measured parameters: a) height of the peak, b) time of onset, c) time to peak, d) FWHM, e) undershoot amplitude.

Figure 2. Between-subject averaged time courses (dark line) for the six different ROIs: cerebellar lobules VIII (CVIII) and V (CV), primary motor cortex (M1), primary sensory cortex (S1), secondary sensory cortex (S2) and supplementary motor area (SMA). The shaded band represents mean standard deviation.

Figure 3. Measures of the HRF parameters (mean with standard deviation) of the CAIPI data (TR=400ms) with significant differences between ROIs for the height, the time of onset and the undershoot. The horizontal bars represent the statistically significant differences (p<.05).

Figure 4. Measures of the HRF parameters (mean with standard deviation) of the downsampled CAIPI data (TR=2000ms) with differences between ROIs for the height only. The horizontal bars represent the statistically significant differences (p<.05).



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