Arterial impulse model for the BOLD response to brief neural activation

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.

The confound of hemodynamic response function variability in human resting-state functional MRI studies

Functional magnetic resonance imaging (fMRI) is an indirect measure of neural activity with the hemodynamic response function (HRF) coupling it with unmeasured neural activity. The HRF, modulated by several non-neural factors, is variable across brain regions, individuals and populations. Yet, a majority of human resting-state fMRI connectivity studies continue to assume a non-variable HRF. In this article, with supportive prior evidence, we argue that HRF variability cannot be ignored as it substantially confounds within-subject connectivity estimates and between-subjects connectivity group differences. We also discuss its clinical relevance with connectivity impairments confounded by HRF aberrations in several disorders. We present limited data on HRF differences between women and men, which resulted in a 15.4% median error in functional connectivity estimates in a group-level comparison. We also discuss the implications of HRF variability for fMRI studies in the spinal cord. There is a need for more dialogue within the community on the HRF confound, and we hope that our article is a catalyst in the process.

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.

Mechanistic model for human brain metabolism and its connection to the neurovascular coupling

The neurovascular and neurometabolic couplings (NVC and NMC) connect cerebral activity, blood flow, and metabolism. This interconnection is used in for instance functional imaging, which analyses the blood-oxygen-dependent (BOLD) signal. The mechanisms underlying the NVC are complex, which warrants a model-based analysis of data. We have previously developed a mechanistically detailed model for the NVC, and others have proposed detailed models for cerebral metabolism. However, existing metabolic models are still not fully utilizing available magnetic resonance spectroscopy (MRS) data and are not connected to detailed models for NVC. Therefore, we herein present a new model that integrates mechanistic modelling of both MRS and BOLD data. The metabolic model covers central metabolism, using a minimal set of interactions, and can describe time-series data for glucose, lactate, aspartate, and glutamate, measured after visual stimuli. Statistical tests confirm that the model can describe both estimation data and predict independent validation data, not used for model training. The interconnected NVC model can simultaneously describe BOLD data and can be used to predict expected metabolic responses in experiments where metabolism has not been measured. This model is a step towards a useful and mechanistically detailed model for cerebral blood flow and metabolism, with potential applications in both basic research and clinical applications.

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.

Real-time tracking of brain oxygen gradients and blood flow during functional activation

Significance

Cerebral metabolic rate of oxygen ( CMRO 2 ) consumption is a key physiological variable that characterizes brain metabolism in a steady state and during functional activation.

Aim

We aim to develop a minimally invasive optical technique for real-time measurement of CMRO 2 concurrently with cerebral blood flow (CBF).

Approach

We used a pair of macromolecular phosphorescent probes with nonoverlapping optical spectra, which were localized in the intra- and extravascular compartments of the brain tissue, thus providing a readout of oxygen gradients between these two compartments. In parallel, we measured CBF using laser speckle contrast imaging.

Results

The method enables computation and tracking of CMRO 2 during functional activation with high temporal resolution ( ∼ 7    Hz ). In contrast to other approaches, our assessment of CMRO 2 does not require measurements of CBF or hemoglobin oxygen saturation.

Conclusions

The independent records of intravascular and extravascular partial pressures of oxygen, CBF, and CMRO 2 provide information about the physiological events that accompany neuronal activation, creating opportunities for dynamic quantification of brain metabolism.

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.

Extracting Dynamical Understanding From Neural-Mass Models of Mouse Cortex

