Model of the transient neurovascular response based on prompt arterial dilation

Global effects of aging on the hemodynamic response function in the human brain

In functional magnetic resonance imaging, the hemodynamic response function (HRF) is a transient, stereotypical response to local changes in cerebral hemodynamics and oxygen metabolism due to briefly (< 4 s) evoked neural activity. Accordingly, the HRF is often used as an impulse response with the assumption of linearity in data analysis. In cognitive aging studies, it has been very common to interpret differences in brain activation as age-related changes in neural activity. Contrary to this assumption, however, evidence has accrued that normal aging may also significantly affect the vasculature, thereby affecting cerebral hemodynamics and metabolism, confounding interpretation of fMRI aging studies. In this study, use was made of a multisensory stimulus to evoke the HRF in ~ 87% of cerebral cortex in cognitively intact adults with ages ranging from 22-75 years. The stimulus evokes both positive and negative HRFs, which were characterized using model-free parameters in native-space coordinates. Results showed significant age trends in HRF parameter distributions in terms of both amplitudes (e.g., peak amplitude and CNR) and temporal dynamics (e.g., full-width-at-half-maximum). This work sets the stage for using HRF methods as a biomarker for age-related pathology.

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 the blood oxygen level dependent hemodynamic response function in human subcortical regions with high spatiotemporal resolution

Subcortical brain regions are absolutely essential for normal human function. These phylogenetically early brain regions play critical roles in human behaviors such as the orientation of attention, arousal, and the modulation of sensory signals to cerebral cortex. Despite the critical health importance of subcortical brain regions, there has been a dearth of research on their neurovascular responses. Blood oxygen level dependent (BOLD) functional MRI (fMRI) experiments can help fill this gap in our understanding. The BOLD hemodynamic response function (HRF) evoked by brief (<4 s) neural activation is crucial for the interpretation of fMRI results because linear analysis between neural activity and the BOLD response relies on the HRF. Moreover, the HRF is a consequence of underlying local blood flow and oxygen metabolism, so characterization of the HRF enables understanding of neurovascular and neurometabolic coupling. We measured the subcortical HRF at 9.4T and 3T with high spatiotemporal resolution using protocols that enabled reliable delineation of HRFs in individual subjects. These results were compared with the HRF in visual cortex. The HRF was faster in subcortical regions than cortical regions at both field strengths. There was no significant undershoot in subcortical areas while there was a significant post-stimulus undershoot that was tightly coupled with its peak amplitude in cortex. The different BOLD temporal dynamics indicate different vascular dynamics and neurometabolic responses between cortex and subcortical nuclei.

Temporal stability of the hemodynamic response function across the majority of human cerebral cortex

The hemodynamic response function (HRF) measured with functional magnetic resonance imaging is generated by vascular and metabolic responses evoked by brief (<4 s) stimuli. It is known that the human HRF varies across cortex, between subjects, with stimulus paradigms, and even between different measurements in the same cortical location. However, our results demonstrate that strong HRFs are remarkably repeatable across sessions separated by time intervals up to 3 months. In this study, a multisensory stimulus was used to activate and measure the HRF across the majority of cortex (>70%, with lesser reliability observed in some areas of prefrontal cortex). HRFs were measured with high spatial resolution (2‐mm voxels) in central gray matter to minimize variations caused by partial‐volume effects. HRF amplitudes and temporal dynamics were highly repeatable across four sessions in 20 subjects. Positive and negative HRFs were consistently observed across sessions and subjects. Negative HRFs were generally weaker and, thus, more variable than positive HRFs. Statistical measurements showed that across‐session variability is highly correlated to the variability across events within a session; these measurements also indicated a normal distribution of variability across cortex. The overall repeatability of the HRFs over long time scales generally supports the long‐term use of event‐related functional magnetic resonance imaging protocols.

Imaging faster neural dynamics with fast fMRI: A need for updated models of the hemodynamic response

Fast fMRI enables the detection of neural dynamics over timescales of hundreds of milliseconds, suggesting it may provide a new avenue for studying subsecond neural processes in the human brain. The magnitudes of these fast fMRI dynamics are far greater than predicted by canonical models of the hemodynamic response. Several studies have established nonlinear properties of the hemodynamic response that have significant implications for fast fMRI. We first review nonlinear properties of the hemodynamic response function that may underlie fast fMRI signals. We then illustrate the breakdown of canonical hemodynamic response models in the context of fast neural dynamics. We will then argue that the canonical hemodynamic response function is not likely to reflect the BOLD response to neuronal activity driven by sparse or naturalistic stimuli or perhaps to spontaneous neuronal fluctuations in the resting state. These properties suggest that fast fMRI is capable of tracking surprisingly fast neuronal dynamics, and we discuss the neuroscientific questions that could be addressed using this approach.

