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

Dissecting Multiparametric Cerebral Hemodynamics using Integrated Ultrafast Ultrasound and Multispectral Photoacoustic Imaging

Understanding brain-wide hemodynamic responses to different stimuli at high spatiotemporal resolutions can help study neuro-disorders and brain functions. However, the existing brain imaging technologies have limited resolution, sensitivity, imaging depth and provide information about only one or two hemodynamic parameters. To address this, we propose a multimodal functional ultrasound and photoacoustic (fUSPA) imaging platform, which integrates ultrafast ultrasound and multispectral photoacoustic imaging methods in a compact head-mountable device, to quantitatively map cerebral blood volume (CBV), cerebral blood flow (CBF), oxygen saturation (SO2) dynamics as well as contrast agent enhanced brain imaging with high spatiotemporal resolutions. After systematic characterization, the fUSPA system was applied to quantitatively study the changes in brain hemodynamics and vascular reactivity at single vessel resolution in response to hypercapnia stimulation. Our results show an overall increase in brain-wide CBV, CBF, and SO2, but regional differences in singular cortical veins and arteries and a reproducible anti-correlation pattern between venous and cortical hemodynamics, demonstrating the capabilities of the fUSPA system for providing multiparametric cerebrovascular information at high-resolution and sensitivity, that can bring insights into the complex mechanisms of neurodiseases.

Arousal state transitions occlude sensory-evoked neurovascular coupling in neonatal mice

In the adult sensory cortex, increases in neural activity elicited by sensory stimulation usually drive vasodilation mediated by neurovascular coupling. However, whether neurovascular coupling is the same in neonatal animals as adults is controversial, as both canonical and inverted responses have been observed. We investigated the nature of neurovascular coupling in unanesthetized neonatal mice using optical imaging, electrophysiology, and BOLD fMRI. We find in neonatal (postnatal day 15, P15) mice, sensory stimulation induces a small increase in blood volume/BOLD signal, often followed by a large decrease in blood volume. An examination of arousal state of the mice revealed that neonatal mice were asleep a substantial fraction of the time, and that stimulation caused the animal to awaken. As cortical blood volume is much higher during REM and NREM sleep than the awake state, awakening occludes any sensory-evoked neurovascular coupling. When neonatal mice are stimulated during an awake period, they showed relatively normal (but slowed) neurovascular coupling, showing that that the typically observed constriction is due to arousal state changes. These result show that sleep-related vascular changes dominate over any sensory-evoked changes, and hemodynamic measures need to be considered in the context of arousal state changes.

Arousal state transitions occlude sensory-evoked neurovascular coupling in neonatal mice

In the adult sensory cortex, increases in neural activity elicited by sensory stimulation usually drives vasodilation mediated by neurovascular coupling. However, whether neurovascular coupling is the same in neonatal animals as adults is controversial, as both canonical and inverted responses have been observed. We investigated the nature of neurovascular coupling in unanesthetized neonatal mice using optical imaging, electrophysiology, and BOLD fMRI. We find in neonatal (postnatal day 15, P15) mice, sensory stimulation induces a small increase in blood volume/BOLD signal, often followed by a large decrease in blood volume. An examination of arousal state of the mice revealed that neonatal mice were asleep a substantial fraction of the time, and that stimulation caused the animal to awaken. As cortical blood volume is much higher during REM and NREM sleep than the awake state, awakening occludes any sensory-evoked neurovascular coupling. When neonatal mice are stimulated during an awake period, they showed relatively normal (but slowed) neurovascular coupling, showing that that the typically observed constriction is due to arousal state changes. These result show that sleep-related vascular changes dominate over any sensory-evoked changes, and hemodynamic measures need to be considered in the context of arousal state changes.

Significance Statement

