Large arteriolar component of oxygen delivery implies a safe margin of oxygen supply to cerebral tissue

Quantitative Optical Imaging of Oxygen in Brain Vasculature

The intimate relationship between neuronal activity and cerebral oxygenation underpins fundamental brain functions like cognition, sensation, and motor control. Optical imaging offers a noninvasive approach to assess brain oxygenation and often serves as an indirect proxy for neuronal activity. However, deciphering neurovascular coupling─the intricate interplay between neuronal activity, blood flow, and oxygen delivery─necessitates independent, high spatial resolution, and high temporal resolution measurements of both microvasculature oxygenation and neuronal activation. This Perspective examines the established optical techniques employed for brain oxygen imaging, specifically functional near-infrared spectroscopy, photoacoustic imaging, optical coherence tomography, and two-photon phosphorescent lifetime microscopy, highlighting their fundamental principles, strengths, and limitations. Several other emerging optical techniques are also introduced. Finally, we discuss key technological challenges and future directions for quantitative optical oxygen imaging, paving the way for a deeper understanding of oxygen metabolism in the brain.

The relationship between ischemic penumbra progression and the oxygen content of cortex microcirculation in acute ischemic stroke

The precise oxygen content thresholds of ischemic deep parenchymal (OCIDP) and that in cortical microcirculation (OCCM), which leads to ischemic penumbra converting into the infarcted core, remain uncertain. This study employed an invasive fiber-optic oxygen meter and a newly developed oxygen-responsive probe called RuA3-Cy5-rtPA (RC-rtPA) based on recombinant tissue-type plasminogen activator (rtPA) to examine the oxygen content thresholds. A mouse model of middle cerebral artery occlusion was generated and animals were randomly divided into a sham, 24-h reperfusion after 3-h ischemia (IR 3-h), and IR 6-h groups, all of which were sacrificed following reperfusion. Stroke severity was evaluated based on the infarction area, neurological symptoms, microcirculation perfusion, and microemboli in microcirculation. OCIDP was characterized based on its extent and distribution, whereas OCCM was measured using RC-rtPA. During ischemia, stroke severity escalation manifested as increasing infarction area, severe neurologic symptoms, and poorer microcirculation perfusion with more microthrombi depositions. OCIDP presented rapid decline following artery occlusion along with a gradual increase in the hypoxic area. Within 3 ​h following ischemia induction, the ischemic tissue that experienced hypoxia could be rescued, and this reversibility would disappear after 6 ​h. Within 6 ​h, OCCM continued to decrease. A significant decrease in oxygen content in cortical venules and cortical parenchyma was observed. These findings assist in establishing the extent of the ischemic penumbra at the microcirculation level and offer a foundation for assessing the ischemic penumbra that could respond positively to reperfusion therapy beyond the typical time window.

Modeling hypoxia-induced radiation resistance and the impact of radiation sources

Hypoxia contributes significantly to resistance in radiotherapy. Our research rigorously examines the influence of microvascular morphology on radiotherapy outcome, specifically focusing on how microvasculature shapes hypoxia within the microenvironment and affects resistance to a standard treatment regimen (30×2GyRBE). Our computational modeling extends to the effects of different radiation sources. For photons and protons, our analysis establishes a clear correlation between hypoxic volume distribution and treatment effectiveness, with vascular density and regularity playing a crucial role in treatment success. On the contrary, carbon ions exhibit distinct effectiveness, even in areas of intense hypoxia and poor vascularization. This finding points to the potential of carbon-based hadron therapy in overcoming hypoxia-induced resistance to RT. Considering that the spatial scale analyzed in this study is closely aligned with that of imaging data voxels, we also address the implications of these findings in a clinical context envisioning the possibility of detecting subvoxel hypoxia.

Neuroinflammation increases oxygen extraction in a mouse model of Alzheimer's disease

BACKGROUND

Neuroinflammation, impaired metabolism, and hypoperfusion are fundamental pathological hallmarks of early Alzheimer's disease (AD). Numerous studies have asserted a close association between neuroinflammation and disrupted cerebral energetics. During AD progression and other neurodegenerative disorders, a persistent state of chronic neuroinflammation reportedly exacerbates cytotoxicity and potentiates neuronal death. Here, we assessed the impact of a neuroinflammatory challenge on metabolic demand and microvascular hemodynamics in the somatosensory cortex of an AD mouse model.

METHODS

We utilized in vivo 2-photon microscopy and the phosphorescent oxygen sensor Oxyphor 2P to measure partial pressure of oxygen (pO2) and capillary red blood cell flux in cortical microvessels of awake mice. Intravascular pO2 and capillary RBC flux measurements were performed in 8-month-old APPswe/PS1dE9 mice and wildtype littermates on days 0, 7, and 14 of a 14-day period of lipopolysaccharide-induced neuroinflammation.

RESULTS

Before the induced inflammatory challenge, AD mice demonstrated reduced metabolic demand but similar capillary red blood cell flux as their wild type counterparts. Neuroinflammation provoked significant reductions in cerebral intravascular oxygen levels and elevated oxygen extraction in both animal groups, without significantly altering red blood cell flux in capillaries.

CONCLUSIONS

This study provides evidence that neuroinflammation alters cerebral oxygen demand at the early stages of AD without substantially altering vascular oxygen supply. The results will guide our understanding of neuroinflammation's influence on neuroimaging biomarkers for early AD diagnosis.

A mathematical model for velocity-selective arterial spin labeling

PURPOSE

To present a theoretical framework that rigorously defines and analyzes key concepts and quantities for velocity selective arterial spin labeling (VSASL).

THEORY AND METHODS

An expression for the VSASL arterial delivery function is derived based on (1) labeling and saturation profiles as a function of velocity and (2) physiologically plausible approximations of changes in acceleration and velocity across the vascular system. The dependence of labeling efficiency on the amplitude and effective bolus width of the arterial delivery function is defined. Factors that affect the effective bolus width are examined, and timing requirements to minimize quantitation errors are derived.

RESULTS

The model predicts that a flow-dependent negative bias in the effective bolus width can occur when velocity selective inversion (VSI) is used for the labeling module and velocity selective saturation (VSS) is used for the vascular crushing module. The bias can be minimized by choosing a nominal labeling cutoff velocity that is lower than the nominal cutoff velocity of the vascular crushing module.

CONCLUSION

The elements of the model are specified in a general fashion such that future advances can be readily integrated. The model can facilitate further efforts to understand and characterize the performance of VSASL and provide critical theoretical insights that can be used to design future experiments and develop novel VSASL approaches.

Comparison of cerebral oxygen extraction fraction using ASE and TRUST methods in patients with sickle cell disease and healthy controls

