Parvalbumin interneuron activity drives fast inhibition-induced vasoconstriction followed by slow substance P-mediated vasodilation.

The role of parvalbumin (PV) interneurons in vascular control is poorly understood. Here, we investigated the hemodynamic responses elicited by optogenetic stimulation of PV interneurons using electrophysiology, functional magnetic resonance imaging (fMRI), wide-field optical imaging (OIS), and pharmacological applications. As a control, forepaw stimulation was used. Stimulation of PV interneurons in the somatosensory cortex evoked a biphasic fMRI response in the photostimulation site and negative fMRI signals in projection regions. Activation of PV neurons engaged two separable neurovascular mechanisms in the stimulation site. First, an early vasoconstrictive response caused by the PV-driven inhibition is sensitive to the brain state affected by anesthesia or wakefulness. Second, a later ultraslow vasodilation lasting a minute is closely dependent on the sum of interneuron multiunit activities, but is not due to increased metabolism, neural or vascular rebound, or increased glial activity. The ultraslow response is mediated by neuropeptide substance P (SP) released from PV neurons under anesthesia, but disappears during wakefulness, suggesting that SP signaling is important for vascular regulation during sleep. Our findings provide a comprehensive perspective about the role of PV neurons in controlling the vascular response.

Relating Pupil Diameter and Blinking to Cortical Activity and Hemodynamics across Arousal States.

Arousal state affects neural activity and vascular dynamics in the cortex, with sleep associated with large changes in the local field potential and increases in cortical blood flow. We investigated the relationship between pupil diameter and blink rate with neural activity and blood volume in the somatosensory cortex in male and female unanesthetized, head-fixed mice. We monitored these variables while the mice were awake, during periods of rapid eye movement (REM), and non-rapid eye movement (NREM) sleep. Pupil diameter was smaller during sleep than in the awake state. Changes in pupil diameter were coherent with both gamma-band power and blood volume in the somatosensory cortex, but the strength and sign of this relationship varied with arousal state. We observed a strong negative correlation between pupil diameter and both gamma-band power and blood volume during periods of awake rest and NREM sleep, although the correlations between pupil diameter and these signals became positive during periods of alertness, active whisking, and REM. Blinking was associated with increases in arousal and decreases in blood volume when the mouse was asleep. Bilateral coherence in gamma-band power and in blood volume dropped following awake blinking, indicating a reset of neural and vascular activity. Using only eye metrics (pupil diameter and eye motion), we could determine the arousal state of the mouse ('Awake,' 'NREM,' 'REM') with >90% accuracy with a 5 s resolution. There is a strong relationship between pupil diameter and hemodynamics signals in mice, reflecting the pronounced effects of arousal on cerebrovascular dynamics.SIGNIFICANCE STATEMENT Determining arousal state is a critical component of any neuroscience experiment. Pupil diameter and blinking are influenced by arousal state, as are hemodynamics signals in the cortex. We investigated the relationship between cortical hemodynamics and pupil diameter and found that pupil diameter was strongly related to the blood volume in the cortex. Mice were more likely to be awake after blinking than before, and blinking resets neural activity. Pupil diameter and eye motion can be used as a reliable, noninvasive indicator of arousal state. As mice transition from wake to sleep and back again over a timescale of seconds, monitoring pupil diameter and eye motion permits the noninvasive detection of sleep events during behavioral or resting-state experiments.

Origins of 1/f-like tissue oxygenation fluctuations in the murine cortex.

The concentration of oxygen in the brain spontaneously fluctuates, and the distribution of power in these fluctuations has a 1/f-like spectra, where the power present at low frequencies of the power spectrum is orders of magnitude higher than at higher frequencies. Though these oscillations have been interpreted as being driven by neural activity, the origin of these 1/f-like oscillations is not well understood. Here, to gain insight of the origin of the 1/f-like oxygen fluctuations, we investigated the dynamics of tissue oxygenation and neural activity in awake behaving mice. We found that oxygen signal recorded from the cortex of mice had 1/f-like spectra. However, band-limited power in the local field potential did not show corresponding 1/f-like fluctuations. When local neural activity was suppressed, the 1/f-like fluctuations in oxygen concentration persisted. Two-photon measurements of erythrocyte spacing fluctuations and mathematical modeling show that stochastic fluctuations in erythrocyte flow could underlie 1/f-like dynamics in oxygenation. These results suggest that the discrete nature of erythrocytes and their irregular flow, rather than fluctuations in neural activity, could drive 1/f-like fluctuations in tissue oxygenation.