In the adult brain, increases in neural activity are often followed by vasodilation, allowing activity to be monitored using optical or magnetic resonance imaging. However, in neonates, sensory stimulation can drive vasoconstriction, whose origin was not understood. We used optical and magnetic resonance imaging approaches to investigate hemodynamics in neonatal mice. We found that sensory-induced vasoconstriction occurred when the mice were asleep, as sleep is associated with dilation of the vasculature of the brain relative to the awake state. The stimulus awakens the mice, causing a constriction due to the arousal state change. Our study shows the importance of monitoring arousal state, particularly when investigating subjects that may sleep, and the dominance arousal effects on brain hemodynamics.

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.

The Influence of Movement on the Cerebrospinal Fluid Pressure of the American Alligator (Alligator mississippiensis)

This study was undertaken to document how the cerebrospinal fluid (CSF) pressure varied during movements and physiological activities. Using surgically implanted pressure catheters; the CSF pressure was recorded from sub-adult American alligators (Alligator mississippiensis) under anesthesia and post-recovery. Pressures were recorded during physiological activities (the cardiac cycle; passive and active ventilation); manual manipulation of the anesthetized animals (foot sweeps; tail oscillations; and body bends); as well as voluntary movements post-recovery (changes in body tone; defensive strikes; and locomotion). The CSF pulsations associated with the cardiac cycle had the lowest mean amplitude (3.7 mm Hg); during active ventilation and defensive strikes; the alligators routinely generated CSF pressure spikes in excess of 100 mm Hg. The recorded CSF pressures appear to be caused by a variety of mechanisms including vascular pressure; fluid inertia; and possible physical displacement of the spinal cord. The results of the study suggest that any model of CSF dynamics or perfusion should incorporate the episodic high-pressure CSF pulsations associated with movement.

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.

Treadmill locomotion in the American alligator (Alligator mississippiensis) produces dynamic changes in intracranial cerebrospinal fluid pressure

To examine the influence of movement on cerebrospinal fluid (CSF) dynamics, intracranial subdural pressure recordings were taken from sub-adult alligators (Alligator mississippiensis) locomoting on a treadmill. Pressure recordings documenting the cardiac, ventilatory, and barostatic influences on the CSF were in good agreement with previous studies. During locomotion the CSF exhibits sinusoidal patterns of pressure change that spanned a mean amplitude of 56 mm Hg, some 16 × the amplitude of the cardiac-linked pulsations. These sinusoidal CSF pulsations were closely linked to the locomotor kinematics, particularly the lateral oscillations of the alligator's head. Data recorded from the freely moving alligators suggest that fluid inertia, body cavity pressures, and likely other factors all influence the CSF pressure. The clear relationship between movement and CSF pressure described in this study suggests that the paucity of studies examining human CSF dynamics during movement should be addressed.

Coupling between cerebrovascular oscillations and CSF flow fluctuations during wakefulness: An fMRI study

It is commonly believed that cerebrospinal fluid (CSF) movement is facilitated by blood vessel wall movements (i.e., hemodynamic oscillations) in the brain. A coherent pattern of low frequency hemodynamic oscillations and CSF movement was recently found during non-rapid eye movement (NREM) sleep via functional MRI. This finding raises other fundamental questions: 1) the explanation of coupling between hemodynamic oscillations and CSF movement from fMRI signals; 2) the existence of the coupling during wakefulness; 3) the direction of CSF movement. In this resting state fMRI study, we proposed a mechanical model to explain the coupling between hemodynamics and CSF movement through the lens of fMRI. Time delays between CSF movement and global hemodynamics were calculated. The observed delays between hemodynamics and CSF movement match those predicted by the model. Moreover, by conducting separate fMRI scans of the brain and neck, we confirmed the low frequency CSF movement at the fourth ventricle is bidirectional. Our finding also demonstrates that CSF movement is facilitated by changes in cerebral blood volume mainly in the low frequency range, even when the individual is awake.

The effects of locomotion on sensory-evoked haemodynamic responses in the cortex of awake mice