Abnormal oxygen extraction fraction (OEF), a putative biomarker of cerebral metabolic stress, may indicate compromised oxygen delivery and ischemic vulnerability in patients with sickle cell disease (SCD). Elevated OEF was observed at the tissue level across the brain using an asymmetric spin echo (ASE) MR method, while variable global OEFs were found from the superior sagittal sinus (SSS) using a T2-relaxation-under-spin-tagging (TRUST) MRI method with different calibration models. In this study, we aimed to compare the average ASE-OEF in the SSS drainage territory and TRUST-OEF in the SSS from the same SCD patients and healthy controls. 74 participants (SCD: N = 49; controls: N = 25) underwent brain MRI. TRUST-OEF was quantified using the Lu-bovine, Bush-HbA and Li-Bush-HbS models. ASE-OEF and TRUST-OEF were significantly associated in healthy controls after controlling for hematocrit using the Lu-bovine or the Bush-HbA model. However, no association was found between ASE-OEF and TRUST-OEF in patients with SCD using either the Bush-HbA or the Li-Bush-HbS model. Plausible explanations include a discordance between spatially volume-averaged oxygenation brain tissue and flow-weighted volume-averaged oxygenation in SSS or sub-optimal calibration in SCD. Further work is needed to refine and validate non-invasive MR OEF measurements in SCD.

Revisiting hemodynamics and blood oxygenation in a microfluidic microvasculature replica

The complexity of microvascular circulation has led to the development of advanced imaging techniques and biomimetic models. This study developed a multifaceted microfluidic-based microdevice as an in vitro model of microvasculature to replicate important geometric and functional features of in vivo perfusion in mice. The microfluidic device consisted of a microchannel for blood perfusion, mirroring the natural hierarchical branching vascular structures found in mice. Additionally, the device incorporated a steady gradient of oxygen (O2) which diffused through the polydimethylsiloxane (PDMS) layer, allowing for dynamic blood oxygenation. The assembled multi-layered microdevice was accompanied by a dual-modal imaging system that combined laser speckle contrast imaging (LSCI) and intrinsic signal optical imaging (ISOI) to visualize full-field blood flow distributions and blood O2 profiles. By closely reproducing in vivo blood perfusion and oxygenation conditions, this microvasculature model, in conjunction with numerical simulation results, can provide quantitative information on physiologically relevant hemodynamics and key O2 transport parameters that are not directly measurable in traditional animal studies.

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.

Neuroinflammation increases oxygen extraction in a mouse model of Alzheimer's disease

Neuroinflammation, impaired metabolism, and hypoperfusion are fundamental pathological hallmarks of early Alzheimer's disease (AD). Numerous studies have asserted a close association between neuroinflammation and disrupted cerebral energetics. During AD progression and other neurodegenerative disorders, a persistent state of chronic neuroinflammation reportedly exacerbates cytotoxicity and potentiates neuronal death. Here, we assessed the impact of a neuroinflammatory challenge on metabolic demand and microvascular hemodynamics in the somatosensory cortex of an AD mouse model. We utilized in vivo 2-photon microscopy and the phosphorescent oxygen sensor Oxyphor 2P to measure partial pressure of oxygen (pO2) and capillary red blood cell flux in cortical microvessels of awake mice. Intravascular pO2 and capillary RBC flux measurements were performed in 8-month-old APPswe/PS1dE9 mice and wildtype littermates on days 0, 7, and 14 of a 14-day period of lipopolysaccaride-induced neuroinflammation. Before the induced inflammatory challenge, AD mice demonstrated reduced metabolic demand but similar capillary red blood cell flux as their wild type counterparts. Neuroinflammation provoked significant reductions in cerebral intravascular oxygen levels and elevated oxygen extraction in both animal groups, without significantly altering red blood cell flux in capillaries. This study provides evidence that neuroinflammation alters cerebral oxygen demand at the early stages of AD without substantially altering vascular oxygen supply. The results will guide our understanding of neuroinflammation's influence on neuroimaging biomarkers for early AD diagnosis.

Cortical microvascular blood flow velocity mapping by combining dynamic light scattering optical coherence tomography and two-photon microscopy

Significance

The accurate large-scale mapping of cerebral microvascular blood flow velocity is crucial for a better understanding of cerebral blood flow (CBF) regulation. Although optical imaging techniques enable both high-resolution microvascular angiography and fast absolute CBF velocity measurements in the mouse cortex, they usually require different imaging techniques with independent system configurations to maximize their performances. Consequently, it is still a challenge to accurately combine functional and morphological measurements to co-register CBF speed distribution from hundreds of microvessels with high-resolution microvascular angiograms.

Aim

We propose a data acquisition and processing framework to co-register a large set of microvascular blood flow velocity measurements from dynamic light scattering optical coherence tomography (DLS-OCT) with the corresponding microvascular angiogram obtained using two-photon microscopy (2PM).

Approach

We used DLS-OCT to first rapidly acquire a large set of microvascular velocities through a sealed cranial window in mice and then to acquire high-resolution microvascular angiograms using 2PM. The acquired data were processed in three steps: (i) 2PM angiogram coregistration with the DLS-OCT angiogram, (ii) 2PM angiogram segmentation and graphing, and (iii) mapping of the CBF velocities to the graph representation of the 2PM angiogram.

Results

We implemented the developed framework on the three datasets acquired from the mice cortices to facilitate the coregistration of the large sets of DLS-OCT flow velocity measurements with 2PM angiograms. We retrieved the distributions of red blood cell velocities in arterioles, venules, and capillaries as a function of the branching order from precapillary arterioles and postcapillary venules from more than 1000 microvascular segments.

Conclusions

The proposed framework may serve as a useful tool for quantitative analysis of large microvascular datasets obtained by OCT and 2PM in studies involving normal brain functioning, progression of various diseases, and numerical modeling of the oxygen advection and diffusion in the realistic microvascular networks.

Deep optoacoustic localization microangiography of ischemic stroke in mice

Super-resolution optoacoustic imaging of microvascular structures deep in mammalian tissues has so far been impeded by strong absorption from densely-packed red blood cells. Here we devised 5 µm biocompatible dichloromethane-based microdroplets exhibiting several orders of magnitude higher optical absorption than red blood cells at near-infrared wavelengths, thus enabling single-particle detection in vivo. We demonstrate non-invasive three-dimensional microangiography of the mouse brain beyond the acoustic diffraction limit (<20 µm resolution). Blood flow velocity quantification in microvascular networks and light fluence mapping was also accomplished. In mice affected by acute ischemic stroke, the multi-parametric multi-scale observations enabled by super-resolution and spectroscopic optoacoustic imaging revealed significant differences in microvascular density, flow and oxygen saturation in ipsi- and contra-lateral brain hemispheres. Given the sensitivity of optoacoustics to functional, metabolic and molecular events in living tissues, the new approach paves the way for non-invasive microscopic observations with unrivaled resolution, contrast and speed.