A Network Architecture for Bidirectional Neurovascular Coupling in Rat Whisker Barrel Cortex

Neurovascular coupling is typically considered as a master-slave relationship between the neurons and the cerebral vessels: the neurons demand energy which the vessels supply in the form of glucose and oxygen. In the recent past, both theoretical and experimental studies have suggested that the neurovascular coupling is a bidirectional system, a loop that includes a feedback signal from the vessels influencing neural firing and plasticity. An integrated model of bidirectionally connected neural network and the vascular network is hence required to understand the relationship between the informational and metabolic aspects of neural dynamics. In this study, we present a computational model of the bidirectional neurovascular system in the whisker barrel cortex and study the effect of such coupling on neural activity and plasticity as manifest in the whisker barrel map formation. In this model, a biologically plausible self-organizing network model of rate coded, dynamic neurons is nourished by a network of vessels modeled using the biophysical properties of blood vessels. The neural layer which is designed to simulate the whisker barrel cortex of rat transmits vasodilatory signals to the vessels. The feedback from the vessels is in the form of available oxygen for oxidative metabolism whose end result is the adenosine triphosphate (ATP) necessary to fuel neural firing. The model captures the effect of the feedback from the vascular network on the neuronal map formation in the whisker barrel model under normal and pathological (Hypoxia and Hypoxia-Ischemia) conditions.

Retinotopic variations of the negative blood-oxygen-level dependent hemodynamic response function in human primary visual cortex

Functional magnetic resonance imaging (fMRI) measures blood-oxygen-level-dependent (BOLD) contrast that is generally assumed to be linearly related to excitatory neural activity. The positive hemodynamic response function (pHRF) is the positive BOLD response (PBR) evoked by a brief neural stimulation; the pHRF is often used as the impulse response for linear analysis of neural excitation. Many fMRI studies have observed a negative BOLD response (NBR) that is often associated with neural suppression. However, the temporal dynamics of the NBR evoked by a brief stimulus, the negative HRF (nHRF), remains unclear. Here, a unilateral visual stimulus was presented in a slow event-related design to elicit both pHRFs in the stimulus representation (SR), and nHRFs elsewhere. The observed nHRFs were not inverted versions of the pHRF previously reported. They were characterized by a stronger initial negative response followed by a significantly later positive peak. In contralateral primary visual cortex (V1), these differences varied with eccentricity from the SR. Similar nHRFs were observed in ipsilateral V1 with less eccentricity variation. Experiments with the blocked version of the same stimulus confirmed that brain regions presenting the unexpected nHRF dynamics correspond to those presenting a strong NBR. These data demonstrated that shift-invariant temporal linearity did not hold for the NBR while confirming that the PBR maintained rough linearity. Modeling indicated that the observed nHRFs can be created by suppression of both blood flow and oxygen metabolism. Critically, the nHRF can be misinterpreted as a pHRF due to their similarity, which could confound linear analysis for event-related fMRI experiments.NEW & NOTEWORTHY We investigate dynamics of the negative hemodynamic response function (nHRF), the negative blood-oxygen-level-dependent (BOLD) response (NBR) evoked by a brief stimulus, in human early visual cortex. Here, we show that the nHRFs are not inverted versions of the corresponding pHRFs. The nHRF has complex dynamics that varied significantly with eccentricity. The results also show shift-invariant temporal linearity does not hold for the NBR.

Dynamics of the cerebral blood flow response to brief neural activity in human visual cortex

The blood oxygen-level dependent (BOLD) functional magnetic resonance imaging (fMRI) signal depends on an interplay of cerebral blood flow (CBF), oxygen metabolism, and cerebral blood volume. Despite wide usage of BOLD fMRI, it is not clear how these physiological components create the BOLD signal. Here, baseline CBF and its dynamics evoked by a brief stimulus (2 s) in human visual cortex were measured at 3T. We found a stereotypical CBF response: immediate increase, rising to a peak a few second after the stimulus, followed by a significant undershoot. The BOLD hemodynamic response function (HRF) was also measured in the same session. Strong correlations between HRF and CBF peak responses indicate that the flow responses evoked by neural activation in nearby gray matter drive the early HRF. Remarkably, peak CBF and HRF were also strongly modulated by baseline perfusion. The CBF undershoot was reliable and significantly correlated with the HRF undershoot. However, late-time dynamics of the HRF and CBF suggest that oxygen metabolism can also contribute to the HRF undershoot. Combined measurement of the CBF and HRF for brief neural activation is a useful tool to understand the temporal dynamics of neurovascular and neurometabolic coupling.

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.