Investigating neurovascular coupling in awake rodents is becoming ever more popular due, in part, to our increasing knowledge of the profound impacts that anaesthesia can have upon brain physiology. Although awake imaging brings with it many advantages, we still do not fully understand how voluntary locomotion during imaging affects sensory-evoked haemodynamic responses. In this study we investigated how evoked haemodynamic responses can be affected by the amount and timing of locomotion. Using an awake imaging set up, we used 2D-Optical Imaging Spectroscopy (2D-OIS) to measure changes in cerebral haemodynamics within the sensory cortex of the brain during either 2 s whisker stimulation or spontaneous (no whisker stimulation) experiments, whilst animals could walk on a spherical treadmill. We show that locomotion alters haemodynamic responses. The amount and timing of locomotion relative to whisker stimulation is important, and can significantly impact sensory-evoked haemodynamic responses. If locomotion occurred before or during whisker stimulation, the amplitude of the stimulus-evoked haemodynamic response was significantly altered. Therefore, monitoring of locomotion during awake imaging is necessary to ensure that conclusions based on comparisons of evoked haemodynamic responses (e.g., between control and disease groups) are not confounded by the effects of locomotion.

Behavioral and physiological monitoring for awake neurovascular coupling experiments: a how-to guide

Significance: Functional brain imaging in awake animal models is a popular and powerful technique that allows the investigation of neurovascular coupling (NVC) under physiological conditions. However, ubiquitous facial and body motions (fidgeting) are prime drivers of spontaneous fluctuations in neural and hemodynamic signals. During periods without movement, animals can rapidly transition into sleep, and the hemodynamic signals tied to arousal state changes can be several times larger than sensory-evoked responses. Given the outsized influence of facial and body motions and arousal signals in neural and hemodynamic signals, it is imperative to detect and monitor these events in experiments with un-anesthetized animals. Aim: To cover the importance of monitoring behavioral state in imaging experiments using un-anesthetized rodents, and describe how to incorporate detailed behavioral and physiological measurements in imaging experiments. Approach: We review the effects of movements and sleep-related signals (heart rate, respiration rate, electromyography, intracranial pressure, whisking, and other body movements) on brain hemodynamics and electrophysiological signals, with a focus on head-fixed experimental setup. We summarize the measurement methods currently used in animal models for detection of those behaviors and arousal changes. We then provide a guide on how to incorporate this measurements with functional brain imaging and electrophysiology measurements. Results: We provide a how-to guide on monitoring and interpreting a variety of physiological signals and their applications to NVC experiments in awake behaving mice. Conclusion: This guide facilitates the application of neuroimaging in awake animal models and provides neuroscientists with a standard approach for monitoring behavior and other associated physiological parameters in head-fixed animals.

Central command-related increases in blood velocity of anterior cerebral artery and prefrontal oxygenation at the onset of voluntary tapping