New brain atlases with high spatial resolution and whole-brain coverage have rapidly advanced our knowledge of the brain's neural architecture, including the systematic variation of excitatory and inhibitory cell densities across the mammalian cortex. But understanding how the brain's microscale physiology shapes brain dynamics at the macroscale has remained a challenge. While physiologically based mathematical models of brain dynamics are well placed to bridge this explanatory gap, their complexity can form a barrier to providing clear mechanistic interpretation of the dynamics they generate. In this work, we develop a neural-mass model of the mouse cortex and show how bifurcation diagrams, which capture local dynamical responses to inputs and their variation across brain regions, can be used to understand the resulting whole-brain dynamics. We show that strong fits to resting-state functional magnetic resonance imaging (fMRI) data can be found in surprisingly simple dynamical regimes—including where all brain regions are confined to a stable fixed point—in which regions are able to respond strongly to variations in their inputs, consistent with direct structural connections providing a strong constraint on functional connectivity in the anesthetized mouse. We also use bifurcation diagrams to show how perturbations to local excitatory and inhibitory coupling strengths across the cortex, constrained by cell-density data, provide spatially dependent constraints on resulting cortical activity, and support a greater diversity of coincident dynamical regimes. Our work illustrates methods for visualizing and interpreting model performance in terms of underlying dynamical mechanisms, an approach that is crucial for building explanatory and physiologically grounded models of the dynamical principles that underpin large-scale brain activity.

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.

FMRI hemodynamic response function (HRF) as a novel marker of brain function: applications for understanding obsessive-compulsive disorder pathology and treatment response

The hemodynamic response function (HRF) represents the transfer function linking neural activity with the functional MRI (fMRI) signal, modeling neurovascular coupling. Since HRF is influenced by non-neural factors, to date it has largely been considered as a confound or has been ignored in many analyses. However, underlying biophysics suggests that the HRF may contain meaningful correlates of neural activity, which might be unavailable through conventional fMRI metrics. Here, we estimated the HRF by performing deconvolution on resting-state fMRI data from a longitudinal sample of 25 healthy controls scanned twice and 44 adults with obsessive-compulsive disorder (OCD) before and after 4-weeks of intensive cognitive-behavioral therapy (CBT). HRF response height, time-to-peak and full-width at half-maximum (FWHM) in OCD were abnormal before treatment and normalized after treatment in regions including the caudate. Pre-treatment HRF predicted treatment outcome (OCD symptom reduction) with 86.4% accuracy, using machine learning. Pre-treatment HRF response height in the caudate head and time-to-peak in the caudate tail were top-predictors of treatment response. Time-to-peak in the caudate tail, a region not typically identified in OCD studies using conventional fMRI activation or connectivity measures, may carry novel importance. Additionally, pre-treatment response height in caudate head predicted post-treatment OCD severity (R = -0.48, P = 0.001), and was associated with treatment-related OCD severity changes (R = -0.44, P = 0.0028), underscoring its relevance. With HRF being a reliable marker sensitive to brain function, OCD pathology, and intervention-related changes, these results could guide future studies towards novel discoveries not possible through conventional fMRI approaches like standard BOLD activation or connectivity.

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.

Developmental patterns of CBF and BOLD responses to visual stimulus

To investigate the developmental changes of cerebral blood flow (CBF) and hemodynamic responses to changing neural activity, we used the arterial spin label (ASL) technique to measure resting CBF and simultaneous CBF / blood-oxygen-level dependent (BOLD) signal changes during visual stimulation in 97 typically developing children and young adults (age 13.35 [6.02, 25.25] (median [min, max]) years old at the first time point). The longitudinal study protocol included three MRIs (2.7 ± 0.06 obtained), one year apart, for each participant. Mixed-effect linear and non-linear statistical models were used to analyze age effects on CBF and BOLD signals. Resting CBF decreased exponentially with age (p = 0.0001) throughout the brain, and developmental trajectories differed across brain lobes. The absolute CBF increase in visual cortex during stimulation was constant over the age range, but the fractional CBF change increased with age (p = 0.0001) and the fractional BOLD signal increased with age (p = 0.0001) correspondingly. These findings suggest that the apparent neural hemodynamic coupling in visual cortex does not change after age six years, but age-related BOLD signal changes continue through adolescence primarily due to the changes with age in resting CBF.

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.