The hemodynamic initial-dip consists of both volumetric and oxymetric changes reflecting localized spiking activity

The initial-dip is a transient decrease frequently observed in functional neuroimaging signals, immediately after stimulus onset, believed to originate from a rise in deoxy-hemoglobin (HbR) caused by local neural activity. It has been shown to be more spatially specific than the hemodynamic response, and is believed to represent focal neuronal activity. However, despite being observed in various neuroimaging modalities (such as fMRI, fNIRS, etc), its origins are disputed, and its precise neuronal correlates are unknown. Here we show that the initial-dip is dominated by a decrease in total-hemoglobin (HbT). We also find a biphasic response in deoxy-Hb (HbR), with an early decrease and later rebound. Both the HbT-dip and HbR-rebound were strongly correlated to highly localized spiking activity. However, HbT decreases were always large enough to counter the spiking-induced increase in HbR. We find that the HbT-dip counters spiking induced HbR increases, imposing an upper-limit to HbR concentration in the capillaries. Building on our results, we explore the possibility of active venule dilation (purging) as a possible mechanism for the HbT dip.

Adaptive spectroscopic visible-light optical coherence tomography for clinical retinal oximetry

BACKGROUND

Retinal oxygen saturation (sO2) provides essential information about the eye's response to pathological changes that can result in vision loss. Visible-light optical coherence tomography (vis-OCT) is a noninvasive tool that has the potential to measure retinal sO2 in a clinical setting. However, its reliability is currently limited by unwanted signals referred to as spectral contaminants (SCs), and a comprehensive strategy to isolate true oxygen-dependent signals from SCs in vis-OCT is lacking.

METHODS

We develop an adaptive spectroscopic vis-OCT (ADS-vis-OCT) technique that can adaptively remove SCs and accurately measure sO2 under the unique conditions of each vessel. We also validate the accuracy of ADS-vis-OCT using ex vivo blood phantoms and assess its repeatability in the retina of healthy volunteers.

RESULTS

In ex vivo blood phantoms, ADS-vis-OCT agrees with a blood gas machine with only a 1% bias in samples with sO2 ranging from 0% to 100%. In the human retina, the root mean squared error between sO2 values in major arteries measured by ADS-vis-OCT and a pulse oximeter is 2.1% across 18 research participants. Additionally, the standard deviations of repeated ADS-vis-OCT measurements of sO2 values in smaller arteries and veins are 2.5% and 2.3%, respectively. Non-adaptive methods do not achieve comparable repeatabilities from healthy volunteers.

CONCLUSIONS

ADS-vis-OCT effectively removes SCs from human images, yielding accurate and repeatable sO2 measurements in retinal arteries and veins with varying diameters. This work could have important implications for the clinical use of vis-OCT to manage eye diseases.

Arteriolar neuropathology in cerebral microvascular disease

Cerebral microvascular disease (MVD) is an important cause of vascular cognitive impairment. MVD is heterogeneous in aetiology, ranging from universal ageing to the sporadic (hypertension, sporadic cerebral amyloid angiopathy [CAA] and chronic kidney disease) and the genetic (e.g., familial CAA, cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy [CADASIL] and cerebral autosomal recessive arteriopathy with subcortical infarcts and leukoencephalopathy [CARASIL]). The brain parenchymal consequences of MVD predominantly consist of lacunar infarcts (lacunes), microinfarcts, white matter disease of ageing and microhaemorrhages. MVD is characterised by substantial arteriolar neuropathology involving ubiquitous vascular smooth muscle cell (SMC) abnormalities. Cerebral MVD is characterised by a wide variety of arteriolar injuries but only a limited number of parenchymal manifestations. We reason that the cerebral arteriole plays a dominant role in the pathogenesis of each type of MVD. Perturbations in signalling and function (i.e., changes in proliferation, apoptosis, phenotypic switch and migration of SMC) are prominent in the pathogenesis of cerebral MVD, making 'cerebral angiomyopathy' an appropriate term to describe the spectrum of pathologic abnormalities. The evidence suggests that the cerebral arteriole acts as both source and mediator of parenchymal injury in MVD.

A quantitative model for human neurovascular coupling with translated mechanisms from animals

Neurons regulate the activity of blood vessels through the neurovascular coupling (NVC). A detailed understanding of the NVC is critical for understanding data from functional imaging techniques of the brain. Many aspects of the NVC have been studied both experimentally and using mathematical models; various combinations of blood volume and flow, local field potential (LFP), hemoglobin level, blood oxygenation level-dependent response (BOLD), and optogenetics have been measured and modeled in rodents, primates, or humans. However, these data have not been brought together into a unified quantitative model. We now present a mathematical model that describes all such data types and that preserves mechanistic behaviors between experiments. For instance, from modeling of optogenetics and microscopy data in mice, we learn cell-specific contributions; the first rapid dilation in the vascular response is caused by NO-interneurons, the main part of the dilation during longer stimuli is caused by pyramidal neurons, and the post-peak undershoot is caused by NPY-interneurons. These insights are translated and preserved in all subsequent analyses, together with other insights regarding hemoglobin dynamics and the LFP/BOLD-interplay, obtained from other experiments on rodents and primates. The model can predict independent validation-data not used for training. By bringing together data with complementary information from different species, we both understand each dataset better, and have a basis for a new type of integrative analysis of human data.

Unique expression of the atypical mitochondrial subunit NDUFA4L2 in cerebral pericytes fine tunes HIF activity in response to hypoxia

A central response to insufficient cerebral oxygen delivery is a profound reprograming of metabolism, which is mainly regulated by the Hypoxia Inducible Factor (HIF). Among other responses, HIF induces the expression of the atypical mitochondrial subunit NDUFA4L2. Surprisingly, NDUFA4L2 is constitutively expressed in the brain in non-hypoxic conditions. Analysis of publicly available single cell transcriptomic (scRNA-seq) data sets coupled with high-resolution multiplexed fluorescence RNA in situ hybridization (RNA F.I.S.H.) revealed that in the murine and human brain NDUFA4L2 is exclusively expressed in mural cells with the highest levels found in pericytes and declining along the arteriole-arterial smooth muscle cell axis. This pattern was mirrored by COX4I2, another atypical mitochondrial subunit. High NDUFA4L2 expression was also observed in human brain pericytes in vitro, decreasing when pericytes are muscularized and further induced by HIF stabilization in a PHD2/PHD3 dependent manner. In vivo, Vhl conditional inactivation in pericyte targeting Ng2-cre transgenic mice dramatically induced NDUFA4L2 expression. Finally NDUFA4L2 inactivation in pericytes increased oxygen consumption and therefore the degree of HIF pathway induction in hypoxia. In conclusion our work reveals that NDUFA4L2 together with COX4I2 is a key hypoxic-induced metabolic marker constitutively expressed in pericytes coupling mitochondrial oxygen consumption and cellular hypoxia response.