The anterior cerebral artery (ACA) supplies blood predominantly to the frontal lobe including the prefrontal cortex. Our laboratory reported that prefrontal oxygenated-hemoglobin concentration (Oxy-Hb) increased before and at exercise onset, as long as exercise is arbitrarily started. Moreover, the increased prefrontal oxygenation seems independent of both exercise intensity and muscle mass. If so, mean blood velocity of the ACA (ACA_BV) should increase with "very light motor effort," concomitantly with the preexercise and initial increase in prefrontal Oxy-Hb. This study aimed to examine the responses in ACA_BV and vascular conductance index (ACA_VCI) of the ACA as well as prefrontal Oxy-Hb during arbitrary or cued finger tapping in 12 subjects, an activity with a Borg scale perceived exertion rating of 7 (median). With arbitrary start, ACA_BV increased at tapping onset (14 ± 9%) via an elevation in ACA_VCI. Likewise, prefrontal Oxy-Hb increased at the onset of tapping with a time course resembling that of ACA_BV. A positive cross correlation between the initial changes in ACA_BV and prefrontal Oxy-Hb was found significant in 67% of subjects, having a time lag of 2 s, whereas a positive linear regression between them was significant in 75% of subjects. When tapping was forced to start by cue, the initial increases in ACA_BV, ACA_VCI, and prefrontal Oxy-Hb were delayed and blunted as compared with an arbitrary start. Thus, active vasodilatation of the ACA vascular bed occurs at tapping onset, as long as tapping is arbitrarily started, and contributes to immediate increases in blood flow and prefrontal oxygenation.NEW & NOTEWORTHY Anterior cerebral artery blood velocity and vascular conductance index along with prefrontal oxygenated-hemoglobin concentration all increased at the onset of finger tapping, peaking immediately after tapping onset, as long as tapping was arbitrarily started. Positive cross correlation and linear regression between the increases in ACA_BV and prefrontal Oxy-Hb were significant in 67%-75% of subjects. Active vasodilatation of the ACA vascular bed occurs with arbitrary tapping onset and contributes to increased ACA_BV and prefrontal oxygenation.

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.

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.

Iliski is a software helping the user to find the relationship between two sets of data, namely transfer functions. Although transfer functions are widely used in many scientific fields to link two signals, their computation can be tricky due to data features such as multisource noise, or to specific shape requirements imposed by the nature of the signals, e.g. in biological data. Iliski offers a user-friendly graphical interface to ease the computation of transfer functions for both experienced and users with no coding skills. It proposes several signal pre-processing methods and allows rapid testing of different computing approaches, either based on deconvolution or on optimization of multi-parametric functions. This article, combined with a User Manual, provides a detailed description of Iliski functionalities and a thorough description of the advantages and drawbacks of each computing method using experimental biological data. In the era of Big Data, scientists strive to find new models for patho-physiological mechanisms, and Iliski fulfils the requirements of rigorous, flexible, and fast data driven hypothesis testing.

Adaptive modulation of brain hemodynamics across stereotyped running episodes

During locomotion, theta and gamma rhythms are essential to ensure timely communication between brain structures. However, their metabolic cost and contribution to neuroimaging signals remain elusive. To finely characterize neurovascular interactions during locomotion, we simultaneously recorded mesoscale brain hemodynamics using functional ultrasound (fUS) and local field potentials (LFP) in numerous brain structures of freely-running overtrained rats. Locomotion events were reliably followed by a surge in blood flow in a sequence involving the retrosplenial cortex, dorsal thalamus, dentate gyrus and CA regions successively, with delays ranging from 0.8 to 1.6 seconds after peak speed. Conversely, primary motor cortex was suppressed and subsequently recruited during reward uptake. Surprisingly, brain hemodynamics were strongly modulated across trials within the same recording session; cortical blood flow sharply decreased after 10-20 runs, while hippocampal responses strongly and linearly increased, particularly in the CA regions. This effect occurred while running speed and theta activity remained constant and was accompanied by an increase in the power of hippocampal, but not cortical, high-frequency oscillations (100-150 Hz). Our findings reveal distinct vascular subnetworks modulated across fast and slow timescales and suggest strong hemodynamic adaptation, despite the repetition of a stereotyped behavior.

Increased prefrontal oxygenation prior to and at the onset of over-ground locomotion in humans

We found using wireless near-infrared spectroscopy that prefrontal oxygenation increased before the onset of arbitrary over-ground walking, whereas the preexercise increase was absent when walking was suddenly started by cue. The difference in prefrontal oxygenation between start modes (considered related to central command) preceded heart rate response variances and demonstrated a positive relationship with the difference in heart rate. The central command-related prefrontal activity may contribute to cardiac adjustment, synchronized with the beginning of over-ground walking.

Neurovascular coupling and bilateral connectivity during NREM and REM sleep

