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

Long-wavelength traveling waves of vasomotion modulate the perfusion of cortex

Brain arterioles are active, multicellular complexes whose diameters oscillate at ∼ 0.1 Hz. We assess the physiological impact and spatiotemporal dynamics of vaso-oscillations in the awake mouse. First, vaso-oscillations in penetrating arterioles, which source blood from pial arterioles to the capillary bed, profoundly impact perfusion throughout neocortex. The modulation in flux during resting-state activity exceeds that of stimulus-induced activity. Second, the change in perfusion through arterioles relative to the change in their diameter is weak. This implies that the capillary bed dominates the hydrodynamic resistance of brain vasculature. Lastly, the phase of vaso-oscillations evolves slowly along arterioles, with a wavelength that exceeds the span of the cortical mantle and sufficient variability to establish functional cortical areas as parcels of uniform phase. The phase-gradient supports traveling waves in either direction along both pial and penetrating arterioles. This implies that waves along penetrating arterioles can mix, but not directionally transport, interstitial fluids.

Cortical spreading depolarization-induced constriction of penetrating arteries can cause watershed ischemia: A potential mechanism for white matter lesions

Periventricular white matter lesions (WMLs) are common MRI findings in migraine with aura (MA). Although hemodynamic disadvantages of vascular supply to this region create vulnerability, the pathophysiological mechanisms causing WMLs are unclear. We hypothesize that prolonged oligemia, a consequence of cortical spreading depolarization (CSD) underlying migraine aura, may lead to ischemia/hypoxia at hemodynamically vulnerable watershed zones fed by long penetrating arteries (PAs). For this, we subjected mice to KCl-triggered single or multiple CSDs. We found that post-CSD oligemia was significantly deeper at medial compared to lateral cortical areas, which induced ischemic/hypoxic changes at watershed areas between the MCA/ACA, PCA/anterior choroidal and at the tip of superficial and deep PAs, as detected by histological and MRI examination of brains 2-4 weeks after CSD. BALB-C mice, in which MCA occlusion causes large infarcts due to deficient collaterals, exhibited more profound CSD-induced oligemia and were more vulnerable compared to Swiss mice such that a single CSD was sufficient to induce ischemic lesions at the tip of PAs. In conclusion, CSD-induced prolonged oligemia has potential to cause ischemic/hypoxic injury at hemodynamically vulnerable brain areas, which may be one of the mechanisms underlying WMLs located at the tip of medullary arteries seen in MA patients.

Infiltration of meningeal macrophages into the Virchow-Robin space after ischemic stroke in rats: Correlation with activated PDGFR-β-positive adventitial fibroblasts

Macrophages play a crucial role in wound healing and fibrosis progression after brain injury. However, a detailed analysis of their initial infiltration and interaction with fibroblasts is yet to be conducted. This study aimed to investigate the possible route for migration of meningeal macrophages into the ischemic brain and whether these macrophages closely interact with neighboring platelet-derived growth factor beta receptor (PDGFR-β)-positive adventitial fibroblasts during this process. A rat model of ischemic stroke induced by middle cerebral artery occlusion (MCAO) was developed. In sham-operated rats, CD206-positive meningeal macrophages were confined to the leptomeninges and the perivascular spaces, and they were not found in the cortical parenchyma. In MCAO rats, the number of CD206-positive meningeal macrophages increased both at the leptomeninges and along the vessels penetrating the cortex 1 day after reperfusion and increased progressively in the extravascular area of the cortical parenchyma by 3 days. Immunoelectron microscopy and correlative light and electron microscopy showed that in the ischemic brain, macrophages were frequently located in the Virchow-Robin space around the penetrating arterioles and ascending venules at the pial surface. This was identified by cells expressing PDGFR-β, a novel biomarker of leptomeningeal cells. Macrophages within penetrating vessels were localized in the perivascular space between smooth muscle cells and PDGFR-β-positive adventitial fibroblasts. In addition, these PDGFR-β-positive fibroblasts showed morphological and molecular characteristics similar to those of leptomeningeal cells: they had large euchromatic nuclei with prominent nucleoli and well-developed rough endoplasmic reticulum; expressed nestin, vimentin, and type I collagen; and were frequently surrounded by collagen fibrils, indicating active collagen synthesis. In conclusion, the perivascular Virchow-Robin space surrounding the penetrating vessels could be an entry route of meningeal macrophages from the subarachnoid space into the ischemic cortical parenchyma, implying that activated PDGFR-β-positive adventitial fibroblasts could be involved in this process.

