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

Increased wiring cost during development is driven by long-range cortical, but not subcortical connections

The brain undergoes a protracted, metabolically expensive maturation process from childhood to adulthood. Therefore, it is crucial to understand how network cost is distributed among different brain systems as the brain matures. To address this issue, here we examined developmental changes in wiring cost and brain network topology using resting-state functional magnetic resonance imaging (rsfMRI) data longitudinally collected in awake rats from the juvenile age to adulthood. We found that the wiring cost increased in the vast majority of cortical connections but decreased in most subcortico-subcortical connections. Importantly, the developmental increase in wiring cost was dominantly driven by long-range cortical, but not subcortical connections, which was consistent with more pronounced increase in network integration in the cortical network. These results collectively indicate that there is a non-uniform distribution of network cost as the brain matures, and network resource is dominantly consumed for the development of the cortex, but not subcortex from the juvenile age to adulthood.

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.

Nitric oxide is fundamental to neurovascular coupling in humans

KEY POINTS

Preclinical models have demonstrated that nitric oxide is a key component of neurovascular coupling; this has yet to be translated to humans. We conducted two separate protocols utilizing intravenous infusion of a nitric oxide synthase inhibitor and isovolumic haemodilution to assess the influence of nitric oxide on neurovascular coupling in humans. Isovolumic haemodilution did not alter neurovascular coupling. Intravenous infusion of a nitric oxide synthase inhibitor reduced the neurovascular coupling response by ∼30%, indicating that nitric oxide is integral to neurovascular coupling in humans.

ABSTRACT

Nitric oxide is a vital neurovascular signalling molecule in preclinical models, yet the mechanisms underlying neurovascular coupling (NVC) in humans have yet to be elucidated. To investigate the contribution of nitric oxide to NVC in humans, we utilized a visual stimulus paradigm to elicit an NVC response in the posterior cerebral circulation. Two distinct mechanistic interventions were conducted on young healthy males: (1) NVC was assessed during intravenous infusion of saline (placebo) and the non-selective competitive nitric oxide synthase inhibitor N^G -monomethyl-l-arginine (l-NMMA, 5 mg kg^-1 bolus & subsequent 50 μg kg^-1 min^-1 maintenance dose; n = 10). The order of infusion was randomized, counterbalanced and single blinded. A subset of participants in this study (n = 4) underwent a separate intervention with phenylephrine infusion to independently consider the influence of blood pressure changes on NVC (0.1-0.6 μg kg^-1 min^-1 constant infusion). (2) NVC was assessed prior to and following isovolumic haemodilution, whereby 20% of whole blood was removed and replaced with 5% human serum albumin to reduce haemoglobin concentration (n = 8). For both protocols, arterial and internal jugular venous blood samples were collected at rest and coupled with volumetric measures of cerebral blood flow (duplex ultrasound) to quantify resting cerebral metabolic parameters. l-NMMA elicited a 30% reduction in the peak (P = 0.01), but not average (P = 0.11), NVC response. Neither phenylephrine nor haemodilution influenced NVC. Nitric oxide signalling is integral to NVC in humans, providing a new direction for research into pharmacological treatment of humans with dementia.

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.

Modeling of vascular space occupancy and BOLD functional MRI from first principles using real microvascular angiograms

PURPOSE

The vascular space occupancy (VASO) is a functional MRI technique for probing cerebral blood volume changes noninvasively, including during neuronal activation in humans. An important consideration when implementing VASO is the BOLD effect in the signal. Assessing the physical origin of this BOLD contamination and the capabilities of correction methods could improve the quantification of cerebral blood volume changes with VASO.

METHODS

Given the heterogeneity of cerebral microvascular architecture, the vascular geometry within an MRI voxel can influence both BOLD and VASO signals. To investigate this effect, 3D high-resolution images of mouse cerebral vasculature measured with two-photon microscopy were used to model BOLD and VASO signals from first principles using Monte Carlo diffusion of water protons. Quantitative plots of VASO together with intravascular and extravascular BOLD signals as a function of TE at B₀ fields 1.5 T to 14 T were obtained.

RESULTS

The BOLD contamination of the VASO response was on the order of 50% for gradient echo and 5% for spin echo at 7 T and TE = 6 ms and significantly increased with TE and B₀ . Two currently used correction schemes were shown to account for most of this contamination and recover accurate relative signal changes, with optimal correction obtained using TEs as short as possible.

CONCLUSION

These results may provide useful information for optimizing sequence parameters in VASO and BOLD functional MRI, leading the way to a wider application of these techniques in healthy and diseased brain.

A quantitative analysis of cell-specific contributions and the role of anesthetics to the neurovascular coupling

