Early and late stimulus-evoked cortical hemodynamic responses provide insight into the neurogenic nature of neurovascular coupling

Cracking and Packing Information about the Features of Expected Rewards in the Orbitofrontal Cortex

The orbitofrontal cortex (OFC) is crucial for tracking various aspects of expected outcomes, thereby helping to guide choices and support learning. Our previous study showed that the effects of reward timing and size on the activity of single units in OFC were dissociable when these attributes were manipulated independently ( Roesch et al., 2006). However, in real-life decision-making scenarios, outcome features often change simultaneously, so here we investigated how OFC neurons in male rats integrate information about the timing and identity (flavor) of reward and respond to changes in these features, according to whether they were changed simultaneously or separately. We found that a substantial number of OFC neurons fired differentially to immediate versus delayed reward and to the different reward flavors. However, contrary to the previous study, selectivity for timing was strongly correlated with selectivity for identity. Taken together with the previous research, these results suggest that when reward features are correlated, OFC tends to "pack" them into unitary constructs, whereas when they are independent, OFC tends to "crack" them into separate constructs. Furthermore, we found that when both reward timing and flavor were changed, reward-responsive OFC neurons showed unique activity patterns preceding and during the omission of an expected reward. Interestingly, this OFC activity is similar and slightly preceded the ventral tegmental area dopamine (VTA DA) activity observed in a previous study ( Takahashi et al., 2023), consistent with the role of OFC in providing predictive information to VTA DA neurons.

Neurovascular coupling and CO2 interrogate distinct vascular regulations

Neurovascular coupling (NVC), which mediates rapid increases in cerebral blood flow in response to neuronal activation, is commonly used to map brain activation or dysfunction. Here we tested the reemerging hypothesis that CO2 generated by neuronal metabolism contributes to NVC. We combined functional ultrasound and two-photon imaging in the mouse barrel cortex to specifically examine the onsets of local changes in vessel diameter, blood flow dynamics, vascular/perivascular/intracellular pH, and intracellular calcium signals along the vascular arbor in response to a short and strong CO2 challenge (10 s, 20%) and whisker stimulation. We report that the brief hypercapnia reversibly acidifies all cells of the arteriole wall and the periarteriolar space 3-4 s prior to the arteriole dilation. During this prolonged lag period, NVC triggered by whisker stimulation is not affected by the acidification of the entire neurovascular unit. As it also persists under condition of continuous inflow of CO2, we conclude that CO2 is not involved in NVC.

The effect of CO2 on the age dependence of neurovascular coupling

Prior studies have identified variable effects of aging on neurovascular coupling (NVC). Carbon dioxide (CO2) affects both cerebral blood velocity (CBv) and NVC, but the effects of age on NVC under different CO2 conditions are unknown. Therefore, we investigated the effects of aging on NVC in different CO2 states during cognitive paradigms. Seventy-eight participants (18-78 yr), with well-controlled comorbidities, underwent continuous recordings of CBv by bilateral insonation of middle (MCA) and posterior (PCA) cerebral arteries (transcranial Doppler), blood pressure, end-tidal CO2, and heart rate during poikilocapnia, hypercapnia (5% CO2 inhalation), and hypocapnia (paced hyperventilation). Neuroactivation via visuospatial (VS) and attention tasks (AT) was used to stimulate NVC. Peak percentage and absolute change in MCAv/PCAv, were compared between CO2 conditions and age groups (≤30, 31-60, and >60 yr). For the VS task, in poikilocapnia, younger adults had a lower NVC response compared with older adults [mean difference (MD): -7.92% (standard deviation (SD): 2.37), P = 0.004], but comparable between younger and middle-aged groups. In hypercapnia, both younger [MD: -4.75% (SD: 1.56), P = 0.009] and middle [MD: -4.58% (SD: 1.69), P = 0.023] age groups had lower NVC responses compared with older adults. Finally, in hypocapnia, both older [MD: 5.92% (SD: 2.21), P = 0.025] and middle [MD: 5.44% (SD: 2.27), P = 0.049] age groups had greater NVC responses, compared with younger adults. In conclusion, the magnitude of NVC response suppression from baseline during hyper- and hypocapnia, did not differ significantly between age groups. However, the middle age group demonstrated a different NVC response while under hypercapnic conditions, compared with hypocapnia.NEW & NOTEWORTHY This study describes the effects of age on neurovascular coupling under altered CO2 conditions. We demonstrated that both hypercapnia and hypocapnia suppress neurovascular coupling (NVC) responses. Furthermore, that middle age exhibits an NVC response comparable with younger adults under hypercapnia, and older adults under hypocapnia.

Laminar CBV and BOLD response-characteristics over time and space in the human primary somatosensory cortex at 7T

Uncovering the cortical representation of the body has been at the core of human brain mapping for decades, with special attention given to the digits. In the last decade, advances in functional magnetic resonance imaging (fMRI) technologies have opened the possibility of noninvasively unraveling the 3rd dimension of digit representations in humans along cortical layers. In laminar fMRI it is common to combine the use of the highly sensitive blood oxygen level dependent (BOLD) contrast with cerebral blood volume sensitive measurements, like vascular space occupancy (VASO), that are more specific to the underlying neuronal populations. However, the spatial and temporal VASO response characteristics across cortical depth to passive stimulation of the digits are still unknown. Therefore, we characterized haemodynamic responses to vibrotactile stimulation of individual digit-tips across cortical depth at 0.75 mm in-plane spatial resolution using BOLD and VASO fMRI at 7T. We could identify digit-specific regions of interest (ROIs) in putative Brodmann area 3b, following the known anatomical organization. In the ROIs, the BOLD response increased towards the cortical surface due to the draining vein effect, while the VASO response was more shifted towards middle cortical layers, likely reflecting bottom-up input from the thalamus, as expected. Interestingly, we also found slightly negative BOLD and VASO responses for non-preferred digits in the ROIs, potentially indicating neuronal surround inhibition. Finally, we explored the temporal signal dynamics for BOLD and VASO as a function of distance from activation peaks resulting from stimulation of contralateral digits. With this analysis, we showed a triphasic response consisting of an initial peak and a subsequent negative deflection during stimulation, followed by a positive post-stimulus response in BOLD and to some extent in VASO. While similar responses were reported with invasive methods in animal models, here we demonstrate a potential neuronal excitation-inhibition mechanism in a center-surround architecture across layers in the human somatosensory cortex. Given that, unlike in animals, human experiments do not rely on anesthesia and can readily implement extensive behavioral testing, obtaining this effect in humans is an important step towards further uncovering the functional significance of the different aspects of the triphasic response.

Whisker-evoked neurovascular coupling is preserved during hypoglycemia in mouse cortical arterioles and capillaries

