Response-mode decomposition of spatio-temporal haemodynamics

The blood oxygen-level dependent (BOLD) response to a neural stimulus is analysed using the transfer function derived from a physiologically based poroelastic model of cortical tissue. The transfer function is decomposed into components that correspond to distinct poles, each related to a response mode with a natural frequency and dispersion relation; together these yield the total BOLD response. The properties of the decomposed components provide a deeper understanding of the nature of the BOLD response, via the components' frequency dependences, spatial and temporal power spectra, and resonances. The transfer function components are then used to separate the BOLD response to a localized impulse stimulus, termed the Green function or spatio-temporal haemodynamic response function, into component responses that are explicitly related to underlying physiological quantities. The analytical results also provide a quantitative tool to calculate the linear BOLD response to an arbitrary neural drive, which is faster to implement than direct Fourier transform methods. The results of this study can be used to interpret functional magnetic resonance imaging data in new ways based on physiology, to enhance deconvolution methods and to design experimental protocols that can selectively enhance or suppress particular responses, to probe specific physiological phenomena.

[1]  S. Rossitti.  Introduction to Functional Magnetic Resonance Imaging, Principles and Techniques , 2002 .

[2]  H. Duvernoy,et al.  Cortical blood vessels of the human brain , 1981, Brain Research Bulletin.

[3]  Cornelius Weiller,et al.  Functional magnetic resonance imaging: A review of methodological aspects and clinical applications , 2003, Journal of magnetic resonance imaging : JMRI.

[4]  Manuel Guizar-Sicairos,et al.  Computation of quasi-discrete Hankel transforms of integer order for propagating optical wave fields. , 2004, Journal of the Optical Society of America. A, Optics, image science, and vision.

[5]  E. Switkes,et al.  Deoxyglucose analysis of retinotopic organization in primate striate cortex. , 1982, Science.

[6]  Karl J. Friston,et al.  Nonlinear Responses in fMRI: The Balloon Model, Volterra Kernels, and Other Hemodynamics , 2000, NeuroImage.

[7]  I. S. Gradshteyn,et al.  Table of Integrals, Series, and Products , 1976 .

[8]  Stanley M. Shinners,et al.  Modern Control System Theory and Design , 1992 .

[9]  R. Buxton,et al.  Modeling the hemodynamic response to brain activation , 2004, NeuroImage.

[10]  Karl J. Friston.  Modalities, Modes, and Models in Functional Neuroimaging , 2009, Science.

[11]  Ravi S. Menon,et al.  Intrinsic signal changes accompanying sensory stimulation: functional brain mapping with magnetic resonance imaging. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[12]  Kazuto Masamoto,et al.  Unveiling astrocytic control of cerebral blood flow with optogenetics , 2015, Scientific Reports.

[13]  Rachid Mansouri,et al.  Approximation of high order integer systems by fractional order reduced-parameters models , 2010, Math. Comput. Model..

[14]  Catie Chang,et al.  Mapping and correction of vascular hemodynamic latency in the BOLD signal , 2008, NeuroImage.

[15]  B. Spehar,et al.  The Foveal Confluence in Human Visual Cortex , 2009, The Journal of Neuroscience.

[16]  R. Buxton,et al.  Dynamics of blood flow and oxygenation changes during brain activation: The balloon model , 1998, Magnetic resonance in medicine.

[17]  Karl J. Friston,et al.  Dynamic causal modelling , 2003, NeuroImage.

[18]  M. Raichle,et al.  The Effects of Changes in PaCO2 Cerebral Blood Volume, Blood Flow, and Vascular Mean Transit Time , 1974, Stroke.

[19]  Alex R. Wade,et al.  Experimental Study of EEG and Bold Responses to Sinusoidal Contrast Modulation , 2007, 2007 IEEE Dallas Engineering in Medicine and Biology Workshop.

[20]  E. DeYoe,et al.  Functional magnetic resonance imaging (FMRI) of the human brain , 1994, Journal of Neuroscience Methods.

[21]  Karl J. Friston,et al.  Comparing hemodynamic models with DCM , 2007, NeuroImage.

[22]  R. Trudnowski,et al.  Specific Gravity of Blood and Plasma at 4 and 37 °C , 1974 .

[23]  P A Robinson,et al.  Spatiotemporal BOLD dynamics from a poroelastic hemodynamic model. , 2010, Journal of theoretical biology.

[24]  Nikolaus Kriegeskorte,et al.  How does an fMRI voxel sample the neuronal activity pattern: Compact-kernel or complex spatiotemporal filter? , 2010, NeuroImage.

[25]  David A. Boas,et al.  A temporal comparison of BOLD, ASL, and NIRS hemodynamic responses to motor stimuli in adult humans , 2006, NeuroImage.

[26]  Essa Yacoub,et al.  Spatio-temporal point-spread function of fMRI signal in human gray matter at 7 Tesla , 2007, NeuroImage.

[27]  P A Robinson,et al.  Spatiotemporal hemodynamic response functions derived from physiology. , 2014, Journal of theoretical biology.

[28]  R. Turner,et al.  Dynamic magnetic resonance imaging of human brain activity during primary sensory stimulation. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[29]  G. Glover,et al.  Retinotopic organization in human visual cortex and the spatial precision of functional MRI. , 1997, Cerebral cortex.

[30]  J. W. Belliveau,et al.  Functional Brain Mapping Using Magnetic Resonance Imaging , 1991, Proceedings of the Annual International Conference of the IEEE Engineering in Medicine and Biology Society Volume 13: 1991.

[31]  Herbert F. Wang.  Theory of Linear Poroelasticity with Applications to Geomechanics and Hydrogeology , 2000 .

[32]  G. Rees,et al.  Predicting the orientation of invisible stimuli from activity in human primary visual cortex , 2005, Nature Neuroscience.

[33]  K. Uğurbil,et al.  Layer-Specific fMRI Reflects Different Neuronal Computations at Different Depths in Human V1 , 2012, PloS one.

[34]  Peter A. Robinson,et al.  BOLD responses to stimuli: Dependence on frequency, stimulus form, amplitude, and repetition rate , 2006, NeuroImage.

[35]  Michael Breakspear,et al.  Hemodynamic Traveling Waves in Human Visual Cortex , 2012, PLoS Comput. Biol..

[36]  Michael Breakspear,et al.  Deconvolution of neural dynamics from fMRI data using a spatiotemporal hemodynamic response function , 2014, NeuroImage.

[37]  Karl J. Friston.  Regulation of rCBF by diffusible signals: An analysis of constraints on diffusion and elimination , 1995 .