Estimating cerebral oxygen metabolism from fMRI with a dynamic multicompartment Windkessel model

Time-Lagged Functional Ultrasound for Multi-Parametric Cerebral Hemodynamic Imaging

We introduce an ultrasound speckle decorrelation-based time-lagged functional ultrasound technique (tl-fUS) for the quantification of the relative changes in cerebral blood flow speed (rCBF [Formula: see text]), cerebral blood volume (rCBV) and cerebral blood flow (rCBF) during functional stimulations. Numerical simulations, phantom validations, and in vivo mouse brain experiments were performed to test the capability of tl-fUS to parse out and quantify the ratio change of these hemodynamic parameters. The blood volume change was found to be more prominent in arterioles compared to venules and the peak blood flow changes were around 2.5 times the peak blood volume change during brain activation, agreeing with previous observations in the literature. The tl-fUS shows the ability of distinguishing the relative changes of rCBFspeed, rCBV, and rCBF, which can inform specific physiological interpretations of the fUS measurements.

A quantitative model for human neurovascular coupling with translated mechanisms from animals

Neurons regulate the activity of blood vessels through the neurovascular coupling (NVC). A detailed understanding of the NVC is critical for understanding data from functional imaging techniques of the brain. Many aspects of the NVC have been studied both experimentally and using mathematical models; various combinations of blood volume and flow, local field potential (LFP), hemoglobin level, blood oxygenation level-dependent response (BOLD), and optogenetics have been measured and modeled in rodents, primates, or humans. However, these data have not been brought together into a unified quantitative model. We now present a mathematical model that describes all such data types and that preserves mechanistic behaviors between experiments. For instance, from modeling of optogenetics and microscopy data in mice, we learn cell-specific contributions; the first rapid dilation in the vascular response is caused by NO-interneurons, the main part of the dilation during longer stimuli is caused by pyramidal neurons, and the post-peak undershoot is caused by NPY-interneurons. These insights are translated and preserved in all subsequent analyses, together with other insights regarding hemoglobin dynamics and the LFP/BOLD-interplay, obtained from other experiments on rodents and primates. The model can predict independent validation-data not used for training. By bringing together data with complementary information from different species, we both understand each dataset better, and have a basis for a new type of integrative analysis of human data.

Characterization of cortical hemodynamic changes following sensory, visual, and speech activation by intraoperative optical imaging utilizing phase-based evaluation methods

Alterations within cerebral hemodynamics are the intrinsic signal source for a wide variety of neuroimaging techniques. Stimulation of specific functions leads due to neurovascular coupling, to changes in regional cerebral blood flow, oxygenation and volume. In this study, we investigated the temporal characteristics of cortical hemodynamic responses following electrical, tactile, visual, and speech activation for different stimulation paradigms using Intraoperative Optical Imaging (IOI). Image datasets from a total of 22 patients that underwent surgical resection of brain tumors were evaluated. The measured reflectance changes at different light wavelength bands, representing alterations in regional cortical blood volume (CBV), and deoxyhemoglobin (HbR) concentration, were assessed by using Fourier-based evaluation methods. We found a decrease of CBV connected to an increase of HbR within the contralateral primary sensory cortex (SI) in patients that were prolonged (30 s/15 s) electrically stimulated. Additionally, we found differences in amplitude as well as localization of activated areas for different stimulation patterns. Contrary to electrical stimulation, prolonged tactile as well as prolonged visual stimulation are provoking increases in CBV within the corresponding activated areas (SI, visual cortex). The processing of the acquired data from awake patients performing speech tasks reveals areas with increased, as well as areas with decreased CBV. The results lead us to the conclusion, that the CBV decreases in connection with HbR increases in SI are associated to processing of nociceptive stimuli and that stimulation type, as well as paradigm have a nonnegligible impact on the temporal characteristics of the following hemodynamic response.

Near-Infrared Optical Spectroscopy In Vivo Distinguishes Subjects with Alzheimer's Disease from Age-Matched Controls

Background:

Medical imaging methods such as PET and MRI aid clinical assessment of Alzheimer’s disease (AD). Less expensive, less technically demanding, and more widely deployable technologies are needed to expand objective screening for diagnosis, treatment, and research. We previously reported brain tissue near-infrared optical spectroscopy (NIR) in vitro indicating the potential to meet this need.

Objective:

To determine whether completely non-invasive, clinical, NIR in vivo can distinguish AD patients from age-matched controls and to show the potential of NIR as a clinical screen and monitor of therapeutic efficacy.

Methods:

NIR spectra were acquired in vivo. Three groups were studied: autopsy-confirmed AD, control and mild cognitive impairment (MCI). A feature selection approach using the first derivative of the intensity normalized spectra was used to discover spectral regions that best distinguished “AD-alone” (i.e., without other significant neuropathology) from controls. The approach was then applied to other autopsy-confirmed AD cases and to clinically diagnosed MCI cases.

Results:

Two regions about 860 and 895 nm completely separate AD patients from controls and differentiate MCI subjects according to the degree of impairment. The 895 nm feature is more important in separating MCI subjects from controls (ratio-of-weights: 1.3); the 860 nm feature is more important for distinguishing MCI from AD (ratio-of-weights: 8.2).

Conclusion:

These results form a proof of the concept that near-infrared spectroscopy can detect and classify diseased and normal human brain in vivo. A clinical trial is needed to determine whether the two features can track disease progression and monitor potential therapeutic interventions.

High spatiotemporal vessel-specific hemodynamic mapping with multi-echo single-vessel fMRI