Could respiration-driven blood oxygen changes modulate neural activity?

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

Depth-Resolved Localization Microangiography in the NIR-II Window

Detailed characterization of microvascular alterations requires high-resolution 3D imaging methods capable of providing both morphological and functional information. Existing optical microscopy tools are routinely used for microangiography, yet offer suboptimal trade-offs between the achievable field of view and spatial resolution with the intense light scattering in biological tissues further limiting the achievable penetration depth. Herein, a new approach for volumetric deep-tissue microangiography based on stereovision combined with super-resolution localization imaging is introduced that overcomes the spatial resolution limits imposed by light diffusion and optical diffraction in wide-field imaging configurations. The method capitalizes on localization and tracking of flowing fluorescent particles in the second near-infrared window (NIR-II, ≈1000-1700 nm), with the third (depth) dimension added by triangulation and stereo-matching of images acquired with two short-wave infrared cameras operating in a dual-view mode. The 3D imaging capability enabled with the proposed method facilitates a detailed visualization of microvascular networks and an accurate blood flow quantification. Experiments performed in tissue-mimicking phantoms demonstrate that high resolution is preserved up to a depth of 4 mm in a turbid medium. Transcranial microangiography of the entire murine cortex and penetrating vessels is further demonstrated at capillary level resolution.

Baseline oxygen consumption decreases with cortical depth

The cerebral cortex is organized in cortical layers that differ in their cellular density, composition, and wiring. Cortical laminar architecture is also readily revealed by staining for cytochrome oxidase-the last enzyme in the respiratory electron transport chain located in the inner mitochondrial membrane. It has been hypothesized that a high-density band of cytochrome oxidase in cortical layer IV reflects higher oxygen consumption under baseline (unstimulated) conditions. Here, we tested the above hypothesis using direct measurements of the partial pressure of O2 (pO2) in cortical tissue by means of 2-photon phosphorescence lifetime microscopy (2PLM). We revisited our previously developed method for extraction of the cerebral metabolic rate of O2 (CMRO2) based on 2-photon pO2 measurements around diving arterioles and applied this method to estimate baseline CMRO2 in awake mice across cortical layers. To our surprise, our results revealed a decrease in baseline CMRO2 from layer I to layer IV. This decrease of CMRO2 with cortical depth was paralleled by an increase in tissue oxygenation. Higher baseline oxygenation and cytochrome density in layer IV may serve as an O2 reserve during surges of neuronal activity or certain metabolically active brain states rather than reflecting baseline energy needs. Our study provides to our knowledge the first quantification of microscopically resolved CMRO2 across cortical layers as a step towards better understanding of brain energy metabolism.

Mapping oxidative metabolism in the human brain with calibrated fMRI in health and disease

Conventional functional MRI (fMRI) with blood-oxygenation level dependent (BOLD) contrast is an important tool for mapping human brain activity non-invasively. Recent interest in quantitative fMRI has renewed the importance of oxidative neuroenergetics as reflected by cerebral metabolic rate of oxygen consumption (CMR_O2) to support brain function. Dynamic CMR_O2 mapping by calibrated fMRI require multi-modal measurements of BOLD signal along with cerebral blood flow (CBF) and/or volume (CBV). In human subjects this "calibration" is typically performed using a gas mixture containing small amounts of carbon dioxide and/or oxygen-enriched medical air, which are thought to produce changes in CBF (and CBV) and BOLD signal with minimal or no CMR_O2 changes. However non-human studies have demonstrated that the "calibration" can also be achieved without gases, revealing good agreement between CMR_O2 changes and underlying neuronal activity (e.g., multi-unit activity and local field potential). Given the simpler set-up of gas-free calibrated fMRI, there is evidence of recent clinical applications for this less intrusive direction. This up-to-date review emphasizes technological advances for such translational gas-free calibrated fMRI experiments, also covering historical progression of the calibrated fMRI field that is impacting neurological and neurodegenerative investigations of the human brain.

Assessment of functional shunting in patients with sickle cell disease

Silent cerebral infarcts (SCI) are common in patients with sickle cell disease (SCD) and are thought to be caused by a mismatch between oxygen delivery and consumption. Functional cerebrovascular shunting is defined as reduced oxygen offloading due to the rapid transit of blood through the capillaries caused by increased flow and has been suggested as a potential mechanism underlying reduced oxygenation and SCI. We investigated the venous arterial spin labeling signal (VS) in the sagittal sinus as a proxy biomarker of cerebral functional shunting, and its association with hemodynamic imaging and hematological laboratory parameters. We included 28 children and 38 adults with SCD, and ten healthy racematched adult controls. VS, cerebral blood flow (CBF), velocity in the brain feeding arteries, oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO2) were measured before and after acetazolamide administration. VS was higher in patients with SCD compared to controls (P<0.01) and was increased after acetazolamide administration in all groups (P<0.01). VS was primarily predicted by CBF (P<0.01), but CBF-corrected VS was also associated with decreased CMRO2 (P<0.01). Additionally, higher disease severity defined by low hemoglobin and increased hemolysis was associated with higher CBF-corrected VS. Finally, CMRO2 was negatively correlated with fetal hemoglobin, and positively correlated with lactate dehydrogenase, which could be explained by changes in oxygen affinity. These findings provide evidence for cerebral functional shunting and encourage future studies investigating the potential link to aberrant capillary exchange in SCD.

Direct vascular contact is a hallmark of cerebral astrocytes

Astrocytes establish extensive networks via gap junctions that allow each astrocyte to connect indirectly to the vasculature. However, the proportion of astrocytes directly associated with blood vessels is unknown. Here, we quantify structural contacts of cortical astrocytes with the vasculature in vivo. We show that all cortical astrocytes are connected to at least one blood vessel. Moreover, astrocytes contact more vessels in deeper cortical layers where vessel density is known to be higher. Further examination of different brain regions reveals that only the hippocampus, which has the lowest vessel density of all investigated brain regions, harbors single astrocytes with no apparent vascular connection. In summary, we show that almost all gray matter astrocytes have direct contact to the vasculature. In addition to the glial network, a direct vascular access may represent a complementary pathway for metabolite uptake and distribution.

Optical measurement of microvascular oxygenation and blood flow responses in awake mouse cortex during functional activation