The neurovascular coupling (NVC) connects neuronal activity to hemodynamic responses in the brain. This connection is the basis for the interpretation of functional magnetic resonance imaging data. Despite the central role of this coupling, we lack detailed knowledge about cell-specific contributions and our knowledge about NVC is mainly based on animal experiments performed during anesthesia. Anesthetics are known to affect neuronal excitability, but how this affects the vessel diameters is not known. Due to the high complexity of NVC data, mathematical modeling is needed for a meaningful analysis. However, neither the relevant neuronal subtypes nor the effects of anesthetics are covered by current models. Here, we present a mathematical model including GABAergic interneurons and pyramidal neurons, as well as the effect of an anesthetic agent. The model is consistent with data from optogenetic experiments from both awake and anesthetized animals, and it correctly predicts data from experiments with different pharmacological modulators. The analysis suggests that no downstream anesthetic effects are necessary if one of the GABAergic interneuron signaling pathways include a Michaelis-Menten expression. This is the first example of a quantitative model that includes both the cell-specific contributions and the effect of an anesthetic agent on the NVC.

EEG spectro-temporal amplitude modulation as a measurement of cortical hemodynamics: an EEG-fNIRS study

Neurovascular coupling provides valuable descriptive information about neural function and communication. In this work, we propose to objectively characterize EEG sub-band modulation in an attempt to compare with local variations of fNIRS hemoglobin concentration. First, full-band EEG signals are decomposed into five well-known frequency sub-bands: delta, theta, alpha, beta, and gamma. The temporal amplitude envelope of each sub-band is then computed via Hilbert transformation. The proposed EEG 'spectro-temporal amplitude modulation' (EEG-AM) feature measures the rate at which each sub-band is modulated. Similarities between EEG-AM features and fNIRS hemoglobin concentration are computed for four neighboring channels over the occipital area during resting-state. Experiments with a database of 29 participants show statistically significant similarities between the total hemoglobin concentration and the alpha band modulating the alpha, beta, and gamma frequencies. These results support the idea that the EEG-AM can carry hemodynamic properties.Clinical relevance- This shows that the EEG spectro-temporal amplitude modulation present similarities with the hemoglobin concentration in co-placed channels.

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.

Fiber photometry for monitoring cerebral oxygen saturation in freely-moving rodents

Hemodynamic parameters, such as tissue oxygen saturation and blood volume fraction, are important markers of brain physiology. They are also widely used surrogate markers of electrophysiological activity. Here, we present a single fiber spectroscopic (SFS) system for monitoring cerebral oxygen saturation in localized, non-line-of-sight brain regions in freely-moving rodents. We adapted the implantation ferrule and patch cable design from commercialized optogenetics and fiber photometry systems, enabling stereotaxic fiber implantation, longitudinal tissue access and measurement from freely-moving animals. The optical system delivers and collects light from the brain through a 200 µm-core-diameter, 0.39NA multimode fiber. We robustly measured oxygen saturation from phantoms with different optical properties mimicking brain tissue. In mice, we demonstrated, for the first time, measurements of oxygen saturation from a highly-localized, targeted brain region over 31 days and continuous measurements from a freely-moving animal for over an hour. These results suggest that single fiber spectroscopy has enormous potential for functional brain monitoring and investigating neurovascular coupling in freely-moving animals. In addition, this technique can potentially be combined with fiber photometry systems to correct for hemodynamic artifacts in the fluorescence detection.

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

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.

An open database of resting-state fMRI in awake rats

Rodent models are essential to translational research in health and disease. Investigation in rodent brain function and organization at the systems level using resting-state functional magnetic resonance imaging (rsfMRI) has become increasingly popular. Due to this rapid progress, publicly shared rodent rsfMRI databases can be of particular interest and importance to the scientific community, as inspired by human neuroscience and psychiatric research that are substantially facilitated by open human neuroimaging datasets. However, such databases in rats are still rare. In this paper, we share an open rsfMRI database acquired in 90 rats with a well-established awake imaging paradigm that avoids anesthesia interference. Both raw and preprocessed data are made publicly available. Procedures in data preprocessing to remove artefacts induced by the scanner, head motion and non-neural physiological noise are described in details. We also showcase inter-regional functional connectivity and functional networks obtained from the database.

The Pharmacokinetics of Medetomidine Administered Subcutaneously during Isoflurane Anaesthesia in Sprague-Dawley Rats