Neurovascular coupling: motive unknown

In the brain, increases in neural activity drive changes in local blood flow via neurovascular coupling. The common explanation for increased blood flow (known as functional hyperemia) is that it supplies the metabolic needs of active neurons. However, there is a large body of evidence that is inconsistent with this idea. Baseline blood flow is adequate to supply oxygen needs even with elevated neural activity. Neurovascular coupling is irregular, absent, or inverted in many brain regions, behavioral states, and conditions. Increases in respiration can increase brain oxygenation without flow changes. Simulations show that given the architecture of the brain vasculature, areas of low blood flow are inescapable and cannot be removed by functional hyperemia. As discussed in this article, potential alternative functions of neurovascular coupling include supplying oxygen for neuromodulator synthesis, brain temperature regulation, signaling to neurons, stabilizing and optimizing the cerebral vascular structure, accommodating the non-Newtonian nature of blood, and driving the production and circulation of cerebrospinal fluid (CSF).

Arterial vasodilation drives convective fluid flow in the brain: a poroelastic model

The movement of fluid into, through, and out of the brain plays an important role in clearing metabolic waste. However, there is controversy regarding the mechanisms driving fluid movement in the fluid-filled paravascular spaces (PVS), and whether the movement of metabolic waste in the brain extracellular space (ECS) is primarily driven by diffusion or convection. The dilation of penetrating arterioles in the brain in response to increases in neural activity (neurovascular coupling) is an attractive candidate for driving fluid circulation, as it drives deformation of the brain tissue and of the PVS around arteries, resulting in fluid movement. We simulated the effects of vasodilation on fluid movement into and out of the brain ECS using a novel poroelastic model of brain tissue. We found that arteriolar dilations could drive convective flow through the ECS radially outward from the arteriole, and that this flow is sensitive to the dynamics of the dilation. Simulations of sleep-like conditions, with larger vasodilations and increased extracellular volume in the brain showed enhanced movement of fluid from the PVS into the ECS. Our simulations suggest that both sensory-evoked and sleep-related arteriolar dilations can drive convective flow of cerebrospinal fluid not just in the PVS, but also into the ECS through the PVS around arterioles.

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.

Fluid transport in the brain

The brain harbors a unique ability to, figuratively speaking, shift its gears. During wakefulness, the brain is geared fully toward processing information and behaving, while homeostatic functions predominate during sleep. The blood-brain barrier establishes a stable environment that is optimal for neuronal function, yet the barrier imposes a physiological problem; transcapillary filtration that forms extracellular fluid in other organs is reduced to a minimum in brain. Consequently, the brain depends on a special fluid [the cerebrospinal fluid (CSF)] that is flushed into brain along the unique perivascular spaces created by astrocytic vascular endfeet. We describe this pathway, coined the term glymphatic system, based on its dependency on astrocytic vascular endfeet and their adluminal expression of aquaporin-4 water channels facing toward CSF-filled perivascular spaces. Glymphatic clearance of potentially harmful metabolic or protein waste products, such as amyloid-β, is primarily active during sleep, when its physiological drivers, the cardiac cycle, respiration, and slow vasomotion, together efficiently propel CSF inflow along periarterial spaces. The brain's extracellular space contains an abundance of proteoglycans and hyaluronan, which provide a low-resistance hydraulic conduit that rapidly can expand and shrink during the sleep-wake cycle. We describe this unique fluid system of the brain, which meets the brain's requisites to maintain homeostasis similar to peripheral organs, considering the blood-brain-barrier and the paths for formation and egress of the CSF.

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.