The cerebral cortex has a number of conserved morphological and functional characteristics across brain regions and species. Among them, the laminar differences in microvascular density and mitochondrial cytochrome c oxidase staining suggest potential laminar variability in the baseline O₂ metabolism and/or laminar variability in both O₂ demand and hemodynamic response. Here, we investigate the laminar profile of stimulus-induced intravascular partial pressure of O₂ (pO2) transients to stimulus-induced neuronal activation in fully awake mice using two-photon phosphorescence lifetime microscopy. Our results demonstrate that stimulus-induced changes in intravascular pO₂ are conserved across cortical layers I-IV, suggesting a tightly controlled neurovascular response to provide adequate O₂ supply across cortical depth. In addition, we observed a larger change in venular O₂ saturation (ΔsO₂) compared to arterioles, a gradual increase in venular ΔsO₂ response towards the cortical surface, and absence of the intravascular "initial dip" previously reported under anesthesia. This study paves the way for quantification of layer-specific cerebral O₂ metabolic responses, facilitating investigation of brain energetics in health and disease and informed interpretation of laminar blood oxygen level dependent functional magnetic resonance imaging signals.

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.

Pericyte Control of Blood Flow Across Microvascular Zones in the Central Nervous System

The vast majority of the brain's vascular length is composed of capillaries, where our understanding of blood flow control remains incomplete. This review synthesizes current knowledge on the control of blood flow across microvascular zones by addressing issues with nomenclature and drawing on new developments from in vivo optical imaging and single-cell transcriptomics. Recent studies have highlighted important distinctions in mural cell morphology, gene expression, and contractile dynamics, which can explain observed differences in response to vasoactive mediators between arteriole, transitional, and capillary zones. Smooth muscle cells of arterioles and ensheathing pericytes of the arteriole-capillary transitional zone control large-scale, rapid changes in blood flow. In contrast, capillary pericytes downstream of the transitional zone act on slower and smaller scales and are involved in establishing resting capillary tone and flow heterogeneity. Many unresolved issues remain, including the vasoactive mediators that activate the different pericyte types in vivo, the role of pericyte-endothelial communication in conducting signals from capillaries to arterioles, and how neurological disease affects these mechanisms.

Neurophotonic tools for microscopic measurements and manipulation: status report

Neurophotonics was launched in 2014 coinciding with the launch of the BRAIN Initiative focused on development of technologies for advancement of neuroscience. For the last seven years, Neurophotonics' agenda has been well aligned with this focus on neurotechnologies featuring new optical methods and tools applicable to brain studies. While the BRAIN Initiative 2.0 is pivoting towards applications of these novel tools in the quest to understand the brain, this status report reviews an extensive and diverse toolkit of novel methods to explore brain function that have emerged from the BRAIN Initiative and related large-scale efforts for measurement and manipulation of brain structure and function. Here, we focus on neurophotonic tools mostly applicable to animal studies. A companion report, scheduled to appear later this year, will cover diffuse optical imaging methods applicable to noninvasive human studies. For each domain, we outline the current state-of-the-art of the respective technologies, identify the areas where innovation is needed, and provide an outlook for the future directions.

Network-driven anomalous transport is a fundamental component of brain microvascular dysfunction

Blood microcirculation supplies neurons with oxygen and nutrients, and contributes to clearing their neurotoxic waste, through a dense capillary network connected to larger tree-like vessels. This complex microvascular architecture results in highly heterogeneous blood flow and travel time distributions, whose origin and consequences on brain pathophysiology are poorly understood. Here, we analyze highly-resolved intracortical blood flow and transport simulations to establish the physical laws governing the macroscopic transport properties in the brain micro-circulation. We show that network-driven anomalous transport leads to the emergence of critical regions, whether hypoxic or with high concentrations of amyloid-β, a waste product centrally involved in Alzheimer's Disease. We develop a Continuous-Time Random Walk theory capturing these dynamics and predicting that such critical regions appear much earlier than anticipated by current empirical models under mild hypoperfusion. These findings provide a framework for understanding and modelling the impact of microvascular dysfunction in brain diseases, including Alzheimer's Disease.

Modelling the effects of cerebral microthrombi on tissue oxygenation and cell death

Thrombectomy, the mechanical removal of a clot, is the most common way to treat ischaemic stroke with large vessel occlusions. However, perfusion cannot always be restored after such an intervention. It has been hypothesised that the absence of reperfusion is at least partially due to the clot fragments that block the downstream vessels. In this paper, we present a new way of quantifying the effects of cerebral microthrombi on oxygen transport to tissue in terms of hypoxia and ischaemia. The oxygen transport was simulated with the Green's function method on physiologically representative microvascular cubes, which was found independent of both microvascular geometry and length scale. The microthrombi occlusions were then simulated in the microvasculature, which were extravasated over time with a new thrombus extravasation model. The tissue hypoxic fraction was fitted as a sigmoidal function of vessel blockage fraction, which was then taken to be a function of time after the formation of microthrombi occlusions. A novel hypoxia-based 3-state cell death model was finally proposed to simulate the hypoxic tissue damage over time. Using the cell death model, the impact of a certain degree of microthrombi occlusions on tissue viability and microinfarct volume can be predicted over time. Quantifying the impact of microthrombi on oxygen transport and tissue death will play an important role in full brain models of ischaemic stroke and thrombectomy.

Functionalizing Collagen with Vessel-Penetrating Two-Photon Phosphorescence Probes: A New In Vivo Strategy to Map Oxygen Concentration in Tumor Microenvironment and Tissue Ischemia

The encapsulation and/or surface modification can stabilize and protect the phosphorescence bio-probes but impede their intravenous delivery across biological barriers. Here, a new class of biocompatible rhenium (ReI ) diimine carbonyl complexes is developed, which can efficaciously permeate normal vessel walls and then functionalize the extravascular collagen matrixes as in situ oxygen sensor. Without protective agents, ReI -diimine complex already exhibits excellent emission yield (34%, λem   = 583 nm) and large two-photon absorption cross-sections (σ2   = 300 GM @ 800 nm) in water (pH 7.4). After extravasation, remarkably, the collagen-bound probes further enhanced their excitation efficiency by increasing the deoxygenated lifetime from 4.0 to 7.5 µs, paving a way to visualize tumor hypoxia and tissue ischemia in vivo. The post-extravasation functionalization of extracellular matrixes demonstrates a new methodology for biomaterial-empowered phosphorescence sensing and imaging.

The Ageing Brain: Investigating the Role of Age in Changes to the Human Cerebral Microvasculature With an in silico Model