Anaesthetic protocols involving the combined use of a sedative agent, medetomidine, and an anaesthetic agent, isoflurane, are increasingly being used in functional magnetic resonance imaging (fMRI) studies of the rodent brain. Despite the popularity of this combination, a standardised protocol for the combined use of medetomidine and isoflurane has not been established, resulting in inconsistencies in the reported use of these drugs. This study investigated the pharmacokinetic detail required to standardise the use of medetomidine and isoflurane in rat brain fMRI studies. Using mass spectrometry, serum concentrations of medetomidine were determined in Sprague-Dawley rats during medetomidine and isoflurane anaesthesia. The serum concentration of medetomidine for administration with 0.5% (vapouriser setting) isoflurane was found to be 14.4 ng/mL (±3.0 ng/mL). The data suggests that a steady state serum concentration of medetomidine when administered with 0.5% (vapouriser setting) isoflurane can be achieved with an initial subcutaneous (SC) dose of 0.12 mg/kg of medetomidine followed by a 0.08 mg/kg/h SC infusion of medetomidine. Consideration of these results for future studies will facilitate standardisation of medetomidine and isoflurane anaesthetic protocols during fMRI data acquisition.

Investigation of robust visual reaction and functional connectivity in the rat brain induced by rocuronium bromide with functional MRI

Functional magnetic resonance imaging (fMRI) has been used extensively to understand the brain function of a wide range of neurological and psychiatric disorders. When applied to animal studies, anesthesia is always used to reduce the movement of the animal and also reduce the impacts on the results of fMRI. Several awake models have been proposed by applying physical animal movement restrictions. However, restraining devices were designed for individual subject which limits the promotion of fMRI in awake animals. Here, a clinical muscle relaxant rocuronium bromide (RB) was introduced to restrain the animal in fMRI scanning time. The fMRI reactions of the animal induced with RB and the other two commonly used anesthesia protocols were investigated. The results of the fMRI showed that there were increased functional connectivity and well-round visual responses in the RB induced state. Furthermore, significant BOLD signal changes were found in the cortex and thalamus regions when the animal revived from isoflurane, which should be essential to further understand the effects of anesthesia on the brain.

Enhanced Cerebral Blood Volume under Normobaric Hyperoxia in the J20-hAPP Mouse Model of Alzheimer's Disease

Early impairments to neurovascular coupling have been proposed to be a key pathogenic factor in the onset and progression of Alzheimer's disease (AD). Studies have shown impaired neurovascular function in several mouse models of AD, including the J20-hAPP mouse. In this study, we aimed to investigate early neurovascular changes using wild-type (WT) controls and J20-hAPP mice at 6 months of age, by measuring cerebral haemodynamics and neural activity to physiological sensory stimulations. A thinned cranial window was prepared to allow access to cortical vasculature and imaged using 2D-optical imaging spectroscopy (2D-OIS). After chronic imaging sessions where the skull was intact, a terminal acute imaging session was performed where an electrode was inserted into the brain to record simultaneous neural activity. We found that cerebral haemodynamic changes were significantly enhanced in J20-hAPP mice compared with controls in response to physiological stimulations, potentially due to the significantly higher neural activity (hyperexcitability) seen in the J20-hAPP mice. Thus, neurovascular coupling remained preserved under a chronic imaging preparation. Further, under hyperoxia, the baseline blood volume and saturation of all vascular compartments in the brains of J20-hAPP mice were substantially enhanced compared to WT controls, but this effect disappeared under normoxic conditions. This study highlights novel findings not previously seen in the J20-hAPP mouse model, and may point towards a potential therapeutic strategy.

Key Aspects of Neurovascular Control Mediated by Specific Populations of Inhibitory Cortical Interneurons

Inhibitory interneurons can evoke vasodilation and vasoconstriction, making them potential cellular drivers of neurovascular coupling. However, the specific regulatory roles played by particular interneuron subpopulations remain unclear. Our purpose was therefore to adopt a cell-specific optogenetic approach to investigate how somatostatin (SST) and neuronal nitric oxide synthase (nNOS)-expressing interneurons might influence the neurovascular relationship. In mice, specific activation of SST- or nNOS-interneurons was sufficient to evoke hemodynamic changes. In the case of nNOS-interneurons, robust hemodynamic changes occurred with minimal changes in neural activity, suggesting that the ability of blood oxygen level dependent functional magnetic resonance imaging (BOLD fMRI) to reliably reflect changes in neuronal activity may be dependent on type of neuron recruited. Conversely, activation of SST-interneurons produced robust changes in evoked neural activity with shallow cortical excitation and pronounced deep layer cortical inhibition. Prolonged activation of SST-interneurons often resulted in an increase in blood volume in the centrally activated area with an accompanying decrease in blood volume in the surrounding brain regions, analogous to the negative BOLD signal. These results demonstrate the role of specific populations of cortical interneurons in the active control of neurovascular function.

Assessment of anesthesia on physiological stability and BOLD signal reliability during visual or acoustic stimulation in the cat

BACKGROUND

Neuroimaging methods including fMRI provide powerful tools to observe whole-brain functional networks. This is particularly powerful in animal models, allowing these networks to be probed using complementary methods. However, most animals must be anesthetized for neuroimaging, giving rise to complications resulting from anesthetic effects on the animal's physiological and neurological functions. For example, an established protocol for feline neuroimaging involves co-administration of ketamine and isoflurane - the latter of which is known to suppress cortical function.