Hypoglycemia is a serious complication of insulin treatment of diabetes that can lead to coma and death. Neurovascular coupling, which mediates increased local blood flow in response to neuronal activity, increases glucose availability to active neurons. This mechanism could be essential for neuronal health during hypoglycemia, when total glucose supplies are low. Previous studies suggest, however, that neurovascular coupling (a transient blood flow increase in response to an increase in neuronal activity) may be reduced during hypoglycemia. Such a reduction in blood flow increase would exacerbate the effects of hypoglycemia, depriving active neurons of glucose. We have reexamined the effects of hypoglycemia on neurovascular coupling by simultaneously monitoring neuronal and vascular responses to whisker stimulation in the awake mouse somatosensory cortex. We find that neurovascular coupling at both penetrating arterioles and at 2nd order capillaries did not change significantly during insulin-induced hypoglycemia compared to euglycemia. In addition, we show that the basal diameter of both arterioles and capillaries increases during hypoglycemia (10.3 and 9.7% increases, respectively). Our results demonstrate that both neurovascular coupling and basal increases in vessel diameter are active mechanisms which help to maintain an adequate supply of glucose to the brain during hypoglycemia.

Characterisation of laminar and vascular spatiotemporal dynamics of CBV and BOLD signals using VASO and ME-GRE at 7T in humans

Interpretation of cortical laminar functional magnetic resonance imaging (fMRI) activity requires detailed knowledge of the spatiotemporal haemodynamic response across vascular compartments due to the well-known vascular biases (e.g. the draining veins). Further complications arise from the spatiotemporal hemodynamic response that differs depending on the duration of stimulation. This information is crucial for future studies using depth-dependent cerebral blood volume (CBV) measurements, which promise higher specificity for the cortical microvasculature than the blood oxygenation level dependent (BOLD) contrast. To date, direct information about CBV dynamics with respect to stimulus duration, cortical depth and vasculature is missing in humans. Therefore, we characterized the cortical depth-dependent CBV-haemodynamic responses across a wide set of stimulus durations with 0.9 mm isotropic spatial and 0.785 seconds effective temporal resolution in humans using slice-selective slab-inversion vascular space occupancy (SS-SI VASO). Additionally, we investigated signal contributions from macrovascular compartments using fine-scale vascular information from multi-echo gradient-echo (ME-GRE) data at 0.35 mm isotropic resolution. In total, this resulted in >7.5h of scanning per participant (n=5). We have three major findings: (I) While we could demonstrate that 1 second stimulation is viable using VASO, more than 12 seconds stimulation provides better CBV responses in terms of specificity to microvasculature, but durations beyond 24 seconds of stimulation may be wasteful for certain applications. (II) We observe that CBV responses show dilation patterns across the cortex. (III) While we found increasingly strong BOLD signal responses in vessel-dominated voxels with longer stimulation durations, we found increasingly strong CBV signal responses in vessel-dominated voxels only until 4 second stimulation durations. After 4 seconds, only the signal from non-vessel dominated voxels kept increasing. This might explain why CBV responses are more specific to the underlying neuronal activity for long stimulus durations.

Layer-fMRI VASO with short stimuli and event-related designs at 7 T

Layers and columns are the dominant processing units in the human (neo)cortex at the mesoscopic scale. While the blood oxygenation dependent (BOLD) signal has a high detection sensitivity, it is biased towards unwanted signals from large draining veins at the cortical surface. The additional fMRI contrast of vascular space occupancy (VASO) has the potential to augment the neuroscientific interpretability of layer-fMRI results by means of capturing complementary information of locally specific changes in cerebral blood volume (CBV). Specifically, VASO is not subject to unwanted sensitivity amplifications of large draining veins. Because of constrained sampling efficiency, it has been mainly applied in combination with efficient block task designs and long trial durations. However, to study cognitive processes in neuroscientific contexts, or probe vascular reactivity, short stimulation periods are often necessary. Here, we developed a VASO acquisition procedure with a short acquisition period and sub-millimeter resolution. During visual event-related stimulation, we show reliable responses in visual cortices within a reasonable number of trials (∼20). Furthermore, the short TR and high spatial specificity of our VASO implementation enabled us to show differences in laminar reactivity and onset times. Finally, we explore the generalizability to a different stimulus modality (somatosensation). With this, we showed that CBV-sensitive VASO provides the means to capture layer-specific haemodynamic responses with high spatio-temporal resolution and is able to be used with event-related paradigms.

Resting-state networks representation of the global phenomena

Resting-state functional magnetic resonance imaging (rsfMRI) has been widely applied to investigate spontaneous neural activity, often based on its macroscopic organization that is termed resting-state networks (RSNs). Although the neurophysiological mechanisms underlying the RSN organization remain largely unknown, accumulating evidence points to a substantial contribution from the global signals to their structured synchronization. This study further explored the phenomenon by taking advantage of the inter- and intra-subject variations of the time delay and correlation coefficient of the signal timeseries in each region using the global mean signal as the reference signal. Consistent with the hypothesis based on the empirical and theoretical findings, the time lag and correlation, which have consistently been proven to represent local hemodynamic status, were shown to organize networks equivalent to RSNs. The results not only provide further evidence that the local hemodynamic status could be the direct source of the RSNs' spatial patterns but also explain how the regional variations in the hemodynamics, combined with the changes in the global events' power spectrum, lead to the observations. While the findings pose challenges to interpretations of rsfMRI studies, they further support the view that rsfMRI can offer detailed information related to global neurophysiological phenomena as well as local hemodynamics that would have great potential as biomarkers.

Concurrent intrinsic optical imaging and fMRI at ultra-high field using magnetic field proof optical components

Intrinsic optical imaging (IOI) is a well established technique to quantify activation-related hemodynamical changes at the surface of the brain, which can be used to investigate the underlying processes of BOLD signal formation. To directly and quantitatively relate IOI and fMRI, simultaneous measurements with the two modalities are necessary. Here, a novel technical solution for a completely in-bore setup is presented, which uses only magnetic field proof components and thus allows concurrent recordings with a quality similar to that obtained in separate experiments. Measurements of the somatosensory cortex of rats with electrical forepaw stimulation were used to verify this approach. The high spatial and temporal resolution of the fMRI data, which is possible due to the high magnetic field of 14.1 T, the use of a point-spread function-based distortion correction and optimized additional anatomical images, allowed accurate colocalization of the images of the two modalities. Accordingly, detailed investigations of the temporal and spatial relationships between the hemodynamic parameters and the fMRI signal, which demonstrate the linear dependence of the BOLD effect on changes in the concentrations of oxygenated and deoxygenated hemoglobin, are possible. Comparisons between the signals emerging from arterial, venous and parenchymal areas are possible and show clearly distinct characteristics. The presented setup allows combining MRI measurements and optical recordings without serious losses in the data quality of either modality. While the proposed combination of fMRI and IOI can help to gain valuable insight into the generation of the BOLD effect, the setup can be easily modified to include different types of optical or MRI measurements.

Validating layer-specific VASO across species