High-resolution fMRI enables noninvasive mapping of the hemodynamic responses from individual penetrating vessels in animal brains. Here, a 2D multi-echo single-vessel fMRI (MESV-fMRI) method has been developed to map the fMRI signal from arterioles and venules with a 100 ms sampling rate at multiple echo times (TE, 3-30 ms) and short acquisition windows (<1 ms). The T₂*-weighted signal shows the increased extravascular effect on venule voxels as a function of TE. In contrast, the arteriole voxels show an increased fMRI signal with earlier onset than venules voxels at the short TE (3 ms) with increased blood inflow and volume effects. MESV-fMRI enables vessel-specific T2* mapping and presents T2*-based fMRI time courses with higher contrast-to-noise ratios (CNRs) than the T2*-weighted fMRI signal at a given TE. The vessel-specific T2* mapping also allows semi-quantitative estimation of the oxygen saturation levels (Y) and their changes (ΔY) at a given blood volume fraction upon neuronal activation. The MESV-fMRI method enables vessel-specific T2* measurements with high spatiotemporal resolution for better modeling of the fMRI signal based on the hemodynamic parameters.

A quantitative analysis of cell-specific contributions and the role of anesthetics to the neurovascular coupling

The neurovascular coupling (NVC) connects neuronal activity to hemodynamic responses in the brain. This connection is the basis for the interpretation of functional magnetic resonance imaging data. Despite the central role of this coupling, we lack detailed knowledge about cell-specific contributions and our knowledge about NVC is mainly based on animal experiments performed during anesthesia. Anesthetics are known to affect neuronal excitability, but how this affects the vessel diameters is not known. Due to the high complexity of NVC data, mathematical modeling is needed for a meaningful analysis. However, neither the relevant neuronal subtypes nor the effects of anesthetics are covered by current models. Here, we present a mathematical model including GABAergic interneurons and pyramidal neurons, as well as the effect of an anesthetic agent. The model is consistent with data from optogenetic experiments from both awake and anesthetized animals, and it correctly predicts data from experiments with different pharmacological modulators. The analysis suggests that no downstream anesthetic effects are necessary if one of the GABAergic interneuron signaling pathways include a Michaelis-Menten expression. This is the first example of a quantitative model that includes both the cell-specific contributions and the effect of an anesthetic agent on the NVC.

Real-Time Dual-Wavelength Time-Resolved Diffuse Optical Tomography System for Functional Brain Imaging Based on Probe-Hosted Silicon Photomultipliers

Near-infrared diffuse optical tomography is a non-invasive photonics-based imaging technology suited to functional brain imaging applications. Recent developments have proved that it is possible to build a compact time-domain diffuse optical tomography system based on silicon photomultipliers (SiPM) detectors. The system presented in this paper was equipped with the same eight SiPM probe-hosted detectors, but was upgraded with six injection fibers to shine the sample at several points. Moreover, an automatic switch was included enabling a complete measurement to be performed in less than one second. Further, the system was provided with a dual-wavelength (670 nm and 820 nm) light source to quantify the oxy- and deoxy-hemoglobin concentration evolution in the tissue. This novel system was challenged against a solid phantom experiment, and two in-vivo tests, namely arm occlusion and motor cortex brain activation. The results show that the tomographic system makes it possible to follow the evolution of brain activation over time with a 1s-resolution.

Dynamic cerebral autoregulation estimates derived from near infrared spectroscopy and transcranial Doppler are similar after correction for transit time and blood flow and blood volume oscillations

We analysed mean arterial blood pressure, cerebral blood flow velocity, oxygenated haemoglobin and deoxygenated haemoglobin signals to estimate dynamic cerebral autoregulation. We compared macrovascular (mean arterial blood pressure-cerebral blood flow velocity) and microvascular (oxygenated haemoglobin-deoxygenated haemoglobin) dynamic cerebral autoregulation estimates during three different conditions: rest, mild hypocapnia and hypercapnia. Microvascular dynamic cerebral autoregulation estimates were created by introducing the constant time lag plus constant phase shift model, which enables correction for transit time, blood flow and blood volume oscillations (TT-BF/BV correction). After TT-BF/BV correction, a significant agreement between mean arterial blood pressure-cerebral blood flow velocity and oxygenated haemoglobin-deoxygenated haemoglobin phase differences in the low frequency band was found during rest (left: intraclass correlation=0.6, median phase difference 29.5° vs. 30.7°, right: intraclass correlation=0.56, median phase difference 32.6° vs. 39.8°) and mild hypocapnia (left: intraclass correlation=0.73, median phase difference 48.6° vs. 43.3°, right: intraclass correlation=0.70, median phase difference 52.1° vs. 61.8°). During hypercapnia, the mean transit time decreased and blood volume oscillations became much more prominent, except for very low frequencies. The transit time related to blood flow oscillations was remarkably stable during all conditions. We conclude that non-invasive microvascular dynamic cerebral autoregulation estimates are similar to macrovascular dynamic cerebral autoregulation estimates, after TT-BF/BV correction is applied. These findings may increase the feasibility of non-invasive continuous autoregulation monitoring and guided therapy in clinical situations.

A review on low-dimensional physics-based models of systemic arteries: application to estimation of central aortic pressure

The physiological processes and mechanisms of an arterial system are complex and subtle. Physics-based models have been proven to be a very useful tool to simulate actual physiological behavior of the arteries. The current physics-based models include high-dimensional models (2D and 3D models) and low-dimensional models (0D, 1D and tube-load models). High-dimensional models can describe the local hemodynamic information of arteries in detail. With regard to an exact model of the whole arterial system, a high-dimensional model is computationally impracticable since the complex geometry, viscosity or elastic properties and complex vectorial output need to be provided. For low-dimensional models, the structure, centerline and viscosity or elastic properties only need to be provided. Therefore, low-dimensional modeling with lower computational costs might be a more applicable approach to represent hemodynamic properties of the entire arterial system and these three types of low-dimensional models have been extensively used in the study of cardiovascular dynamics. In recent decades, application of physics-based models to estimate central aortic pressure has attracted increasing interest. However, to our best knowledge, there has been few review paper about reconstruction of central aortic pressure using these physics-based models. In this paper, three types of low-dimensional physical models (0D, 1D and tube-load models) of systemic arteries are reviewed, the application of three types of models on estimation of central aortic pressure is taken as an example to discuss their advantages and disadvantages, and the proper choice of models for specific researches and applications are advised.