Characterization of the hemodynamic response function in white matter tracts for event-related fMRI

Accurate estimates of the BOLD hemodynamic response function (HRF) are crucial for the interpretation and analysis of event-related functional MRI data. To date, however, there have been no comprehensive measurements of the HRF in white matter (WM) despite increasing evidence that BOLD signals in WM change after a stimulus. We performed an event-related cognitive task (Stroop color-word interference) to measure the HRF in selected human WM pathways. The task was chosen in order to produce robust, distributed centers of activity throughout the cortex. To measure the HRF in WM, fiber tracts were reconstructed between each pair of activated cortical areas. We observed clear task-specific HRFs with reduced magnitudes, delayed onsets and prolonged initial dips in WM tracts compared with activated grey matter, thus calling for significant changes to current standard models for accurately characterizing the HRFs in WM and for modifications of standard methods of analysis of functional imaging data.

Existence of Initial Dip for BCI: An Illusion or Reality

A tight coupling between the neuronal activity and the cerebral blood flow (CBF) is the motivation of many hemodynamic response (HR)-based neuroimaging modalities. The increase in neuronal activity causes the increase in CBF that is indirectly measured by HR modalities. Upon functional stimulation, the HR is mainly categorized in three durations: (i) initial dip, (ii) conventional HR (i.e., positive increase in HR caused by an increase in the CBF), and (iii) undershoot. The initial dip is a change in oxygenation prior to any subsequent increase in CBF and spatially more specific to the site of neuronal activity. Despite additional evidence from various HR modalities on the presence of initial dip in human and animal species (i.e., cat, rat, and monkey); the existence/occurrence of an initial dip in HR is still under debate. This article reviews the existence and elusive nature of the initial dip duration of HR in intrinsic signal optical imaging (ISOI), functional magnetic resonance imaging (fMRI), and functional near-infrared spectroscopy (fNIRS). The advent of initial dip and its elusiveness factors in ISOI and fMRI studies are briefly discussed. Furthermore, the detection of initial dip and its role in brain-computer interface using fNIRS is examined in detail. The best possible application for the initial dip utilization and its future implications using fNIRS are provided.

Characterization of the hemodynamic response function across the majority of human cerebral cortex

A brief (<4 s) period of neural activation evokes a stereotypical sequence of vascular and metabolic events to create the hemodynamic response function (HRF) measured using functional magnetic resonance imaging (fMRI). Linear analysis of fMRI data requires that the HRF be treated as an impulse response, so the character and temporal stability of the HRF are critical issues. Here, a simple audiovisual stimulus combined with a fast-paced task was used to evoke a strong HRF across a majority, ∼77%, of cortex during a single scanning session. High spatiotemporal resolution (2-mm voxels, 1.25-s acquisition time) was used to focus HRF measurements specifically on the gray matter for whole brain. The majority of activated cortex responds with positive HRFs, while ∼27% responds with negative (inverted) HRFs. Spatial patterns of the HRF response amplitudes were found to be similar across subjects. Timing of the initial positive lobe of the HRF was relatively stable across the cortical surface with a mean of 6.1 ± 0.6 s across subjects, yet small but significant timing variations were also evident in specific regions of cortex. The results provide guidance for linear analysis of fMRI data. More importantly, this method provides a means to quantify neurovascular function across most of the brain, with potential clinical utility for the diagnosis of brain pathologies such as traumatic brain injury.

Alterations in neurovascular coupling following acute traumatic brain injury