Cerebral blood volume (CBV) has been shown to be a robust and important physiological parameter for quantitative interpretation of functional (f)MRI, capable of delivering highly localized mapping of neural activity. Indeed, with recent advances in ultra-high-field (≥7T) MRI hardware and associated sequence libraries, it has become possible to capture non-invasive CBV weighted fMRI signals across cortical layers. One of the most widely used approaches to achieve this (in humans) is through vascular-space-occupancy (VASO) fMRI. Unfortunately, the exact contrast mechanisms of layer-dependent VASO fMRI have not been validated for human fMRI and thus interpretation of such data is confounded. Here we validate the signal source of layer-dependent SS-SI VASO fMRI using multi-modal imaging in a rat model in response to neuronal activation (somatosensory cortex) and respiratory challenge (hypercapnia). In particular VASO derived CBV measures are directly compared to concurrent measures of total haemoglobin changes from high resolution intrinsic optical imaging spectroscopy (OIS). Quantified cortical layer profiling is demonstrated to be in agreement between VASO and contrast enhanced fMRI (using monocrystalline iron oxide nanoparticles, MION). Responses show high spatial localisation to layers of cortical processing independent of confounding large draining veins which can hamper BOLD fMRI studies, (depending on slice positioning). Thus, a cross species comparison is enabled using VASO as a common measure. We find increased VASO based CBV reactivity (3.1 ± 1.2 fold increase) in humans compared to rats. Together, our findings confirm that the VASO contrast is indeed a reliable estimate of layer-specific CBV changes. This validation study increases the neuronal interpretability of human layer-dependent VASO fMRI as an appropriate method in neuroscience application studies, in which the presence of large draining intracortical and pial veins limits neuroscientific inference with BOLD fMRI.

Origin of the Time Lag Phenomenon and the Global Signal in Resting-State fMRI

The global mean signal of resting-state fMRI (rs-fMRI) shows a characteristic spatiotemporal pattern that is closely related to the pattern of vascular perfusion. Although being increasingly adopted in the mapping of the flow of neural activity, the mechanism that gives rise to the BOLD signal time lag remains controversial. In the present study, we compared the time lag of the global mean signal with those of the local network components obtained by applying temporal independent component analysis to the resting-state fMRI data, as well as by using simultaneous wide-field visual stimulation, and demonstrated that the time lag patterns are highly similar across all types of data. These results suggest that the time lag of the rs-fMRI signal reflects the local variance of the hemodynamic responses rather than the arrival or transit time of the stimulus, whether the trigger is neuronal or non-neuronal in origin as long as it is mediated by local hemodynamic responses. Examinations of the internal carotid artery signal further confirmed that the arterial signal is tightly inversely coupled with the global mean signal in accordance with previous studies, presumably reflecting the blood flow or blood pressure changes that are occurring almost simultaneously in the internal carotid artery and the cerebral pial/capillary arteries, within the low-frequency component in human rs-fMRI.

Neurovascular coupling preserved in a chronic mouse model of Alzheimer's disease: Methodology is critical

Impaired neurovascular coupling has been suggested as an early pathogenic factor in Alzheimer's disease (AD), which could serve as an early biomarker of cerebral pathology. We have established an anaesthetic regime to allow repeated measurements of neurovascular function over three months in the J20 mouse model of AD (J20-AD) and wild-type (WT) controls. Animals were 9-12 months old at the start of the experiment. Mice were chronically prepared with a cranial window through which 2-Dimensional optical imaging spectroscopy (2D-OIS) was used to generate functional maps of the cerebral blood volume and saturation changes evoked by whisker stimulation and vascular reactivity challenges. Unexpectedly, the hemodynamic responses were largely preserved in the J20-AD group. This result failed to confirm previous investigations using the J20-AD model. However, a final acute electrophysiology and 2D-OIS experiment was performed to measure both neural and hemodynamic responses concurrently. In this experiment, previously reported deficits in neurovascular coupling in the J20-AD model were observed. This suggests that J20-AD mice may be more susceptible to the physiologically stressing conditions of an acute experimental procedure compared to WT animals. These results therefore highlight the importance of experimental procedure when determining the characteristics of animal models of human disease.

BOLD fMRI and hemodynamic responses to somatosensory stimulation in anesthetized mice: spontaneous breathing vs. mechanical ventilation

Mouse functional MRI (fMRI) has been of great interest due to the abundance of transgenic models. Due to a mouse's small size, spontaneous breathing has often been used. Because the vascular physiology affecting fMRI might not be controlled normally, its effects on functional responses were investigated with optical intrinsic signal (OIS) imaging and 9.4 T BOLD fMRI. Three conditions were tested in C57BL/6 mice: spontaneous breathing under ketamine and xylazine anesthesia (KX), mechanical ventilation under KX, and mechanical ventilation under isoflurane. Spontaneous breathing under KX induced an average pCO₂ of 83 mmHg, whereas a mechanical ventilation condition achieved a pCO₂ of 37-41 mmHg within a physiological range. The baseline diameter of arterial and venous vessels was only 7%-9% larger with spontaneous breathing than with mechanical ventilation under KX, but it was much smaller than that in normocapnic isoflurane-anesthetized mice. Three major functional studies were performed. First, CBV-weighted OIS and arterial dilations to 4-second forepaw stimulation were rapid and larger at normocapnia than hypercapnia under KX, but very small under isoflurane. Second, CBV-weighted OIS and arterial dilations by vasodilator acetazolamide were measured for investigating vascular reactivity and were larger in the normocapnic condition than in the hypercapnic condition under KX. Third, evoked OIS and BOLD fMRI responses in the contralateral mouse somatosensory cortex to 20-second forepaw stimulation were faster and larger in the mechanical ventilation than spontaneous breathing. BOLD fMRI peaked at the end of the 20-second stimulation under hypercapnic spontaneous breathing, and at ~9 seconds under mechanical ventilation. The peak amplitude of BOLD fMRI was 2.2% at hypercapnia and ~3.4% at normocapnia. Overall, spontaneous breathing induces sluggish reduced hemodynamic and fMRI responses, but it is still viable for KX anesthesia due to its simplicity, noninvasiveness, and well-localized BOLD activity in the somatosensory cortex.

Interactions between stimuli-evoked cortical activity and spontaneous low frequency oscillations measured with neuronal calcium

Spontaneous brain activity has been widely used to map brain connectivity. The interactions between task-evoked brain responses and the spontaneous cortical oscillations, especially within the low frequency range of ~0.1 ​Hz, are not fully understood. Trial-to-trial variabilities in brain's response to sensory stimuli and the ability for brain to detect under noisy conditions suggest an appreciable impact of the brain state. Using a multimodality imaging platform, we simultaneously imaged neuronal Ca²⁺ and cerebral hemodynamics at baseline and in response to single-pulse forepaw stimuli in rat's somatosensory cortex. The high sensitivity of this system enables detection of responses to very weak and strong stimuli and real time determination of low frequency oscillations without averaging. Results show that the ongoing neuronal oscillations inversely modulate Ca²⁺ transients evoked by sensory stimuli. High intensity stimuli reset the spontaneous neuronal oscillations to an unpreferable excitability following the stimulus. Cerebral hemodynamic responses also inversely interact with the spontaneous hemodynamic oscillations, correlating with the neuronal Ca²⁺ transient changes. The results reveal competing interactions between spontaneous oscillations and stimulation-evoked brain activities in somatosensory cortex and the resultant hemodynamics.

Directional sensitivity of dynamic cerebral autoregulation in squat-stand maneuvers