Gradual Not Sudden Change: Multiple Sites of Functional Transition Across the Microvascular Bed

In understanding the role of the neurovascular unit as both a biomarker and target for disease interventions, it is vital to appreciate how the function of different components of this unit change along the vascular tree. The cells of the neurovascular unit together perform an array of vital functions, protecting the brain from circulating toxins and infection, while providing nutrients and clearing away waste products. To do so, the brain's microvasculature dilates to direct energy substrates to active neurons, regulates access to circulating immune cells, and promotes angiogenesis in response to decreased blood supply, as well as pulsating to help clear waste products and maintain the oxygen supply. Different parts of the cerebrovascular tree contribute differently to various aspects of these functions, and previously, it has been assumed that there are discrete types of vessel along the vascular network that mediate different functions. Another option, however, is that the multiple transitions in function that occur across the vascular network do so at many locations, such that vascular function changes gradually, rather than in sharp steps between clearly distinct vessel types. Here, by reference to new data as well as by reviewing historical and recent literature, we argue that this latter scenario is likely the case and that vascular function gradually changes across the network without clear transition points between arteriole, precapillary arteriole and capillary. This is because classically localized functions are in fact performed by wide swathes of the vasculature, and different functional markers start and stop being expressed at different points along the vascular tree. Furthermore, vascular branch points show alterations in their mural cell morphology that suggest functional specializations irrespective of their position within the network. Together this work emphasizes the need for studies to consider where transitions of different functions occur, and the importance of defining these locations, in order to better understand the vascular network and how to target it to treat disease.

Cortical depth-dependent modeling of visual hemodynamic responses

A physiologically based three-dimensional (3D) hemodynamic model is developed to predict the experimentally observed blood oxygen level dependent (BOLD) responses versus the cortical depth induced by visual stimuli. Prior 2D approximations are relaxed in order to analyze 3D blood flow dynamics as a function of cortical depth. Comparison of the predictions with experimental data for evoked stimuli demonstrates that the full 3D model performs at least as well as previous approaches while remaining parsimonious. In particular, the 3D model requires significantly fewer assumptions and model parameters than previous models such that there is no longer need to define depth-specific parameter values for spatial spreading, peak amplitude, and hemodynamic velocity.

Sexual dimorphism in the cerebrovascular network: Brain MRI shows lower arterial density in women

BACKGROUND AND PURPOSE

Accumulating evidence suggests that there is a sexual dimorphism in brain health, with women exhibiting greater disability following strokes of comparable size and having a higher prevalence of cognitive impairment later in life. Despite the critical implication of the cerebrovascular architecture in brain perfusion and brain health, it remains unclear whether structural differences in vessel density exist across the sexes.

METHODS

In this study, we used high-density MRI imaging to characterize the intracerebral arterial and venous density of 28 (14 women) sex-matched healthy young volunteers in vivo. Using an in-house vessel segmentation algorithm, we quantified and compared these vascular features across the cortical and subcortical deep gray matter, white matter, and periventricular white matter.

RESULTS

We found that, on average, women have reduced intracerebral arterial density in comparison to men (F 2.34 ± 0.48%, M 2.67 ± 0.39%; p<.05). This difference was most pronounced in the subcortical deep gray matter (F 1.78 ± 0.53%, M 2.38 ± 0.82%; p<.05) and periventricular white matter (F 0.68 ± 0.15%, M 1.14 ± 0.33%; p<.0005), indicating a potential sex-specific vulnerability to hypoperfusion in areas critical to core cerebral functions. In contrast, venous density did not exhibit a significant difference between sexes.

CONCLUSIONS

While this research remains exploratory, it raises important pathophysiological considerations for brain health, adverse cerebrovascular events, and dementia across the sexes. Our findings also highlight the need to take into account sex differences when investigating cerebral characteristics in humans.

Imaging faster neural dynamics with fast fMRI: A need for updated models of the hemodynamic response