Following acute traumatic brain injury (TBI), timely transport to a hospital can significantly improve the prognosis for recovery. There is, however, a dearth of quantitative biomarkers for brain injury that can be rapidly acquired and interpreted in active, field environments in which TBIs are frequently incurred. We explored potential functional indicators for TBI that can be noninvasively obtained through portable detection modalities, namely optical and electrophysiological approaches. By combining diffuse correlation spectroscopy with colocalized electrophysiological measurements in a mouse model of TBI, we observed concomitant alterations in sensory-evoked cerebral blood flow (CBF) and electrical potentials following controlled cortical impact. Injury acutely reduced the peak amplitude of both electrophysiological and CBF responses, which mostly recovered to baseline values within 30 min, and intertrial variability for these parameters was also acutely altered. Notably, the postinjury dynamics of the CBF overshoot and undershoot amplitudes differed significantly; whereas the amplitude of the initial peak of stimulus-evoked CBF recovered relatively rapidly, the ensuing undershoot did not appear to recover within 30 min of injury. Additionally, acute injury induced apparent low-frequency oscillatory behavior in CBF ([Formula: see text]). Histological assessment indicated that these physiological alterations were not associated with any major, persisting anatomical changes. Several time-domain features of the blood flow and electrophysiological responses showed strong correlations in recovery kinetics. Overall, our results reveal an array of stereotyped, injury-induced alterations in electrophysiological and hemodynamic responses that can be rapidly obtained using a combination of portable detection techniques.

The hemodynamic signal as a first-order low-pass temporal filter: Evidence and implications for neuroimaging studies

Neuronal activation triggers local changes in blood flow and hemoglobin oxygenation. These hemodynamic signals can be recorded through functional magnetic resonance imaging or intrinsic optical imaging, and allows inferring neural activity in response to stimuli. These techniques are widely used to uncover functional brain architectures. However, their accuracy suffers from distortions inherent to hemodynamic responses and noise. The analysis of these signals currently relies on models of impulse hemodynamic responses to brief stimuli. Here, in order to infer precise functional architectures, we focused on integrated signals associated to the dynamic response of functional maps. To this end, we recorded orientation and direction maps in cat primary visual cortex and compared two protocols: the conventional episodic stimulation technique and a continuous, periodic stimulation paradigm. Conventional methods show that the dynamics of activation and deactivation of the functional maps follows a linear first-order differential equation representing a low-pass filter. Comparison with the periodic stimulation methods confirmed this observation: the phase shifts and magnitude attenuations extracted at various frequencies were consistent with a low-pass filter with a 5s time constant. This dynamics presumably reflects the variations in deoxyhemoglobin mediated by arterial dilations. This dynamics open new avenues in the analysis of neuroimaging data that differs from common methods based on the hemodynamic response function. In particular, we demonstrate that inverting this first-order low-pass filter minimized the distortions of the signal and enabled a much faster and accurate reconstruction of functional maps.

Reliability of the depth-dependent high-resolution BOLD hemodynamic response in human visual cortex and vicinity

Functional magnetic resonance imaging (fMRI) often relies on a hemodynamic response function (HRF), the stereotypical blood oxygen level dependent (BOLD) response elicited by a brief (<4s) stimulus. Early measurements of the HRF used coarse spatial resolutions (≥3mm voxels) that would generally include contributions from white matter, gray matter, and the extra-pial compartment (the space between the pial surface and skull including pial blood vessels) within each voxel. To resolve these contributions, high-resolution fMRI (0.9-mm voxels) was performed at 3T in early visual cortex and its apposed white-matter and extra-pial compartments. The results characterized the depth dependence of the HRF and its reliability during nine fMRI sessions. Significant HRFs were observed in white-matter and extra-pial compartments as well as in gray matter. White-matter HRFs were faster and weaker than in the gray matter, while extra-pial HRFs were comparatively slower and stronger. Depth trends of the HRF peak amplitude were stable throughout a broad depth range that included all three compartments for each session. Across sessions, however, the depth trend of HRF peak amplitudes was stable only in the white matter and deep-intermediate gray matter, while there were strong session-to-session variations in the superficial gray matter and the extra-pial compartment. Thus, high-resolution fMRI can resolve significant and dynamically distinct HRFs in gray matter, white matter, and extra-pial compartments.

Mechanistic Mathematical Modeling Tests Hypotheses of the Neurovascular Coupling in fMRI