Dynamic cerebral autoregulation (CA), the transient response of cerebral blood flow (CBF) to rapid changes in arterial blood pressure (BP), is usually modeled as a linear mechanism. We tested the hypothesis that dynamic CA can display nonlinear behavior resulting from differential efficiency dependent on the direction of BP changes. Cerebral blood velocity (CBV) (transcranial Doppler), heart rate (HR) (three-lead ECG), continuous BP (Finometer), and end-tidal CO2 (capnograph) were measured in 10 healthy young subjects during 15 squat-stand maneuvers (SSM) with a frequency of 0.05 Hz. The protocol was repeated with a median (interquartile range) of 44 (35–64) days apart. Dynamic CA was assessed with the autoregulation index (ARI) obtained from CBV step responses estimated with an autoregressive moving-average model. Mean BP, HR, and CBV were different (all P < 0.001) between squat and stand, regardless of visits. ARI showed a strong interaction ( P < 0.001) of SSM with the progression of transients; in general, the mean ARI was higher for the squat phase compared with standing. The changes in ARI were partially explained by concomitant changes in CBV ( P = 0.023) and pulse pressure ( P < 0.001), but there was no evidence that ARI differed between visits ( P = 0.277). These results demonstrate that dynamic CA is dependent on the direction of BP change, but further work is needed to confirm if this finding can be generalized to other physiological conditions and also to assess its dependency on age, sex and pathology.

The use of intensity-based Doppler variance method for single vessel response to functional neurovascular activation

This study presents 1 use of optical coherence tomography (OCT) angiography technique to examine neurovascular coupling effect. Repeated B-scans OCT recording is performed on the rat somatosensory cortex with cranial window preparation while its contralateral forepaw is electrically stimulated to activate the neurons in rest. We use an intensity-based Doppler variance (IBDV) algorithm mapped cerebral blood vessels in the cortex, and the temporal alteration in blood perfusion during neurovascular activation is analyzed using the proposed IBDV quantitative parameters. By using principal component analysis-based Fuzzy C Means clustering method, the stimulus-evoked vasomotion patterns were classified into 3 categories. We found that the response time of small vessels (resting diameter 14.9 ±6.6 μm), middle vessels (resting diameter 21.1 ±7.9 μm) and large vessels (resting diameter 50.7 ±6.5 μm) to achieve 5% change of vascular dilation after stimulation was 1.5, 2 and 5.5 seconds, respectively. Approximately 5% peak change of relative blood flow (RBF) in both small and middle vessels was observed. The large vessels react slowly and their responses nearly 4 seconds delayed, but no significant change in RBF of the large vessels was seen.

Techniques for blood volume fMRI with VASO: From low-resolution mapping towards sub-millimeter layer-dependent applications

Quantitative cerebral blood volume (CBV) fMRI has the potential to overcome several specific limitations of BOLD fMRI. It provides direct physiological interpretability and promises superior localization specificity in applications of sub-millimeter resolution fMRI applications at ultra-high magnetic fields (7T and higher). Non-invasive CBV fMRI using VASO (vascular space occupancy), however, is inherently limited with respect to its data acquisition efficiency, restricting its imaging coverage and achievable spatial and temporal resolution. This limitation may be reduced with recent advanced acceleration and reconstruction strategies that allow two-dimensional acceleration, such as in simultaneous multi-slice (SMS) 2D-EPI or 3D-EPI in combination with CAIPIRINHA field-of-view shifting. In this study, we sought to determine the functional sensitivity and specificity of these readout strategies with VASO over a broad range of spatial resolutions; spanning from low spatial resolution (3mm) whole-cortex to sub-millimeter (0.75mm) slab-of-cortex (for cortical layer-dependent applications). In the thermal-noise-dominated regime of sub-millimeter resolutions, 3D-EPI-VASO provides higher temporal stability and sensitivity to detect changes in CBV compared to 2D-EPI-VASO. In this regime, 3D-EPI-VASO unveils task activation located in the cortical laminae with little contamination from surface veins, in contrast to the cortical surface weighting of GE-BOLD fMRI. In the physiological-noise-dominated regime of lower resolutions, however, 2D-SMS-VASO shows superior performance compared to 3D-EPI-VASO. Due to its superior sensitivity at a layer-dependent level, 3D-EPI VASO promises to play an important role in future neuroscientific applications of layer-dependent fMRI.

Wide-field optical mapping of neural activity and brain haemodynamics: considerations and novel approaches

Although modern techniques such as two-photon microscopy can now provide cellular-level three-dimensional imaging of the intact living brain, the speed and fields of view of these techniques remain limited. Conversely, two-dimensional wide-field optical mapping (WFOM), a simpler technique that uses a camera to observe large areas of the exposed cortex under visible light, can detect changes in both neural activity and haemodynamics at very high speeds. Although WFOM may not provide single-neuron or capillary-level resolution, it is an attractive and accessible approach to imaging large areas of the brain in awake, behaving mammals at speeds fast enough to observe widespread neural firing events, as well as their dynamic coupling to haemodynamics. Although such wide-field optical imaging techniques have a long history, the advent of genetically encoded fluorophores that can report neural activity with high sensitivity, as well as modern technologies such as light emitting diodes and sensitive and high-speed digital cameras have driven renewed interest in WFOM. To facilitate the wider adoption and standardization of WFOM approaches for neuroscience and neurovascular coupling research, we provide here an overview of the basic principles of WFOM, considerations for implementation of wide-field fluorescence imaging of neural activity, spectroscopic analysis and interpretation of results.

This article is part of the themed issue ‘Interpreting BOLD: a dialogue between cognitive and cellular neuroscience’.

Cerebral haemodynamic response to somatosensory stimulation in neonatal lambs

KEY POINTS

Cerebral haemodynamic response to neural stimulation has been extensively studied in adults, but little is known about cerebral haemodynamic response in the fetal and neonatal brain. The present study describes the cerebral haemodynamic response measured by near infrared spectroscopy to somatosensory stimulation in newborn lambs, in comparison to recent findings in fetal sheep. The cerebral haemodynamic responses in the newborn lamb brain can involve an increase in oxyhaemoglobin (oxyHb), or a decrease of oxyHb suggestive of reduced perfusion and oxygenation. Positive correlations between changes in oxyHb and mean arterial blood pressure were found in newborn but not fetal sheep, which suggests the result is unlikely to be due to immature autoregulation alone. In contrast to adult studies, hypercapnia increased the changes in cerebral blood flow and oxyHb in most of the lambs in response to somatosensory stimulation.

ABSTRACT