Fast fMRI enables the detection of neural dynamics over timescales of hundreds of milliseconds, suggesting it may provide a new avenue for studying subsecond neural processes in the human brain. The magnitudes of these fast fMRI dynamics are far greater than predicted by canonical models of the hemodynamic response. Several studies have established nonlinear properties of the hemodynamic response that have significant implications for fast fMRI. We first review nonlinear properties of the hemodynamic response function that may underlie fast fMRI signals. We then illustrate the breakdown of canonical hemodynamic response models in the context of fast neural dynamics. We will then argue that the canonical hemodynamic response function is not likely to reflect the BOLD response to neuronal activity driven by sparse or naturalistic stimuli or perhaps to spontaneous neuronal fluctuations in the resting state. These properties suggest that fast fMRI is capable of tracking surprisingly fast neuronal dynamics, and we discuss the neuroscientific questions that could be addressed using this approach.

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.

Hidden Topological Structure of Flow Network Functionality

The ability to reroute and control flow is vital to the function of venation networks across a wide range of organisms. By modifying individual edges in these networks, either by adjusting edge conductances or creating and destroying edges, organisms robustly control the propagation of inputs to perform specific tasks. However, a fundamental disconnect exists between the structure and function: networks with different local architectures can perform the same functions. Here, we answer the question of how changes at the level of individual edges collectively create functionality at the scale of an entire network. Using persistent homology, we analyze networks tuned to perform complex tasks. We find that the responses of such networks encode a hidden topological structure composed of sectors of nearly uniform pressure. Although these sectors are not apparent in the underlying network structure, they correlate strongly with the tuned function. The connectivity of these sectors, rather than that of individual nodes, provides a quantitative relationship between structure and function in flow networks.

Cerebral amyloid angiopathy aggravates perivascular clearance impairment in an Alzheimer's disease mouse model

Cerebral amyloid angiopathy (CAA), defined as the accumulation of amyloid-beta (Aβ) on the vascular wall, is a major pathology of Alzheimer’s disease (AD) and has been thought to be caused by the failure of Aβ clearance. Although two types of perivascular clearance mechanisms, intramural periarterial drainage (IPAD) and the perivascular cerebrospinal fluid (CSF) influx, have been identified, the exact contribution of CAA on perivascular clearance is still not well understood. In this study, we investigated the effect of CAA on the structure and function of perivascular clearance systems in the APP/PS1 transgenic mouse model. To investigate the pathological changes accompanied by CAA progression, the key elements of perivascular clearance such as the perivascular basement membrane, vascular smooth muscle cells (vSMCs), and vascular pulsation were evaluated in middle-aged (7–9 months) and old-aged (19–21 months) mice using in vivo imaging and immunofluorescence staining. Changes in IPAD and perivascular CSF influx were identified by ex vivo fluorescence imaging after dextran injection into the parenchyma or cisterna magna. Amyloid deposition on the vascular wall disrupted the integrity and morphology of the arterial basement membrane. With CAA progression, vascular pulsation was augmented, and conversely, vSMC coverage was decreased. These pathological changes were more pronounced in the surface arteries with earlier amyloid accumulation than in penetrating arteries. IPAD and perivascular CSF influx were impaired in the middle-aged APP/PS1 mice and further aggravated in old age, showing severely impaired tracer influx and efflux patterns. Reduced clearance was also observed in old wild-type mice without changing the tracer distribution pattern in the influx and efflux pathway. These findings suggest that CAA is not merely a consequence of perivascular clearance impairment, but rather a contributor to this process, causing changes in arterial function and structure and increasing AD severity.

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.

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.

Non-Invasive Ultrasound Detection of Cerebrovascular Changes in a Mouse Model of Traumatic Brain Injury