Functional magnetic resonance imaging (fMRI) measures brain activity by detecting the blood-oxygen-level dependent (BOLD) response to neural activity. The BOLD response depends on the neurovascular coupling, which connects cerebral blood flow, cerebral blood volume, and deoxyhemoglobin level to neuronal activity. The exact mechanisms behind this neurovascular coupling are not yet fully investigated. There are at least three different ways in which these mechanisms are being discussed. Firstly, mathematical models involving the so-called Balloon model describes the relation between oxygen metabolism, cerebral blood volume, and cerebral blood flow. However, the Balloon model does not describe cellular and biochemical mechanisms. Secondly, the metabolic feedback hypothesis, which is based on experimental findings on metabolism associated with brain activation, and thirdly, the neurotransmitter feed-forward hypothesis which describes intracellular pathways leading to vasoactive substance release. Both the metabolic feedback and the neurotransmitter feed-forward hypotheses have been extensively studied, but only experimentally. These two hypotheses have never been implemented as mathematical models. Here we investigate these two hypotheses by mechanistic mathematical modeling using a systems biology approach; these methods have been used in biological research for many years but never been applied to the BOLD response in fMRI. In the current work, model structures describing the metabolic feedback and the neurotransmitter feed-forward hypotheses were applied to measured BOLD responses in the visual cortex of 12 healthy volunteers. Evaluating each hypothesis separately shows that neither hypothesis alone can describe the data in a biologically plausible way. However, by adding metabolism to the neurotransmitter feed-forward model structure, we obtained a new model structure which is able to fit the estimation data and successfully predict new, independent validation data. These results open the door to a new type of fMRI analysis that more accurately reflects the true neuronal activity.

Functional magnetic resonance imaging (fMRI) is a widely used technique for measuring brain activity. However, the signal registered by fMRI is not a direct measurement of the neuronal activity in the brain, but it is influenced by the interplay between the metabolism, blood flow and blood volume in the active area. This signal is called the blood-oxygen-level dependent (BOLD) response and occurs when the blood supply to the active area increases in response to neuronal activity. The mechanisms that the cells use to influence the blood supply are not fully known, and therefore it is difficult to know the true neuronal signalling only from inspection of the fMRI signal. In this article, we present a new mathematical model built on the physiological mechanisms thought to underlie the BOLD response. We could successfully fit the model to data and predict the activity caused by new stimuli. By using the validated model we investigated physiological mechanisms that cause different parts of the BOLD response.

Investigating Human Neurovascular Coupling Using Functional Neuroimaging: A Critical Review of Dynamic Models

The mechanisms that link a transient neural activity to the corresponding increase of cerebral blood flow (CBF) are termed neurovascular coupling (NVC). They are possibly impaired at early stages of small vessel or neurodegenerative diseases. Investigation of NVC in humans has been made possible with the development of various neuroimaging techniques based on variations of local hemodynamics during neural activity. Specific dynamic models are currently used for interpreting these data that can include biophysical parameters related to NVC. After a brief review of the current knowledge about possible mechanisms acting in NVC we selected seven models with explicit integration of NVC found in the literature. All these models were described using the same procedure. We compared their physiological assumptions, mathematical formalism, and validation. In particular, we pointed out their strong differences in terms of complexity. Finally, we discussed their validity and their potential applications. These models may provide key information to investigate various aspects of NVC in human pathology.

Arterial impulse model for the BOLD response to brief neural activation

The blood oxygen level dependent (BOLD) signal evoked by brief neural stimulation, the hemodynamic response function (HRF), is a critical feature of neurovascular coupling. The HRF is directly related to local transient changes in oxygen supplied by cerebral blood flow (CBF) and oxygen demand, the cerebral metabolic rate of oxygen (CMRO2). Previous efforts to explain the HRF have relied upon the hypothesis that CBF produces a non-linear venous dilation within the cortical parenchyma. Instead, the observed dynamics correspond to prompt arterial dilation without venous volume change. This work develops an alternative biomechanical model for the BOLD response based on the hypothesis that prompt upstream dilation creates an arterial flow impulse amenable to linear description. This flow model is coupled to a continuum description of oxygen transport. Measurements using high-resolution fMRI demonstrate the efficacy of the model. The model predicts substantial spatial variations of the oxygen saturation along the length of capillaries and veins, and fits the varied gamut of measured HRFs by the combined effects of corresponding CBF and CMRO2 responses. Three interesting relationships among the hemodynamic parameters are predicted. First, there is an offset linear correlation with approximately unity slope between CBF and CMRO2 responses. Second, the HRF undershoot is strongly correlated to the corresponding CBF undershoot. Third, late-time-CMRO2 response can contribute to a slow recovery to baseline, lengthening the HRF undershoot. The model provides a powerful mathematical framework to understand the dynamics of neurovascular and neurometabolic responses that form the BOLD HRF.

Mechanical restriction of intracortical vessel dilation by brain tissue sculpts the hemodynamic response

