Nonlinear Coupling between Evoked rCBF and BOLD Signals: A Simulation Study of Hemodynamic Responses

The aim of this work was to investigate the dependence of BOLD responses on different patterns of stimulus input/neuronal changes. In an earlier report, we described an input-state-output model that combined (i) the Balloon/Windkessel model of nonlinear coupling between rCBF and BOLD signals, and (ii) a linear model of how regional flow changes with synaptic activity. In the present investigation, the input-state-output model was used to explore the dependence of simulated PET (rCBF) and fMRI (BOLD) signals on various parameters pertaining to experimental design. Biophysical simulations were used to estimate rCBF and BOLD responses as functions of (a) a prior stimulus, (b) epoch length (for a fixed SOA), (c) SOA (for a fixed number of events), and (d) stimulus amplitude. We also addressed the notion that a single neuronal response may differ, in terms of the relative contributions of early and late neural components, and investigated the effect of (e) the relative size of the late or "endogenous" neural component. We were interested in the estimated average rCBF and BOLD responses per stimulus or event, not in the statistical efficiency with which these responses are detected. The BOLD response was underestimated relative to rCBF with a preceding stimulus, increasing epoch length, and increasing SOA. Furthermore, the BOLD response showed some highly nonlinear behaviour when varying stimulus amplitude, suggesting some form of hemodynamic "rectification." Finally, the BOLD response was underestimated in the context of large late neuronal components. The difference between rCBF and BOLD is attributed to the nonlinear transduction of rCBF to BOLD signal. Our simulations support the idea that varying parameters that specify the experimental design may have differential effects in PET and fMRI. Moreover, they show that fMRI can be asymmetric in its ability to detect deactivations relative to activations when an absolute baseline is stipulated. Finally, our simulations suggest that relative insensitivity to BOLD signal in specific regions, such as the temporal lobe, may be partly explained by higher cognitive functions eliciting a relatively large late endogenous neuronal component.

[1]  M. Kutas,et al.  Reading senseless sentences: brain potentials reflect semantic incongruity. , 1980, Science.

[2]  B. Rosen,et al.  Evidence of a Cerebrovascular Postarteriole Windkessel with Delayed Compliance , 1999, Journal of cerebral blood flow and metabolism : official journal of the International Society of Cerebral Blood Flow and Metabolism.

[3]  T. Allison,et al.  Electrophysiological studies of human face perception. III: Effects of top-down processing on face-specific potentials. , 1999, Cerebral cortex.

[4]  A. Dale,et al.  Functional-Anatomic Correlates of Object Priming in Humans Revealed by Rapid Presentation Event-Related fMRI , 1998, Neuron.

[5]  J. Haxby,et al.  fMRI study of face perception and memory using random stimulus sequences. , 1998, Journal of neurophysiology.

[6]  Karl J. Friston,et al.  Quantitative Comparison of Functional Magnetic Resonance Imaging with Positron Emission Tomography Using a Force-Related Paradigm , 1996, NeuroImage.

[7]  Karl J. Friston,et al.  Nonlinear event‐related responses in fMRI , 1998, Magnetic resonance in medicine.

[8]  G. McCarthy,et al.  Evidence for a Refractory Period in the Hemodynamic Response to Visual Stimuli as Measured by MRI , 2000, NeuroImage.

[9]  Functional Contrast in Event-Related fMRI: Interstimulus Interval Dependency and Blocked Design Comparison , 1998, NeuroImage.

[10]  M. Rugg,et al.  Electrophysiology of Mind. , 1996 .

[11]  R Turner,et al.  The Effect of Inter-Stimulus Interval on Signal Response in fMRI , 1998, NeuroImage.

[12]  Karl J. Friston,et al.  The Effects of Presentation Rate During Word and Pseudoword Reading: A Comparison of PET and fMRI , 2000, Journal of Cognitive Neuroscience.

[13]  Karl J. Friston,et al.  Regional response differences within the human auditory cortex when listening to words , 1992, Neuroscience Letters.

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

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

[16]  Karl J. Friston,et al.  Regionally Specific Sensitivity Differences in fMRI and PET: Where Do They Come From? , 2000, NeuroImage.

[17]  D. Noll,et al.  Nonlinear Aspects of the BOLD Response in Functional MRI , 1998, NeuroImage.

[18]  S. Petersen,et al.  Characterizing the Hemodynamic Response: Effects of Presentation Rate, Sampling Procedure, and the Possibility of Ordering Brain Activity Based on Relative Timing , 2000, NeuroImage.

[19]  A. Dale,et al.  Selective averaging of rapidly presented individual trials using fMRI , 1997, Human brain mapping.

[20]  J. Talairach,et al.  Co-Planar Stereotaxic Atlas of the Human Brain: 3-Dimensional Proportional System: An Approach to Cerebral Imaging , 1988 .

[21]  A W Toga,et al.  Refractory periods observed by intrinsic signal and fluorescent dye imaging. , 1998, Journal of neurophysiology.