Traumatic brain injury (TBI) can induce changes in vascular architecture. Although ultrasound metrics such as pulsatility index (PI) are sensitive to changes in hemodynamic resistance downstream from major arteries, these metrics depend on features unrelated to vessel architecture, such as blood pressure and heart rate. In contrast, input impedance and reflection coefficient that are derived from wave reflection theory seek to minimize the effects of altered cardiac output or heart rate. In this article, we investigate the use of ultrasound to assess changes in vascular impedance and wave reflection in the common carotid arteries of mice exposed to a controlled cortical impact. Focusing on the first harmonics of the reflected waves, the impedance phase was increased ipsilaterally in impacted mice compared with shams, whereas the magnitude of the impedance was unchanged. In contrast, PI was reduced bilaterally. Interestingly, PI and the first harmonic magnitude of input impedance in the carotid artery were correlated on the contralateral but not ipsilateral side. We investigated the use of these metrics to classify mice as sham or TBI, finding an area under the receiver operating characteristic curve ipsilaterally of 0.792 (confidence interval [CI]: 0.648-0.936) for correct classification with first harmonic impedance magnitude and phase as predictors and 0.716 (CI: 0.553-0.879) using carotid artery PI and diameter as predictors. Overall, the findings support the use of wave reflection analysis as a more specific measure of vascular changes following TBI and motivate the translation of this approach for monitoring vascular changes in humans affected by TBI.

Power spectrum of resting-state blood-oxygen-level-dependent signal

Hemodynamic modeling is used to explore the origin, predict, and analyze the power spectrum of the resting-state blood-oxygen-level-dependent (BOLD) signal measured by functional magnetic resonance imaging (fMRI), which has been reported to have a power-law form, i.e., P(f)∝f^{-s}, where P(f) is the power, f is the frequency, and s>0 is the power-law exponent. However, current fMRI experimental paradigms have limited acquisition durations, affecting the spectral resolution of fMRI data at the low-frequency regime. Here, the claimed power-law spectrum is investigated by using a recent hemodynamic model to analytically derive the BOLD power spectrum, with parameters that are related to neurophysiology. The theoretical results show that, for all realistic parameter combinations, the BOLD power spectrum is flat at f≲0.01Hz, has a weak resonance originating from intrinsic oscillations of vasodilatory response, and becomes a power law for high frequencies, all of which is in agreement with an empirical data set that describes the spectrum of one subject and brain region. However, the results are contrary to studies reporting a pure power-law spectrum at f≲0.2Hz. The discrepancy is attributed largely to data averaging employed by current approaches that averages together important properties of the BOLD power spectrum, such as its resonance, that biases the spectrum to only show a power law. Data averaging also reduces the high-frequency power-law exponent relative to individual cases. Overall, this work demonstrates how the model can reproduce BOLD dynamics and further analyze its low-frequency behavior. Moreover, it also uses the model to explain the impact of procedures, such as data averaging, on the reported features of the BOLD power spectrum.

Measurement and visualization of stimulus-evoked tissue dynamics in mouse barrel cortex using phase-sensitive optical coherence tomography

We describe a method to measure tissue dynamics in mouse barrel cortex during functional activation via phase-sensitive optical coherence tomography (PhS-OCT). The method measures the phase changes in OCT signals, which are induced by the tissue volume change, upon which to localize the activated tissue region. Phase unwrapping, compensation and normalization are applied to increase the dynamic range of the OCT phase detection. To guide the OCT scanning, intrinsic optical signal imaging (IOSI) system equipped with a green light laser source (532 nm) is integrated with the PhS-OCT system to provide a full field time-lapsed images of the reflectance that is used to identify the transversal 2D localized tissue response in the mouse brain. The OCT results show a localized decrease in the OCT phase signal in the activated region of the mouse brain tissue. The decrease in the phase signal may be originated from the brain tissue compression caused by the vasodilatation in the activated region. The activated region revealed in the cross-sectional OCT image is consistent with that identified by the IOSI imaging, indicating the phase change in the OCT signals may associate with the changes in the corresponding hemodynamics. In vivo localized tissue dynamics in the barrel cortex at depth during whisker stimulation is observed and monitored in this study.

Robust blind spectral unmixing for fluorescence microscopy using unsupervised learning