Validation and optimization of hypercapnic-calibrated fMRI from oxygen-sensitive two-photon microscopy

Hypercapnic-calibrated fMRI allows the estimation of the relative changes in the cerebral metabolic rate of oxygen (rCMRO2) 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'.

Temporal Cortex Activation to Audiovisual Speech in Normal-Hearing and Cochlear Implant Users Measured with Functional Near-Infrared Spectroscopy

BACKGROUND

Speech understanding may rely not only on auditory, but also on visual information. Non-invasive functional neuroimaging techniques can expose the neural processes underlying the integration of multisensory processes required for speech understanding in humans. Nevertheless, noise (from functional MRI, fMRI) limits the usefulness in auditory experiments, and electromagnetic artifacts caused by electronic implants worn by subjects can severely distort the scans (EEG, fMRI). Therefore, we assessed audio-visual activation of temporal cortex with a silent, optical neuroimaging technique: functional near-infrared spectroscopy (fNIRS).

METHODS

We studied temporal cortical activation as represented by concentration changes of oxy- and deoxy-hemoglobin in four, easy-to-apply fNIRS optical channels of 33 normal-hearing adult subjects and five post-lingually deaf cochlear implant (CI) users in response to supra-threshold unisensory auditory and visual, as well as to congruent auditory-visual speech stimuli.

RESULTS

Activation effects were not visible from single fNIRS channels. However, by discounting physiological noise through reference channel subtraction (RCS), auditory, visual and audiovisual (AV) speech stimuli evoked concentration changes for all sensory modalities in both cohorts (p < 0.001). Auditory stimulation evoked larger concentration changes than visual stimuli (p < 0.001). A saturation effect was observed for the AV condition.

CONCLUSIONS

Physiological, systemic noise can be removed from fNIRS signals by RCS. The observed multisensory enhancement of an auditory cortical channel can be plausibly described by a simple addition of the auditory and visual signals with saturation.

Nonlinear extension of a hemodynamic linear model for coherent hemodynamics spectroscopy

In this work, we are proposing an extension of a recent hemodynamic model (Fantini, 2014a), which was developed within the framework of a novel approach to the study of tissue hemodynamics, named coherent hemodynamics spectroscopy (CHS). The previous hemodynamic model, from a signal processing viewpoint, treats the tissue microvasculature as a linear time-invariant system, and considers changes of blood volume, capillary blood flow velocity and the rate of oxygen diffusion as inputs, and the changes of oxy-, deoxy-, and total hemoglobin concentrations (measured in near infrared spectroscopy) as outputs. The model has been used also as a forward solver in an inversion procedure to retrieve quantitative parameters that assess physiological and biological processes such as microcirculation, cerebral autoregulation, tissue metabolic rate of oxygen, and oxygen extraction fraction. Within the assumption of "small" capillary blood flow velocity oscillations the model showed that the capillary and venous compartments "respond" to this input as low pass filters, characterized by two distinct impulse response functions. In this work, we do not make the assumption of "small" perturbations of capillary blood flow velocity by solving without approximations the partial differential equation that governs the spatio-temporal behavior of hemoglobin saturation in capillary and venous blood. Preliminary comparison between the linear time-invariant model and the extended model (here identified as nonlinear model) are shown for the relevant parameters measured in CHS as a function of the oscillation frequency (CHS spectra). We have found that for capillary blood flow velocity oscillations with amplitudes up to 10% of the baseline value (which reflect typical scenarios in CHS), the discrepancies between CHS spectra obtained with the linear and nonlinear models are negligible. For larger oscillations (~50%) the linear and nonlinear models yield CHS spectra with differences within typical experimental errors, but further investigation is needed to assess the effect of these differences. Flow oscillations larger than 10-20% are not typically induced in CHS; therefore, the results presented in this work indicate that a linear hemodynamic model, combined with a method to elicit controlled hemodynamic oscillations (as done for CHS), is appropriate for the quantitative assessment of cerebral microcirculation.

Hematocrit distribution and tissue oxygenation in large microcirculatory networks

OBJECTIVE

Oxygen tension in the brain is controlled by the microcirculatory supply of RBC, but the effect of non-Newtonian blood flow rheology on tissue oxygenation is not well characterized. This study assesses different biphasic blood flow models for predicting tissue oxygen tension as a function of microcirculatory hemodynamics.

METHODS

Two existing plasma-skimming laws are compared against measured RBC distributions in rat and hamster microcirculatory networks. A novel biphasic blood flow model is introduced. The computational models predict tissue oxygenation in the mesentery, cremaster muscle, and the human secondary cortex.

RESULTS

This investigation shows deficiencies in prior models, including inconsistent plasma-skimming trends and insufficient oxygen perfusion due to the high prevalence (33%) of RBC-free microvessels. Our novel method yields physiologically sound RBC distributions and tissue oxygen tensions within one standard deviation of experimental measurements.

CONCLUSIONS

A simple, novel biphasic blood flow model is introduced with equal or better predictive power when applied to historic raw data sets. It can overcome limitations of prior models pertaining to trifurcations, anastomoses, and loops. This new plasma-skimming law eases the computations of bulk blood flow and hematocrit fields in large microcirculatory networks and converges faster than prior procedures.

Dynamic causal modelling for functional near-infrared spectroscopy