Voxelized simulation of cerebral oxygen perfusion elucidates hypoxia in aged mouse cortex

Departures of normal blood flow and metabolite distribution from the cerebral microvasculature into neuronal tissue have been implicated with age-related neurodegeneration. Mathematical models informed by spatially and temporally distributed neuroimage data are becoming instrumental for reconstructing a coherent picture of normal and pathological oxygen delivery throughout the brain. Unfortunately, current mathematical models of cerebral blood flow and oxygen exchange become excessively large in size. They further suffer from boundary effects due to incomplete or physiologically inaccurate computational domains, numerical instabilities due to enormous length scale differences, and convergence problems associated with condition number deterioration at fine mesh resolutions. Our proposed simple finite volume discretization scheme for blood and oxygen microperfusion simulations does not require expensive mesh generation leading to the critical benefit that it drastically reduces matrix size and bandwidth of the coupled oxygen transfer problem. The compact problem formulation yields rapid and stable convergence. Moreover, boundary effects can effectively be suppressed by generating very large replica of the cortical microcirculation in silico using an image-based cerebrovascular network synthesis algorithm, so that boundaries of the perfusion simulations are far removed from the regions of interest. Massive simulations over sizeable portions of the cortex with feature resolution down to the micron scale become tractable with even modest computer resources. The feasibility and accuracy of the novel method is demonstrated and validated with in vivo oxygen perfusion data in cohorts of young and aged mice. Our oxygen exchange simulations quantify steep gradients near penetrating blood vessels and point towards pathological changes that might cause neurodegeneration in aged brains. This research aims to explain mechanistic interactions between anatomical structures and how they might change in diseases or with age. Rigorous quantification of age-related changes is of significant interest because it might aide in the search for imaging biomarkers for dementia and Alzheimer’s disease.

Brain function critically depends on the maintenance of physiological blood supply and metabolism in the cortex. Disturbances to adequate perfusion have been linked to age-related neurodegeneration. However, the precise correlation between age-related hemodynamic changes and the resulting decline in oxygen delivery is not well understood and has not been quantified. Therefore, we introduce a new compact, and therefore highly scalable, computational method for predicting the physiological relationship between hemodynamics and cortical oxygen perfusion for large sections of the cortical microcirculation. We demonstrate the novel mesh generation-free (MGF), multi-scale simulation approach through realistic in vivo case studies of cortical microperfusion in the mouse brain. We further validate mechanistic correlations and a quantitative relationship between blood flow and brain oxygenation using experimental data from cohorts of young, middle aged and old mouse brains. Our computational approach overcomes size and performance limitations of previous unstructured meshing techniques to enable the prediction of oxygen tension with a spatial resolution of least two orders of magnitude higher than previously possible. Our simulation results support the hypothesis that structural changes in the microvasculature induce hypoxic pockets in the aged brain that are absent in the healthy, young mouse.

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.

Axial variation of deoxyhemoglobin density as a source of the low-frequency time lag structure in blood oxygenation level-dependent signals

Perfusion-related information is reportedly embedded in the low-frequency component of a blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal. The blood-propagation pattern through the cerebral vascular tree is detected as an interregional lag variation of spontaneous low-frequency oscillations (sLFOs). Mapping of this lag, or phase, has been implicitly treated as a projection of the vascular tree structure onto real space. While accumulating evidence supports the biological significance of this signal component, the physiological basis of the “perfusion lag structure,” a requirement for an integrative resting-state fMRI-signal model, is lacking. In this study, we conducted analyses furthering the hypothesis that the sLFO is not only largely of systemic origin, but also essentially intrinsic to blood, and hence behaves as a virtual tracer. By summing the small fluctuations of instantaneous phase differences between adjacent vascular regions, a velocity response to respiratory challenges was detected. Regarding the relationship to neurovascular coupling, the removal of the whole lag structure, which can be considered as an optimized global-signal regression, resulted in a reduction of inter-individual variance while preserving the fMRI response. Examination of the T2* and S₀, or non-BOLD, components of the fMRI signal revealed that the lag structure is deoxyhemoglobin dependent, while paradoxically presenting a signal-magnitude reduction in the venous side of the cerebral vasculature. These findings provide insight into the origin of BOLD sLFOs, suggesting that they are highly intrinsic to the circulating blood.