Ageing causes extensive structural changes to the human cerebral microvasculature, which have a significant effect on capillary bed perfusion and oxygen transport. Current models of brain capillary networks in the literature focus on healthy adult brains and do not capture the effects of ageing, which is critical when studying neurodegenerative diseases. This study builds upon a statistically accurate model of the human cerebral microvasculature based on ex-vivo morphological data. This model is adapted for “healthy” ageing using in-vivo measurements from mice at three distinct age groups—young, middle-aged, and old. From this new model, blood and molecular exchange parameters are calculated such as permeability and surface-area-to-volume ratio, and compared across the three age groups. The ability to alter the model vessel-by-vessel is used to create a continuous gradient of ageing. It was found that surface-area-to-volume ratio reduced in old age by 6% and permeability by 24% from middle-age to old age, and variability within the networks also increased with age. The ageing gradient indicated a threshold in the ageing process around 75 years old, after which small changes have an amplified effect on blood flow properties. This gradient enables comparison of studies measuring cerebral properties at discrete points in time. The response of middle aged and old aged capillary beds to micro-emboli showed a lower robustness of the old age capillary bed to vessel occlusion. As the brain ages, there is thus increased vulnerability of the microvasculature—with a “tipping point” beyond which further remodeling of the microvasculature has exaggerated effects on the brain. When developing in-silico models of the brain, age is a very important consideration to accurately assess risk factors for cognitive decline and isolate early biomarkers of microvascular health.

Reduced global cerebral oxygen metabolic rate in sickle cell disease and chronic anemias

Anemia is the most common blood disorder in the world. In patients with chronic anemia, such as sickle cell disease or major thalassemia, cerebral blood flow increases to compensate for decreased oxygen content. However, the effects of chronic anemia on oxygen extraction fraction (OEF) and cerebral metabolic rate of oxygen (CMRO₂ ) are less well understood. In this study, we examined 47 sickle-cell anemia subjects (age 21.7 ± 7.1, female 45%), 27 non-sickle anemic subjects (age 25.0 ± 10.4, female 52%) and 44 healthy controls (age 26.4 ± 10.6, female 71%) using MRI metrics of brain oxygenation and flow. Phase contrast MRI was used to measure resting cerebral blood flow, while T₂ -relaxation-under-spin-tagging (TRUST) MRI with disease appropriate calibrations were used to measure OEF and CMRO₂ . We observed that patients with sickle cell disease and other chronic anemias have decreased OEF and CMRO₂ (respectively 27.4 ± 4.1% and 3.39 ± 0.71 ml O₂ /100 g/min in sickle cell disease, 30.8 ± 5.2% and 3.53 ± 0.64 ml O₂ /100 g/min in other anemias) compared to controls (36.7 ± 6.0% and 4.00 ± 0.65 ml O₂ /100 g/min). Impaired CMRO₂ was proportional to the degree of anemia severity. We further demonstrate striking concordance of the present work with pooled historical data from patients having broad etiologies for their anemia. The reduced cerebral oxygen extraction and metabolism are consistent with emerging data demonstrating increased non-nutritive flow, or physiological shunting, in sickle cell disease patients.

Neurovascular coupling and oxygenation are decreased in hippocampus compared to neocortex because of microvascular differences

The hippocampus is essential for spatial and episodic memory but is damaged early in Alzheimer's disease and is very sensitive to hypoxia. Understanding how it regulates its oxygen supply is therefore key for designing interventions to preserve its function. However, studies of neurovascular function in the hippocampus in vivo have been limited by its relative inaccessibility. Here we compared hippocampal and visual cortical neurovascular function in awake mice, using two photon imaging of individual neurons and vessels and measures of regional blood flow and haemoglobin oxygenation. We show that blood flow, blood oxygenation and neurovascular coupling were decreased in the hippocampus compared to neocortex, because of differences in both the vascular network and pericyte and endothelial cell function. Modelling oxygen diffusion indicates that these features of the hippocampal vasculature may restrict oxygen availability and could explain its sensitivity to damage during neurological conditions, including Alzheimer's disease, where the brain's energy supply is decreased.

The severity of microstrokes depends on local vascular topology and baseline perfusion

Cortical microinfarcts are linked to pathologies like cerebral amyloid angiopathy and dementia. Despite their relevance for disease progression, microinfarcts often remain undetected and the smallest scale of blood flow disturbance has not yet been identified. We employed blood flow simulations in realistic microvascular networks from the mouse cortex to quantify the impact of single-capillary occlusions. Our simulations reveal that the severity of a microstroke is strongly affected by the local vascular topology and the baseline flow rate in the occluded capillary. The largest changes in perfusion are observed in capillaries with two inflows and two outflows. This specific topological configuration only occurs with a frequency of 8%. The majority of capillaries have one inflow and one outflow and is likely designed to efficiently supply oxygen and nutrients. Taken together, microstrokes bear potential to induce a cascade of local disturbances in the surrounding tissue, which might accumulate and impair energy supply locally.

Effects of Genotypes and Treatment on Oxygenscan Parameters in Sickle Cell Disease

(1) Background: The aim of the present study was to compare oxygen gradient ektacytometry parameters between sickle cell patients of different genotypes (SS, SC, and S/β⁺) or under different treatments (hydroxyurea or chronic red blood cell exchange). (2) Methods: Oxygen gradient ektacytometry was performed in 167 adults and children at steady state. In addition, five SS patients had oxygenscan measurements at steady state and during an acute complication requiring hospitalization. (3) Results: Red blood cell (RBC) deformability upon deoxygenation (EImin) and in normoxia (EImax) was increased, and the susceptibility of RBC to sickle upon deoxygenation was decreased in SC patients when compared to untreated SS patients older than 5 years old. SS patients under chronic red blood cell exchange had higher EImin and EImax and lower susceptibility of RBC to sickle upon deoxygenation compared to untreated SS patients, SS patients younger than 5 years old, and hydroxyurea-treated SS and SC patients. The susceptibility of RBC to sickle upon deoxygenation was increased in the five SS patients during acute complication compared to steady state, although the difference between steady state and acute complication was variable from one patient to another. (4) Conclusions: The present study demonstrates that oxygen gradient ektacytometry parameters are affected by sickle cell disease (SCD) genotype and treatment.

Simulation of oxygen transport and estimation of tissue perfusion in extensive microvascular networks: Application to cerebral cortex

Advanced imaging techniques have made available extensive three-dimensional microvascular network structures. Simulation of oxygen transport by such networks requires information on blood flow rates and oxygen levels in vessels crossing boundaries of the imaged region, which is difficult to obtain experimentally. Here, a computational method is presented for estimating blood flow rates, oxygen levels, tissue perfusion and oxygen extraction, based on incomplete boundary conditions. Flow rates in all segments are estimated using a previously published method. Vessels crossing the region boundary are classified as arterioles, capillaries or venules. Oxygen levels in inflowing capillaries are assigned based on values in outflowing capillaries, and similarly for venules. Convective and diffusive oxygen transport is simulated. Contributions of each vessel to perfusion are computed in proportion to the decline in oxygen concentration along that vessel. For a vascular network in the mouse cerebral cortex, predicted tissue oxygen levels show a broad distribution, with 99% of tissue in the range of 20 to 80 mmHg under reference conditions, and steep gradients near arterioles. Perfusion and extraction estimates are consistent with experimental values. A 30% reduction in perfusion or a 30% increase in oxygen demand, relative to reference levels, is predicted to result in tissue hypoxia.