Functional near-infrared spectroscopy (fNIRS) is an emerging technique for measuring changes in cerebral hemoglobin concentration via optical absorption changes. Although there is great interest in using fNIRS to study brain connectivity, current methods are unable to infer the directionality of neuronal connections. In this paper, we apply Dynamic Causal Modelling (DCM) to fNIRS data. Specifically, we present a generative model of how observed fNIRS data are caused by interactions among hidden neuronal states. Inversion of this generative model, using an established Bayesian framework (variational Laplace), then enables inference about changes in directed connectivity at the neuronal level. Using experimental data acquired during motor imagery and motor execution tasks, we show that directed (i.e., effective) connectivity from the supplementary motor area to the primary motor cortex is negatively modulated by motor imagery, and this suppressive influence causes reduced activity in the primary motor cortex during motor imagery. These results are consistent with findings of previous functional magnetic resonance imaging (fMRI) studies, suggesting that the proposed method enables one to infer directed interactions in the brain mediated by neuronal dynamics from measurements of optical density changes.

Resting-state functional magnetic resonance imaging: review of neurosurgical applications

Recent research in brain imaging has highlighted the role of different neural networks in the resting state (ie, no task) in which the brain displays spontaneous low-frequency neuronal oscillations. These can be indirectly measured with resting-state functional magnetic resonance imaging, and functional connectivity can be inferred as the spatiotemporal correlations of this signal. This technique has proliferated in recent years and has allowed the noninvasive investigation of large-scale, distributed functional networks. In this review, we give a brief overview of resting-state networks and examine the use of resting-state functional magnetic resonance imaging in neurosurgical contexts, specifically with respect to neurooncology, epilepsy surgery, and deep brain stimulation. We discuss the advantages and disadvantages compared with task-based functional magnetic resonance imaging, the limitations of resting-state functional magnetic resonance imaging, and the emerging directions of this relatively new technology.

Social neuroscience and hyperscanning techniques: past, present and future

This paper reviews the published literature on the hyperscanning methodologies using hemodynamic or neuro-electric modalities. In particular, we describe how different brain recording devices have been employed in different experimental paradigms to gain information about the subtle nature of human interactions. This review also included papers based on single-subject recordings in which a correlation was found between the activities of different (non-simultaneously recorded) participants in the experiment. The descriptions begin with the methodological issues related to the simultaneous measurements and the descriptions of the results generated by such approaches will follow. Finally, a discussion of the possible future uses of such new approaches to explore human social interactions will be presented.

Further improvement in reducing superficial contamination in NIRS using double short separation measurements

Near-Infrared Spectroscopy (NIRS) allows the recovery of the evoked hemodynamic response to brain activation. In adult human populations, the NIRS signal is strongly contaminated by systemic interference occurring in the superficial layers of the head. An approach to overcome this difficulty is to use additional NIRS measurements with short optode separations to measure the systemic hemodynamic fluctuations occurring in the superficial layers. These measurements can then be used as regressors in the post-experiment analysis to remove the systemic contamination and isolate the brain signal. In our previous work, we showed that the systemic interference measured in NIRS is heterogeneous across the surface of the scalp. As a consequence, the short separation measurement used in the regression procedure must be located close to the standard NIRS channel from which the evoked hemodynamic response of the brain is to be recovered. Here, we demonstrate that using two short separation measurements, one at the source optode and one at the detector optode, further increases the performance of the short separation regression method compared to using a single short separation measurement. While a single short separation channel produces an average reduction in noise of 33% for HbO, using a short separation channel at both source and detector reduces noise by 59% compared to the standard method using a general linear model (GLM) without short separation. For HbR, noise reduction of 3% is achieved using a single short separation and this number goes to 47% when two short separations are used. Our work emphasizes the importance of integrating short separation measurements both at the source and at the detector optode of the standard channels from which the hemodynamic response is to be recovered. While the implementation of short separation sources presents some difficulties experimentally, the improvement in noise reduction is significant enough to justify the practical challenges.

Further improvement in reducing superficial contamination in NIRS using double short separation measurements

Near-Infrared Spectroscopy (NIRS) allows the recovery of the evoked hemodynamic response to brain activation. In adult human populations, the NIRS signal is strongly contaminated by systemic interference occurring in the superficial layers of the head. An approach to overcome this difficulty is to use additional NIRS measurements with short optode separations to measure the systemic hemodynamic fluctuations occurring in the superficial layers. These measurements can then be used as regressors in the post-experiment analysis to remove the systemic contamination and isolate the brain signal. In our previous work, we showed that the systemic interference measured in NIRS is heterogeneous across the surface of the scalp. As a consequence, the short separation measurement used in the regression procedure must be located close to the standard NIRS channel from which the evoked hemodynamic response of the brain is to be recovered. Here, we demonstrate that using two short separation measurements, one at the source optode and one at the detector optode, further increases the performance of the short separation regression method compared to using a single short separation measurement. While a single short separation channel produces an average reduction in noise of 33% for HbO, using a short separation channel at both source and detector reduces noise by 59% compared to the standard method using a general linear model (GLM) without short separation. For HbR, noise reduction of 3% is achieved using a single short separation and this number goes to 47% when two short separations are used. Our work emphasizes the importance of integrating short separation measurements both at the source and at the detector optode of the standard channels from which the hemodynamic response is to be recovered. While the implementation of short separation sources presents some difficulties experimentally, the improvement in noise reduction is significant enough to justify the practical challenges.

Calibrating the BOLD signal during a motor task using an extended fusion model incorporating DOT, BOLD and ASL data