To understand how arousal state impacts cerebral hemodynamics and neurovascular coupling, we monitored neural activity, behavior, and hemodynamic signals in un-anesthetized, head-fixed mice. Mice frequently fell asleep during imaging, and these sleep events were interspersed with periods of wake. During both NREM and REM sleep, mice showed large increases in cerebral blood volume ([HbT]) and arteriole diameter relative to the awake state, two to five times larger than those evoked by sensory stimulation. During NREM, the amplitude of bilateral low-frequency oscillations in [HbT] increased markedly, and coherency between neural activity and hemodynamic signals was higher than the awake resting and REM states. Bilateral correlations in neural activity and [HbT] were highest during NREM, and lowest in the awake state. Hemodynamic signals in the cortex are strongly modulated by arousal state, and changes during sleep are substantially larger than sensory-evoked responses.

nNOS-expressing interneurons control basal and behaviorally evoked arterial dilation in somatosensory cortex of mice

Cortical neural activity is coupled to local arterial diameter and blood flow. However, which neurons control the dynamics of cerebral arteries is not well understood. We dissected the cellular mechanisms controlling the basal diameter and evoked dilation in cortical arteries in awake, head-fixed mice. Locomotion drove robust arterial dilation, increases in gamma band power in the local field potential (LFP), and increases calcium signals in pyramidal and neuronal nitric oxide synthase (nNOS)-expressing neurons. Chemogenetic or pharmocological modulation of overall neural activity up or down caused corresponding increases or decreases in basal arterial diameter. Modulation of pyramidal neuron activity alone had little effect on basal or evoked arterial dilation, despite pronounced changes in the LFP. Modulation of the activity of nNOS-expressing neurons drove changes in the basal and evoked arterial diameter without corresponding changes in population neural activity.

Ultra-slow Oscillations in fMRI and Resting-State Connectivity: Neuronal and Vascular Contributions and Technical Confounds

Ultra-slow, ∼0.1-Hz variations in the oxygenation level of brain blood are widely used as an fMRI-based surrogate of "resting-state" neuronal activity. The temporal correlations among these fluctuations across the brain are interpreted as "functional connections" for maps and neurological diagnostics. Ultra-slow variations in oxygenation follow a cascade. First, they closely track changes in arteriole diameter. Second, interpretable functional connections arise when the ultra-slow changes in amplitude of γ-band neuronal oscillations, which are shared across even far-flung but synaptically connected brain regions, entrain the ∼0.1-Hz vasomotor oscillation in diameter of local arterioles. Significant confounds to estimates of functional connectivity arise from residual vasomotor activity as well as arteriole dynamics driven by self-generated movements and subcortical common modulatory inputs. Last, methodological limitations of fMRI can lead to spurious functional connections. The neuronal generator of ultra-slow variations in γ-band amplitude, including that associated with self-generated movements, remains an open issue.

Functional hyperemia drives fluid exchange in the paravascular space

The brain lacks a conventional lymphatic system to remove metabolic waste. It has been proposed that directional fluid movement through the arteriolar paravascular space (PVS) promotes metabolite clearance. We performed simulations to examine if arteriolar pulsations and dilations can drive directional CSF flow in the PVS and found that arteriolar wall movements do not drive directional CSF flow. We propose an alternative method of metabolite clearance from the PVS, namely fluid exchange between the PVS and the subarachnoid space (SAS). In simulations with compliant brain tissue, arteriolar pulsations did not drive appreciable fluid exchange between the PVS and the SAS. However, when the arteriole dilated, as seen during functional hyperemia, there was a marked exchange of fluid. Simulations suggest that functional hyperemia may serve to increase metabolite clearance from the PVS. We measured blood vessels and brain tissue displacement simultaneously in awake, head-fixed mice using two-photon microscopy. These measurements showed that brain deforms in response to pressure changes in PVS, consistent with our simulations. Our results show that the deformability of the brain tissue needs to be accounted for when studying fluid flow and metabolite transport.

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.