NEW METHOD

Here, we compare this established protocol to alfaxalone, a single-agent anesthetic for functional neuroimaging. We first compare the two in a controlled environment to assess relative safety and to measure physiological stability over an extended time window. We then compare patterns of auditory and visually-evoked activity measured at 7  T to assess mean signal strength and between-subjects signal variability.

RESULTS IN COMPARISON WITH EXISTING METHODS

We show that alfaxalone results in more stable respiratory rates over the 120 min testing period, with evidence of smaller between-measurements variability within this time window, when compared to ketamine plus isoflurane. Moreover, we demonstrate that both agents evoke similar mean BOLD signals across animals, but that alfaxalone elicits more consistent BOLD activity in response to sound stimuli across all ROIs observed.

CONCLUSIONS

Alfaxalone is observed to be more physiologically stable, evoking a more consistent BOLD signal across animals than the co-administration of ketamine and isoflurane. Thus, an alfaxalone-based protocol may represent a better approach for neuroimaging in animal models requiring anesthesia.

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.

Global brain signal in awake rats

Although often used as a nuisance in resting-state functional magnetic resonance imaging (rsfMRI), the global brain signal in humans and anesthetized animals has important neural basis. However, our knowledge of the global signal in awake rodents is sparse. To bridge this gap, we systematically analyzed rsfMRI data acquired with a conventional single-echo (SE) echo planar imaging (EPI) sequence in awake rats. The spatial pattern of rsfMRI frames during peaks of the global signal exhibited prominent co-activations in the thalamo-cortical and hippocampo-cortical networks, as well as in the basal forebrain, hinting that these neural networks might contribute to the global brain signal in awake rodents. To validate this concept, we acquired rsfMRI data using a multi-echo (ME) EPI sequence and removed non-neural components in the rsfMRI signal. Consistent co-activation patterns were obtained in extensively de-noised ME-rsfMRI data, corroborating the finding from SE-rsfMRI data. Furthermore, during rsfMRI experiments, we simultaneously recorded neural spiking activities in the hippocampus using GCaMP-based fiber photometry. The hippocampal calcium activity exhibited significant correspondence with the global rsfMRI signal. These data collectively suggest that the global rsfMRI signal contains significant neural components that involve coordinated activities in the thalamo-cortical and hippocampo-cortical networks. These results provide important insight into the neural substrate of the global brain signal in awake rodents.

Temporal stability of fMRI in medetomidine-anesthetized rats

Medetomidine has become a popular choice for anesthetizing rats during long-lasting sessions of blood-oxygen-level dependent (BOLD) functional magnetic resonance imaging (fMRI). Despite this, it has not yet been thoroughly established how commonly reported fMRI readouts evolve over several hours of medetomidine anesthesia and how they are affected by the precise timing, dose, and route of administration. We used four different protocols of medetomidine administration to anesthetize rats for up to six hours and repeatedly evaluated somatosensory stimulus-evoked BOLD responses and resting state functional connectivity. We found that the temporal evolution of fMRI readouts strongly depended on the method of administration. Intravenous administration of a medetomidine bolus (0.05 mg/kg), combined with a subsequent continuous infusion (0.1 mg/kg/h), led to temporally stable measures of stimulus-evoked activity and functional connectivity throughout the anesthesia. Deviating from the above protocol-by omitting the bolus, lowering the medetomidine dose, or using the subcutaneous route-compromised the stability of these measures in the initial two-hour period. We conclude that both an appropriate protocol of medetomidine administration and a suitable timing of fMRI experiments are crucial for obtaining consistent results. These factors should be considered for the design and interpretation of future rat fMRI studies.

Neurovascular Coupling under Chronic Stress Is Modified by Altered GABAergic Interneuron Activity

Neurovascular coupling (NVC), the interaction between neural activity and vascular response, ensures normal brain function by maintaining brain homeostasis. We previously reported altered cerebrovascular responses during functional hyperemia in chronically stressed animals. However, the underlying neuronal-level changes associated with those hemodynamic changes remained unclear. Here, using in vivo and ex vivo experiments, we investigate the neuronal origins of altered NVC dynamics under chronic stress conditions in adult male mice. Stimulus-evoked hemodynamic and neural responses, especially beta and gamma-band local field potential activity, were significantly lower in chronically stressed animals, and the NVC relationship, itself, had changed. Further, using acute brain slices, we discovered that the underlying cause of this change was dysfunction of neuronal nitric oxide synthase (nNOS)-mediated vascular responses. Using FISH to check the mRNA expression of several GABAergic subtypes, we confirmed that only nNOS mRNA was significantly decreased in chronically stressed mice. Ultimately, chronic stress impairs NVC by diminishing nNOS-mediated vasodilation responses to local neural activity. Overall, these findings provide useful information in understanding NVC dynamics in the healthy brain. More importantly, this study reveals that impaired nNOS-mediated NVC function may be a contributory factor in the progression of stress-related diseases.