Multimodal imaging improves the accuracy of the localization and the quantification of brain activation when measuring different manifestations of the hemodynamic response associated with cerebral activity. In this study, we incorporated cerebral blood flow (CBF) changes measured with arterial spin labeling (ASL), Diffuse Optical Tomography (DOT) and blood oxygen level-dependent (BOLD) recordings to reconstruct changes in oxy- (ΔHbO(2)) and deoxyhemoglobin (ΔHbR). Using the Grubb relation between relative changes in CBF and cerebral blood volume (CBV), we incorporated the ASL measurement as a prior to the total hemoglobin concentration change (ΔHbT). We applied this ASL fusion model to both synthetic data and experimental multimodal recordings during a 2-s finger-tapping task. Our results show that the new approach is very powerful in estimating ΔHbO(2) and ΔHbR with high spatial and quantitative accuracy. Moreover, our approach allows the computation of baseline total hemoglobin concentration (HbT(0)) as well as of the BOLD calibration factor M on a single subject basis. We obtained an average HbT(0) of 71 μM, an average M value of 0.18 and an average increase of 13% in cerebral metabolic rate of oxygen (CMRO(2)), all of which are in agreement with values previously reported in the literature. Our method yields an independent measurement of M, which provides an alternative measurement to validate the hypercapnic calibration of the BOLD signal.

Evolution of the dynamic changes in functional cerebral oxidative metabolism from tissue mitochondria to blood oxygen

The dynamic properties of the cerebral metabolic rate of oxygen consumption (CMR(O2)) during changes in brain activity remain unclear. Therefore, the spatial and temporal evolution of functional increases in CMR(O2) was investigated in the rat somato-sensory cortex during forelimb stimulation under a suppressed blood flow response condition. Temporally, stimulation elicited a fast increase in tissue mitochondria CMR(O2) described by a time constant of ~1 second measured using flavoprotein autofluorescence imaging. CMR(O2)-driven changes in the tissue oxygen tension measured using an oxygen electrode and blood oxygenation measured using optical imaging of intrinsic signal followed; however, these changes were slow with time constants of ~5 and ~10 seconds, respectively. This slow change in CMR(O2)-driven blood oxygenation partly explains the commonly observed post-stimulus blood oxygen level-dependent (BOLD) undershoot. Spatially, the changes in mitochondria CMR(O2) were similar to the changes in blood oxygenation. Finally, the increases in CMR(O2) were well correlated with the evoked multi-unit spiking activity. These findings show that dynamic CMR(O2) calculations made using only blood oxygenation data (e.g., BOLD functional magnetic resonance imaging (fMRI)) do not directly reflect the temporal changes in the tissue's mitochondria metabolic rate; however, the findings presented can bridge the gap between the changes in cellular oxidative rate and blood oxygenation.

Quantitative measurement of cerebral physiology using respiratory-calibrated MRI

Functional magnetic resonance imaging typically measures signal increases arising from changes in the transverse relaxation rate over small regions of the brain and associates these with local changes in cerebral blood flow, blood volume and oxygen metabolism. Recent developments in pulse sequences and image analysis methods have improved the specificity of the measurements by focussing on changes in blood flow or changes in blood volume alone. However, FMRI is still unable to match the physiological information obtainable from positron emission tomography (PET), which is capable of quantitative measurements of blood flow and volume, and can indirectly measure resting metabolism. The disadvantages of PET are its cost, its availability, its poor spatial resolution and its use of ionising radiation. The MRI techniques introduced here address some of these limitations and provide physiological data comparable with PET measurements. We present an 18-minute MRI protocol that produces multi-slice whole-brain coverage and yields quantitative images of resting cerebral blood flow, cerebral blood volume, oxygen extraction fraction, CMRO(2), arterial arrival time and cerebrovascular reactivity of the human brain in the absence of any specific functional task. The technique uses a combined hyperoxia and hypercapnia paradigm with a modified arterial spin labelling sequence.

Nonlinear hemodynamic responses in human epilepsy: a multimodal analysis with fNIRS-EEG and fMRI-EEG

Functional magnetic resonance imaging (fMRI) combined with electroencephalography (fMRI-EEG) is a neuroimaging technique based on the blood oxygenation level dependent (BOLD) signal which has been shown to be useful in the study of epilepsy for the localization of the epileptogenic focus. Functional near-infrared spectroscopy (fNIRS) combined with EEG (fNIRS-EEG) is another imaging technique based on the measurement of oxygenated and deoxygenated hemoglobin with complementary clinical potential in epilepsy, for continuous patient monitoring, language lateralization, and focus localization. In this work fMRI-EEG and fNIRS-EEG are used to quantify nonlinear hemodynamic responses in three cases of human refractory focal epilepsy, by using the Volterra kernel expansion up to second order. Prior to analyzing real data, extensive simulations are carried out to show that nonlinearities are estimable. The Volterra methodology is then applied to multimodal data recorded from 3 epileptic patients selected for their frequent spiking activity. Care is taken to account for variability of hemodynamic responses due to other causes than Volterra nonlinearities. Statistically significant nonlinearities are observed for all patients and all modalities. Good concordance between fNIRS and fMRI is found for both the amplitude of the Volterra responses, and, with limitations, in the localization of the epileptic focus and regions of inverted responses (negative BOLD signals). In one patient, Volterra nonlinearities allowed epileptic focus identification with fMRI, while analyses without nonlinearities failed to see it. In simulations when nonlinearities were included, analysis without Volterra nonlinearities performed poorly. These two observations suggest routinely checking for nonlinearities in functional imaging of patients presenting with frequent spikes.

The utility of near-infrared spectroscopy in the regression of low-frequency physiological noise from functional magnetic resonance imaging data