Understanding the spatial dynamics of dilation in the cerebral vasculature is essential for deciphering the vascular basis of hemodynamic signals in the brain. We used two-photon microscopy to image neural activity and vascular dynamics in the somatosensory cortex of awake behaving mice during voluntary locomotion. Arterial dilations within the histologically-defined forelimb/hindlimb (FL/HL) representation were larger than arterial dilations in the somatosensory cortex immediately outside the FL/HL representation, demonstrating that the vascular response during natural behaviors was spatially localized. Surprisingly, we found that locomotion drove dilations in surface vessels that were nearly three times the amplitude of intracortical vessel dilations. The smaller dilations of the intracortical arterioles were not due to saturation of dilation. Anatomical imaging revealed that, unlike surface vessels, intracortical vessels were tightly enclosed by brain tissue. A mathematical model showed that mechanical restriction by the brain tissue surrounding intracortical vessels could account for the reduced amplitude of intracortical vessel dilation relative to surface vessels. Thus, under normal conditions, the mechanical properties of the brain may play an important role in sculpting the laminar differences of hemodynamic responses.

Coupling mechanism and significance of the BOLD signal: a status report

Functional magnetic resonance imaging (fMRI) provides a unique view of the working human mind. The blood-oxygen-level-dependent (BOLD) signal, detected in fMRI, reflects changes in deoxyhemoglobin driven by localized changes in brain blood flow and blood oxygenation, which are coupled to underlying neuronal activity by a process termed neurovascular coupling. Over the past 10 years, a range of cellular mechanisms, including astrocytes, pericytes, and interneurons, have been proposed to play a role in functional neurovascular coupling. However, the field remains conflicted over the relative importance of each process, while key spatiotemporal features of BOLD response remain unexplained. Here, we review current candidate neurovascular coupling mechanisms and propose that previously overlooked involvement of the vascular endothelium may provide a more complete picture of how blood flow is controlled in the brain. We also explore the possibility and consequences of conditions in which neurovascular coupling may be altered, including during postnatal development, pathological states, and aging, noting relevance to both stimulus-evoked and resting-state fMRI studies.

Modeling the role of osmotic forces in the cerebrovascular response to CO2

Increases in blood osmolarity have been shown to exert a vasodilatory effect on cerebral and other vasculature, with accompanying increases in blood flow. It has also been shown that, through an influence on blood concentration of the bicarbonate ion and pH, changes in blood levels of CO2 can alter blood osmolarity sufficiently to have an impact on vessel diameter. We propose here that this phenomenon plays a previously unappreciated role in CO2-mediated vasodilation, and present a biophysical model of osmotically driven vasodilation. Our model, which is based on literature data describing CO2-dependent changes in blood osmolarity and hydraulic conductivity (Lp) of the blood-brain barrier, is used to predict the change in cerebral blood flow (CBF) associated with osmotic forces arising from a specific hypercapnic challenge. Modeled changes were then compared with actual CBF changes determined using arterial spin-labeling (ASL) MRI. For changes in the arterial partial pressure of CO2 (PaCO2) of 20 mmHg, our model predicted increases of 80% from baseline CBF with a temporal evolution that was comparable to the measured hemodynamic responses. Our modeling results suggest that osmotic forces could play a significant role in the cerebrovascular response to CO2.

Quantitative separation of arterial and venous cerebral blood volume increases during voluntary locomotion

Voluntary locomotion is accompanied by large increases in cortical activity and localized increases in cerebral blood volume (CBV). We sought to quantitatively determine the spatial and temporal dynamics of voluntary locomotion-evoked cerebral hemodynamic changes. We measured single vessel dilations using two-photon microscopy and cortex-wide changes in CBV-related signal using intrinsic optical signal (IOS) imaging in head-fixed mice freely locomoting on a spherical treadmill. During bouts of locomotion, arteries dilated rapidly, while veins distended slightly and recovered slowly. The dynamics of diameter changes of both vessel types could be captured using a simple linear convolution model. Using these single vessel measurements, we developed a novel analysis approach to separate out spatially and temporally distinct arterial and venous components of the location-specific hemodynamic response functions (HRF) for IOS. The HRF of each pixel of was well fit by a sum of a fast arterial and a slow venous component. The HRFs of pixels in the limb representations of somatosensory cortex had a large arterial contribution, while in the frontal cortex the arterial contribution to the HRF was negligible. The venous contribution was much less localized, and was substantial in the frontal cortex. The spatial pattern and amplitude of these HRFs in response to locomotion in the cortex were robust across imaging sessions. Separating the more localized arterial component from the diffuse venous signals will be useful for dealing with the dynamic signals generated by naturalistic stimuli.