Due to the overlapping emission spectra of fluorophores, fluorescence microscopy images often have bleed-through problems, leading to a false positive detection. This problem is almost unavoidable when the samples are labeled with three or more fluorophores, and the situation is complicated even further when imaged under a multiphoton microscope. Several methods have been developed and commonly used by biologists for fluorescence microscopy spectral unmixing, such as linear unmixing, non-negative matrix factorization, deconvolution, and principal component analysis. However, they either require pre-knowledge of emission spectra or restrict the number of fluorophores to be the same as detection channels, which highly limits the real-world applications of those spectral unmixing methods. In this paper, we developed a robust and flexible spectral unmixing method: Learning Unsupervised Means of Spectra (LUMoS), which uses an unsupervised machine learning clustering method to learn individual fluorophores’ spectral signatures from mixed images, and blindly separate channels without restrictions on the number of fluorophores that can be imaged. This method highly expands the hardware capability of two-photon microscopy to simultaneously image more fluorophores than is possible with instrumentation alone. Experimental and simulated results demonstrated the robustness of LUMoS in multi-channel separations of two-photon microscopy images. We also extended the application of this method to background/autofluorescence removal and colocalization analysis. Lastly, we integrated this tool into ImageJ to offer an easy to use spectral unmixing tool for fluorescence imaging. LUMoS allows us to gain a higher spectral resolution and obtain a cleaner image without the need to upgrade the imaging hardware capabilities.

Limits of multifunctionality in tunable networks

Nature is rife with networks that are functionally optimized to propagate inputs to perform specific tasks. Whether via genetic evolution or dynamic adaptation, many networks create functionality by locally tuning interactions between nodes. Here we explore this behavior in two contexts: strain propagation in mechanical networks and pressure redistribution in flow networks. By adding and removing links, we are able to optimize both types of networks to perform specific functions. We define a single function as a tuned response of a single "target" link when another, predetermined part of the network is activated. Using network structures generated via such optimization, we investigate how many simultaneous functions such networks can be programed to fulfill. We find that both flow and mechanical networks display qualitatively similar phase transitions in the number of targets that can be tuned, along with the same robust finite-size scaling behavior. We discuss how these properties can be understood in the context of constraint-satisfaction problems.

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 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.

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

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

Cerebral vascular structure in the motor cortex of adult mice is stable and is not altered by voluntary exercise

The cerebral vasculature provides blood flow throughout the brain, and local changes in blood flow are regulated to match the metabolic demands of the active brain regions. This neurovascular coupling is mediated by real-time changes in vessel diameter and depends on the underlying vascular network structure. Neurovascular structure is configured during development by genetic and activity-dependent factors. In adulthood, it can be altered by experiences such as prolonged hypoxia, sensory deprivation and seizure. Here, we have sought to determine whether exercise could alter cerebral vascular structure in the adult mouse. We performed repeated in vivo two-photon imaging in the motor cortex of adult transgenic mice expressing membrane-anchored green fluorescent protein in endothelial cells (tyrosine endothelial kinase 2 receptor (Tie2)-Cre:mTmG). This strategy allows for high-resolution imaging of the vessel walls throughout the lifespan. Vascular structure, as measured by capillary branch point number and position, segment diameter and length remained stable over a time scale of months as did pericyte number and position. Furthermore, we compared the vascular structure before, during, and after periods of voluntary wheel running and found no alterations in these same parameters. In both running and control mice, we observed a low rate of capillary segment subtraction. Interestingly, these rare subtraction events preferentially remove short vascular loops.

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.

3D morphological analysis of the mouse cerebral vasculature: Comparison of in vivo and ex vivo methods