Near-infrared spectroscopy (NIRS) signals have been shown to correlate with resting-state BOLD-fMRI data across the whole brain volume, particularly at frequencies below 0.1Hz. While the physiological origins of this correlation remain unclear, its existence may have a practical application in minimizing the background physiological noise present in BOLD-fMRI recordings. We performed simultaneous, resting-state fMRI and 28-channel NIRS in seven adult subjects in order to assess the utility of NIRS signals in the regression of physiological noise from fMRI data. We calculated the variance of the residual error in a general linear model of the baseline fMRI signal, and the reduction of this variance achieved by including NIRS signals in the model. In addition, we introduced a sequence of simulated hemodynamic response functions (HRFs) into the resting-state fMRI data of each subject in order to quantify the effectiveness of NIRS signals in optimizing the recovery of that HRF. For comparison, these calculations were also performed using a pulse and respiration RETROICOR model. Our results show that the use of 10 or more NIRS channels can reduce variance in the residual error by as much as 36% on average across the whole cortex. However the same number of low-pass filtered white noise regressors is shown to produce a reduction of 19%. The RETROICOR model obtained a variance reduction of 6.4%. Our HRF simulation showed that the mean-squared error (MSE) between the recovered and true HRFs is reduced by 21% on average when 10 NIRS channels are applied and by introducing an optimized time lag between the NIRS and fMRI time series, a single NIRS channel can provide an average MSE reduction of 14%. The RETROICOR model did not provide a significant change in MSE. By each of the metrics calculated, NIRS recording is shown to be of significant benefit to the regression of low-frequency physiological noise from fMRI data.

Quantification of the cortical contribution to the NIRS signal over the motor cortex using concurrent NIRS-fMRI measurements

Near-Infrared Spectroscopy (NIRS) measures the functional hemodynamic response occurring at the surface of the cortex. Large pial veins are located above the surface of the cerebral cortex. Following activation, these veins exhibit oxygenation changes but their volume likely stays constant. The back-reflection geometry of the NIRS measurement renders the signal very sensitive to these superficial pial veins. As such, the measured NIRS signal contains contributions from both the cortical region as well as the pial vasculature. In this work, the cortical contribution to the NIRS signal was investigated using (1) Monte Carlo simulations over a realistic geometry constructed from anatomical and vascular MRI and (2) multimodal NIRS-BOLD recordings during motor stimulation. A good agreement was found between the simulations and the modeling analysis of in vivo measurements. Our results suggest that the cortical contribution to the deoxyhemoglobin signal change (ΔHbR) is equal to 16-22% of the cortical contribution to the total hemoglobin signal change (ΔHbT). Similarly, the cortical contribution of the oxyhemoglobin signal change (ΔHbO) is equal to 73-79% of the cortical contribution to the ΔHbT signal. These results suggest that ΔHbT is far less sensitive to pial vein contamination and therefore, it is likely that the ΔHbT signal provides better spatial specificity and should be used instead of ΔHbO or ΔHbR to map cerebral activity with NIRS. While different stimuli will result in different pial vein contributions, our finger tapping results do reveal the importance of considering the pial contribution.

NIRS-based hyperscanning reveals increased interpersonal coherence in superior frontal cortex during cooperation

We used Near-Infrared Spectroscopy (NIRS) to simultaneously measure brain activity in two people while they played a computer-based cooperation game side by side. Inter-brain activity coherence was calculated between the two participants. We found that the coherence between signals generated by participants' right superior frontal cortices increased during cooperation, but not during competition. Increased coherence was also associated with better cooperation performance. To our knowledge, this work represents the first use of a single NIRS instrument for simultaneous measurements of brain activity in two people. This study demonstrates the use of NIRS-based hyperscanning in studies of social interaction in a naturalistic environment.

Improved recovery of the hemodynamic response in diffuse optical imaging using short optode separations and state-space modeling

Diffuse optical imaging (DOI) allows the recovery of the hemodynamic response associated with evoked brain activity. The signal is contaminated with systemic physiological interference which occurs in the superficial layers of the head as well as in the brain tissue. The back-reflection geometry of the measurement makes the DOI signal strongly contaminated by systemic interference occurring in the superficial layers. A recent development has been the use of signals from small source-detector separation (1cm) optodes as regressors. Since those additional measurements are mainly sensitive to superficial layers in adult humans, they help in removing the systemic interference present in longer separation measurements (3 cm). Encouraged by those findings, we developed a dynamic estimation procedure to remove global interference using small optode separations and to estimate simultaneously the hemodynamic response. The algorithm was tested by recovering a simulated synthetic hemodynamic response added over baseline DOI data acquired from 6 human subjects at rest. The performance of the algorithm was quantified by the Pearson R(2) coefficient and the mean square error (MSE) between the recovered and the simulated hemodynamic responses. Our dynamic estimator was also compared with a static estimator and the traditional adaptive filtering method. We observed a significant improvement (two-tailed paired t-test, p<0.05) in both HbO and HbR recovery using our Kalman filter dynamic estimator compared to the traditional adaptive filter, the static estimator and the standard GLM technique.

A quantitative comparison of NIRS and fMRI across multiple cognitive tasks

Near infrared spectroscopy (NIRS) is an increasingly popular technology for studying brain function. NIRS presents several advantages relative to functional magnetic resonance imaging (fMRI), such as measurement of concentration changes in both oxygenated and deoxygenated hemoglobin, finer temporal resolution, and ease of administration, as well as disadvantages, most prominently inferior spatial resolution and decreased signal-to-noise ratio (SNR). While fMRI has become the gold standard for in vivo imaging of the human brain, in practice NIRS is a more convenient and less expensive technology than fMRI. It is therefore of interest to many researchers how NIRS compares to fMRI in studies of brain function. In the present study we scanned participants with simultaneous NIRS and fMRI on a battery of cognitive tasks, placing NIRS probes over both frontal and parietal brain regions. We performed detailed comparisons of the signals in both temporal and spatial domains. We found that NIRS signals have significantly weaker SNR, but are nonetheless often highly correlated with fMRI measurements. Both SNR and the distance between the scalp and the brain contributed to variability in the NIRS/fMRI correlations. In the spatial domain, we found that a photon path forming an ellipse between the NIRS emitter and detector correlated most strongly with the BOLD response. Taken together these findings suggest that, while NIRS can be an appropriate substitute for fMRI for studying brain activity related to cognitive tasks, care should be taken when designing studies with NIRS to ensure that: 1) the spatial resolution is adequate for answering the question of interest and 2) the design accounts for weaker SNR, especially in brain regions more distal from the scalp.