Laplacian Flow Dynamics on Geometric Graphs for Anatomical Modeling of Cerebrovascular Networks

Generating computational anatomical models of cerebrovascular networks is vital for improving clinical practice and understanding brain oxygen transport. This is achieved by extracting graph-based representations based on pre-mapping of vascular structures. Recent graphing methods can provide smooth vessels trajectories and well-connected vascular topology. However, they require water-tight surface meshes as inputs. Furthermore, adding vessels radii information on their graph compartments restricts their alignment along vascular centerlines. Here, we propose a novel graphing scheme that works with relaxed input requirements and intrinsically captures vessel radii information. The proposed approach is based on deforming geometric graphs constructed within vascular boundaries. Under a laplacian optimization framework, we assign affinity weights on the initial geometry that drives its iterative contraction toward vessels centerlines. We present a mechanism to decimate graph structure at each run and a convergence criterion to stop the process. A refinement technique is then introduced to obtain final vascular models. Our implementation is available on https://github.com/Damseh/VascularGraph. We benchmarked our results with that obtained using other efficient and state-of-the-art graphing schemes, validating on both synthetic and real angiograms acquired with different imaging modalities. The experiments indicate that the proposed scheme produces the lowest geometric and topological error rates on various angiograms. Furthermore, it surpasses other techniques in providing representative models that capture all anatomical aspects of vascular structures.

Sural Nerve Perfusion in Mice

Peripheral nerve function is metabolically demanding and nerve energy failure has been implicated in the onset and development of diabetic peripheral neuropathy and neuropathic pain conditions. Distal peripheral nerve oxygen supply relies on the distribution of red blood cells (RBCs) in just a few, nearby capillary-sized vessels and is therefore technically challenging to characterize. We developed an approach to characterize distal sural nerve hemodynamics in anesthetized, adult male mice using in vivo two-photon laser scanning microscopy. Our results show that RBC velocities in mouse sural nerve vessels are higher than those previously measured in mouse brain, and are sensitive to hindlimb temperatures. Nerve blood flow, measured as RBC flux, however, was similar to that of mouse brain and unaffected by local temperature. Power spectral density analysis of fluctuations in RBC velocities over short time intervals suggest that the technique is sufficiently sensitive and robust to detect subtle flow oscillations over time scales from 0.1 to tens of seconds. We conclude that in vivo two-photon laser scanning microscopy provides a suitable approach to study peripheral nerve hemodynamics in mice, and that local temperature control is important during such measurements.

Three-Photon Spectroscopy of Porphyrins

Recent advances in laser technology have made three-photon (3P) microscopy a real possibility, raising interest in the phenomenon of 3P absorption (3PA). Understanding 3PA of organic chromophores is especially important in view of those imaging applications that rely on exogenous probes, whose optical properties can be manipulated and optimized. Here, we present measurements and theoretical analysis of the degenerate 3PA spectra of several phosphorescent metalloporphyrins, which are used in the construction of biological oxygen probes. The effective 3PA cross sections (σ⁽³⁾) of these porphyrins near 1700 nm, a new promising biological optical window, were found to be on the order of 1000 GM3 (1 GM3 = 10^-83 cm⁶ s²), therefore being among the highest values reported to date for organic chromophores. To interpret our data, we developed a qualitative four-state model specific for porphyrins and used it in conjunction with quantitative analysis based on the time-dependent density functional theory (TDDFT)/a posteriori Tamm-Dancoff approximation (ATDA)/sum-over-states (SOS) formalism. The analysis revealed that B (Soret) state plays a key role in the enhancement of 3PA of porphyrins in the Q band region, while the low-lying two-photon (2P)-allowed gerade states interfere negatively and diminish the 3PA strength. This study features the first systematic examination of 3PA properties of porphyrins, suggesting ways to improve their performance and optimize them for imaging and other biomedical applications.

High spatiotemporal vessel-specific hemodynamic mapping with multi-echo single-vessel fMRI

High-resolution fMRI enables noninvasive mapping of the hemodynamic responses from individual penetrating vessels in animal brains. Here, a 2D multi-echo single-vessel fMRI (MESV-fMRI) method has been developed to map the fMRI signal from arterioles and venules with a 100 ms sampling rate at multiple echo times (TE, 3-30 ms) and short acquisition windows (<1 ms). The T₂*-weighted signal shows the increased extravascular effect on venule voxels as a function of TE. In contrast, the arteriole voxels show an increased fMRI signal with earlier onset than venules voxels at the short TE (3 ms) with increased blood inflow and volume effects. MESV-fMRI enables vessel-specific T2* mapping and presents T2*-based fMRI time courses with higher contrast-to-noise ratios (CNRs) than the T2*-weighted fMRI signal at a given TE. The vessel-specific T2* mapping also allows semi-quantitative estimation of the oxygen saturation levels (Y) and their changes (ΔY) at a given blood volume fraction upon neuronal activation. The MESV-fMRI method enables vessel-specific T2* measurements with high spatiotemporal resolution for better modeling of the fMRI signal based on the hemodynamic parameters.

The mass transfer coefficient for oxygen transport from blood to tissue in cerebral cortex

The functioning of cerebral cortex depends on adequate tissue oxygenation. MRI-based techniques allow estimation of blood oxygen levels, tissue perfusion, and oxygen consumption rate (CMRO₂), but do not directly measure partial pressure of oxygen (PO₂) in tissue. To address the estimation of tissue PO₂, the oxygen mass transfer coefficient (KTO₂) is here defined as the CMRO₂ divided by the difference in spatially averaged PO₂ between blood and tissue, and is estimated by analyzing Krogh-cylinder type models. Resistance to radial diffusion of oxygen from microvessels to tissue is distributed within vessels and in the extravascular tissue. The value of KTO₂ is shown to depend strongly on vascular length density and also on microvessel tube hematocrits and diameters, but to be insensitive to blood flow rate and to transient changes in flow or oxygen consumption. Estimated values of KTO₂ are higher than implied by previous studies, implying smaller declines in PO₂ from blood to tissue. Average tissue PO₂ can be estimated from MRI-based measurements as average blood PO₂ minus the product of KTO₂ and CMRO₂. For oxygen consumption rates and vascular densities typical of mouse cortex, the predicted difference between average blood and tissue PO₂ is about 10 mmHg.

Deep tissue volumetric optoacoustic tracking of individual circulating tumor cells in an intracardially perfused mouse model