The neurovascular coupling response has been defined for the adult brain, but in the neonate non-invasive measurement of local cerebral perfusion using near infrared spectroscopy or blood oxygen level-dependent functional magnetic resonance imaging have yielded variable and inconsistent results, including negative responses suggesting decreased perfusion and localized tissue tissue hypoxia. Also, the impact of permissive hypercapnia (P aC O2 > 50 mmHg) in the management of neonates on cerebrovascular responses to somatosensory input is unknown. Using near infrared spectroscopy to measure changes in cerebral oxy- and deoxyhaemoglobin (ΔoxyHb, ΔdeoxyHb) in eight anaesthetized newborn lambs, we studied the cerebral haemodynamic functional response to left median nerve stimulation using stimulus trains of 1.8, 4.8 and 7.8 s. Stimulation always produced a somatosensory evoked response, and superficial cortical perfusion measured by laser Doppler flowmetry predominantly increased following median nerve stimulation. However, with 1.8 s stimulation, oxyHb responses in the contralateral hemisphere were either positive (i.e. increased oxyHb), negative, or absent; and with 4.8 and 7.8 s stimulations, both positive and negative responses were observed. Hypercapnia increased baseline oxyHb and total Hb consistent with cerebral vasodilatation, and six of seven lambs tested showed increased Δtotal Hb responses after the 7.8 s stimulation, among which four lambs also showed increased ΔoxyHb responses. In two of three lambs, the negative ΔoxyHb response became a positive pattern during hypercapnia. These results show that instead of functional hyperaemia, somatosensory stimulation can evoke negative (decreased oxyHb, total Hb) functional responses in the neonatal brain suggestive of decreased local perfusion and vasoconstriction, and that hypercapnia produces both baseline hyperperfusion and increased functional hyperaemia.

Vasodilation Mechanism of Cerebral Microvessels Induced by Neural Activation under High Baseline Cerebral Blood Flow Level Results from Hypercapnia in Awake Mice

OBJECTIVE

We investigated the effects of the baseline CBF level at resting state on neurovascular coupling.

METHODS

Diameters of arterioles, capillaries, and venulas in awake mouse brain were measured by a two-photon microscope. Vasodilation in each of the cerebral vessels was caused by three experimental conditions: (1) sensory stimulation, (2) 5% CO2 inhalation (hypercapnia), (3) simultaneous exposure to sensory stimulation and 5% CO2 inhalation. CBF and CBV were also measured by a microscope and a CCD camera.

RESULTS

Increases in CBF and CBV were observed under all experimental conditions. After the increases in CBF and CBV due to hypercapnia, additional increases in CBF and CBV occurred during sensory stimulation. Diameter changes in arterioles were significantly larger than those in capillaries and venulas under both sensory stimulation and 5% CO2 inhalation. Additional vasodilation from sensory stimulation was observed under hypercapnia. The diameter change in each vessel type during sensory stimulation was maintained under simultaneous exposure to sensory stimulation and hypercapnia.

CONCLUSIONS

The diameter change of cerebral vessels during neural activation is reproducible regardless of whether baseline CBF has increased or not. Our finding directly demonstrates the concept of uncoupling between energy consumption and energy supply during cortical activation.

The neurogenesis of P1 and N1: A concurrent EEG/LFP study

It is generally recognised that event related potentials (ERPs) of electroencephalogram (EEG) primarily reflect summed post-synaptic activity of the local pyramidal neural population(s). However, it is still not understood how the positive and negative deflections (e.g. P1, N1 etc) observed in ERP recordings are related to the underlying excitatory and inhibitory post-synaptic activity. We investigated the neurogenesis of P1 and N1 in ERPs by pharmacologically manipulating inhibitory post-synaptic activity in the somatosensory cortex of rodent, and concurrently recording EEG and local field potentials (LFPs). We found that the P1 wave in the ERP and LFP of the supragranular layers is determined solely by the excitatory post-synaptic activity of the local pyramidal neural population, as is the initial segment of the N1 wave across cortical depth. The later part of the N1 wave was modulated by inhibitory post-synaptic activity, with its peak and the pulse width increasing as inhibition was reduced. These findings suggest that the temporal delay of inhibition with respect to excitation observed in intracellular recordings is also reflected in extracellular field potentials (FPs), resulting in a temporal window during which only excitatory post-synaptic activity and leak channel activity are recorded in the ERP and evoked LFP time series. Based on these findings, we provide clarification on the interpretation of P1 and N1 in terms of the excitatory and inhibitory post-synaptic activities of the local pyramidal neural population(s).

Hemodynamic changes in cortical sensorimotor systems following hand and orofacial motor tasks and pulsed pneumotactile stimulation

We performed a functional near-infrared spectroscopy (fNIRS) study of the evoked hemodynamic responses seen in hand and face sensorimotor cortical representations during (1) active motor tasks and (2) pulsed pneumotactile stimulation. Contralateral fNIRS measurements were performed on 22 healthy adult participants using a block paradigm that consisted of repetitive right hand and right oral angle somatosensory stimulation using a pulsed pneumotactile array stimulator, and repetitive right-hand grip compression and bilabial compressions on strain gages. Results revealed significant oxyhemoglobin (HbO) modulation across stimulus conditions in corresponding somatotopic cortical regions. Of the 22 participants, 86% exhibited a decreased HbO response during at least one of the stimulus conditions, which may be indicative of cortical steal, or hypo-oxygenation occurring in channels adjacent to the primary areas of activation. Across all conditions, 56% of participants' HbO responses were positive and 44% were negative. Hemodynamic responses most likely differed across hand and face motor and somatosensory cortical regions due to differences in regional arterial/venous anatomy, cortical vascular beds, extent and orientation of somatotopy, task dynamics, and mechanoreceptor typing in hand and face. The combination of optical imaging and task conditions allowed for non-invasive examination of hemodynamic changes in somatosensory and motor cortices using natural, pneumatic stimulation of glabrous hand and hairy skin of the lower face and functionally relevant and measurable motor tasks involving the same structures.

A novel method for classifying cortical state to identify the accompanying changes in cerebral hemodynamics

•

We classified brain state using a vector-based categorisation of neural frequencies.

•

Changes in cerebral blood volume (CBV) were observed when brain state altered.

•

During these state alterations, changes in blood oxygenation were also found.

•

State dependent haemodynamic changes could affect blood based brain imaging.

Background

Many brain imaging techniques interpret the haemodynamic response as an indirect indicator of underlying neural activity. However, a challenge when interpreting this blood based signal is how changes in brain state may affect both baseline and stimulus evoked haemodynamics.

New method

We developed an Automatic Brain State Classifier (ABSC), validated on data from anaesthetised rodents. It uses vectorised information obtained from the windowed spectral frequency power of the Local Field Potential. Current state is then classified by comparing this vectorised information against that calculated from state specific training datasets.

Results

The ABSC identified two user defined brain states (synchronised and desynchronised), with high accuracy (∼90%). Baseline haemodynamics were found to be significantly different in the two identified states. During state defined periods of elevated baseline haemodynamics we found significant decreases in evoked haemodynamic responses to somatosensory stimuli.

Comparison to existing methods

State classification – The ABSC (∼90%) demonstrated greater accuracy than clustering (∼66%) or ‘power threshold’ (∼64%) methods of comparison.

Haemodynamic averaging – Our novel approach of selectively averaging stimulus evoked haemodynamic trials by brain state yields higher quality data than creating a single average from all trials.

Conclusions

The ABSC can account for some of the commonly observed trial-to-trial variability in haemodynamic responses which arises from changes in cortical state. This variability might otherwise be incorrectly attributed to alternative interpretations. A greater understanding of the effects of cortical state on haemodynamic changes could be used to inform techniques such as general linear modelling (GLM), commonly used in fMRI.

Vascular autorescaling of fMRI (VasA fMRI) improves sensitivity of population studies: A pilot study

