Cerebral microvascular network geometry changes in response to functional stimulation.

The cortical microvessels are organized in an intricate, hierarchical, three-dimensional network. Superimposed on this anatomical complexity is the highly complicated signaling that drives the focal blood flow adjustments following a rise in the activity of surrounding neurons. The microvascular response to neuronal activation remains incompletely understood. We developed a custom two photon fluorescence microscopy acquisition and analysis to obtain 3D maps of neuronal activation-induced changes in the geometry of the microvascular network of the primary somatosensory cortex of anesthetized rats. An automated, model-based tracking algorithm was employed to reconstruct the 3D microvascular topology and represent it as a graph. The changes in the geometry of this network were then tracked, over time, in the course of electrical stimulation of the contralateral forepaw. Both dilatory and constrictory responses were observed across the network. Early dilatory and late constrictory responses propagated from deeper to more superficial cortical layers while the response of the vertices that showed initial constriction followed by later dilation spread from cortical surface toward increasing cortical depths. Overall, larger caliber adjustments were observed deeper inside the cortex. This work yields the first characterization of the spatiotemporal pattern of geometric changes on the level of the cortical microvascular network as a whole and provides the basis for bottom-up modeling of the hemodynamically-weighted neuroimaging signals.

Quantification of blood flow and volume in arterioles and venules of the rat cerebral cortex using functional micro-ultrasound.

Relative cerebral blood volume (rCBV), relative cerebral blood flow (rCBF), and blood flow speed are key parameters that characterize cerebral hemodynamics. We used contrast-enhanced functional micro-ultrasound (fMUS) imaging employing a disruption-replenishment imaging sequence to quantify these hemodynamic parameters in the anesthetized rat brain. The method has a spatial resolution of about 100 µm in-plane and around 600 µm through-plane, which is comparable to fMRI, and it has a superior temporal resolution of 40 ms per frame. We found no significant difference in rCBV of cortical and subcortical gray matter (0.89 ± 0.08 and 0.61 ± 0.09 times the brain-average value, respectively). The rCBV was significantly higher in the vascular regions on the pial surface (3.89 ± 0.71) and in the area of major vessels in the subcortical gray matter (2.02 ± 0.31). Parametric images of rCBV, rCBF, and blood flow speed demonstrate spatial heterogeneity of these parameters on the 100 µm scale. Segmentation of the cortex in arteriolar and venular-dominated regions identified through color Doppler imaging showed that rCBV is higher and flow speed is lower in venules than in arterioles. Finally, we show that the dependence of rCBV on rCBF was significantly different in cortical versus subcortical gray matter: the exponent α in the power law relation rCBV=s·rCBF(α) was 0.37 ± 0.13 in cortical and 0.75 ± 0.16 in subcortical gray matter. This work demonstrates that functional micro-ultrasound imaging affords quantification of hemodynamic parameters in the anesthetized rodent brain. This modality is a promising tool for neuroscientists studying these parameters in rodent models of diseases with a cerebrovascular component, such as stroke, neurodegeneration, and venous collagenosis. It is of particular import for studying conditions that selectively affect arteriolar versus venular compartments.

Robust quantification of microvascular transit times via linear dynamical systems using two-photon fluorescence microscopy data.

Vascular transit time is an important indicator of microcirculatory health. We present a second-order-plus-dead-time (SOPDT) model for robust estimation of kinetic parameters characterizing microvascular bolus passage using two-photon fluorescence microscopy (2PFM) in anesthetized rats receiving somatosensory stimulation. This methodology enables quantification of transit time, time-to-peak, overshoot, and rate of bolus passage through the microvascular network. The overall transit time during stimulation, of 2.2 ± 0.1 seconds, was shorter (P ~ 0.0008) than that at rest (2.7 ± 0.2 seconds). When compared with conventional γ-variate modeling, the SOPDT modeling yielded better quality of fit both at rest (P<0.0001) and on activation (P<0.001).