Functional imaging evidence for task-induced deactivation and disconnection of a major default mode network hub in the mouse brain

Our report first validates an imaging technique, functional ultrasound, an experimental approach providing significantly higher sensitivity and spatiotemporal resolution than fMRI, for the study of sensory stimulation-induced changes in activation and connectivity patterns in the awake mouse brain. Next, by using this approach in lightly sedated conscious mice, we provide functional evidence supporting an evolutionary preserved function of the default mode network, a set of highly connected brain areas, which is characteristically impaired in major neuropsychiatric diseases in humans. Uncovering the existence and function of the default mode network in mice, the main preclinical model organism for neuropsychiatric diseases, opens new avenues of high translational relevance for brain research and drug development.

The default mode network (DMN) has been defined in functional brain imaging studies as a set of highly connected brain areas, which are active during wakeful rest and inactivated during task-based stimulation. DMN function is characteristically impaired in major neuropsychiatric diseases, emphasizing its interest for translational research. However, in the mouse, a major preclinical rodent model, there is still no functional imaging evidence supporting DMN deactivation and deconnection during high-demanding cognitive/sensory tasks. Here we have developed functional ultrasound (fUS) imaging to properly visualize both activation levels and functional connectivity patterns, in head-restrained awake and behaving mice, and investigated their modulation during a sensory-task, whisker stimulation. We identified reproducible and highly symmetric resting-state networks, with overall connectivity strength directly proportional to the wakefulness level of the animal. We show that unilateral whisker stimulation leads to the expected activation of the contralateral barrel cortex in lightly sedated mice, while interhemispheric inhibition reduces activity in the ipsilateral barrel cortex. Whisker stimulation also leads to elevated bilateral connectivity in the hippocampus. Importantly, in addition to functional changes in these major hubs of tactile information processing, whisker stimulation during genuine awake resting-state periods leads to highly specific reductions both in activation and interhemispheric correlation within the restrosplenial cortex, a major hub of the DMN. These results validate an imaging technique for the study of activation and connectivity in the lightly sedated awake mouse brain and provide evidence supporting an evolutionary preserved function of the DMN, putatively improving translational relevance of preclinical models of neuropsychiatric diseases.

Separation of hemodynamic signals from GCaMP fluorescence measured with wide-field imaging

Wide-field calcium imaging is often used to measure brain dynamics in behaving mice. With a large field of view and a high sampling rate, wide-field imaging can monitor activity from several distant cortical areas simultaneously, revealing cortical interactions. Interpretation of wide-field images is complicated, however, by the absorption of light by hemoglobin, which can substantially affect the measured fluorescence. One approach to separating hemodynamics and calcium signals is to use multiwavelength backscatter recordings to measure light absorption by hemoglobin. Following this approach, we develop a spatially detailed regression-based method to estimate hemodynamics. This Spatial Model is based on a linear form of the Beer–Lambert relationship but is fit at every pixel in the image and does not rely on the estimation of physical parameters. In awake mice of three transgenic lines, the Spatial Model offers improved separation of hemodynamics and changes in GCaMP fluorescence. The improvement is pronounced near blood vessels and, in contrast with the Beer–Lambert equations, can remove vascular artifacts along the sagittal midline and in general permits more accurate fluorescence-based determination of neuronal activity across the cortex.

NEW & NOTEWORTHY This paper addresses a well-known and strong source of contamination in wide-field calcium-imaging data: hemodynamics. To guide researchers toward the best method to separate calcium signals from hemodynamics, we compare the performance of several methods in three commonly used mouse lines and present a novel regression model that outperforms the other techniques we consider.

Cerebral oxygenation during locomotion is modulated by respiration

