Along-tract quantification of resting-state BOLD hemodynamic response functions in white matter

Detailed knowledge of the BOLD hemodynamic response function (HRF) is crucial for accurate analysis and interpretation of functional MRI data. Considerable efforts have been made to characterize the HRF in gray matter (GM) but much less is known about BOLD effects in white matter (WM). However, recent reports have demonstrated reliable detection and analyses of WM BOLD signals after stimulation and in a resting state. WM and GM differ in energy requirements and blood flow, so neurovascular couplings may well be different. We aimed to derive a comprehensive characterization of the HRF in WM across a population, including accurate measurements of its shape and its variation along and between WM pathways, using resting-state fMRI acquisitions. Our results show that the HRF is significantly different between WM and GM. Features of the HRF, such as a prominent initial dip, show strong relationships with features of the tissue microstructure derived from diffusion imaging, and these relationships differ between WM and GM, consistent with BOLD signal fluctuations reflecting different energy demands and differences in neurovascular coupling between tissues of different composition. We also show that the HRF varies significantly along WM pathways, and is different between different WM pathways. Thus, much like in GM, changes in flow and/or oxygenation are different for different parts of the WM. These features of the HRF in WM are especially relevant for interpretation of the biophysical basis of BOLD effects in WM.