Quantitative estimates of stimulation-induced perfusion response using two-photon fluorescence microscopy of cortical microvascular networks.

Functional hyperemia, or the increase in focal perfusion elicited by neuronal activation, is one of the primary functions of the neurovascular unit and a hallmark of healthy brain functioning. While much is known about the hemodynamics on the millimeter to tenths of millimeter-scale accessible by MRI, there is a paucity of quantitative data on the micrometer-scale changes in perfusion in response to functional stimulation. We present a novel methodology for quantification of perfusion and intravascular flow across the 3D microvascular network in the rat somatosensory cortex using two-photon fluorescence microscopy (2PFM). For approximately 96% of responding microvessels in the forelimb representation of the primary somatosensory cortex, brief (~2s) forepaw stimulation resulted in an increase of perfusion 20±4% (mean±sem). The perfusion levels associated with the remaining 4% of the responding microvessels decreased 10±9% upon stimulation. Vessels irrigating regions of lower vascular density were found to exhibit higher flow (p<0.02), supporting the notion that local vascular morphology and hemodynamics reflect the metabolic needs of the surrounding parenchyma. High dispersion (~77%) in perfusion levels suggests high spatial variation in tissue susceptibility to hypoxia. The current methodology enables quantification of absolute perfusion associated with individual vessels of the cortical microvascular bed and its changes in response to functional stimulation and thereby provides an important tool for studying the cellular mechanisms of functional hyperemia, the spatial specificity of perfusion response to functional stimulation, and, broadly, the micrometer-scale relationship between vascular morphology and function in health and disease.

A constrained independent component analysis technique for artery-vein separation of two-photon laser scanning microscopy images of the cerebral microvasculature.

Understanding brain hemodynamics as well as the coupling between microvascular hemodynamics and neural activity is important in pathophysiology of cerebral microvasculature. When local increases in neuronal activity occur, the blood volume changes in the surrounding brain vasculature. Dynamic contrast enhanced imaging (DCE) is a powerful technique that quantifies these changes in the blood flow by repeatedly imaging the vasculature over time. Separating artery, vein and capillaries in the images and extracting their intensity-time curves from the DCE image sequence is an important first step in understanding vascular function. A constrained independent component analysis (ICA) technique is developed to analyze the two photon laser scanning microscopy (2PLSM) images of rat brain microvasculature, where a bolus of fluorescent dye is administered to the vascular system as the contrast agent. A priori information inferred from the gamma variate model of cerebral microvasculature is incorporated with the data driven technique in temporal and spatial domains using two constraints. The constraints are: no independent component (IC) is allowed to have negative contribution in forming the images (positivity constraint) and the component curves follow a gamma variate function (model fitting constraint). Experimental and simulation studies are conducted to demonstrate the improved performance of the proposed constrained ICA (CICA) technique over the most commonly used classical ICA algorithm (fast-ICA) in providing physiologically meaningful ICs and its ability to separate the model following factors from other factors are shown. The efficiency of CICA in handling noise is compared to model based techniques. Its capability in providing improved separation between artery, vein and capillaries compared to the other two techniques is also demonstrated.

Functional micro-ultrasound imaging of rodent cerebral hemodynamics.

Healthy cerebral microcirculation is crucial to neuronal functioning. We present a new method to investigate microvascular hemodynamics in living rodent brain through a focal cranial window based on high-frequency ultrasound imaging. The method has a temporal resolution of 40ms, and a 100µm in-plane and 600µm through-plane spatial resolution. We use a commercially available high-frequency ultrasound imaging system to quantify changes in the relative cerebral blood volume (CBV) by measuring the scattered signal intensity from an ultrasound contrast agent circulating in the vasculature. Generalized linear model analysis is then used to produce effect size and significance maps of changes in cerebral blood volume upon electrical stimulation of the forepaw. We observe larger CBV increases in the forelimb representation of the primary somatosensory cortex than in the deep gray matter with stimuli as short as 2s (5.1 ± 1.3% vs. 3.3 ± 0.6%). We also investigate the temporal evolution of the blood volume changes in cortical and subcortical gray matter, pial vessels and subcortical major vessels, and show shorter response onset times in the parenchymal regions than in the neighboring large vessels (1.6 ± 1.0s vs. 2.6 ± 1.3s in the cortex for a 10 second stimulus protocol). This method, which we termed functional micro-ultrasound imaging or fMUS, is a novel, highly accessible, and cost-effective way of imaging rodent brain microvascular topology and hemodynamics in vivo at 100micron resolution over a 1-by-1cm field of view with 10s-100s frames per second that opens up a new set of questions regarding brain function in preclinical models of health and disease.