More than BOLD: Dual-spin populations create functional contrast

PURPOSE

Functional MRI contrast has generally been associated with changes in transverse relaxivity caused by blood oxygen concentration, the so-called blood oxygen level dependent contrast. However, this interpretation of fMRI contrast has been called into question by several recent experiments at high spatial resolution. Experiments were conducted to examine contrast dependencies that cannot be explained only by differences in relaxivity in a single-spin population.

METHODS

Measurements of functional signal and contrast were obtained in human early visual cortex during a high-contrast visual stimulation over a large range of TEs and for several flip angles. Small voxels (1.5 mm) were used to restrict the measurements to cortical gray matter in early visual areas identified using retinotopic mapping procedures.

RESULTS

Measurements were consistent with models that include 2 spin populations. The dominant population has a relatively short transverse lifetime that is strongly modulated by activation. However, functional contrast is also affected by volume changes between this short-lived population and the longer-lived population.

CONCLUSION

Some of the previously observed "nonclassical" behaviors of functional contrast can be explained by these interacting dual-spin populations.

Vascular and neural basis of the BOLD signal

Neural activity in the brain is usually coupled to increases in local cerebral blood flow, leading to the increase in oxygenation that generates the BOLD fMRI signal. Recent work has begun to elucidate the vascular and neural mechanisms underlying the BOLD signal. The dilatory response is distributed throughout the vascular network. Arteries actively dilate within a second following neural activity increases, while venous distensions are passive and have a time course that last tens of seconds. Vasodilation, and thus local blood flow, is controlled by the activity of both neurons and astrocytes via multiple different pathways. The relationship between sensory-driven neural activity and the vascular dynamics in sensory areas are well-captured with a linear convolution model. However, depending on the behavioral state or brain region, the coupling between neural activity and hemodynamic signals can be weak or even inverted.

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.

Simulations of blood as a suspension predicts a depth dependent hematocrit in the circulation throughout the cerebral cortex

Recent advances in modeling oxygen supply to cortical brain tissue have begun to elucidate the functional mechanisms of neurovascular coupling. While the principal mechanisms of blood flow regulation after neuronal firing are generally known, mechanistic hemodynamic simulations cannot yet pinpoint the exact spatial and temporal coordination between the network of arteries, arterioles, capillaries and veins for the entire brain. Because of the potential significance of blood flow and oxygen supply simulations for illuminating spatiotemporal regulation inside the cortical microanatomy, there is a need to create mathematical models of the entire cerebral circulation with realistic anatomical detail. Our hypothesis is that an anatomically accurate reconstruction of the cerebrocirculatory architecture will inform about possible regulatory mechanisms of the neurovascular interface. In this article, we introduce large-scale networks of the murine cerebral circulation spanning the Circle of Willis, main cerebral arteries connected to the pial network down to the microcirculation in the capillary bed. Several multiscale models were generated from state-of-the-art neuroimaging data. Using a vascular network construction algorithm, the entire circulation of the middle cerebral artery was synthesized. Blood flow simulations indicate a consistent trend of higher hematocrit in deeper cortical layers, while surface layers with shorter vascular path lengths seem to carry comparatively lower red blood cell (RBC) concentrations. Moreover, the variability of RBC flux decreases with cortical depth. These results support the notion that plasma skimming serves a self-regulating function for maintaining uniform oxygen perfusion to neurons irrespective of their location in the blood supply hierarchy. Our computations also demonstrate the practicality of simulating blood flow for large portions of the mouse brain with existing computer resources. The efficient simulation of blood flow throughout the entire middle cerebral artery (MCA) territory is a promising milestone towards the final aim of predicting blood flow patterns for the entire brain.