A random-effects model for group-level analysis of diffuse optical brain imaging

Diffuse optical imaging is a non-invasive technique for measuring changes in blood oxygenation in the brain. This technique is based on the temporally and spatially resolved recording of optical absorption in tissue within the near-infrared range of light. Optical imaging can be used to study functional brain activity similar to functional MRI. However, group level comparisons of brain activity from diffuse optical data are difficult due to registration of optical sensors between subjects. In addition, optical signals are sensitive to inter-subject differences in cranial anatomy and the specific arrangement of optical sensors relative to the underlying functional region. These factors can give rise to partial volume errors and loss of sensitivity and therefore must be accounted for in combining data from multiple subjects. In this work, we describe an image reconstruction approach using a parametric Bayesian model that simultaneously reconstructs group-level images of brain activity in the context of a random-effects analysis. Using this model, we demonstrate that localization accuracy and the statistical effects size of group-level reconstructions can be improved when compared to individualized reconstructions. In this model, we use the Restricted Maximum Likelihood (ReML) method to optimize a Bayesian random-effects model.

A quantitative comparison of NIRS and fMRI across multiple cognitive tasks

Near infrared spectroscopy (NIRS) is an increasingly popular technology for studying brain function. NIRS presents several advantages relative to functional magnetic resonance imaging (fMRI), such as measurement of concentration changes in both oxygenated and deoxygenated hemoglobin, finer temporal resolution, and ease of administration, as well as disadvantages, most prominently inferior spatial resolution and decreased signal-to-noise ratio (SNR). While fMRI has become the gold standard for in vivo imaging of the human brain, in practice NIRS is a more convenient and less expensive technology than fMRI. It is therefore of interest to many researchers how NIRS compares to fMRI in studies of brain function. In the present study we scanned participants with simultaneous NIRS and fMRI on a battery of cognitive tasks, placing NIRS probes over both frontal and parietal brain regions. We performed detailed comparisons of the signals in both temporal and spatial domains. We found that NIRS signals have significantly weaker SNR, but are nonetheless often highly correlated with fMRI measurements. Both SNR and the distance between the scalp and the brain contributed to variability in the NIRS/fMRI correlations. In the spatial domain, we found that a photon path forming an ellipse between the NIRS emitter and detector correlated most strongly with the BOLD response. Taken together these findings suggest that, while NIRS can be an appropriate substitute for fMRI for studying brain activity related to cognitive tasks, care should be taken when designing studies with NIRS to ensure that: 1) the spatial resolution is adequate for answering the question of interest and 2) the design accounts for weaker SNR, especially in brain regions more distal from the scalp.

Hierarchical Bayesian regularization of reconstructions for diffuse optical tomography using multiple priors

Diffuse optical tomography (DOT) is a non-invasive brain imaging technique that uses low-levels of near-infrared light to measure optical absorption changes due to regional blood flow and blood oxygen saturation in the brain. By arranging light sources and detectors in a grid over the surface of the scalp, DOT studies attempt to spatially localize changes in oxy- and deoxy-hemoglobin in the brain that result from evoked brain activity during functional experiments. However, the reconstruction of accurate spatial images of hemoglobin changes from DOT data is an ill-posed linearized inverse problem, which requires model regularization to yield appropriate solutions. In this work, we describe and demonstrate the application of a parametric restricted maximum likelihood method (ReML) to incorporate multiple statistical priors into the recovery of optical images. This work is based on similar methods that have been applied to the inverse problem for magnetoencephalography (MEG). Herein, we discuss the adaptation of this model to DOT and demonstrate that this approach provides a means to objectively incorporate reconstruction constraints and demonstrate this approach through a series of simulated numerical examples.

Neuroimaging with near-infrared spectroscopy demonstrates speech-evoked activity in the auditory cortex of deaf children following cochlear implantation

Cochlear implants (CI) are commonly used to treat deafness in young children. While many factors influence the ability of a deaf child who is hearing through a CI to develop speech and language skills, an important factor is that the CI has to stimulate the auditory cortex. Obtaining behavioral measurements from young children with CIs can often be unreliable. While a variety of noninvasive techniques can be used for detecting cortical activity in response to auditory stimuli, many have critical limitations when applied to the pediatric CI population. We tested the ability of near-infrared spectroscopy (NIRS) to detect cortical responses to speech stimuli in pediatric CI users. Neuronal activity leads to changes in blood oxy- and deoxy-hemoglobin concentrations that can be detected by measuring the transmission of near-infrared light through the tissue. To verify the efficacy of NIRS, we first compared auditory cortex responses measured with NIRS and fMRI in normal-hearing adults. We then examined four different participant cohorts with NIRS alone. Speech-evoked cortical activity was observed in 100% of normal-hearing adults (11 of 11), 82% of normal-hearing children (9 of 11), 78% of deaf children who have used a CI > 4 months (28 of 36), and 78% of deaf children who completed NIRS testing on the day of CI initial activation (7 of 9). Therefore, NIRS can measure cortical responses in pediatric CI users, and has the potential to be a powerful adjunct to current CI assessment tools.

Diffuse Optics for Tissue Monitoring and Tomography

This review describes the diffusion model for light transport in tissues and the medical applications of diffuse light. Diffuse optics is particularly useful for measurement of tissue hemodynamics, wherein quantitative assessment of oxy- and deoxy-hemoglobin concentrations and blood flow are desired. The theoretical basis for near-infrared or diffuse optical spectroscopy (NIRS or DOS, respectively) is developed, and the basic elements of diffuse optical tomography (DOT) are outlined. We also discuss diffuse correlation spectroscopy (DCS), a technique whereby temporal correlation functions of diffusing light are transported through tissue and are used to measure blood flow. Essential instrumentation is described, and representative brain and breast functional imaging and monitoring results illustrate the workings of these new tissue diagnostics.