Could respiration-driven blood oxygen changes modulate neural activity?

Oxygen is critical for neural metabolism, but under most physiological conditions, oxygen levels in the brain are far more than are required. Oxygen levels can be dynamically increased by increases in respiration rate that are tied to the arousal state of the brain and cognition, and not necessarily linked to exertion by the body. Why these changes in respiration occur when oxygen is already adequate has been a long-standing puzzle. In humans, performance on cognitive tasks can be affected by very high or very low oxygen levels, but whether the physiological changes in blood oxygenation produced by respiration have an appreciable effect is an open question. Oxygen has direct effects on potassium channels, increases the degradation rate of nitric oxide, and is rate limiting for the synthesis of some neuromodulators. We discuss whether oxygenation changes due to respiration contribute to neural dynamics associated with attention and arousal.

High-frequency neuronal signal better explains multi-phase BOLD response.

Visual stimulation-evoked blood-oxygen-level dependent (BOLD) responses can exhibit more complex temporal dynamics than a simple monophasic response. For instance, BOLD responses sometimes include a phase of positive response followed by a phase of post-stimulus undershoot. Whether the BOLD response during these phases reflects the underlying neuronal signal fluctuations or is contributed by non-neuronal physiological factors remains elusive. When presenting blocks of sustained (i.e. DC) light ON-OFF stimulations to unanesthetized rats, we observed that the response following a decrease in illumination (i.e. OFF stimulation-evoked BOLD response) in the visual cortices displayed reproducible multiple phases, including an initial positive BOLD response, followed by an undershoot and then an overshoot before the next ON trial. This multi-phase BOLD response did not result from the entrainment of the periodic stimulation structure. When we measured the neural correlates of these responses, we found that the high-frequency band from the LFP power (300 - 3000 Hz, multi-unit activity (MUA)), but not the power in the gamma band (30 - 100 Hz) exhibited the same multiphasic dynamics as the BOLD signal. This study suggests that the post-stimulus phases of the BOLD response can be better explained by the high-frequency neuronal signal.

Brain-wide ongoing activity is responsible for significant cross-trial BOLD variability.

A notorious issue of task-based functional magnetic resonance imaging (fMRI) is its large cross-trial variability. To quantitatively characterize this variability, the blood oxygenation level-dependent (BOLD) signal can be modeled as a linear summation of a stimulation-relevant and an ongoing (i.e. stimulation-irrelevant) component. However, systematic investigation on the spatiotemporal features of the ongoing BOLD component and how these features affect the BOLD response is still lacking. Here we measured fMRI responses to light onsets and light offsets in awake rats. The neuronal response was simultaneously recorded with calcium-based fiber photometry. We established that between-region BOLD signals were highly correlated brain-wide at zero time lag, including regions that did not respond to visual stimulation, suggesting that the ongoing activity co-fluctuates across the brain. Removing this ongoing activity reduced cross-trial variability of the BOLD response by ~30% and increased its coherence with the Ca2+ signal. Additionally, the negative ongoing BOLD activity sometimes dominated over the stimulation-driven response and contributed to the post-stimulation BOLD undershoot. These results suggest that brain-wide ongoing activity is responsible for significant cross-trial BOLD variability, and this component can be reliably quantified and removed to improve the reliability of fMRI response. Importantly, this method can be generalized to virtually all fMRI experiments without changing stimulation paradigms.

Photoacoustic imaging for microcirculation.

Microcirculation facilitates the blood-tissue exchange of nutrients and regulates blood perfusion. It is, therefore, essential in maintaining tissue health. Aberrations in microcirculation are potentially indicative of underlying cardiovascular and metabolic pathologies. Thus, quantitative information about it is of great clinical relevance. Photoacoustic imaging (PAI) is a capable technique that relies on the generation of imaging contrast via the absorption of light and can image at micron-scale resolution. PAI is especially desirable to map microvasculature as hemoglobin strongly absorbs light and can generate a photoacoustic signal. This paper reviews the current state of the art for imaging microvascular networks using photoacoustic imaging. We further describe how quantitative information about blood dynamics such as the total hemoglobin concentration, oxygen saturation, and blood flow rate is obtained using PAI. We also discuss its importance in understanding key pathophysiological processes in neurovascular, cardiovascular, ophthalmic, and cancer research fields. We then discuss the current challenges and limitations of PAI and the approaches that can help overcome these limitations. Finally, we provide the reader with an overview of future trends in the field of PAI for imaging microcirculation.

