Louis Gagnon

Hypercapnic-calibrated fMRI allows the estimation of the relative changes in the cerebral metabolic rate of oxygen (rCMRO ₂ ) from combined BOLD and arterial spin labelling measurements during a functional task, and promises to permit more quantitative analyses of brain activity patterns. The estimation relies on a macroscopic model of the BOLD effect that balances oxygen delivery and consumption to predict haemoglobin oxygenation and the BOLD signal. The accuracy of calibrated fMRI approaches has not been firmly established, which is limiting their broader adoption. We use our recently developed microscopic vascular anatomical network model in mice as a ground truth simulator to test the accuracy of macroscopic, lumped-parameter BOLD models. In particular, we investigate the original Davis model and a more recent heuristic simplification. We find that these macroscopic models are inaccurate using the originally defined parameters, but that the accuracy can be significantly improved by redefining the model parameters to take on new values. In particular, we find that the parameter α that relates cerebral blood-volume changes to cerebral blood-flow changes is significantly smaller than typically assumed and that the optimal value changes with magnetic field strength. The results are encouraging in that they support the use of simple BOLD models to quantify BOLD signals, but further work is needed to understand the physiological interpretation of the redefined model parameters. This article is part of the themed issue ‘Interpreting BOLD: a dialogue between cognitive and cellular neuroscience’.

The computational properties of the human brain arise from an intricate interplay between billions of neurons connected in complex networks. However, our ability to study these networks in healthy human brain is limited by the necessity to use non-invasive technologies. This is in contrast to animal models where a rich, detailed view of cellular-level brain function with cell-type-specific molecular identity has become available due to recent advances in microscopic optical imaging and genetics. Thus, a central challenge facing neuroscience today is leveraging these mechanistic insights from animal studies to accurately draw physiological inferences from non-invasive signals in humans. On the essential path towards this goal is the development of a detailed ‘bottom-up’ forward model bridging neuronal activity at the level of cell-type-specific populations to non-invasive imaging signals. The general idea is that specific neuronal cell types have identifiable signatures in the way they drive changes in cerebral blood flow, cerebral metabolic rate of O ₂ (measurable with quantitative functional Magnetic Resonance Imaging), and electrical currents/potentials (measurable with magneto/electroencephalography). This forward model would then provide the ‘ground truth’ for the development of new tools for tackling the inverse problem—estimation of neuronal activity from multimodal non-invasive imaging data. This article is part of the themed issue ‘Interpreting BOLD: a dialogue between cognitive and cellular neuroscience’.

Oxygen is delivered to brain tissue by a dense network of microvessels, which actively control cerebral blood flow (CBF) through vasodilation and contraction in response to changing levels of neural activity. Understanding these network-level processes is immediately relevant for (1) interpretation of functional Magnetic Resonance Imaging (fMRI) signals, and (2) investigation of neurological diseases in which a deterioration of neurovascular and neuro-metabolic physiology contributes to motor and cognitive decline. Experimental data on the structure, flow and oxygen levels of microvascular networks are needed, together with theoretical methods to integrate this information and predict physiologically relevant properties that are not directly measurable. Recent progress in optical imaging technologies for high-resolution in vivo measurement of the cerebral microvascular architecture, blood flow, and oxygenation enables construction of detailed computational models of cerebral hemodynamics and oxygen transport based on realistic three-dimensional microvascular networks. In this article, we review state-of-the-art optical microscopy technologies for quantitative in vivo imaging of cerebral microvascular structure, blood flow and oxygenation, and theoretical methods that utilize such data to generate spatially resolved models for blood flow and oxygen transport. These “bottom-up” models are essential for the understanding of the processes governing brain oxygenation in normal and disease states and for eventual translation of the lessons learned from animal studies to humans.

The blood oxygenation level-dependent (BOLD) contrast is widely used in functional magnetic resonance imaging (fMRI) studies aimed at investigating neuronal activity. However, the BOLD signal reflects changes in blood volume and oxygenation rather than neuronal activity per se. Therefore, understanding the transformation of microscopic vascular behavior into macroscopic BOLD signals is at the foundation of physiologically informed noninvasive neuroimaging. Here, we use oxygen-sensitive two-photon microscopy to measure the BOLD-relevant microvascular physiology occurring within a typical rodent fMRI voxel and predict the BOLD signal from first principles using those measurements. The predictive power of the approach is illustrated by quantifying variations in the BOLD signal induced by the morphological folding of the human cortex. This framework is then used to quantify the contribution of individual vascular compartments and other factors to the BOLD signal for different magnet strengths and pulse sequences.

Near-infrared spectroscopy (NIRS) is susceptible to signal artifacts caused by relative motion between NIRS optical fibers and the scalp. These artifacts can be very damaging to the utility of functional NIRS, particularly in challenging subject groups where motion can be unavoidable. A number of approaches to the removal of motion artifacts from NIRS data have been suggested. In this paper we systematically compare the utility of a variety of published NIRS motion correction techniques using a simulated functional activation signal added to 20 real NIRS datasets which contain motion artifacts. Principle component analysis, spline interpolation, wavelet analysis, and Kalman filtering approaches are compared to one another and to standard approaches using the accuracy of the recovered, simulated hemodynamic response function (HRF). Each of the four motion correction techniques we tested yields a significant reduction in the mean-squared error (MSE) and significant increase in the contrast-to-noise ratio (CNR) of the recovered HRF when compared to no correction and compared to a process of rejecting motion-contaminated trials. Spline interpolation produces the largest average reduction in MSE (55%) while wavelet analysis produces the highest average increase in CNR (39%). On the basis of this analysis, we recommend the routine application of motion correction techniques (particularly spline interpolation or wavelet analysis) to minimize the impact of motion artifacts on functional NIRS data.