Interpreting oxygenation-based neuroimaging signals: the importance and the challenge of understanding brain oxygen metabolism

Functional magnetic resonance imaging is widely used to map patterns of brain activation based on blood oxygenation level dependent (BOLD) signal changes associated with changes in neural activity. However, because oxygenation changes depend on the relative changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen (CMRO(2)), a quantitative interpretation of BOLD signals, and also other functional neuroimaging signals related to blood or tissue oxygenation, is fundamentally limited until we better understand brain oxygen metabolism and how it is related to blood flow. However, the positive side of the complexity of oxygenation signals is that when combined with dynamic CBF measurements they potentially provide the best tool currently available for investigating the dynamics of CMRO(2). This review focuses on the problem of interpreting oxygenation-based signals, the challenges involved in measuring CMRO(2) in general, and what is needed to put oxygenation-based estimates of CMRO(2) on a firm foundation. The importance of developing a solid theoretical framework is emphasized, both as an essential tool for analyzing oxygenation-based multimodal measurements, and also potentially as a way to better understand the physiological phenomena themselves. The existing data, integrated within a simple theoretical framework of O(2) transport, suggests the hypothesis that an important functional role of the mismatch of CBF and CMRO(2) changes with neural activation is to prevent a fall of tissue pO(2). Future directions for better understanding brain oxygen metabolism are discussed.

Quantification of CMRO(2) without hypercapnia using simultaneous near-infrared spectroscopy and fMRI measurements

Estimation of the cerebral metabolic rate of oxygen (CMRO(2)) and cerebral blood flow (CBF) is important to investigate the neurovascular coupling and physiological components in blood oxygenation level-dependent (BOLD) signals quantitatively. Although there are methods that can determine CMRO(2) changes using functional MRI (fMRI) or near-infrared spectroscopy (NIRS), current approaches require a separate hypercapnia calibration process and have the potential to incur bias in many assumed model parameters. In this paper, a novel method to estimate CMRO(2) without hypercapnia is described using simultaneous measurements of NIRS and fMRI. Specifically, an optimization framework is proposed that minimizes the differences between the two forms of the relative CMRO(2)-CBF coupling ratio from BOLD and NIRS biophysical models, from which hypercapnia calibration and model parameters are readily estimated. Based on the new methods, we found that group average CBF, CMRO(2), cerebral blood volume (CBV), and BOLD changes within activation of the primary motor cortex during a finger tapping task increased by 39.5 +/- 21.4%, 18.4 +/- 8.7%, 12.9 +/- 6.7%, and 0.5 +/- 0.2%, respectively. The group average estimated flow-metabolism coupling ratio was 2.38 +/- 0.65 and the hypercapnia parameter was 7.7 +/- 1.7%. These values are within the range of values reported from other literatures. Furthermore, the activation maps from CBF and CMRO(2) were well localized on the primary motor cortex, which is the main target region of the finger tapping task.

HomER: a review of time-series analysis methods for near-infrared spectroscopy of the brain

Near-infrared spectroscopy (NIRS) is a noninvasive neuroimaging tool for studying evoked hemodynamic changes within the brain. By this technique, changes in the optical absorption of light are recorded over time and are used to estimate the functionally evoked changes in cerebral oxyhemoglobin and deoxyhemoglobin concentrations that result from local cerebral vascular and oxygen metabolic effects during brain activity. Over the past three decades this technology has continued to grow, and today NIRS studies have found many niche applications in the fields of psychology, physiology, and cerebral pathology. The growing popularity of this technique is in part associated with a lower cost and increased portability of NIRS equipment when compared with other imaging modalities, such as functional magnetic resonance imaging and positron emission tomography. With this increasing number of applications, new techniques for the processing, analysis, and interpretation of NIRS data are continually being developed. We review some of the time-series and functional analysis techniques that are currently used in NIRS studies, we describe the practical implementation of various signal processing techniques for removing physiological, instrumental, and motion-artifact noise from optical data, and we discuss the unique aspects of NIRS analysis in comparison with other brain imaging modalities. These methods are described within the context of the MATLAB-based graphical user interface program, HomER, which we have developed and distributed to facilitate the processing of optical functional brain data.

Oxygen advection and diffusion in a three- dimensional vascular anatomical network

There is an increasing need for quantitative and computationally affordable models for analyzing tissue metabolism and hemodynamics in microvascular networks. In this work, we develop a hybrid model to solve for the time-varying oxygen advection-diffusion equation in the vessels and tissue. To obtain a three-dimensional temporal evolution of tissue oxygen concentration for realistic complex vessel networks, we used a graph-based advection model combined with a finite-element based diffusion model and an implicit time-advancing scheme. We validated this algorithm for both static and dynamic conditions. We also applied it to a complex vascular network obtained from a rodent somatosensory cortex. Qualitative agreement was found with in-vivo experiments.

Oxygen advection and diffusion in a three- dimensional vascular anatomical network

There is an increasing need for quantitative and computationally affordable models for analyzing tissue metabolism and hemodynamics in microvascular networks. In this work, we develop a hybrid model to solve for the time-varying oxygen advection-diffusion equation in the vessels and tissue. To obtain a three-dimensional temporal evolution of tissue oxygen concentration for realistic complex vessel networks, we used a graph-based advection model combined with a finite-element based diffusion model and an implicit time-advancing scheme. We validated this algorithm for both static and dynamic conditions. We also applied it to a complex vascular network obtained from a rodent somatosensory cortex. Qualitative agreement was found with in-vivo experiments.