The blood oxygenation level-dependent (BOLD) signal is widely used for functional magnetic resonance imaging (fMRI) of brain function in health and disease. The statistical power of fMRI group studies is significantly hampered by high inter-subject variance due to differences in baseline vascular physiology. Several methods have been proposed to account for physiological vascularization differences between subjects and hence improve the sensitivity in group studies. However, these methods require the acquisition of additional reference scans (such as a full resting-state fMRI session or ASL-based calibrated BOLD). We present a vascular autorescaling (VasA) method, which does not require any additional reference scans. VasA is based on the observation that slow oscillations (<0.1Hz) in arterial blood CO2 levels occur naturally due to changes in respiration patterns. These oscillations yield fMRI signal changes whose amplitudes reflect the blood oxygenation levels and underlying local vascularization and vascular responsivity. VasA estimates proxies of the amplitude of these CO2-driven oscillations directly from the residuals of task-related fMRI data without the need for reference scans. The estimates are used to scale the amplitude of task-related fMRI responses, to account for vascular differences. The VasA maps compared well to cerebrovascular reactivity (CVR) maps and cerebral blood volume maps based on vascular space occupancy (VASO) measurements in four volunteers, speaking to the physiological vascular basis of VasA. VasA was validated in a wide variety of tasks in 138 volunteers. VasA increased t-scores by up to 30% in specific brain areas such as the visual cortex. The number of activated voxels was increased by up to 200% in brain areas such as the orbital frontal cortex while still controlling the nominal false-positive rate. VasA fMRI outperformed previously proposed rescaling approaches based on resting-state fMRI data and can be readily applied to any task-related fMRI data set, even retrospectively.

Long-latency reductions in gamma power predict hemodynamic changes that underlie the negative BOLD signal

Studies that use prolonged periods of sensory stimulation report associations between regional reductions in neural activity and negative blood oxygenation level-dependent (BOLD) signaling. However, the neural generators of the negative BOLD response remain to be characterized. Here, we use single-impulse electrical stimulation of the whisker pad in the anesthetized rat to identify components of the neural response that are related to "negative" hemodynamic changes in the brain. Laminar multiunit activity and local field potential recordings of neural activity were performed concurrently with two-dimensional optical imaging spectroscopy measuring hemodynamic changes. Repeated measurements over multiple stimulation trials revealed significant variations in neural responses across session and animal datasets. Within this variation, we found robust long-latency decreases (300 and 2000 ms after stimulus presentation) in gamma-band power (30-80 Hz) in the middle-superficial cortical layers in regions surrounding the activated whisker barrel cortex. This reduction in gamma frequency activity was associated with corresponding decreases in the hemodynamic responses that drive the negative BOLD signal. These findings suggest a close relationship between BOLD responses and neural events that operate over time scales that outlast the initiating sensory stimulus, and provide important insights into the neurophysiological basis of negative neuroimaging signals.

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.

Low-frequency calcium oscillations accompany deoxyhemoglobin oscillations in rat somatosensory cortex

Spontaneous low-frequency oscillations (LFOs) of blood-oxygen-level-dependent (BOLD) signals are used to map brain functional connectivity with functional MRI, but their source is not well understood. Here we used optical imaging to assess whether LFOs from vascular signals covary with oscillatory intracellular calcium (Ca(2+)i) and with local field potentials in the rat's somatosensory cortex. We observed that the frequency of Ca(2+)i oscillations in tissue (∼0.07 Hz) was similar to the LFOs of deoxyhemoglobin (HbR) and oxyhemoglobin (HbO2) in both large blood vessels and capillaries. The HbR and HbO2 fluctuations within tissue correlated with Ca(2+)i oscillations with a lag time of ∼5-6 s. The Ca(2+)i and hemoglobin oscillations were insensitive to hypercapnia. In contrast, cerebral-blood-flow velocity (CBFv) in arteries and veins fluctuated at a higher frequency (∼0.12 Hz) and was sensitive to hypercapnia. However, in parenchymal tissue, CBFv oscillated with peaks at both ∼0.06 Hz and ∼0.12 Hz. Although the higher-frequency CBFv oscillation (∼0.12 Hz) was decreased by hypercapnia, its lower-frequency component (∼0.06 Hz) was not. The sensitivity of the higher CBFV oscillations to hypercapnia, which triggers blood vessel vasodilation, suggests its dependence on vascular effects that are distinct from the LFOs detected in HbR, HbO2, Ca(2+)i, and the lower-frequency tissue CBFv, which were insensitive to hypercapnia. Hemodynamic LFOs correlated both with Ca(2+)i and neuronal firing (local field potentials), indicating that they directly reflect neuronal activity (perhaps also glial). These findings show that HbR fluctuations (basis of BOLD oscillations) are linked to oscillatory cellular activity and detectable throughout the vascular tree (arteries, capillaries, and veins).

Contralateral dissociation between neural activity and cerebral blood volume during recurrent acute focal neocortical seizures

Objective

Whether epileptic events disrupt normal neurovascular coupling mechanisms locally or remotely is unclear. We sought to investigate neurovascular coupling in an acute model of focal neocortical epilepsy, both within the seizure onset zone and in contralateral homotopic cortex.

Methods

Neurovascular coupling in both ipsilateral and contralateral vibrissal cortices of the urethane-anesthetized rat were examined during recurrent 4-aminopyridine (4-AP, 15 mm, 1 μl) induced focal seizures. Local field potential (LFP) and multiunit spiking activity (MUA) were recorded via two bilaterally implanted 16-channel microelectrodes. Concurrent two-dimensional optical imaging spectroscopy was used to produce spatiotemporal maps of cerebral blood volume (CBV).

Results

Recurrent acute seizures in right vibrissal cortex (RVC) produced robust ipsilateral increases in LFP and MUA activity, most prominently in layer 5, that were nonlinearly correlated to local increases in CBV. In contrast, contralateral left vibrissal cortex (LVC) exhibited relatively smaller nonlaminar specific increases in neural activity coupled with a decrease in CBV, suggestive of dissociation between neural and hemodynamic responses.

Significance

These findings provide insights into the impact of epileptic events on the neurovascular unit, and have important implications both for the interpretation of perfusion-based imaging signals in the disorder and understanding the widespread effects of epilepsy.

A critical role for the vascular endothelium in functional neurovascular coupling in the brain

Background

The functional modulation of blood flow in the brain is critical for brain health and is the basis of contrast in functional magnetic resonance imaging. There is evident coupling between increases in neuronal activity and increases in local blood flow; however, many aspects of this neurovascular coupling remain unexplained by current models. Based on the rapid dilation of distant pial arteries during cortical functional hyperemia, we hypothesized that endothelial signaling may play a key role in the long‐range propagation of vasodilation during functional hyperemia in the brain. Although well characterized in the peripheral vasculature, endothelial involvement in functional neurovascular coupling has not been demonstrated.

Methods and Results

We combined in vivo exposed‐cortex multispectral optical intrinsic signal imaging (MS‐OISI) with a novel in vivo implementation of the light‐dye technique to record the cortical hemodynamic response to somatosensory stimulus in rats before and after spatially selective endothelial disruption. We demonstrate that discrete interruption of endothelial signaling halts propagation of stimulus‐evoked vasodilation in pial arteries, and that wide‐field endothelial disruption in pial arteries significantly attenuates the hemodynamic response to stimulus, particularly the early, rapid increase and peak in hyperemia.