In the brain, increased neural activity is correlated with increases of cerebral blood flow and tissue oxygenation. However, how cerebral oxygen dynamics are controlled in the behaving animal remains unclear. We investigated to what extent cerebral oxygenation varies during locomotion. We measured oxygen levels in the cortex of awake, head-fixed mice during locomotion using polarography, spectroscopy, and two-photon phosphorescence lifetime measurements of oxygen sensors. We find that locomotion significantly and globally increases cerebral oxygenation, specifically in areas involved in locomotion, as well as in the frontal cortex and the olfactory bulb. The oxygenation increase persists when neural activity and functional hyperemia are blocked, occurred both in the tissue and in arteries feeding the brain, and is tightly correlated with respiration rate and the phase of respiration cycle. Thus, breathing rate is a key modulator of cerebral oxygenation and should be monitored during hemodynamic imaging, such as in BOLD fMRI.

Vascular and neural basis of the BOLD signal

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

Anatomical basis and physiological role of cerebrospinal fluid transport through the murine cribriform plate

Cerebrospinal fluid (CSF) flows through the brain, transporting chemical signals and removing waste. CSF production in the brain is balanced by a constant outflow of CSF, the anatomical basis of which is poorly understood. Here, we characterized the anatomy and physiological function of the CSF outflow pathway along the olfactory sensory nerves through the cribriform plate, and into the nasal epithelia. Chemical ablation of olfactory sensory nerves greatly reduced outflow of CSF through the cribriform plate. The reduction in CSF outflow did not cause an increase in intracranial pressure (ICP), consistent with an alteration in the pattern of CSF drainage or production. Our results suggest that damage to olfactory sensory neurons (such as from air pollution) could contribute to altered CSF turnover and flow, providing a potential mechanism for neurological diseases.

Twitches, Blinks, and Fidgets: Important Generators of Ongoing Neural Activity

Animals and humans continuously engage in small, spontaneous motor actions, such as blinking, whisking, and postural adjustments ("fidgeting"). These movements are accompanied by changes in neural activity in sensory and motor regions of the brain. The frequency of these motions varies in time, is affected by sensory stimuli, arousal levels, and pathology. These fidgeting behaviors can be entrained by sensory stimuli. Fidgeting behaviors will cause distributed, bilateral functional activation in the 0.01 to 0.1 Hz frequency range that will show up in functional magnetic resonance imaging and wide-field calcium neuroimaging studies, and will contribute to the observed functional connectivity among brain regions. However, despite the large potential of these behaviors to drive brain-wide activity, these fidget-like behaviors are rarely monitored. We argue that studies of spontaneous and evoked brain dynamics in awake animals and humans should closely monitor these fidgeting behaviors. Differences in these fidgeting behaviors due to arousal or pathology will "contaminate" ongoing neural activity, and lead to apparent differences in functional connectivity. Monitoring and accounting for the brain-wide activations by these behaviors is essential during experiments to differentiate fidget-driven activity from internally driven neural dynamics.

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.

The pial vasculature of the mouse develops according to a sensory-independent program

The cerebral vasculature is organized to supply the brain’s metabolic needs. Sensory deprivation during the early postnatal period causes altered neural activity and lower metabolic demand. Neural activity is instructional for some aspects of vascular development, and deprivation causes changes in capillary density in the deprived brain region. However, it is not known if the pial arteriole network, which contains many leptomeningeal anastomoses (LMAs) that endow the network with redundancy against occlusions, is also affected by sensory deprivation. We quantified the effects of early-life sensory deprivation via whisker plucking on the densities of LMAs and penetrating arterioles (PAs) in anatomically-identified primary sensory regions (vibrissae cortex, forelimb/hindlimb cortex, visual cortex and auditory cortex) in mice. We found that the densities of penetrating arterioles were the same across cortical regions, though the hindlimb representation had a higher density of LMAs than other sensory regions. We found that the densities of PAs and LMAs, as well as quantitative measures of network topology, were not affected by sensory deprivation. Our results show that the postnatal development of the pial arterial network is robust to sensory deprivation.

Vision and Locomotion Shape the Interactions between Neuron Types in Mouse Visual Cortex