Ex vivo 2-photon fluorescence microscopy (2PFM) with optical clearing enables vascular imaging deep into tissue. However, optical clearing may also produce spherical aberrations if the objective lens is not index-matched to the clearing material, while the perfusion, clearing, and fixation procedure may alter vascular morphology. We compared in vivo and ex vivo 2PFM in mice, focusing on apparent differences in microvascular signal and morphology. Following in vivo imaging, the mice (four total) were perfused with a fluorescent gel and their brains fructose-cleared. The brain regions imaged in vivo were imaged ex vivo. Vessels were segmented in both images using an automated tracing algorithm that accounts for the spatially varying PSF in the ex vivo images. This spatial variance is induced by spherical aberrations caused by imaging fructose-cleared tissue with a water-immersion objective. Alignment of the ex vivo image to the in vivo image through a non-linear warping algorithm enabled comparison of apparent vessel diameter, as well as differences in signal. Shrinkage varied as a function of diameter, with capillaries rendered smaller ex vivo by 13%, while penetrating vessels shrunk by 34%. The pial vasculature attenuated in vivo microvascular signal by 40% 300 μm below the tissue surface, but this effect was absent ex vivo. On the whole, ex vivo imaging was found to be valuable for studying deep cortical vasculature.

Influence of cerebral blood vessel movements on the position of perivascular synapses

Synaptic activity is regulated and limited by blood flow, which is controlled by blood vessel dilation and contraction. Traditionally, the study of neurovascular coupling has mainly focused on energy consumption and oxygen delivery. However, the mechanical changes that blood vessel movements induce in the surrounding tissue have not been considered. We have modeled the mechanical changes that movements of blood vessels cause in neighboring synapses. Our simulations indicate that synaptic densities increase or decrease during vascular dilation and contraction, respectively, near the blood vessel walls. This phenomenon may alter the concentration of neurotransmitters and vasoactive substances in the immediate vicinity of the vessel wall and thus may have an influence on local blood flow.

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.

Mechanistic Mathematical Modeling Tests Hypotheses of the Neurovascular Coupling in fMRI

Functional magnetic resonance imaging (fMRI) measures brain activity by detecting the blood-oxygen-level dependent (BOLD) response to neural activity. The BOLD response depends on the neurovascular coupling, which connects cerebral blood flow, cerebral blood volume, and deoxyhemoglobin level to neuronal activity. The exact mechanisms behind this neurovascular coupling are not yet fully investigated. There are at least three different ways in which these mechanisms are being discussed. Firstly, mathematical models involving the so-called Balloon model describes the relation between oxygen metabolism, cerebral blood volume, and cerebral blood flow. However, the Balloon model does not describe cellular and biochemical mechanisms. Secondly, the metabolic feedback hypothesis, which is based on experimental findings on metabolism associated with brain activation, and thirdly, the neurotransmitter feed-forward hypothesis which describes intracellular pathways leading to vasoactive substance release. Both the metabolic feedback and the neurotransmitter feed-forward hypotheses have been extensively studied, but only experimentally. These two hypotheses have never been implemented as mathematical models. Here we investigate these two hypotheses by mechanistic mathematical modeling using a systems biology approach; these methods have been used in biological research for many years but never been applied to the BOLD response in fMRI. In the current work, model structures describing the metabolic feedback and the neurotransmitter feed-forward hypotheses were applied to measured BOLD responses in the visual cortex of 12 healthy volunteers. Evaluating each hypothesis separately shows that neither hypothesis alone can describe the data in a biologically plausible way. However, by adding metabolism to the neurotransmitter feed-forward model structure, we obtained a new model structure which is able to fit the estimation data and successfully predict new, independent validation data. These results open the door to a new type of fMRI analysis that more accurately reflects the true neuronal activity.

Functional magnetic resonance imaging (fMRI) is a widely used technique for measuring brain activity. However, the signal registered by fMRI is not a direct measurement of the neuronal activity in the brain, but it is influenced by the interplay between the metabolism, blood flow and blood volume in the active area. This signal is called the blood-oxygen-level dependent (BOLD) response and occurs when the blood supply to the active area increases in response to neuronal activity. The mechanisms that the cells use to influence the blood supply are not fully known, and therefore it is difficult to know the true neuronal signalling only from inspection of the fMRI signal. In this article, we present a new mathematical model built on the physiological mechanisms thought to underlie the BOLD response. We could successfully fit the model to data and predict the activity caused by new stimuli. By using the validated model we investigated physiological mechanisms that cause different parts of the BOLD response.