Conclusions

Involvement of endothelial pathways in functional neurovascular coupling provides new explanations for the spatial and temporal features of the hemodynamic response to stimulus and could explain previous results that were interpreted as evidence for astrocyte‐mediated control of functional hyperemia. Our results unify many aspects of blood flow regulation in the brain and body and prompt new investigation of direct links between systemic cardiovascular disease and neural deficits.

Investigation of the neurovascular coupling in positive and negative BOLD responses in human brain at 7 T

Decreases in stimulus-dependent blood oxygenation level dependent (BOLD) signal and their underlying neurovascular origins have recently gained considerable interest. In this study a multi-echo, BOLD-corrected vascular space occupancy (VASO) functional magnetic resonance imaging (fMRI) technique was used to investigate neurovascular responses during stimuli that elicit positive and negative BOLD responses in human brain at 7 T. Stimulus-induced BOLD, cerebral blood volume (CBV), and cerebral blood flow (CBF) changes were measured and analyzed in 'arterial' and 'venous' blood compartments in macro- and microvasculature. We found that the overall interplay of mean CBV, CBF and BOLD responses is similar for tasks inducing positive and negative BOLD responses. Some aspects of the neurovascular coupling however, such as the temporal response, cortical depth dependence, and the weighting between 'arterial' and 'venous' contributions, are significantly different for the different task conditions. Namely, while for excitatory tasks the BOLD response peaks at the cortical surface, and the CBV change is similar in cortex and pial vasculature, inhibitory tasks are associated with a maximum negative BOLD response in deeper layers, with CBV showing strong constriction of surface arteries and a faster return to baseline. The different interplays of CBV, CBF and BOLD during excitatory and inhibitory responses suggests different underlying hemodynamic mechanisms.

Coupling between gamma-band power and cerebral blood volume during recurrent acute neocortical seizures

Characterization of neural and hemodynamic biomarkers of epileptic activity that can be measured using non-invasive techniques is fundamental to the accurate identification of the epileptogenic zone (EZ) in the clinical setting. Recently, oscillations at gamma-band frequencies and above (>30 Hz) have been suggested to provide valuable localizing information of the EZ and track cortical activation associated with epileptogenic processes. Although a tight coupling between gamma-band activity and hemodynamic-based signals has been consistently demonstrated in non-pathological conditions, very little is known about whether such a relationship is maintained in epilepsy and the laminar etiology of these signals. Confirmation of this relationship may elucidate the underpinnings of perfusion-based signals in epilepsy and the potential value of localizing the EZ using hemodynamic correlates of pathological rhythms. Here, we use concurrent multi-depth electrophysiology and 2-dimensional optical imaging spectroscopy to examine the coupling between multi-band neural activity and cerebral blood volume (CBV) during recurrent acute focal neocortical seizures in the urethane-anesthetized rat. We show a powerful correlation between gamma-band power (25-90 Hz) and CBV across cortical laminae, in particular layer 5, and a close association between gamma measures and multi-unit activity (MUA). Our findings provide insights into the laminar electrophysiological basis of perfusion-based imaging signals in the epileptic state and may have implications for further research using non-invasive multi-modal techniques to localize epileptogenic tissue.

Influence of CO2 on neurovascular coupling: interaction with dynamic cerebral autoregulation and cerebrovascular reactivity

PaCO₂ affects cerebral blood flow (CBF) and its regulatory mechanisms, but the interaction between neurovascular coupling (NVC), cerebral autoregulation (CA), and cerebrovascular reactivity to CO₂ (CVR), in response to hypercapnia, is not known. Recordings of cerebral blood flow velocity (CBFv), blood pressure (BP), heart rate, and end‐tidal CO₂ (EtCO₂) were performed in 18 subjects during normocapnia and 5% CO₂ inhalation while performing a passive motor paradigm. Together with BP and EtCO₂, a gate signal to represent the effect of stimulation was used as input to a multivariate autoregressive‐moving average model to calculate their separate effects on CBFv. Hypercapnia led to a depression of dynamic CA at rest and during stimulation in both hemispheres (P <0.02) as well as impairment of the NVC response, particularly in the ipsilateral hemisphere (P <0.01). Neither hypercapnia nor the passive motor stimulation influenced CVR. Dynamic CA was not influenced by the motor paradigm during normocapnia. The CBFv step responses to each individual input (BP, EtCO₂, stimulation) allowed identification of the influences of hypercapnia and neuromotor stimulation on CA, CVR, and NVC, which have not been previously described, and also confirmed the depressing effects of hypercapnia on CA and NVC. The stability of CVR during these maneuvers and the lack of influence of stimulation on dynamic CA are novel findings which deserve further investigation. Dynamic multivariate modeling can identify the complex interplay between different CBF regulatory mechanisms and should be recommended for studies involving similar interactions, such as the effects of exercise or posture on cerebral hemodynamics.

The influence of hypercapnia on dynamic cerebral autoregulation (CA), CO₂ vasoreactivity (CVR), and neurovascular coupling (NVC) was described based on a single recording during motor stimulation coupled to a new multivariate modeling approach. Hypercapnia led to a depression of CA at rest and during stimulation in both hemispheres as well as impairment of the NVC response. Neither hypercapnia nor the passive motor stimulation influenced CVR. Dynamic CA was not influenced by the motor paradigm during normocapnia. The stability of CVR during these maneuvers and the lack of influence of stimulation on dynamic CA are novel findings which deserve further investigation.

Device for simultaneous positron emission tomography, laser speckle imaging and RGB reflectometry: validation and application to cortical spreading depression and brain ischemia in rats

Brain function critically relies on the supply with energy substrates (oxygen and glucose) via blood flow. Alterations in energy demand as during neuronal activation induce dynamic changes in substrate fluxes and blood flow. To study the complex system that regulates cerebral metabolism requires the combination of methods for the simultaneous assessment of multiple parameters. We developed a multimodal imaging device to combine positron emission tomography (PET) with laser speckle imaging (LSI) and RGB reflectometry (RGBR). Depending on the radiotracer, PET provides 3-dimensional quantitative information of specific molecular processes, while LSI and RGBR measure cerebral blood flow (CBF) and hemoglobin oxygenation at high temporal and spatial resolution. We first tested the functional capability of each modality within our system and showed that interference between the modalities is negligible. We then cross-calibrated the system by simultaneously measuring absolute CBF using (15)O-H2O PET (CBF(PET)) and the inverse correlation time (ICT), the LSI surrogate for CBF. ICT and CBF(PET) correlated in multiple measurements in individuals as well as across different animals (R(2)=0.87, n=44 measurements) indicating that ICT can be used for absolute quantitative assessment of CBF. To demonstrate the potential of the combined system, we applied it to cortical spreading depression (CSD), a wave of transient cellular depolarization that served here as a model system for neurovascular and neurometabolic coupling. We analyzed time courses of hemoglobin oxygenation and CBF alterations coupled to CSD, and simultaneously measured regional uptake of (18)F-2-fluoro-2-deoxy-D-glucose ((18)F-FDG) used as a radiotracer for regional glucose metabolism, in response to a single CSD and to a cluster of CSD waves. With this unique combination, we characterized the changes in cerebral metabolic rate of oxygen (CMRO2) in real-time and showed a correlation between (18)F-FDG uptake and the number of CSD waves that passed the local tissue. Finally, we examined CSD spontaneously occurring during focal ischemia also referred to as peri-infarct depolarization (PID). In the vicinity of the ischemic territory, we observed PIDs that were characterized by reduced CMRO2 and increased oxygen extraction fraction (OEF), indicating a limitation of oxygen supply. Simultaneously measured PET showed an increased (18)F-FDG uptake in these regions. Our combined system proved to be a novel tool for the simultaneous study of dynamic spatiotemporal alterations of cortical blood flow, oxygen metabolism and glucose consumption under normal and pathologic conditions.