SIGNIFICANCE STATEMENT The correlation between neuronal activity and cerebral vascular dynamics is defined as neurovascular coupling (NVC), which plays an important role for meeting the metabolic demands of the brain. However, the impact of chronic stress, which is a contributory factor of many cerebrovascular diseases, on NVC is poorly understood. We therefore investigated the effects of chronic stress on impaired neurovascular response to sensory stimulation and their underlying mechanisms. Multimodal approaches, from in vivo hemodynamic imaging and electrophysiology to ex vivo vascular imaging with pharmacological treatment, patch-clamp recording, FISH, and immunohistochemistry revealed that chronic stress-induced dysfunction of nNOS-expressing interneurons contributes to NVC impairment. These findings will provide useful information to understand the role of nNOS interneurons in NVC in normal and pathological conditions.

Propagations of spontaneous brain activity in awake rats

Slow propagations of spontaneous brain activity have been reported in multiple species. However, systematical investigation of the organization of such brain activity is still lacking. In this study, we analyzed propagations of spontaneous brain activity using a reference library of characteristic resting-state functional connectivity (RSFC) patterns in awake rodents. We found that transitions through multiple distinct RSFC patterns were reproducible not only in transition sequences but also in transition time delays. In addition, the organization of these transitions and their spatiotemporal dynamic patterns were revealed using a graphical model. We further identified prominent brain regions involved in these transitions. These results provide a comprehensive framework of brainwide propagations of spontaneous activity in awake rats. This study also offers a new tool to study the spatiotemporal dynamics of activity in the resting brain.

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.

Differential coupling between subcortical calcium and BOLD signals during evoked and resting state through simultaneous calcium fiber photometry and fMRI

Task based and resting state fMRI has been widely utilized to study brain functions. As the foundation of fMRI, the underlying neural basis of the BOLD signal has been extensively studied, but the detailed mechanism remains elusive, particularly during the resting state. To examine the neurovascular coupling, it is important to simultaneously record neural and vascular signals. Here we developed a novel setup of camera based, scalable simultaneous calcium fiber photometry and fMRI in rats. Using this setup, we recorded calcium signals of superior colliculus (SC) and lateral geniculate nucleus (LGN) and fMRI simultaneously during visual stimulation and the resting state. Our results revealed robust, region-specific coupling between calcium and BOLD signals in the task state and weaker, whole brain correlation in the resting state. Interestingly, the spatial specificity of such correlation in the resting state was improved upon regression of white matter, ventricle signals and global signals in fMRI data. Overall, our results suggest differential coupling of calcium and BOLD signals for subcortical regions between evoked and resting states, and the coupling relationship in the resting state was related with resting state BOLD preprocessing strategies.

Individual variability in behavior and functional networks predicts vulnerability using an animal model of PTSD

Only a minority of individuals experiencing trauma subsequently develop post-traumatic stress disorder (PTSD). However, whether differences in vulnerability to PTSD result from a predisposition or trauma exposure remains unclear. A major challenge in differentiating these possibilities is that clinical studies focus on individuals already exposed to trauma without pre-trauma conditions. Here, using the predator scent model of PTSD in rats and a longitudinal design, we measure pre-trauma brain-wide neural circuit functional connectivity, behavioral and corticosterone responses to trauma exposure, and post-trauma anxiety. Freezing during predator scent exposure correlates with functional connectivity in a set of neural circuits, indicating pre-existing circuit function can predispose animals to differential fearful responses to threats. Counterintuitively, rats with lower freezing show more avoidance of the predator scent, a prolonged corticosterone response, and higher anxiety long after exposure. This study provides a framework of pre-existing circuit function that determines threat responses, which might directly relate to PTSD-like behaviors.

A Brief History of Simulation Neuroscience