Widespread metastasis is the major cause of death from melanoma and other types of cancer. At present, the dynamic aspects of the metastatic cascade remain enigmatic. The feasibility to track circulating melanoma cells deep within living intact organisms can greatly impact our knowledge on tumor metastasis, but existing imaging approaches lack the sensitivity, spatio-temporal resolution or penetration depth to capture flowing tumor cells over large fields of view within optically-opaque biological tissues. Vast progress with the development of optoacoustic tomography technologies has recently enabled two- and three-dimensional imaging at unprecedented frame rates in the order of hundreds of Hertz, effectively mapping up to a million image voxels within a single volumetric snapshot. Herein, we employ volumetric optoacoustic tomography for real-time visualization of passage and trapping of individual B16 melanoma cells in the whole mouse brain. Detection of individual circulating melanoma cells was facilitated by substituting blood with an artificial cerebrospinal fluid that removes the strong absorption background in the optoacoustic images. The approach can provide new opportunities for studying trafficking and accumulation of metastatic melanoma cells in different organs.

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.

Automated Analysis of Brain Microvasculature: From Segmentation to Anatomical Modeling

Two-photon microscopy (TPM) can provide a detailed microscopic information of cerebrovascular structures. Extracting anatomical vascular models from TPM angiograms remains a tedious task due to image degeneration associated with TPM acquisitions and the complexity of microvascular networks. Here, we propose a fully automated pipeline capable of providing useful anatomical models of vascular structures captured with TPM. In the proposed method, we segment blood vessels using a fully convolutional neural network and employ the resulting binary labels to create an initial geometric graph enclosed within vessels boundaries. The initial geometry is then decimated and refined to form graphed curve skeletons that can retain both the vascular shape and its topology. We validate the proposed method on 3D realistic TPM angiographies and compare our results with that obtained through manual annotations.

Retinal capillary oximetry with visible light optical coherence tomography

Assessing oxygen saturation (sO₂) remains challenging but is nonetheless necessary for understanding retinal metabolism. We and others previously achieved oximetry on major retinal vessels and measured the total retinal oxygen metabolic rate in rats using visible-light optical coherence tomography. Here we extend oximetry measurements to capillaries and investigate all three retinal vascular plexuses by amplifying and extracting the spectroscopic signal from each capillary segment under the guidance of optical coherence tomography (OCT) angiography. Using this approach, we measured capillary sO₂ in the retinal circulation in rats, demonstrated reproducibility of the results, validated the measurements in superficial capillaries with known perfusion pathways, and determined sO₂ responses to hypoxia and hyperoxia in the different retinal capillary beds. OCT capillary oximetry has the potential to provide new insights into the retinal circulation in the normal eye as well as in retinal vascular diseases.

Longitudinal optical coherence tomography imaging of tissue repair and microvasculature regeneration and function after targeted cerebral ischemia

SIGNIFICANCE

Understanding how the brain recovers from cerebral tissue and vascular damage after an ischemic event can help develop new therapeutic strategies for the treatment of stroke.

AIM

We investigated cerebral tissue repair and microvasculature regeneration and function after a targeted ischemic stroke.

APPROACH

Following photothrombosis occlusion of microvasculature, chronic optical coherence tomography (OCT)-based angiography was used to track ischemic tissue repair and microvasculature regeneration at three different cortical depths and up to 28 days in awake animals. Capillary network orientation analysis was performed to study the structural pattern of newly formed microvasculature. Based on the time-resolved OCT-angiography, we also investigated capillary stalling, which is likely related to ischemic stroke-induced inflammation.

RESULTS

Deeper cerebral tissue was found to have a larger ischemic area than shallower regions at any time point during the course of poststroke recovery, which suggests that cerebral tissue located deep in the cortex is more vulnerable. Regenerated microvasculature had a highly organized pattern at all cortical depths with a higher degree of structural reorganization in deeper regions. Additionally, capillary stalling event analysis revealed that cerebral ischemia augmented stalling events considerably.

CONCLUSION

Longitudinal OCT angiography reveals that regenerated capillary network has a highly directional pattern and an increased density and incidence of capillary stalling event.

Acoustic-spectrum-compensated photoacoustic microscopy

Photoacoustic microscopy (PAM) can label-free image oxy- and deoxy-hemoglobin (${{\rm HbO}_2}$HbO₂ and Hb) concentrations in vivo, providing useful information for metabolic researches and diagnostic applications. Conventional PAM assumes a linear relationship between the photoacoustic amplitude and the absorption coefficient. However, many factors, including absorber size, laser pulse width, and frequency response of the ultrasound transducer, may affect the measured acoustic spectrum and the shape of the temporal photoacoustic signal. The ultrasound transducer may weigh the blood vessels differently according to their diameters. In addition, the pulse width also affects the photoacoustic signal amplitude. These factors may cause inaccurate measurement of Hb and ${{\rm HbO}_2}$HbO₂ concentrations. To address this issue, we develop an acoustic-spectrum-compensated optical-resolution PAM (OR-PAM) that corrects the nonuniform acoustic spectrum and makes the quantitative results to be independent of the vessel diameter and pulse width. In dual-wavelength OR-PAM, we demonstrate that the acoustic spectrum compensation can improve the accuracy of oxygen saturation imaging by $\sim{15}\% $∼15%.

NIR-II fluorescence microscopic imaging of cortical vasculature in non-human primates

Vasculature architecture of the brain can provide revealing information about mental and neurological function and disease. Fluorescence imaging in the second near-infrared (NIR-II) regime with less light scattering is a more promising method for detecting cortical vessels than traditional visible and NIR-I modes. Methods: Clinically approved dye indocyanine green (ICG) was used for NIR-II fluorescence imaging. Here, for the first time, we developed two NIR-II fluorescence microscopy systems for brain vasculature imaging in macaque monkeys. The first is a wide-field microscope with high temporal resolution for measuring blood flow velocity and cardiac impulse period, while the second is a high spatial resolution confocal microscope producing three-dimensional maps of the cortical microvascular network. Both were designed with flexibility to image various cortical locations on the head. Results: Here, ICG was proved to have high brightness in NIR-II region and an 8-fold QY increase in serum than in water. We achieved cerebrovascular functional imaging of monkey with high temporal resolution (25 frames/second) with wide-field microscope. The blood flow velocity of capillaries can be precisely calculated and the cardiac impulse period can be monitored as well. In vivo structural imaging of cerebrovasculature was accomplished with both high spatial lateral resolution (~8 µm) and high signal to background ratio (SBR). Vivid 3D reconstructed NIR-II fluorescence confocal microscopic images up to depth of 470 μm were also realized. Conclusion: This work comprises an important advance towards studies of neurovascular coupling, stroke, and other diseases relevant to neurovascular health in humans.