Awake mouse brain photoacoustic and optical imaging through a transparent ultrasound cranial window.

Optical resolution photoacoustic microscopy (OR-PAM) can map the cerebral vasculature at capillary-level resolution. However, the OR-PAM setup's bulky imaging head makes awake mouse brain imaging challenging and inhibits its integration with other optical neuroimaging modalities. Moreover, the glass cranial windows used for optical microscopy are unsuitable for OR-PAM due to the acoustic impedance mismatch between the glass plate and the tissue. To overcome these challenges, we propose a lithium niobate based transparent ultrasound transducer (TUT) as a cranial window on a thinned mouse skull. The TUT cranial window simplifies the imaging head considerably due to its dual functionality as an optical window and ultrasound transducer. The window remains stable for six weeks, with no noticeable inflammation and minimal bone regrowth. The TUT window's potential is demonstrated by imaging the awake mouse cerebral vasculature using OR-PAM, intrinsic optical signal imaging, and two-photon microscopy. The TUT cranial window can potentially also be used for ultrasound stimulation and simultaneous multimodal imaging of the awake mouse brain.

Iliski, a software for robust calculation of transfer functions.

Understanding the relationships between biological processes is paramount to unravel pathophysiological mechanisms. These relationships can be modeled with Transfer Functions (TFs), with no need of a priori hypotheses as to the shape of the transfer function. Here we present Iliski, a software dedicated to TFs computation between two signals. It includes different pre-treatment routines and TF computation processes: deconvolution, deterministic and non-deterministic optimization algorithms that are adapted to disparate datasets. We apply Iliski to data on neurovascular coupling, an ensemble of cellular mechanisms that link neuronal activity to local changes of blood flow, highlighting the software benefits and caveats in the computation and evaluation of TFs. We also propose a workflow that will help users to choose the best computation according to the dataset. Iliski is available under the open-source license CC BY 4.0 on GitHub (https://github.com/alike-aydin/Iliski) and can be used on the most common operating systems, either within the MATLAB environment, or as a standalone application.

Spatial and temporal patterns of nitric oxide diffusion and degradation drive emergent cerebrovascular dynamics.

Nitric oxide (NO) is a gaseous signaling molecule that plays an important role in neurovascular coupling. NO produced by neurons diffuses into the smooth muscle surrounding cerebral arterioles, driving vasodilation. However, the rate of NO degradation in hemoglobin is orders of magnitude higher than in brain tissue, though how this might impact NO signaling dynamics is not completely understood. We used simulations to investigate how the spatial and temporal patterns of NO generation and degradation impacted dilation of a penetrating arteriole in cortex. We found that the spatial location of NO production and the size of the vessel both played an important role in determining its responsiveness to NO. The much higher rate of NO degradation and scavenging of NO in the blood relative to the tissue drove emergent vascular dynamics. Large vasodilation events could be followed by post-stimulus constrictions driven by the increased degradation of NO by the blood, and vasomotion-like 0.1-0.3 Hz oscillations could also be generated. We found that these dynamics could be enhanced by elevation of free hemoglobin in the plasma, which occurs in diseases such as malaria and sickle cell anemia, or following blood transfusions. Finally, we show that changes in blood flow during hypoxia or hyperoxia could be explained by altered NO degradation in the parenchyma. Our simulations suggest that many common vascular dynamics may be emergent phenomena generated by NO degradation by the blood or parenchyma.

Effects of Voluntary Locomotion and Calcitonin Gene-Related Peptide on the Dynamics of Single Dural Vessels in Awake Mice.

The dura mater is a vascularized membrane surrounding the brain and is heavily innervated by sensory nerves. Our knowledge of the dural vasculature has been limited to pathological conditions, such as headaches, but little is known about the dural blood flow regulation during behavior. To better understand the dynamics of dural vessels during behavior, we used two-photon laser scanning microscopy (2PLSM) to measure the diameter changes of single dural and pial vessels in the awake mouse during voluntary locomotion. Surprisingly, we found that voluntary locomotion drove the constriction of dural vessels, and the dynamics of these constrictions could be captured with a linear convolution model. Dural vessel constrictions did not mirror the large increases in intracranial pressure (ICP) during locomotion, indicating that dural vessel constriction was not caused passively by compression. To study how behaviorally driven dynamics of dural vessels might be altered in pathological states, we injected the vasodilator calcitonin gene-related peptide (CGRP), which induces headache in humans. CGRP dilated dural, but not pial, vessels and significantly reduced spontaneous locomotion but did not block locomotion-induced constrictions in dural vessels. Sumatriptan, a drug commonly used to treat headaches, blocked the vascular and behavioral the effects of CGRP. These findings suggest that, in the awake animal, the diameters of dural vessels are regulated dynamically during behavior and during drug-induced pathological states.

Time to wake up: Studying neurovascular coupling and brain-wide circuit function in the un-anesthetized animal.

Functional magnetic resonance imaging (fMRI) has allowed the noninvasive study of task-based and resting-state brain dynamics in humans by inferring neural activity from blood-oxygenation-level dependent (BOLD) signal changes. An accurate interpretation of the hemodynamic changes that underlie fMRI signals depends on the understanding of the quantitative relationship between changes in neural activity and changes in cerebral blood flow, oxygenation and volume. While there has been extensive study of neurovascular coupling in anesthetized animal models, anesthesia causes large disruptions of brain metabolism, neural responsiveness and cardiovascular function. Here, we review work showing that neurovascular coupling and brain circuit function in the awake animal are profoundly different from those in the anesthetized state. We argue that the time is right to study neurovascular coupling and brain circuit function in the awake animal to bridge the physiological mechanisms that underlie animal and human neuroimaging signals, and to interpret them in light of underlying neural mechanisms. Lastly, we discuss recent experimental innovations that have enabled the study of neurovascular coupling and brain-wide circuit function in un-anesthetized and behaving animal models.

Neurovascular coupling and decoupling in the cortex during voluntary locomotion.

Hemodynamic signals are widely used to infer neural activity in the brain. We tested the hypothesis that hemodynamic signals faithfully report neural activity during voluntary behaviors by measuring cerebral blood volume (CBV) and neural activity in the somatosensory cortex and frontal cortex of head-fixed mice during locomotion. Locomotion induced a large and robust increase in firing rate and gamma-band (40-100 Hz) power in the local field potential in the limb representations in somatosensory cortex, and was accompanied by increases in CBV, demonstrating that hemodynamic signals are coupled with neural activity in this region. However, in the frontal cortex, CBV did not change during locomotion, but firing rate and gamma-band power both increased, indicating a decoupling of neural activity from the hemodynamic signal. These results show that hemodynamic signals are not faithful indicators of the mean neural activity in the frontal cortex during locomotion; thus, the results from fMRI and other hemodynamic imaging methodologies for studying neural processes must be interpreted with caution.

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

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

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

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

Venous cerebral blood volume increase during voluntary locomotion reflects cardiovascular changes.

Understanding how changes in the cardiovascular system contribute to cerebral blood flow (CBF) and volume (CBV) increases is critical for interpreting hemodynamic signals. Here we investigated how systemic cardiovascular changes affect the cortical hemodynamic response during voluntary locomotion. In the mouse, voluntary locomotion drives an increase in cortical CBF and arterial CBV that is localized to the forelimb/hindlimb representation in the somatosensory cortex, as well as a diffuse venous CBV increase. To determine if the heart rate increases that accompany locomotion contribute to locomotion-induced CBV and CBF increases, we occluded heart rate increases with the muscarinic cholinergic receptor antagonist glycopyrrolate, and reduced heart rate with the ß1-adrenergic receptor antagonist atenolol. We quantified the effects of these cardiovascular manipulations on CBV and CBF dynamics by comparing the hemodynamic response functions (HRF) to locomotion across these conditions. Neither the CBF HRF nor the arterial component of the CBV HRF was significantly affected by pharmacological disruption of the heart rate. In contrast, the amplitude and spatial extent of the venous component of the CBV HRF were decreased by atenolol. These results suggest that the increase in venous CBV during locomotion was partially driven by peripheral cardiovascular changes, whereas CBF and arterial CBV increases associated with locomotion reflect central processes.

Brief anesthesia, but not voluntary locomotion, significantly alters cortical temperature.

Changes in brain temperature can alter electrical properties of neurons and cause changes in behavior. However, it is not well understood how behaviors, like locomotion, or experimental manipulations, like anesthesia, alter brain temperature. We implanted thermocouples in sensorimotor cortex of mice to understand how cortical temperature was affected by locomotion, as well as by brief and prolonged anesthesia. Voluntary locomotion induced small (∼ 0.1 °C) but reliable increases in cortical temperature that could be described using a linear convolution model. In contrast, brief (90-s) exposure to isoflurane anesthesia depressed cortical temperature by ∼ 2 °C, which lasted for up to 30 min after the cessation of anesthesia. Cortical temperature decreases were not accompanied by a concomitant decrease in the γ-band local field potential power, multiunit firing rate, or locomotion behavior, which all returned to baseline within a few minutes after the cessation of anesthesia. In anesthetized animals where core body temperature was kept constant, cortical temperature was still > 1 °C lower than in the awake animal. Thermocouples implanted in the subcortex showed similar temperature changes under anesthesia, suggesting these responses occur throughout the brain. Two-photon microscopy of individual blood vessel dynamics following brief isoflurane exposure revealed a large increase in vessel diameter that ceased before the brain temperature significantly decreased, indicating cerebral heat loss was not due to increased cerebral blood vessel dilation. These data should be considered in experimental designs recording in anesthetized preparations, computational models relating temperature and neural activity, and awake-behaving methods that require brief anesthesia before experimental procedures.

Robust and fragile aspects of cortical blood flow in relation to the underlying angioarchitecture.

We review the organizational principles of the cortical vasculature and the underlying patterns of blood flow under normal conditions and in response to occlusion of single vessels. The cortex is sourced by a two-dimensional network of pial arterioles that feeds a three-dimensional network of subsurface microvessels in close proximity to neurons and glia. Blood flow within the surface and subsurface networks is largely insensitive to occlusion of a single vessel within either network. However, the penetrating arterioles that connect the pial network to the subsurface network are bottlenecks to flow; occlusion of even a single penetrating arteriole results in the death of a 500 µm diameter cylinder of cortical tissue despite the potential for collateral flow through microvessels. This pattern of flow is consistent with that calculated from a full reconstruction of the angioarchitecture. Conceptually, collateral flow is insufficient to compensate for the occlusion of a penetrating arteriole because penetrating venules act as shunts of blood that flows through collaterals. Future directions that stem from the analysis of the angioarchitecture concern cellular-level issues, in particular the regulation of blood flow within the subsurface microvascular network, and system-level issues, in particular the role of penetrating arteriole occlusions in human cognitive impairment.

Fluctuating and sensory-induced vasodynamics in rodent cortex extend arteriole capacity.

Neural activity in the brain is followed by localized changes in blood flow and volume. We address the relative change in volume for arteriole vs. venous blood within primary vibrissa cortex of awake, head-fixed mice. Two-photon laser-scanning microscopy was used to measure spontaneous and sensory evoked changes in flow and volume at the level of single vessels. We find that arterioles exhibit slow (<1 Hz) spontaneous increases in their diameter, as well as pronounced dilation in response to both punctate and prolonged stimulation of the contralateral vibrissae. In contrast, venules dilate only in response to prolonged stimulation. We conclude that stimulation that occurs on the time scale of natural stimuli leads to a net increase in the reservoir of arteriole blood. Thus, a "bagpipe" model that highlights arteriole dilation should augment the current "balloon" model of venous distension in the interpretation of fMRI images.

Multiparametric Brain Hemodynamics Imaging Using a Combined Ultrafast Ultrasound and Photoacoustic System.

Studying brain-wide hemodynamic responses to different stimuli at high spatiotemporal resolutions can help gain new insights into the mechanisms of neuro- diseases and -disorders. Nonetheless, this task is challenging, primarily due to the complexity of neurovascular coupling, which encompasses interdependent hemodynamic parameters including cerebral blood volume (CBV), cerebral blood flow (CBF), and cerebral oxygen saturation (SO2). The current brain imaging technologies exhibit inherent limitations in resolution, sensitivity, and imaging depth, restricting their capacity to comprehensively capture the intricacies of cerebral functions. To address this, a multimodal functional ultrasound and photoacoustic (fUSPA) imaging platform is reported, which integrates ultrafast ultrasound and multispectral photoacoustic imaging methods in a compact head-mountable device, to quantitatively map individual dynamics of CBV, CBF, and SO2 as well as contrast agent enhanced brain imaging at high spatiotemporal resolutions. Following systematic characterization, the fUSPA system is applied to study brain-wide cerebrovascular reactivity (CVR) at single-vessel resolution via relative changes in CBV, CBF, and SO2 in response to hypercapnia stimulation. These results show that cortical veins and arteries exhibit differences in CVR in the stimulated state and consistent anti-correlation in CBV oscillations during the resting state, demonstrating the multiparametric fUSPA system's unique capabilities in investigating complex mechanisms of brain functions.