Our knowledge of the brain has evolved over millennia in philosophical, experimental and theoretical phases. We suggest that the next phase is simulation neuroscience. The main drivers of simulation neuroscience are big data generated at multiple levels of brain organization and the need to integrate these data to trace the causal chain of interactions within and across all these levels. Simulation neuroscience is currently the only methodology for systematically approaching the multiscale brain. In this review, we attempt to reconstruct the deep historical paths leading to simulation neuroscience, from the first observations of the nerve cell to modern efforts to digitally reconstruct and simulate the brain. Neuroscience began with the identification of the neuron as the fundamental unit of brain structure and function and has evolved towards understanding the role of each cell type in the brain, how brain cells are connected to each other, and how the seemingly infinite networks they form give rise to the vast diversity of brain functions. Neuronal mapping is evolving from subjective descriptions of cell types towards objective classes, subclasses and types. Connectivity mapping is evolving from loose topographic maps between brain regions towards dense anatomical and physiological maps of connections between individual genetically distinct neurons. Functional mapping is evolving from psychological and behavioral stereotypes towards a map of behaviors emerging from structural and functional connectomes. We show how industrialization of neuroscience and the resulting large disconnected datasets are generating demand for integrative neuroscience, how the scale of neuronal and connectivity maps is driving digital atlasing and digital reconstruction to piece together the multiple levels of brain organization, and how the complexity of the interactions between molecules, neurons, microcircuits and brain regions is driving brain simulation to understand the interactions in the multiscale brain.

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.

Gastric stimulation drives fast BOLD responses of neural origin

Functional magnetic resonance imaging (fMRI) is commonly thought to be too slow to capture any neural dynamics faster than 0.1 Hz. However, recent findings demonstrate the feasibility of detecting fMRI activity at higher frequencies beyond 0.2 Hz. The origin, reliability, and generalizability of fast fMRI responses are still under debate and await confirmation through animal experiments with fMRI and invasive electrophysiology. Here, we acquired single-echo and multi-echo fMRI, as well as local field potentials, from anesthetized rat brains given gastric electrical stimulation modulated at 0.2, 0.4 and 0.8 Hz. Such gastric stimuli could drive widespread fMRI responses at corresponding frequencies from the somatosensory and cingulate cortices. Such fast fMRI responses were linearly dependent on echo times and thus indicative of blood oxygenation level dependent nature (BOLD). Local field potentials recorded during the same gastric stimuli revealed transient and phase-locked broadband neural responses, preceding the fMRI responses by as short as 0.5 s. Taken together, these results suggest that gastric stimulation can drive widespread and rapid fMRI responses of BOLD and neural origin, lending support to the feasibility of using fMRI to detect rapid changes in neural activity up to 0.8 Hz under visceral stimulation.

Standardization of Small Animal Imaging-Current Status and Future Prospects

The benefit of small animal imaging is directly linked to the validity and reliability of the collected data. If the data (regardless of the modality used) are not reproducible and/or reliable, then the outcome of the data is rather questionable. Therefore, standardization of the use of small animal imaging equipment, as well as of animal handling in general, is of paramount importance. In a recent paper, guidance for efficient small animal imaging quality control was offered and discussed, among others, the use of phantoms in setting up a quality control program (Osborne et al. 2016). The same phantoms can be used to standardize image quality parameters for multi-center studies or multi-scanners within center studies. In animal experiments, the additional complexity due to animal handling needs to be addressed to ensure standardized imaging procedures. In this review, we will address the current status of standardization in preclinical imaging, as well as potential benefits from increased levels of standardization.

Mapping stress networks using functional magnetic resonance imaging in awake animals

The neurobiology of stress is studied through behavioral neuroscience, endocrinology, neuronal morphology and neurophysiology. There is a shift in focus toward progressive changes throughout stress paradigms and individual susceptibility to stress that requires methods that allow for longitudinal study design and study of individual differences in stress response. Functional magnetic resonance imaging (fMRI), with the advantages of noninvasiveness and a large field of view, can be used for functionally mapping brain-wide regions and circuits critical to the stress response, making it suitable for longitudinal studies and understanding individual variability of short-term and long-term consequences of stress exposure. In addition, fMRI can be applied to both animals and humans, which is highly valuable in translating findings across species and examining whether the physiology and neural circuits involved in the stress response are conserved in mammals. However, compared to human fMRI studies, there are a number of factors that are essential for the success of fMRI studies in animals. This review discussed the use of fMRI in animal studies of stress. It reviewed advantages, challenges and technical considerations of the animal fMRI methodology as well as recent literature of stress studies using fMRI in animals. It also highlighted the development of combining fMRI with other methods and the future potential of fMRI in animal studies of stress. We conclude that animal fMRI studies, with their flexibility, low cost and short time frame compared to human studies, are crucial to advancing our understanding of the neurobiology of stress.

Awake Rat Brain Functional Magnetic Resonance Imaging Using Standard Radio Frequency Coils and a 3D Printed Restraint Kit