The brain’s astonishing cognitive capacity depends on the coordination between neurons and the cerebral circulation, a system known as the neurovascular unit. The spatial and temporal coupling between these two networks is the object of intense research. However, the concise anatomical description of the cerebral circulation has so far been intractable. This paper introduces a methodology for the in silico creation of realistic models for the entire cerebral circulation. This innovation incorporates topological data from several neuroimaging modalities covering three lengths scales as input into a computer algorithm, which assembles anatomically accurate circulatory networks. When simulating blood flow as red blood cells suspended in plasma for experimental and synthetic cortical network models, we discovered that red blood cells tend to be more concentrated in deeper layers of the cortex compared to the surface. RBC fluxes are more homogenous in deeper layers. The phenomenon of depth dependent red blood cell supply supports the notion that the intricate architecture of the cortical microcirculation serves a self-regulating function to maintain uniform oxygen perfusion to neurons. We also demonstrate the practicality of predicting blood flow patterns for the entire brain with existing computer power.

Stimulus-dependent hemodynamic response timing across the human subcortical-cortical visual pathway identified through high spatiotemporal resolution 7T fMRI

Recent developments in fMRI acquisition techniques now enable fast sampling with whole-brain coverage, suggesting fMRI can be used to track changes in neural activity at increasingly rapid timescales. When images are acquired at fast rates, the limiting factor for fMRI temporal resolution is the speed of the hemodynamic response. Given that HRFs may vary substantially in subcortical structures, characterizing the speed of subcortical hemodynamic responses, and how the hemodynamic response shape changes with stimulus duration (i.e. the hemodynamic nonlinearity), is needed for designing and interpreting fast fMRI studies of these regions. We studied the temporal properties and nonlinearities of the hemodynamic response function (HRF) across the human subcortical visual system, imaging superior colliculus (SC), lateral geniculate nucleus of the thalamus (LGN) and primary visual cortex (V1) with high spatiotemporal resolution 7 Tesla fMRI. By presenting stimuli of varying durations, we mapped the timing and nonlinearity of hemodynamic responses in these structures at high spatiotemporal resolution. We found that the hemodynamic response is consistently faster and narrower in subcortical structures than in cortex. However, the nonlinearity in LGN is similar to that in cortex, with shorter duration stimuli eliciting larger and faster responses than would have been predicted by a linear model. Using oscillatory visual stimuli, we tested the frequency response in LGN and found that its BOLD response tracked high-frequency (0.5 Hz) oscillations. The LGN response magnitudes were comparable to V1, allowing oscillatory BOLD signals to be detected in LGN despite the small size of this structure. These results suggest that the increase in the speed and amplitude of the hemodynamic response when neural activity is brief may be the key physiological driver of fast fMRI signals, enabling detection of high-frequency oscillations with fMRI. We conclude that subcortical visual structures exhibit fast and nonlinear hemodynamic responses, and that these dynamics enable detection of fast BOLD signals even within small deep brain structures when imaging is performed at ultra-high field.

Simulating laminar neuroimaging data for a visual delayed match-to-sample task