The effects of focal epileptic activity on regional sensory-evoked neurovascular coupling and postictal modulation of bilateral sensory processing

While it is known that cortical sensory dysfunction may occur in focal neocortical epilepsy, it is unknown whether sensory-evoked neurovascular coupling is also disrupted during epileptiform activity. Addressing this open question may help to elucidate both the effects of focal neocortical epilepsy on sensory responses and the neurovascular characteristics of epileptogenic regions in sensory cortex. We therefore examined bilateral sensory-evoked neurovascular responses before, during, and after 4-aminopyridine (4-AP, 15 mmol/L, 1 μL) induced focal neocortical seizures in right vibrissal cortex of the rat. Stimulation consisted of electrical pulse trains (16 seconds, 5 Hz, 1.2 mA) presented to the mystacial pad. Consequent current-source density neural responses and epileptic activity in both cortices and across laminae were recorded via two 16-channel microelectrodes bilaterally implanted in vibrissal cortices. Concurrent two-dimensional optical imaging spectroscopy was used to produce spatiotemporal maps of total, oxy-, and deoxy-hemoglobin concentration. Compared with control, sensory-evoked neurovascular coupling was altered during ictal activity, but conserved postictally in both ipsilateral and contralateral vibrissal cortices, despite neurovascular responses being significantly reduced in the former, and enhanced in the latter. Our results provide insights into sensory-evoked neurovascular dynamics and coupling in epilepsy, and may have implications for the localization of epileptogenic foci and neighboring eloquent cortex.

Using carbogen for calibrated fMRI at 7Tesla: comparison of direct and modelled estimation of the M parameter

Task-evoked changes in cerebral oxygen metabolism can be measured using calibrated functional Magnetic Resonance Imaging (fMRI). This technique requires the use of breathing manipulations such as hypercapnia, hyperoxia or a combination of both to determine a calibration factor M. The M-value is usually obtained by extrapolating the BOLD signal measured during the gas manipulation to its upper theoretical physiological limit using a biophysical model. However, a recently introduced technique uses a combination of increased inspired concentrations of O2 and CO2 to saturate the BOLD signal completely. In this study, we used this BOLD saturation technique to measure M directly at 7Tesla (T). Simultaneous carbogen-7 (7% CO2 in 93% O2) inhalation and visuo-motor task performance were used to elevate venous oxygen saturation in visual and motor areas close to their maximum, and the BOLD signal measured during this manipulation was used as an estimate of M. As accurate estimation of M is crucial for estimation of valid oxidative metabolism values, these directly estimated M-values were assessed and compared with M-values obtained via extrapolation modelling using the generalized calibration model (GCM) on the same dataset. Average M-values measured using both methods were 10.4±3.9% (modelled) and 7.5±2.2% (direct) for a visual-related ROI, and 11.3±5.2% (modelled) and 8.1±2.6% (direct) for a motor-related ROI. Results from this study suggest that, for the CO2 concentration used here, modelling is necessary for the accurate estimation of the M parameter. Neither gas inhalation alone, nor gas inhalation combined with a visuo-motor task, was sufficient to completely saturate venous blood in most subjects. Calibrated fMRI studies should therefore rely on existing models for gas inhalation-based calibration of the BOLD signal.

Does hypercapnia-induced impairment of cerebral autoregulation affect neurovascular coupling? A functional TCD study

Neurovascular coupling (NVC) and dynamic cerebral autoregulation (dCA) are both impaired in the acute phase of ischemic stroke, but their reciprocal interactions are difficult to predict. To clarify these aspects, the present study explored NVC in a healthy volunteer population during a surrogate state of impaired dCA induced by hypercapnia. This study aimed to test whether hypercapnia leads to a depression of NVC through an impairment of dCA. Continuous recordings of middle cerebral arteries cerebral blood flow velocity (CBFv), blood pressure (BP), heart rate, and end-tidal CO2 were performed in 19 right-handed subjects (aged >45 yr) before, during, and after 60 s of a passive paradigm during normocapnia and hypercapnia. The CBFv response was broken down into subcomponents describing the relative contributions of BP (VBP), critical closing pressure (VCrCP), and resistance area product (VRAP). VRAP reflects myogenic activity in response to BP changes, whereas VCrCP is more indicative of metabolic control. The results revealed that hypercapnia significantly affected NVC, with significant reductions in the relative contribution of VCrCP to the paradigm-induced increase in CBFv. The present study suggests that hypercapnia impairs both dCA and NVC, probably acting through an impairment of the metabolic component of CBF control.

Is optical imaging spectroscopy a viable measurement technique for the investigation of the negative BOLD phenomenon? A concurrent optical imaging spectroscopy and fMRI study at high field (7 T)

Traditionally functional magnetic resonance imaging (fMRI) has been used to map activity in the human brain by measuring increases in the Blood Oxygenation Level Dependent (BOLD) signal. Often accompanying positive BOLD fMRI signal changes are sustained negative signal changes. Previous studies investigating the neurovascular coupling mechanisms of the negative BOLD phenomenon have used concurrent 2D-optical imaging spectroscopy (2D-OIS) and electrophysiology (Boorman et al., 2010). These experiments suggested that the negative BOLD signal in response to whisker stimulation was a result of an increase in deoxy-haemoglobin and reduced multi-unit activity in the deep cortical layers. However, Boorman et al. (2010) did not measure the BOLD and haemodynamic response concurrently and so could not quantitatively compare either the spatial maps or the 2D-OIS and fMRI time series directly. Furthermore their study utilised a homogeneous tissue model in which is predominantly sensitive to haemodynamic changes in more superficial layers.

Here we test whether the 2D-OIS technique is appropriate for studies of negative BOLD. We used concurrent fMRI with 2D-OIS techniques for the investigation of the haemodynamics underlying the negative BOLD at 7 Tesla. We investigated whether optical methods could be used to accurately map and measure the negative BOLD phenomenon by using 2D-OIS haemodynamic data to derive predictions from a biophysical model of BOLD signal changes. We showed that despite the deep cortical origin of the negative BOLD response, if an appropriate heterogeneous tissue model is used in the spectroscopic analysis then 2D-OIS can be used to investigate the negative BOLD phenomenon.