Functional magnetic resonance imaging (fMRI) is a powerful noninvasive tool for studying spontaneous resting state functional connectivity (RSFC) in laboratory animals. Brain function can be significantly affected by generally used anesthetics, however, rendering the need for awake imaging. Only a few different awake animal habituation protocols have been presented, and there is a critical need for practical and improved low-stress techniques. Here we demonstrate a novel restraint approach for awake rat RSFC studies. Our custom-made 3D printed restraint kit is compatible with a standard Bruker Biospin MRI rat bed, rat brain receiver coil, and volume transmitter coil. We also implemented a progressive habituation protocol aiming to minimize the stress experienced by the rats, and compared RSFC between awake, lightly sedated, and isoflurane-anesthetized rats. Our results demonstrated that the 3D printed restraint kit was suitable for RSFC studies of awake rats. During the short 4-day habituation period, the plasma corticosterone concentration, movement, and heart rate, which were measured as stress indicators, decreased significantly, indicating adaptation to the restraint protocol. Additionally, 10 days after the awake MRI session, rats exhibited no signs of depression or anxiety based on open-field and sucrose preference behavioral tests. The RSFC data revealed significant changes in the thalamo-cortical and cortico-cortical networks between the awake, lightly sedated, and anesthetized groups, emphasizing the need for awake imaging. The present work demonstrates the feasibility of our custom-made 3D printed restraint kit. Using this kit, we found that isoflurane markedly affected brain connectivity compared with that in awake rats, and that the effect was less pronounced, but still significant, when light isoflurane sedation was used instead.

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.

Development of brain-wide connectivity architecture in awake rats

Childhood and adolescence are both critical developmental periods, evidenced by complex neurophysiological changes the brain undergoes and high occurrence rates of neuropsychiatric disorders during these periods. Despite substantial progress in elucidating the developmental trajectories of individual neural circuits, our knowledge of developmental changes of whole-brain connectivity architecture in animals is sparse. To fill this gap, here we longitudinally acquired rsfMRI data in awake rats during five developmental stages from juvenile to adulthood. We found that the maturation timelines of brain circuits were heterogeneous and system specific. Functional connectivity (FC) tended to decrease in subcortical circuits, but increase in cortical circuits during development. In addition, the developing brain exhibited hemispheric functional specialization, evidenced by reduced inter-hemispheric FC between homotopic regions, and lower similarity of region-to-region FC patterns between the two hemispheres. Finally, we showed that whole-brain network development was characterized by reduced clustering (i.e. local communication) but increased integration (distant communication). Taken together, the present study has systematically characterized the development of brain-wide connectivity architecture from juvenile to adulthood in awake rats. It also serves as a critical reference point for understanding circuit- and network-level changes in animal models of brain development-related disorders. Furthermore, FC data during brain development in awake rodents contain high translational value and can shed light onto comparative neuroanatomy.

Delineating the Macroscale Areal Organization of the Macaque Cortex In Vivo

Complementing long-standing traditions centered on histology, fMRI approaches are rapidly maturing in delineating brain areal organization at the macroscale. The non-human primate (NHP) provides the opportunity to overcome critical barriers in translational research. Here, we establish the data requirements for achieving reproducible and internally valid parcellations in individuals. We demonstrate that functional boundaries serve as a functional fingerprint of the individual animals and can be achieved under anesthesia or awake conditions (rest, naturalistic viewing), though differences between awake and anesthetized states precluded the detection of individual differences across states. Comparison of awake and anesthetized states suggested a more nuanced picture of changes in connectivity for higher-order association areas, as well as visual and motor cortex. These results establish feasibility and data requirements for the generation of reproducible individual-specific parcellations in NHPs, provide insights into the impact of scan state, and motivate efforts toward harmonizing protocols.

Fiber-optic implant for simultaneous fluorescence-based calcium recordings and BOLD fMRI in mice

Despite the growing popularity of blood oxygen level-dependent (BOLD) functional MRI (fMRI), understanding of its underlying principles is still limited. This protocol describes a technique for simultaneous measurement of neural activity using fluorescent calcium indicators together with the corresponding hemodynamic BOLD fMRI response in the mouse brain. Our early work using small-molecule fluorophores in rats gave encouraging results but was limited to acute measurements using synthetic dyes. Our latest procedure combines fMRI with optical detection of cell-type-specific virally delivered GCaMP6, a genetically encoded calcium indicator (GECI). GCaMP6 fluorescence, which increases upon calcium binding, is collected by a chronically implanted optical fiber, allowing longitudinal studies in mice. The chronic implant, placed horizontally on the skull, has an angulated tip that reflects light into the brain and is connected via fiber optics to a remote optical setup. The technique allows access to the neocortex and does not require adaptations of commercial MRI hardware. The hybrid approach permits fiber-optic calcium recordings with simultaneous artifact-free BOLD fMRI with full brain coverage and 1-s temporal resolution using standard gradient-echo echo-planar imaging (GE-EPI) sequences. The method provides robust, cell-type-specific readouts to link neural activity to BOLD signals, as emonstrated for task-free ('resting-state') conditions and in response to hind-paw stimulation. These results highlight the power of fiber photometry combined with fMRI, which we aim to further advance in this protocol. The approach can be easily adapted to study other molecular processes using suitable fluorescent indicators.