Two-photon fluorescence microscopy of cerebral hemodynamics.

Under physiological conditions, neuronal activity is tightly coupled to hemodynamics of the surrounding microvessels. Conversely, most brain diseases are associated with a disturbance in neurovascular coupling. Measuring the hemodynamic response of the microvascular network to neuronal stimulation in vivo involves two major challenges: maintaining a stable systemic physiological state in the animal and imaging with sufficient temporal and spatial resolution to capture the hemodynamic changes across the three-dimensional cortical microvascular network. Two-photon fluorescence microscopy allows imaging of intact cortex in situ at micrometer spatial and microsecond temporal resolution. However, this modality necessitates focal opening of the skull because of its high scattering of light waves. This protocol describes in detail the requisite surgical preparation and physiological maintenance of an adolescent rat with a closed cranial window. Two-photon fluorescence laser scanning microscopy is then used to image the hemodynamic response of the microvasculature of the primary somatosensory cortex to electrical stimulation of the forepaw, with the end goal of quantitative analysis of brain hemodynamics.

In vivo imaging of cerebral hemodynamics using high-frequency micro-ultrasound.

Assessment of cerebral vascular response is important in neuroscience research. Some imaging modalities that are commonly used to detect flow and/or vessel diameter changes in the brain include magnetic resonance imaging, positron emission tomography, and optical intrinsic signal imaging. Ultrasound has not typically been used to assess neurovascular response but recent advances in the technology have led to the development of micro-ultrasound systems with significant potential for this application. The state of the art in high frequency (15-50 MHz) micro-ultrasound is based on linear arrays specifically designed for small animal imaging. These systems can achieve axial resolution ranging from 30 to 200 microm. They are capable of quantifying brain hemodynamics in terms of red blood cell (RBC) velocity, flow, and vascular density in real time, up to 35 mm below the cortical surface, and can achieve temporal resolution of up to 1000 frames per second. This protocol describes imaging of the rat brain using various ultrasound imaging modes (power Doppler, color Doppler, pulsed-wave Doppler, and nonlinear contrast-enhanced imaging) to assess the state of cerebral microcirculation.

Temperature change near microbubbles within a capillary network during focused ultrasound.

Preformed gas bubbles can increase energy absorption from an ultrasound beam and therefore they have been proposed for an enhancer of ultrasound treatments. Although tissue temperature measurements performed in vivo using invasive thermocouple probes and MRI thermometry have demonstrated increased tissue temperature, the microscopic temperature distribution has not been investigated so far. In this study the transfer of heat between bubbles and tissue during focused ultrasound was simulated. Microbubble oscillations were simulated within a rat cortical microvascular network reconstructed from in vivo dual-photon microscopy images and the power density of these oscillations was used as an input term in the Pennes bioheat transfer equation. The temperature solution from the bioheat transfer equation was mapped onto vascular data to produce a three-dimensional temperature map. The results showed high temperatures near the bubbles and slow temperature rise in the tissue. Heating was shown to increase with increasing bubble frequency and insonation pressure, and showed a frequency-dependent peak. The goal of this research is to characterize the effect of various parameters on bubble-enhanced therapeutic ultrasound to allow better treatment planning. These results show that the induced temperature elevations have nonuniformities which may have a significant impact on the bio-effects of the exposure.