Abstract:Conventionally, analysis of functional MRI (fMRI) data relies on available information about the experimental paradigm to establish hypothesized models of brain activity. However, this information can be inaccurate, incomplete or unavailable in multiple scenarios such as resting-state, naturalistic paradigms or clinical conditions. In these cases, blind estimates of neuronal-related activity can be obtained with paradigm-free analysis methods such as hemodynamic deconvolution. Yet, current formulations of the hemodynamic deconvolution problem have three important limitations: (1) their efficacy strongly depends on the appropriate selection of regularization parameters, (2) being univariate, they do not take advantage of the information present across the brain, and (3) they do not provide any measure of statistical certainty associated with each detected event. Here we propose a novel approach that addresses all these limitations. Specifically, we introduce multivariate sparse paradigm free mapping (Mv-SPFM), a novel hemodynamic deconvolution algorithm that operates at the whole brain level and adds spatial information via a mixed-norm regularization term over all voxels. Additionally, Mv-SPFM employs a stability selection procedure that removes the need to select regularization parameters and also lets us obtain an estimate of the true probability of having a neuronal-related BOLD event at each voxel and time-point based on the area under the curve (AUC) of the stability paths. Besides, we present a formulation tailored for multi-echo fMRI acquisitions (MvME-SPFM), which allows us to better isolate fluctuations of BOLD origin on the basis of their linear dependence with the echo time (TE) and to assign physiologically interpretable units (i.e., changes in the apparent transverse relaxation ΔR2∗) to the resulting deconvolved events. Remarkably, we demonstrate that Mv-SPFM achieves comparable performance even when using a single-echo formulation. We demonstrate that this algorithm outperforms existing state-of-the-art deconvolution approaches, and shows higher spatial and temporal agreement with the activation maps and BOLD signals obtained with a standard model-based linear regression approach, even at the level of individual neuronal events. Furthermore, we show that by employing stability selection, the performance of the algorithm depends less on the selection of temporal and spatial regularization parameters λ and ρ. Consequently, the proposed algorithm provides more reliable estimates of neuronal-related activity, here in terms of ΔR2∗, for the study of the dynamics of brain activity when no information about the timings of the BOLD events is available. This algorithm will be made publicly available as part of the splora Python package.

Flexible Amortized Variational Inference in qBOLD MRI

Exploiting Magnetic Resonance Angiography Imaging Improves Model Estimation Of Bold Signal

A frequency-domain machine learning method for dual-calibrated fMRI mapping of oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen consumption (CMRO2)

Quantitative Evaluation of Activation State in Functional Brain Imaging

Validation of a new 3D quantitative BOLD based cerebral oxygen extraction mapping

Learning Large Number of Local Statistical Models Via Variational Bayesian Inference for Brain Voxel Classification in Magnetic Resonance Images

Nonlinear estimation of BOLD signals with the aid of cerebral blood volume imaging

Oxygen extraction fraction mapping with multi-parametric quantitative BOLD MRI: Reduced transverse relaxation bias using 3D-GraSE imaging

Multi-Echo Complex Quantitative Susceptibility Mapping and Quantitative Blood Oxygen Level-Dependent Magnitude (mcQSM + qBOLD or mcQQ) for Oxygen Extraction Fraction (OEF) Mapping

Blood oxygenation level-dependent (BOLD)-based techniques for the quantification of brain hemodynamic and metabolic properties - theoretical models and experimental approaches

Sensitivity assessment of QSM+qBOLD (or QQ) in detecting elevated oxygen extraction fraction (OEF) in physiological change

Whole-brain multivariate hemodynamic deconvolution for functional MRI with stability selection

Bayesian Modelling Approaches for Breath-Hold Induced Cerebrovascular Reactivity

Exploiting Mr Venography Segmentation for the Accurate Model Estimation of Bold Signal

qMRI Diffuser: Quantitative T1 Mapping of the Brain using a Denoising Diffusion Probabilistic Model

Comparison of compartmental analytical Blood-Oxygen-Level-Dependent functional Magnetic Resonance Imaging models against Monte Carlo simulations performed over cortical micro-angiograms

Modelling the Distribution of 3D Brain MRI Using a 2D Slice VAE

Neural Orientation Distribution Fields for Estimation and Uncertainty Quantification in Diffusion MRI

Multi-Parameter Molecular MRI Quantification using Physics-Informed Self-Supervised Learning

Modified test statistics by inter-voxel variance shrinkage with an application to f MRI.

Measurement of Cerebral Oxygen Extraction Fraction Using Quantitative BOLD Approach: A Review