Temporal transitions of spontaneous brain activity

Spontaneous brain activity, typically investigated using resting-state fMRI (rsfMRI), provides a measure of inter-areal resting-state functional connectivity (RSFC). Although it has been established that RSFC is non-stationary, previous dynamic rsfMRI studies mainly focused on revealing the spatial characteristics of dynamic RSFC patterns, but the temporal relationship between these RSFC patterns remains elusive. Here we investigated the temporal organization of characteristic RSFC patterns in awake rats and humans. We found that transitions between RSFC patterns were not random but followed specific sequential orders. The organization of RSFC pattern transitions was further analyzed using graph theory, and pivotal RSFC patterns in transitions were identified. This study has demonstrated that spontaneous brain activity is not only nonrandom spatially, but also nonrandom temporally, and this feature is well conserved between rodents and humans. These results offer new insights into understanding the spatiotemporal dynamics of spontaneous activity in the mammalian brain.

Concussion, microvascular injury, and early tauopathy in young athletes after impact head injury and an impact concussion mouse model

The mechanisms underpinning concussion, traumatic brain injury, and chronic traumatic encephalopathy, and the relationships between these disorders, are poorly understood. We examined post-mortem brains from teenage athletes in the acute-subacute period after mild closed-head impact injury and found astrocytosis, myelinated axonopathy, microvascular injury, perivascular neuroinflammation, and phosphorylated tau protein pathology. To investigate causal mechanisms, we developed a mouse model of lateral closed-head impact injury that uses momentum transfer to induce traumatic head acceleration. Unanaesthetized mice subjected to unilateral impact exhibited abrupt onset, transient course, and rapid resolution of a concussion-like syndrome characterized by altered arousal, contralateral hemiparesis, truncal ataxia, locomotor and balance impairments, and neurobehavioural deficits. Experimental impact injury was associated with axonopathy, blood-brain barrier disruption, astrocytosis, microgliosis (with activation of triggering receptor expressed on myeloid cells, TREM2), monocyte infiltration, and phosphorylated tauopathy in cerebral cortex ipsilateral and subjacent to impact. Phosphorylated tauopathy was detected in ipsilateral axons by 24 h, bilateral axons and soma by 2 weeks, and distant cortex bilaterally at 5.5 months post-injury. Impact pathologies co-localized with serum albumin extravasation in the brain that was diagnostically detectable in living mice by dynamic contrast-enhanced MRI. These pathologies were also accompanied by early, persistent, and bilateral impairment in axonal conduction velocity in the hippocampus and defective long-term potentiation of synaptic neurotransmission in the medial prefrontal cortex, brain regions distant from acute brain injury. Surprisingly, acute neurobehavioural deficits at the time of injury did not correlate with blood-brain barrier disruption, microgliosis, neuroinflammation, phosphorylated tauopathy, or electrophysiological dysfunction. Furthermore, concussion-like deficits were observed after impact injury, but not after blast exposure under experimental conditions matched for head kinematics. Computational modelling showed that impact injury generated focal point loading on the head and seven-fold greater peak shear stress in the brain compared to blast exposure. Moreover, intracerebral shear stress peaked before onset of gross head motion. By comparison, blast induced distributed force loading on the head and diffuse, lower magnitude shear stress in the brain. We conclude that force loading mechanics at the time of injury shape acute neurobehavioural responses, structural brain damage, and neuropathological sequelae triggered by neurotrauma. These results indicate that closed-head impact injuries, independent of concussive signs, can induce traumatic brain injury as well as early pathologies and functional sequelae associated with chronic traumatic encephalopathy. These results also shed light on the origins of concussion and relationship to traumatic brain injury and its aftermath.awx350media15713427811001.

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.

The Neurovascular Unit Coming of Age: A Journey through Neurovascular Coupling in Health and Disease

The concept of the neurovascular unit (NVU), formalized at the 2001 Stroke Progress Review Group meeting of the National Institute of Neurological Disorders and Stroke, emphasizes the intimate relationship between the brain and its vessels. Since then, the NVU has attracted the interest of the neuroscience community, resulting in considerable advances in the field. Here the current state of knowledge of the NVU will be assessed, focusing on one of its most vital roles: the coupling between neural activity and blood flow. The evidence supports a conceptual shift in the mechanisms of neurovascular coupling, from a unidimensional process involving neuronal-astrocytic signaling to local blood vessels to a multidimensional one in which mediators released from multiple cells engage distinct signaling pathways and effector systems across the entire cerebrovascular network in a highly orchestrated manner. The recently appreciated NVU dysfunction in neurodegenerative diseases, although still poorly understood, supports emerging concepts that maintaining neurovascular health promotes brain health.