Quantifying hemodynamic refractory bold effects in normal subjects at the single-subject level using an inverse logit fitting procedure

Statistical modeling of time-dependent fMRI activation effects

Functional magnetic resonance imaging (fMRI) activation detection within stimulus-based experimental paradigms is conventionally based on the assumption that activation effects remain constant over time. This assumption neglects the fact that the strength of activation may vary, for example, due to habituation processes or changing attention. Neither the functional form of time variation can be retrieved nor short-lasting effects can be detected by conventional methods. In this work, a new dynamic approach is proposed that allows to estimate time-varying effect profiles and hemodynamic response functions in event-related fMRI paradigms. To this end, we incorporate the time-varying coefficient methodology into the fMRI general regression framework. Inference is based on a voxelwise penalized least squares procedure. We assess the strength of activation and corresponding time variation on the basis of pointwise confidence intervals on a voxel level. Additionally, spatial clusters of effect curves are presented. Results of the analysis of an active oddball experiment show that activation effects deviating from a constant trend coexist with time-varying effects that exhibit different types of shapes, such as linear, (inversely) U-shaped or fluctuating forms. In a comparison to conventional approaches, like classical SPM, we observe that time-constant methods are rather insensitive to detect temporary effects, because these do not emerge when aggregated across the entire experiment. Hence, it is recommended to base activation detection analyses not merely on time-constant procedures but to include flexible time-varying effects that harbour valuable information on individual response patterns.

fMRI hemodynamics accurately reflects neuronal timing in the human brain measured by MEG

Neuronal activation sequence information is essential for understanding brain functions. Extracting such timing information from blood oxygenation level dependent (BOLD) fMRI is confounded by interregional neurovascular differences and poorly understood relations between BOLD and electrophysiological response delays. Here, we recorded whole-head BOLD fMRI at 100 ms resolution and magnetoencephalography (MEG) during a visuomotor reaction-time task. Both methods detected the same activation sequence across five regions, from visual towards motor cortices, with linearly correlated interregional BOLD and MEG response delays. The smallest significant interregional BOLD delay was 100 ms; all delays ≥400 ms were significant. Switching the order of external events reversed the sequence of BOLD activations, indicating that interregional neurovascular differences did not confound the results. This may open new avenues for using fMRI to follow rapid activation sequences in the brain.