Cortical computation arises from the interaction of multiple neuronal types, including pyramidal (Pyr) cells and interneurons expressing Sst, Vip, or Pvalb. To study the circuit underlying such interactions, we imaged these four types of cells in mouse primary visual cortex (V1). Our recordings in darkness were consistent with a "disinhibitory" model in which locomotion activates Vip cells, thus inhibiting Sst cells and disinhibiting Pyr cells. However, the disinhibitory model failed when visual stimuli were present: locomotion increased Sst cell responses to large stimuli and Vip cell responses to small stimuli. A recurrent network model successfully predicted each cell type's activity from the measured activity of other types. Capturing the effects of locomotion, however, required allowing it to increase feedforward synaptic weights and modulate recurrent weights. This network model summarizes interneuron interactions and suggests that locomotion may alter cortical computation by changing effective synaptic connectivity.

Weak correlations between hemodynamic signals and ongoing neural activity during the resting state

Spontaneous fluctuations in hemodynamic signals in the absence of a task or overt stimulation are used to infer neural activity. We tested this coupling by simultaneously measuring neural activity and changes in cerebral blood volume (CBV) in the somatosensory cortex of awake, head-fixed mice during periods of true rest, and during whisker stimulation and volitional whisking. Here we show that neurovascular coupling was similar across states, and large spontaneous CBV changes in the absence of sensory input were driven by volitional whisker and body movements. Hemodynamic signals during periods of rest were weakly correlated with neural activity. Spontaneous fluctuations in CBV and vessel diameter persisted when local neural spiking and glutamatergic input was blocked, and during blockade of noradrenergic receptors, suggesting a non-neuronal origin for spontaneous CBV fluctuations. Spontaneous hemodynamic signals reflect a combination of behavior, local neural activity, and putatively non-neural processes.

Chronic imaging through "transparent skull" in mice

Growing interest in long-term visualization of cortical structure and function requires methods that allow observation of an intact cortex in longitudinal imaging studies. Here we describe a detailed protocol for the "transparent skull" (TS) preparation based on skull clearing with cyanoacrylate, which is applicable for long-term imaging through the intact skull in mice. We characterized the properties of the TS in imaging of intrinsic optical signals and compared them with the more conventional cranial window preparation. Our results show that TS is less invasive, maintains stabile transparency for at least two months, and compares favorably to data obtained from the conventional cranial window. We applied this method to experiments showing that a four-week treatment with the antidepressant fluoxetine combined with one week of monocular deprivation induced a shift in ocular dominance in the mouse visual cortex, confirming that fluoxetine treatment restores critical-period-like plasticity. Our results demonstrate that the TS preparation could become a useful method for long-term visualization of the living mouse brain.

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

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

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.

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.

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.

A method for longitudinal, transcranial imaging of blood flow and remodeling of the cerebral vasculature in postnatal mice

In the weeks following birth, both the brain and the vascular network that supplies it undergo dramatic alteration. While studies of the postnatal evolution of the pial vasculature and blood flow through its vessels have been previously done histologically or acutely, here we describe a neonatal reinforced thin‐skull preparation for longitudinally imaging the development of the pial vasculature in mice using two‐photon laser scanning microscopy. Starting with mice as young as postnatal day 2 (P2), we are able to chronically image cortical areas >1 mm², repeatedly for several consecutive days, allowing us to observe the remodeling of the pial arterial and venous networks. We used this method to measure blood velocity in individual vessels over multiple days, and show that blood flow through individual pial venules was correlated with subsequent diameter changes. This preparation allows the longitudinal imaging of the developing mammalian cerebral vascular network and its physiology.

We developed a technique to longitudinally image blood vessels in the neonatal mouse cortex transcranially using two‐photon microscopy. The blood vessels on the surface of the brain undergo substantial pruning after birth. Blood flow through a vessel was correlated with the subsequent diameter change of the vessel.