BOLD signal changes can oppose oxygen metabolism across the human cortex
BOLD signal changes can oppose oxygen metabolism across the human cortex
Using multiparametric quantitative fMRI (mqBOLD and pCASL) across tasks, the authors report that ~40% of voxels with significant BOLD responses demonstrate 'discordant' oxygen metabolism—CMRO2 changes opposite to the BOLD sign—particularly in the default mode network; these discordant voxels regulate oxygen demand via OEF changes, while concordant voxels depend mainly on CBF. The findings challenge the canonical hemodynamic interpretation of ∆BOLD and argue that quantitative fMRI more reliably captures absolute and relative changes in neuronal activity.
Community Review
This is a community review of an externally published paper. The original authors retain all rights to their work. TOE-Share provides independent AI analysis — full content is available at the original source linked below.
AI Review Rating
Composite of the review dimensions below, on a 0–5 scale.
This is a rigorous, well-executed study that makes a substantive empirical contribution to the understanding of neurovascular coupling across the human cortex. Using a thoughtfully designed multimodal protocol combining mqBOLD and pCASL in the same imaging session, the authors demonstrate that approximately 40% of cortical voxels with significant BOLD responses exhibit oxygen metabolism changes that oppose the BOLD signal direction. This finding, replicated across tasks, baselines, and an independent cohort, directly challenges the universality of the canonical hemodynamic response model—a challenge that is theoretically grounded, empirically supported, and appropriately caveated.
Particularly compelling is the conjunction voxel analysis, which demonstrates that the same voxels can show canonical responses for positive BOLD during a memory task while showing discordant responses for negative BOLD during a calculation task. This within-voxel design elegantly eliminates voxel-specific artifact explanations. The additional finding that baseline OEF predicts hemodynamic response type—with discordant voxels showing lower OEF and higher CBV, suggesting a larger oxygen buffer—provides a coherent mechanistic account and is independently validated using external QSM datasets. The empirical confirmation of Davis model predictions at the voxel level, using the authors' own data to derive model parameters, is a valuable methodological contribution to the calibrated fMRI literature.
The study has some limitations that are largely acknowledged: the temporal resolution mismatch between BOLD and qfMRI acquisitions could in principle introduce subtle mismatches in the hemodynamic state being characterized, and the neurophysiological basis of OEF-dominated coupling remains speculative. The use of a BOLD-informed surrogate for R2' (rather than fully quantitative T2* mapping) in the main voxel-level analysis is methodologically reasonable but means the CMRO2 estimates in Figure 3 are semi-quantitative and partially depend on the BOLD signal itself—a circularity that warrants clear labeling. These are minor issues relative to the overall strength of the work. This paper represents a meaningful advance for the field of quantitative fMRI and for the interpretation of functional neuroimaging data more broadly.
This review was generated by AI for research and educational purposes. It is not a substitute for formal peer review. All analyses are advisory; publication decisions are based on numerical score thresholds.
This work departs from mainstream consensus physics in the following ways. These are not penalties - they are informational flags that highlight where the author proposes alternative interpretations of physical phenomena. The scores below evaluate rigor, not orthodoxy.
- ◈Challenges the canonical hemodynamic response assumption that positive BOLD reliably indicates increased neuronal activity and negative BOLD indicates decreased activity across the whole cortex.
- ◈Demonstrates that approximately 40% of significant BOLD responses are 'discordant'—a substantially higher proportion than previously quantified in the human cortex using in vivo qfMRI.
- ◈Argues that OEF changes, rather than CBF changes, are the primary mechanism regulating oxygen demand in a substantial fraction (~22%) of task-responsive cortical voxels, departing from the CBF-centric canonical model.
- ◈Shows that negative BOLD responses in the default mode network during cognitive tasks are not reliably associated with decreased CMRO2, contrary to common neuroscientific interpretation of DMN deactivation.
- ◈Empirically validates Davis model parameters at the voxel level across the whole cortex, extending a model previously validated only in specific sensory regions or under hypercapnia conditions.
The paper is internally highly consistent. The central claim—that ~40% of voxels with significant BOLD responses show discordant CMRO2 changes—is supported through multiple complementary analyses: PLS statistics, GLM validation, conjunction voxel analysis, Davis model comparison, and a replication study. The mechanistic explanation (discordant voxels rely on OEF rather than CBF changes) follows logically from Fick's principle and is empirically supported. Claims about baseline OEF predicting hemodynamic response type are coherently connected to the observed task-related findings. The 'mixed voxel' analysis elegantly rules out voxel-specific confounds by demonstrating that the same voxels can show concordant responses in one condition and discordant in another.
The core equations are well-established and correctly applied: Fick's principle for CMRO2, the mqBOLD relation for OEF (R2'/c×sCBV), the BOLD-informed R2' estimation, and the Davis model. The BOLD-informed surrogate for R2' (Equation 4) is a reasonable approximation, and the authors appropriately compare it against fully quantitative estimates in a sensitivity analysis (Supplementary Fig. 5). The Davis model parameter derivation from their own data (α=0.38, M=11.2%) is a meaningful cross-validation. One minor concern is the 25% upscaling correction applied to CBF for background suppression pulses—while referenced, the propagation of this correction into CMRO2 uncertainty is not explicitly quantified in the main text. The CBV normalization to a WM reference value of 2.5% is standard but introduces assumptions not fully discussed. These are minor issues in an otherwise mathematically sound framework.
The work is strongly falsifiable and in fact directly falsifies aspects of the canonical neurovascular coupling assumption. Specific quantitative predictions are tested: the Davis model's prediction of discordant BOLD is empirically confirmed. The finding that baseline OEF predicts hemodynamic response type is a testable hypothesis that could be examined with independent datasets. The use of external QSM datasets to validate lower susceptibility and higher venous density in discordant voxels represents prospective validation. The claim that discordant voxels are distributed uniformly across BOLD amplitude quartiles (no amplitude bias) is a specific, checkable prediction. The entire framework could be falsified if quantitative measurements with PET-MRI showed concordant CMRO2 in these voxels, or if higher-resolution CBF measurements revealed CBF decreases in the negative BOLD regions.
The paper is well-written and logically structured. The introduction clearly motivates the departure from canonical assumptions. Figure 1's explanation of the hemodynamic model and the mqBOLD pipeline is pedagogically effective. The terminology of 'concordant' and 'discordant' voxels is consistently applied throughout. The distinction between BOLD-informed and fully quantitative CMRO2 estimation could be more prominently flagged in the main text (it is relegated to Methods and a supplementary figure). The description of the conjunction voxel analysis is clear and the logic for ruling out confounds is well-articulated. Some figure captions are dense with statistical notation that may require careful reading, but this is appropriate for the technical audience.
This paper makes a genuinely important empirical contribution. While the theoretical possibility of discordant BOLD responses was predicted by the Davis model and suggested by prior animal and limited human studies, no prior work has systematically quantified the prevalence (~40%) and spatial distribution of discordant voxels across the whole human cortex using in-session multimodal qfMRI. The identification of two mechanistically distinct hemodynamic response modes—CBF-dominated versus OEF-dominated—with baseline OEF as a predictor, is novel and provides a physiological framework for understanding regional heterogeneity. The validation of decades-old Davis model parameters empirically from voxel-level data is itself a noteworthy contribution. The 'mixed voxel' finding (same voxel concordant for one task, discordant for another) is a particularly elegant and original demonstration.
The paper is substantially complete within its stated scope. Multiple control analyses address the major alternative explanations: partial volume effects, signal drift, sensitivity limitations, and voxel-specific confounds. The replication study with harmonized voxel sizes is important and strengthens the claims. Areas that could be more fully developed: (1) The neurophysiological interpretation of OEF-mediated CMRO2 regulation is discussed speculatively (astrocytic activity, E/I balance, neuromodulation) but not empirically tested—this is acknowledged as future work; (2) The restriction to healthy young adults (mean age ~32y) and the question of generalizability to clinical populations is acknowledged but not empirically addressed; (3) White matter is appropriately excluded but subcortical gray matter receives limited discussion; (4) The temporal resolution mismatch between BOLD (30-s blocks) and qfMRI (~6-min blocks) is a design limitation that is acknowledged but whose potential impact on the concordance/discordance classification warrants more explicit discussion.
Key Equations (3)
mqBOLD relation that converts R2' and scaled CBV into oxygen extraction fraction (OEF); c is a constant depending on susceptibility, hematocrit and field strength.
Fick's principle used to compute cerebral metabolic rate of oxygen (CMRO2) from OEF, cerebral blood flow (CBF) and arterial oxygen content (CaO2).
Davis (calibrated fMRI) model relating fractional BOLD signal change to relative changes in CBF and CMRO2 via parameters M, α and β.
Other Equations (3)
Reversible transverse relaxation rate (R2') derived from T2* and T2 maps; reflects voxel deoxyhemoglobin content.
Semiquantitative approximation used to estimate task-induced R2' changes from measured BOLD signal change (ΔS/S0), baseline R2' and echo time TE.
Definition of percent BOLD signal change (task T versus baseline B) used in the study.
Testable Predictions (4)
Approximately 40% of cortical voxels with significant task-related BOLD changes exhibit 'discordant' oxygen metabolism, i.e., ∆CMRO2 with sign opposite to ∆BOLD.
Falsifiable if: Independent datasets or replication studies using mqBOLD/pCASL or PET-calibrated measures show a substantially different proportion (e.g., <15% or not statistically different from zero) of discordant voxels across comparable tasks and cortical coverage.
Discordant voxels accommodate task-induced metabolic demand primarily through changes in OEF rather than changes in CBF (∆OEF explains the largest share of ∆CMRO2 variance in discordant voxels).
Falsifiable if: Voxel-wise regression in independent data indicates ∆CBF, not ∆OEF, explains the majority (>50–60%) of ∆CMRO2 variance in the same voxels, or ∆OEF is not a significant predictor after controlling for confounds.
Baseline OEF predicts a voxel's hemodynamic response mode: voxels with lower baseline OEF are more likely to show discordant BOLD–CMRO2 relationships, whereas voxels with higher baseline OEF tend to show canonical concordant coupling.
Falsifiable if: No statistically significant relationship is found between baseline OEF and subsequent classification of voxels as concordant/discordant in independent cohorts or datasets, or the direction of association is reversed.
Negative ∆BOLD does not reliably indicate decreased oxygen metabolism across cortex; negative BOLD areas can show zero or increased ∆CMRO2.
Falsifiable if: Multiple independent, well-powered studies consistently demonstrate that negative ∆BOLD is accompanied by reliably negative ∆CMRO2 across tasks and cortical regions.
Tags & Keywords
Keywords: quantitative fMRI, BOLD signal, cerebral metabolic rate of oxygen (CMRO2), oxygen extraction fraction (OEF), cerebral blood flow (CBF), mqBOLD, Davis model, default mode network (DMN)
Full content is available at the original source:
www.nature.com/articles/s41593-025-02132-9You Might Also Find Interesting
Semantically similar papers and frameworks on TOE-Share