Invasive electrophysiological and neuroanatomical studies in nonhuman mammalian experimental preparations have helped elucidate the lamina (layer) dependence of neural computations and interregional connections. Noninvasive functional neuroimaging can, in principle, resolve cortical laminae (layers), and thus provide insight into human neural computations and interregional connections. However human neuroimaging data are noisy and difficult to interpret; biologically realistic simulations can aid experimental interpretation by relating the neuroimaging data to simulated neural activity. We illustrate the potential of laminar neuroimaging by upgrading an existing large-scale, multiregion neural model that simulates a visual delayed match-to-sample (DMS) task. The new laminar-based neural unit incorporates spiny stellate, pyramidal, and inhibitory neural populations which are divided among supragranular, granular, and infragranular laminae (layers). We simulated neural activity which is translated into local field potential-like data used to simulate conventional and laminar fMRI activity. We implemented the laminar connectivity schemes proposed by Felleman and Van Essen (Cerebral Cortex, 1991) for interregional connections. The hemodynamic model that we employ is a modified version of one due to Heinzle et al. (Neuroimage, 2016) that incorporates the effects of draining veins. We show that the laminar version of the model replicates the findings of the existing model. The laminar model shows the finer structure in fMRI activity and functional connectivity. Laminar differences in the magnitude of neural activities are a prominent finding; these are also visible in the simulated fMRI. We illustrate differences between task and control conditions in the fMRI signal, and demonstrate differences in interregional laminar functional connectivity that reflect the underlying connectivity scheme. These results indicate that multi-layer computational models can aid in interpreting layer-specific fMRI, and suggest that increased use of laminar fMRI could provide unique and fundamental insights to human neuroscience.

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.

Post-stimulus fMRI and EEG responses: Evidence for a neuronal origin hypothesised to be inhibitory

Post-stimulus undershoots, negative responses following cessation of stimulation, are widely observed in functional magnetic resonance (fMRI) blood oxygenation level dependent (BOLD) data. However, the debate surrounding whether the origin of this response phase is neuronal or vascular, and whether it provides functionally relevant information, that is additional to what is contained in the primary response, means that undershoots are widely overlooked. We simultaneously recorded electroencephalography (EEG), BOLD and cerebral blood-flow (CBF) [obtained from arterial spin labelled (ASL) fMRI] fMRI responses to hemifield checkerboard stimulation to test the potential neural origin of the fMRI post-stimulus undershoot. The post-stimulus BOLD and CBF signal amplitudes in both contralateral and ipsilateral visual cortex depended on the post-stimulus power of the occipital 8-13Hz (alpha) EEG neuronal activity, such that trials with highest EEG power showed largest fMRI undershoots in contralateral visual cortex. This correlation in post-stimulus EEG-fMRI responses was not predicted by the primary response amplitude. In the contralateral visual cortex we observed a decrease in both cerebral rate of oxygen metabolism (CMRO₂) and CBF during the post-stimulus phase. In addition, the coupling ratio (n) between CMRO₂ and CBF was significantly lower during the positive contralateral primary response phase compared with the post-stimulus phase and we propose that this reflects an altered balance of excitatory and inhibitory neuronal activity. Together our data provide strong evidence that the post-stimulus phase of the BOLD response has a neural origin which reflects, at least partially, an uncoupling of the neuronal responses driving the primary and post-stimulus responses, explaining the uncoupling of the signals measured in the two response phases. We suggest our results are consistent with inhibitory processes driving the post-stimulus EEG and fMRI responses. We therefore propose that new methods are required to model the post-stimulus and primary responses independently, enabling separate investigation of response phases in cognitive function and neurological disease.

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.

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.

Using an achiasmic human visual system to quantify the relationship between the fMRI BOLD signal and neural response

Achiasma in humans causes gross mis-wiring of the retinal-fugal projection, resulting in overlapped cortical representations of left and right visual hemifields. We show that in areas V1-V3 this overlap is due to two co-located but non-interacting populations of neurons, each with a receptive field serving only one hemifield. Importantly, the two populations share the same local vascular control, resulting in a unique organization useful for quantifying the relationship between neural and fMRI BOLD responses without direct measurement of neural activity. Specifically, we can non-invasively double local neural responses by stimulating both neuronal populations with identical stimuli presented symmetrically across the vertical meridian to both visual hemifields, versus one population by stimulating in one hemifield. Measurements from a series of such doubling experiments show that the amplitude of BOLD response is proportional to approximately 0.5 power of the underlying neural response. Reanalyzing published data shows that this inferred relationship is general.