Vascular density and distribution in neocortex

Distal activity patterns shape the spatial specificity of neurovascular coupling

Neurovascular coupling links brain activity to local changes in blood flow, forming the basis for non-invasive brain mapping. Using multiscale imaging, we investigated how vascular activity spatially relates to neuronal activity elicited by single whiskers across different columns and layers of mouse cortex. Here we show that mesoscopic hemodynamic signals quantitatively reflect neuronal activity across space but are composed of a highly heterogeneous pattern of responses across individual vessel segments that is poorly predicted by local neuronal activity. Rather, this heterogeneity is dependent on vessel directionality, specifically in thalamocortical input layer 4, where capillaries respond preferentially to neuronal activity patterns along their downstream perfusion domain. Thus, capillaries fine-tune blood flow based on distant activity and encode laminar-specific activity patterns. These findings imply that vascular anatomy sets a resolution limit on functional imaging signals, where individual blood vessels inaccurately report neuronal activity in their immediate vicinity but, instead, integrate activity patterns along the vascular arbor.

Sensorless adaptive optics in the second near-infrared window for deep vascular imaging in vivo

We have experimentally validated the use of sensorless adaptive optics (AO) to enhance laser scanning confocal microscopy in the second near-infrared (NIR II) spectral range, termed as AO-NIR II confocal microscopy. This approach harnesses a NIR II fluorophore, excited by an 808 nm wavelength and emitting beyond 1000 nm, to visualize intricate structures in deep brain tissues with the intact skull. By leveraging the reduced scattering and aberrations in the NIR II spectrum, we successfully captured a three-dimensional (3D) vascular structure map extending 310 µm beneath the skull. AO typically boosts the fluorescence signal by approximately 2-3 times, leading to a superior contrast and diminished smearing effects. Consequently, small blood vessels at various depths can be clearly visualized, which might otherwise remain undetectable without AO corrections.

Mapping vascular network architecture in primate brain using ferumoxytol-weighted laminar MRI

Mapping the vascular organization of the brain is of great importance across various domains of basic neuroimaging research, diagnostic radiology, and neurology. However, the intricate task of precisely mapping vasculature across brain regions and cortical layers presents formidable challenges, resulting in a limited understanding of neurometabolic factors influencing the brain's microvasculature. Addressing this gap, our study investigates whole-brain vascular volume using ferumoxytol-weighted laminar-resolution multi-echo gradient-echo imaging in macaque monkeys. We validate the results with published data for vascular densities and compare them with cytoarchitecture, neuron and synaptic densities. The ferumoxytol-induced change in transverse relaxation rate (ΔR2*), an indirect proxy measure of cerebral blood volume (CBV), was mapped onto twelve equivolumetric laminar cortical surfaces. Our findings reveal that CBV varies 3-fold across the brain, with the highest vascular volume observed in the inferior colliculus and lowest in the corpus callosum. In the cerebral cortex, CBV is notably high in early primary sensory areas and low in association areas responsible for higher cognitive functions. Classification of CBV into distinct groups unveils extensive replication of translaminar vascular network motifs, suggesting distinct computational energy supply requirements in areas with varying cytoarchitecture types. Regionally, baseline R2* and CBV exhibit positive correlations with neuron density and negative correlations with receptor densities. Adjusting image resolution based on the critical sampling frequency of penetrating cortical vessels, allows us to delineate approximately 30% of the arterial-venous vessels. Collectively, these results mark significant methodological and conceptual advancements, contributing to the refinement of cerebrovascular MRI. Furthermore, our study establishes a linkage between neurometabolic factors and the vascular network architecture in the primate brain.

The contribution of the vascular architecture and cerebrovascular reactivity to the BOLD signal formation across cortical depth

Assessment of neuronal activity using blood oxygenation level-dependent (BOLD) is confounded by how the cerebrovascular architecture modulates hemodynamic responses. To understand brain function at the laminar level, it is crucial to distinguish neuronal signal contributions from those determined by the cortical vascular organization. Therefore, our aim was to investigate the purely vascular contribution in the BOLD signal by using vasoactive stimuli and compare that with neuronal-induced BOLD responses from a visual task. To do so, we estimated the hemodynamic response function (HRF) across cortical depth following brief visual stimulations under different conditions using ultrahigh-field (7 Tesla) functional (f)MRI. We acquired gradient-echo (GE)-echo-planar-imaging (EPI) BOLD, containing contributions from all vessel sizes, and spin-echo (SE)-EPI BOLD for which signal changes predominately originate from microvessels, to distinguish signal weighting from different vascular compartments. Non-neuronal hemodynamic changes were induced by hypercapnia and hyperoxia to estimate cerebrovascular reactivity and venous cerebral blood volume ( C B V v O 2 ). Results show that increases in GE HRF amplitude from deeper to superficial layers coincided with increased macrovascular C B V v O 2 . C B V v O 2 -normalized GE-HRF amplitudes yielded similar cortical depth profiles as SE, thereby possibly improving specificity to neuronal activation. For GE BOLD, faster onset time and shorter time-to-peak were observed toward the deeper layers. Hypercapnia reduced the amplitude of visual stimulus-induced signal responses as denoted by lower GE-HRF amplitudes and longer time-to-peak. In contrast, the SE-HRF amplitude was unaffected by hypercapnia, suggesting that these responses reflect predominantly neurovascular processes that are less contaminated by macrovascular signal contributions.

On the influence of the vascular architecture on Gradient Echo and Spin Echo BOLD fMRI signals across cortical depth: a simulation approach based on realistic 3D vascular networks

GE-BOLD contrast stands out as the predominant technique in functional MRI experiments for its high sensitivity and straightforward implementation. GE-BOLD exhibits rather similar sensitivity to vessels independent of their size at submillimeter resolution studies like those examining cortical columns and laminae. However, the presence of nonspecific macrovascular contributions poses a challenge to accurately isolate neuronal activity. SE-BOLD increases specificity towards small vessels, thereby enhancing its specificity to neuronal activity, due to the effective suppression of extravascular contributions caused by macrovessels with its refocusing pulse. However, even SE-BOLD measurements may not completely remove these macrovascular contributions. By simulating hemodynamic signals across cortical depth, we gain insights into vascular contributions to the laminar BOLD signal. In this study, we employed four realistic 3D vascular models to simulate oxygen saturation states in various vascular compartments, aiming to characterize both intravascular and extravascular contributions to GE and SE signals, and corresponding BOLD signal changes, across cortical depth at 7T. Simulations suggest that SE-BOLD cannot completely reduce the macrovascular contribution near the pial surface. Simulations also show that both the specificity and signal amplitude of BOLD signals at 7T depend on the spatial arrangement of large vessels throughout cortical depth and on the pial surface.

A fully synthetic three-dimensional human cerebrovascular model based on histological characteristics to investigate the hemodynamic fingerprint of the layer BOLD fMRI signal formation

Recent advances in functional magnetic resonance imaging (fMRI) at ultra-high field (≥7 tesla), novel hardware, and data analysis methods have enabled detailed research on neurovascular function, such as cortical layer-specific activity, in both human and nonhuman species. A widely used fMRI technique relies on the blood oxygen level-dependent (BOLD) signal. BOLD fMRI offers insights into brain function by measuring local changes in cerebral blood volume, cerebral blood flow, and oxygen metabolism induced by increased neuronal activity. Despite its potential, interpreting BOLD fMRI data is challenging as it is only an indirect measurement of neuronal activity. Computational modeling can help interpret BOLD data by simulating the BOLD signal formation. Current developments have focused on realistic 3D vascular models based on rodent data to understand the spatial and temporal BOLD characteristics. While such rodent-based vascular models highlight the impact of the angioarchitecture on the BOLD signal amplitude, anatomical differences between the rodent and human vasculature necessitate the development of human-specific models. Therefore, a computational framework integrating human cortical vasculature, hemodynamic changes, and biophysical properties is essential. Here, we present a novel computational approach: a three-dimensional VAscular MOdel based on Statistics (3D VAMOS), enabling the investigation of the hemodynamic fingerprint of the BOLD signal within a model encompassing a fully synthetic human 3D cortical vasculature and hemodynamics. Our algorithm generates microvascular and macrovascular architectures based on morphological and topological features from the literature on human cortical vasculature. By simulating specific oxygen saturation states and biophysical interactions, our framework characterizes the intravascular and extravascular signal contributions across cortical depth and voxel-wise levels for gradient-echo and spin-echo readouts. Thereby, the 3D VAMOS computational framework demonstrates that using human characteristics significantly affects the BOLD fingerprint, making it an essential step in understanding the fundamental underpinnings of layer-specific fMRI experiments.

Oxygen and the Spark of Human Brain Evolution: Complex Interactions of Metabolism and Cortical Expansion across Development and Evolution

Scientific theories on the functioning and dysfunction of the human brain require an understanding of its development-before and after birth and through maturation to adulthood-and its evolution. Here we bring together several accounts of human brain evolution by focusing on the central role of oxygen and brain metabolism. We argue that evolutionary expansion of human transmodal association cortices exceeded the capacity of oxygen delivery by the vascular system, which led these brain tissues to rely on nonoxidative glycolysis for additional energy supply. We draw a link between the resulting lower oxygen tension and its effect on cytoarchitecture, which we posit as a key driver of genetic developmental programs for the human brain-favoring lower intracortical myelination and the presence of biosynthetic materials for synapse turnover. Across biological and temporal scales, this protracted capacity for neural plasticity sets the conditions for cognitive flexibility and ongoing learning, supporting complex group dynamics and intergenerational learning that in turn enabled improved nutrition to fuel the metabolic costs of further cortical expansion. Our proposed model delineates explicit mechanistic links among metabolism, molecular and cellular brain heterogeneity, and behavior, which may lead toward a clearer understanding of brain development and its disorders.

Morphofunctional and histological patterns of blood vessels in the superficial cervicofacial musculoaponeurotic system in midlateral face regions

OBJECTIVE

The superficial cervicofacial musculoaponeurotic system (SMAS) is a complex network formed by mimic muscles and conjunctive tissue of the superficial fascia of the face.This study aimed to introduce new anatomofunctional data on the importance of the trans-SMAS distribution pattern of the skin microperfusion of the face and to underline the role of SMAS in maintaining the homeostasis of the vascular network that crosses it. Considering the fibrous and muscular matrix of the SMAS, using COLIII and MyoH2 antibodies, together with endothelial immunohistochemistry(IHC)intercellular adhesion molecule 2 marker, we determined the correlation of these structures and their interaction.

METHODS

This study included 33donors of SMAS tissues, which have been stained withregular hematoxylin and eosin (HE), and three different IHC markers have been used (collagen III, muscular tissue, and blood vessels). The samples were collected from parotid, masseteric, jugal, and zygomatic regions. Magnetic resonance angiography was used to identify the main vascular sources of the midlateral regions of the face of another 47 patients.

RESULTS

Significant differences in topographic arrangement, density, and relations of the microsopic vasculature were observed between each of the four regions. Major differences were identified between the role of SMAS in each of these regions, from the parotid capsule to masseteric fascia, transition mobile part, and attaching manners in the zygomatic subunit.

CONCLUSIONS

Blood vessel topography must be related with the surrounding conjunctive and muscular tissue, especially regarding facial SMAS. Intrinsic relations between these three components of the SMAS and nervous fibers can provide us important hints on the functionality of the whole system.

A Closer Look at the Perivascular Unit in the Development of Enlarged Perivascular Spaces in Obesity, Metabolic Syndrome, and Type 2 Diabetes Mellitus

The recently described perivascular unit (PVU) resides immediately adjacent to the true capillary neurovascular unit (NVU) in the postcapillary venule and contains the normal-benign perivascular spaces (PVS) and pathological enlarged perivascular spaces (EPVS). The PVS are important in that they have recently been identified to be the construct and the conduit responsible for the delivery of metabolic waste from the interstitial fluid to the ventricular cerebrospinal fluid for disposal into the systemic circulation, termed the glymphatic system. Importantly, the outermost boundary of the PVS is lined by protoplasmic perivascular astrocyte endfeet (pvACef) that communicate with regional neurons. As compared to the well-recognized and described neurovascular unit (NVU) and NVU coupling, the PVU is less well understood and remains an emerging concept. The primary focus of this narrative review is to compare the similarities and differences between these two units and discuss each of their structural and functional relationships and how they relate not only to brain homeostasis but also how they may relate to the development of multiple clinical neurological disease states and specifically how they may relate to obesity, metabolic syndrome, and type 2 diabetes mellitus. Additionally, the concept and importance of a perisynaptic astrocyte coupling to the neuronal synapses with pre- and postsynaptic neurons will also be considered as a perisynaptic unit to provide for the creation of the information transfer in the brain via synaptic transmission and brain homeostasis. Multiple electron microscopic images and illustrations will be utilized in order to help explain these complex units.

The impact of lipopolysaccharide on cerebrovascular function and cognition resulting from obesity-induced gut dysbiosis

Obesity is a worldwide epidemic coinciding with a concomitant increase in the incidence of neurodegenerative diseases, particularly dementia. Obesity is characterised by increased adiposity, chronic low-grade systemic inflammation, and oxidative stress, which promote endothelial dysfunction. Endothelial dysfunction reduces cerebrovascular function leading to reduced cerebral blood flow and, eventually, cognitive decline, thus predisposing to a neurodegenerative disease. Obesity is also characterised by gut dysbiosis and a subsequent increase in the lipopolysaccharide which increasingly activates toll-like receptor 4 (TLR4) and further promotes chronic low-grade systemic inflammation. This also disrupts the crosstalk within the gut-brain axis, thus influencing the functions of the central nervous system, including cognition. However, the mechanisms by which obesity-related increases in oxidative stress, inflammation and endothelial dysfunction are driven by, or associated with, increased systemic lipopolysaccharide leading to reduced cerebrovascular function and cognition, beyond normal ageing, have not been elucidated. Hence, this review examines how increased concentrations of lipopolysaccharide and the subsequent increased TLR4 activation observed in obesity exacerbate the development of obesity-induced reductions in cerebrovascular function and cognition.

Distance from main arteries influences microstructural and functional brain tissue characteristics

Given the substantial dependence of neurons on continuous supply of energy, the distribution of major cerebral arteries opens a question whether the distance from the main supply arteries constitutes a modulating factor for the microstructural and functional properties of brain tissue. To tackle this question, multimodal MRI acquisitions of 102 healthy volunteers over the full adult age span were utilised. Relaxation along a fictitious field in the rotating frame of rank n = 4 (RAFF4), adiabatic T1ρ, T2ρ,  and intracellular volume fraction (fICVF) derived from diffusion-weighted imaging were implemented to quantify microstructural (cellularity, myelin density, iron concentration) tissue characteristics and degree centrality and fractional amplitude of low-frequency fluctuations to probe for functional metrics. Inverse correlation of arterial distance with robust homogeneity was detected for T1ρ, T2ρ and RAFF4 for cortical grey matter and white matter, showing substantial complex microstructural differences between brain tissue close and farther from main arterial trunks. Albeit with wider variability, functional metrics pointed to increased connectivity and neuronal activity in areas farther from main arteries. Surprisingly, multiple of these microstructural and functional distance-based gradients diminished with higher age, pointing to uniformization of brain tissue with ageing. All in all, this pilot study provides a novel insight on brain regionalisation based on artery distance, which merits further investigation to validate its biological underpinnings.

Fetal head-down posture may explain the rapid brain evolution in humans and other primates: An interpretative review

Evolutionary cerebrovascular consequences of upside-down postural verticality of the anthropoid fetus have been largely overlooked in the literature. This working hypothesis-based report provides a literature interpretation from an aspect that the rapid evolution of the human brain has been promoted by fetal head-down position due to maternal upright and semi-upright posture. Habitual vertical torso posture is a feature not only of humans, but also of monkeys and non-human apes that spend considerable time in a sitting position. Consequently, the head-down position of the fetus may have caused physiological craniovascular hypertension that stimulated expansion of the intracranial vessels and acted as an epigenetic physiological stress, which enhanced neurogenesis and eventually, along with other selective pressures, led to the progressive growth of the anthropoid brain and its organization. This article collaterally opens a new insight into the conundrum of high cephalopelvic proportions (i.e., the tight fit between the pelvic birth canal and fetal head) in phylogenetically distant lineages of monkeys, lesser apes, and humans. Low cephalopelvic proportions in non-human great apes could be accounted for by their energetically efficient horizontal nest-sleeping and consequently by their larger body mass compared to monkeys and lesser apes that sleep upright. One can further hypothesize that brain size varies in anthropoids according to the degree of exposure of the fetus to postural verticality. The supporting evidence for this postulation includes a finding that in fossil hominins cerebral blood flow rate increased faster than brain volume. This testable hypothesis opens a perspective for research on fetal postural cerebral hemodynamics.

A mechanistic computational framework to investigate the hemodynamic fingerprint of the blood oxygenation level-dependent signal

Blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) is one of the most used imaging techniques to map brain activity or to obtain clinical information about human cortical vasculature, in both healthy and disease conditions. Nevertheless, BOLD fMRI is an indirect measurement of brain functioning triggered by neurovascular coupling. The origin of the BOLD signal is quite complex, and the signal formation thus depends, among other factors, on the topology of the cortical vasculature and the associated hemodynamic changes. To understand the hemodynamic evolution of the BOLD signal response in humans, it is beneficial to have a computational framework available that virtually resembles the human cortical vasculature, and simulates hemodynamic changes and corresponding MRI signal changes via interactions of intrinsic biophysical and magnetic properties of the tissues. To this end, we have developed a mechanistic computational framework that simulates the hemodynamic fingerprint of the BOLD signal based on a statistically defined, three-dimensional, vascular model that approaches the human cortical vascular architecture. The microvasculature is approximated through a Voronoi tessellation method and the macrovasculature is adapted from two-photon microscopy mice data. Using this computational framework, we simulated hemodynamic changes-cerebral blood flow, cerebral blood volume, and blood oxygen saturation-induced by virtual arterial dilation. Then we computed local magnetic field disturbances generated by the vascular topology and the corresponding blood oxygen saturation changes. This mechanistic computational framework also considers the intrinsic biophysical and magnetic properties of nearby tissue, such as water diffusion and relaxation properties, resulting in a dynamic BOLD signal response. The proposed mechanistic computational framework provides an integrated biophysical model that can offer better insights regarding the spatial and temporal properties of the BOLD signal changes.

Neuronal functional connectivity is impaired in a layer dependent manner near the chronically implanted microelectrodes

Objective

This study aims to reveal longitudinal changes in functional network connectivity within and across different brain structures near the chronically implanted microelectrode. While it is well established that the foreign-body response (FBR) contributes to the gradual decline of the signals recorded from brain implants over time, how does the FBR impact affect the functional stability of neural circuits near implanted Brain-Computer Interfaces (BCIs) remains unknown. This research aims to illuminate how the chronic FBR can alter local neural circuit function and the implications for BCI decoders.

Approach

This study utilized multisite Michigan-style microelectrodes that span all cortical layers and the hippocampal CA1 region to collect spontaneous and visually-evoked electrophysiological activity. Alterations in neuronal activity near the microelectrode were tested assessing cross-frequency synchronization of LFP and spike entrainment to LFP oscillatory activity throughout 16 weeks after microelectrode implantation.

Main Results

The study found that cortical layer 4, the input-receiving layer, maintained activity over the implantation time. However, layers 2/3 rapidly experienced severe impairment, leading to a loss of proper intralaminar connectivity in the downstream output layers 5/6. Furthermore, the impairment of interlaminar connectivity near the microelectrode was unidirectional, showing decreased connectivity from Layers 2/3 to Layers 5/6 but not the reverse direction. In the hippocampus, CA1 neurons gradually became unable to properly entrain to the surrounding LFP oscillations.

Significance

This study provides a detailed characterization of network connectivity dysfunction over long-term microelectrode implantation periods. This new knowledge could contribute to the development of targeted therapeutic strategies aimed at improving the health of the tissue surrounding brain implants and potentially inform engineering of adaptive decoders as the FBR progresses. Our study's understanding of the dynamic changes in the functional network over time opens the door to developing interventions for improving the long-term stability and performance of intracortical microelectrodes.

The role of leptomeningeal collaterals in redistributing blood flow during stroke

Leptomeningeal collaterals (LMCs) connect the main cerebral arteries and provide alternative pathways for blood flow during ischaemic stroke. This is beneficial for reducing infarct size and reperfusion success after treatment. However, a better understanding of how LMCs affect blood flow distribution is indispensable to improve therapeutic strategies. Here, we present a novel in silico approach that incorporates case-specific in vivo data into a computational model to simulate blood flow in large semi-realistic microvascular networks from two different mouse strains, characterised by having many and almost no LMCs between middle and anterior cerebral artery (MCA, ACA) territories. This framework is unique because our simulations are directly aligned with in vivo data. Moreover, it allows us to analyse perfusion characteristics quantitatively across all vessel types and for networks with no, few and many LMCs. We show that the occlusion of the MCA directly caused a redistribution of blood that was characterised by increased flow in LMCs. Interestingly, the improved perfusion of MCA-sided microvessels after dilating LMCs came at the cost of a reduced blood supply in other brain areas. This effect was enhanced in regions close to the watershed line and when the number of LMCs was increased. Additional dilations of surface and penetrating arteries after stroke improved perfusion across the entire vasculature and partially recovered flow in the obstructed region, especially in networks with many LMCs, which further underlines the role of LMCs during stroke.

The Venular Side of Cerebral Amyloid Angiopathy: Proof of Concept of a Neglected Issue

Small vessel diseases (SVD) is an umbrella term including several entities affecting small arteries, arterioles, capillaries, and venules in the brain. One of the most relevant and prevalent SVDs is cerebral amyloid angiopathy (CAA), whose pathological hallmark is the deposition of amyloid fragments in the walls of small cortical and leptomeningeal vessels. CAA frequently coexists with Alzheimer's Disease (AD), and both are associated with cerebrovascular events, cognitive impairment, and dementia. CAA and AD share pathophysiological, histopathological and neuroimaging issues. The venular involvement in both diseases has been neglected, although both animal models and human histopathological studies found a deposition of amyloid beta in cortical venules. This review aimed to summarize the available information about venular involvement in CAA, starting from the biological level with the putative pathomechanisms of cerebral damage, passing through the definition of the peculiar angioarchitecture of the human cortex with the functional organization and consequences of cortical arteriolar and venular occlusion, and ending to the hypothesized links between cortical venular involvement and the main neuroimaging markers of the disease.

Microvascular Specificity of Spin Echo BOLD fMRI: Impact of EPI Echo Train Length

A spatially specific fMRI acquisition requires specificity to the microvasculature that serves active neuronal sites. Macrovascular contributions will reduce the microvascular specificity but can be reduced by using spin echo (SE) sequences that use a π pulse to refocus static field inhomogeneities near large veins. The microvascular specificity of a SE-echo planar imaging (SE-EPI) scan depends on the echo train length (ETL)-duration, but the dependence is not well-characterized in humans at 7T. To determine how microvascular-specific SE-EPI BOLD is in humans at 7T, we developed a Monte Carlo voxel model that computes the signal of a proton ensemble residing in a vasculature subjected to a SE-EPI pulse sequence. We characterized the ETL-duration dependence of the microvascular specificity by simulating the BOLD signal as a function of ETL, the range adhering to experimentally realistic readouts. We performed a validation experiment for our simulation observations, in which we acquired a set of SE-EPI BOLD time series with varying ETL during a hyperoxic gas challenge. Both our simulations and measurements show an increase in macrovascular contamination as a function of ETL, with an increase of 30% according to our simulation and 60% according to our validation experiment between the shortest and longest ETL durations (23.1 - 49.7 ms). We conclude that the microvascular specificity decreases heavily with increasing ETL-durations. We recommend reducing the ETL-duration as much as possible to minimize macrovascular contamination in SE-EPI BOLD experiments. We additionally recommend scanning at high resolutions to minimize partial volume effects with CSF. CSF voxels show a large BOLD response, which can be attributed to both the presence of large veins (high blood volume) and molecular oxygen-induced T1-shortening (significant in a hyperoxia experiment). The magnified BOLD signal in a GM-CSF partial volume voxel reduces the desired microvascular specificity and, therefore, will hinder the interpretation of functional MRI activation patterns.

Spatially heterogeneous structure-function coupling in haemodynamic and electromagnetic brain networks

The relationship between structural and functional connectivity in the brain is a key question in connectomics. Here we quantify patterns of structure-function coupling across the neocortex, by comparing structural connectivity estimated using diffusion MRI with functional connectivity estimated using both neurophysiological (MEG-based) and haemodynamic (fMRI-based) recordings. We find that structure-function coupling is heterogeneous across brain regions and frequency bands. The link between structural and functional connectivity is generally stronger in multiple MEG frequency bands compared to resting state fMRI. Structure-function coupling is greater in slower and intermediate frequency bands compared to faster frequency bands. We also find that structure-function coupling systematically follows the archetypal sensorimotor-association hierarchy, as well as patterns of laminar differentiation, peaking in granular layer IV. Finally, structure-function coupling is better explained using structure-informed inter-regional communication metrics than using structural connectivity alone. Collectively, these results place neurophysiological and haemodynamic structure-function relationships in a common frame of reference and provide a starting point for a multi-modal understanding of structure-function coupling in the brain.

Stretch-activated ion channel TMEM63B associates with developmental and epileptic encephalopathies and progressive neurodegeneration

By converting physical forces into electrical signals or triggering intracellular cascades, stretch-activated ion channels allow the cell to respond to osmotic and mechanical stress. Knowledge of the pathophysiological mechanisms underlying associations of stretch-activated ion channels with human disease is limited. Here, we describe 17 unrelated individuals with severe early-onset developmental and epileptic encephalopathy (DEE), intellectual disability, and severe motor and cortical visual impairment associated with progressive neurodegenerative brain changes carrying ten distinct heterozygous variants of TMEM63B, encoding for a highly conserved stretch-activated ion channel. The variants occurred de novo in 16/17 individuals for whom parental DNA was available and either missense, including the recurrent p.Val44Met in 7/17 individuals, or in-frame, all affecting conserved residues located in transmembrane regions of the protein. In 12 individuals, hematological abnormalities co-occurred, such as macrocytosis and hemolysis, requiring blood transfusions in some. We modeled six variants (p.Val44Met, p.Arg433His, p.Thr481Asn, p.Gly580Ser, p.Arg660Thr, and p.Phe697Leu), each affecting a distinct transmembrane domain of the channel, in transfected Neuro2a cells and demonstrated inward leak cation currents across the mutated channel even in isotonic conditions, while the response to hypo-osmotic challenge was impaired, as were the Ca2+ transients generated under hypo-osmotic stimulation. Ectopic expression of the p.Val44Met and p.Gly580Cys variants in Drosophila resulted in early death. TMEM63B-associated DEE represents a recognizable clinicopathological entity in which altered cation conductivity results in a severe neurological phenotype with progressive brain damage and early-onset epilepsy associated with hematological abnormalities in most individuals.

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 change of attention network functions and physiological adaptation during high-altitude hypoxia and reoxygenation

Previous studies have not reached a definitive conclusion regarding the effect of high-altitude hypoxia and reoxygenation on attention. To clarify the influence of altitude and exposure time on attention and the relations between physiological activity and attention, we conducted a longitudinal study to track attention network functions in 26 college students. The scores on the attention network test and physiological data, including heart rate, percutaneous arterial oxygen saturation (SpO₂), blood pressure, and vital capacity in pulmonary function measurement, were collected at five time-points: two weeks before arriving at high altitude (baseline), within 3 days after arriving at high altitude (HA3), 21 days after arriving at high altitude (HA21), 7 days after returning to sea level (POST7) and 30 days after returning to sea level (POST30). The alerting scores at POST30 were significantly higher than those at baseline, HA3 and HA21; the orienting scores at HA3 were lower than those at POST7 and POST30; the executive control scores at POST7 were significantly lower than those at baseline, HA3, HA21, and POST30; and the executive control scores at HA3 were significantly higher than those at POST30. The change in SpO₂ during high-altitude acclimatization (from HA3 to HA21) was positively correlated with the orienting score at HA21. Vital capacity changes during acute deacclimatization positively correlated with orienting scores at POST7. Attention network functions at the behavioral level did not decline after acute hypoxia exposure compared with baseline. Attention network functions after returning to sea level were improved compared with those during acute hypoxia; additionally, alerting and executive function scores were improved compared with those at baseline. Thus, the speed of physiological adaptation could facilitate the recovery of orienting function during acclimatization and deacclimatization.

Single-microvessel occlusion produces lamina-specific microvascular flow vasodynamics and signs of neurodegenerative change

Recent studies have highlighted the importance of understanding the architecture and function of microvasculature, and dysfunction of these microvessels may underlie neurodegenerative disease. Here, we utilize a high-precision ultrafast laser-induced photothrombosis (PLP) method to occlude single capillaries and then quantitatively study the effects on vasodynamics and surrounding neurons. Analysis of the microvascular architecture and hemodynamics after single-capillary occlusion reveals distinct changes upstream vs. downstream branches, which shows rapid regional flow redistribution and local downstream blood-brain barrier (BBB) leakage. Focal ischemia via capillary occlusions surrounding labeled target neurons induces dramatic and rapid lamina-specific changes in neuronal dendritic architecture. Further, we find that micro-occlusion at two different depths within the same vascular arbor results in distinct effects on flow profiles in layers 2/3 vs layer 4. The current results reveal laminar-scale regulation distinctions in microinfarct response and raise the possibility that relatively greater impacts on microvascular function contribute to cognitive decline in neurodegenerative disease.

Capsaicin: A Potential Treatment to Improve Cerebrovascular Function and Cognition in Obesity and Ageing

Impaired cognition is the primary symptom of dementia, which can lead to functional disability and reduced quality of life among an increasingly ageing population. Ageing is associated with increased oxidative stress, chronic low-grade systemic inflammation, and endothelial dysfunction, which reduces cerebrovascular function leading to cognitive decline. Chronic low-grade systemic inflammatory conditions, such as obesity, exacerbate this decline beyond normal ageing and predispose individuals to neurodegenerative diseases, such as dementia. Capsaicin, the major pungent molecule of chilli, has recently demonstrated improvements in cognition in animal models via activation of the transient receptor potential vanilloid channel 1 (TRPV1). Capsaicin-induced TRPV1 activation reduces adiposity, chronic low-grade systemic inflammation, and oxidative stress, as well as improves endothelial function, all of which are associated with cerebrovascular function and cognition. This review examines the current literature on capsaicin and Capsimax, a capsaicin supplement associated with reduced gastrointestinal irritation compared to capsaicin. Acute and chronic capsaicin treatment can improve cognition in animals. However, studies adequately assessing the effects of capsaicin on cerebrovascular function, and cognition in humans do not exist. Capsimax may be a potentially safe therapeutic intervention for future clinical trials testing the effects of capsaicin on cerebrovascular function and cognition.

Estimates of the permeability of extra-cellular pathways through the astrocyte endfoot sheath

BACKGROUND

Astrocyte endfoot processes are believed to cover all micro-vessels in the brain cortex and may play a significant role in fluid and substance transport into and out of the brain parenchyma. Detailed fluid mechanical models of diffusive and advective transport in the brain are promising tools to investigate theories of transport.

METHODS

We derive theoretical estimates of astrocyte endfoot sheath permeability for advective and diffusive transport and its variation in microvascular networks from mouse brain cortex. The networks are based on recently published experimental data and generated endfoot patterns are based on Voronoi tessellations of the perivascular surface. We estimate corrections for projection errors in previously published data.

RESULTS

We provide structural-functional relationships between vessel radius and resistance that can be directly used in flow and transport simulations. We estimate endfoot sheath filtration coefficients in the range [Formula: see text] to [Formula: see text], diffusion membrane coefficients for small solutes in the range [Formula: see text] to [Formula: see text], and gap area fractions in the range 0.2-0.6%, based on a inter-endfoot gap width of 20 nm.

CONCLUSIONS

The astrocyte endfoot sheath surrounding microvessels forms a secondary barrier to extra-cellular transport, separating the extra-cellular space of the parenchyma and the perivascular space outside the endothelial layer. The filtration and membrane diffusion coefficients of the endfoot sheath are estimated to be an order of magnitude lower than those of the extra-cellular matrix while being two orders of magnitude higher than those of the vessel wall.

Endothelial Cells are Heterogeneous in Different Brain Regions and are Dramatically Altered in Alzheimer's Disease

Vascular endothelial cells play an important role in maintaining brain health, but their contribution to Alzheimer's disease (AD) is obscured by limited understanding of the cellular heterogeneity in normal aged brain and in disease. To address this, we performed single nucleus RNAseq on tissue from 32 AD and non-AD donors each with five cortical regions: entorhinal cortex, inferior temporal gyrus, prefrontal cortex, visual association cortex and primary visual cortex. Analysis of 51,586 endothelial cells revealed unique gene expression patterns across the five regions in non-AD donors. Alzheimer's brain endothelial cells were characterized by upregulated protein folding genes and distinct transcriptomic differences in response to amyloid beta plaques and cerebral amyloid angiopathy (CAA). This dataset demonstrates previously unrecognized regional heterogeneity in the endothelial cell transcriptome in both aged non-AD and AD brain.

Significance Statement

In this work, we show that vascular endothelial cells collected from five different brain regions display surprising variability in gene expression. In the presence of Alzheimer's disease pathology, endothelial cell gene expression is dramatically altered with clear differences in regional and temporal changes. These findings help explain why certain brain regions appear to differ in susceptibility to disease-related vascular remodeling events that may impact blood flow.

A unifying model for discordant and concordant results in human neuroimaging studies of facial viewpoint selectivity

Our ability to recognize faces regardless of viewpoint is a key property of the primate visual system. Traditional theories hold that facial viewpoint is represented by view-selective mechanisms at early visual processing stages and that representations become increasingly tolerant to viewpoint changes in higher-level visual areas. Newer theories, based on single-neuron monkey electrophysiological recordings, suggest an additional intermediate processing stage invariant to mirror-symmetric face views. Consistent with traditional theories, human studies combining neuroimaging and multivariate pattern analysis (MVPA) methods have provided evidence of view-selectivity in early visual cortex. However, contradictory results have been reported in higher-level visual areas concerning the existence in humans of mirror-symmetrically tuned representations. We believe these results reflect low-level stimulus confounds and data analysis choices. To probe for low-level confounds, we analyzed images from two popular face databases. Analyses of mean image luminance and contrast revealed biases across face views described by even polynomials-i.e., mirror-symmetric. To explain major trends across human neuroimaging studies of viewpoint selectivity, we constructed a network model that incorporates three biological constraints: cortical magnification, convergent feedforward projections, and interhemispheric connections. Given the identified low-level biases, we show that a gradual increase of interhemispheric connections across network layers is sufficient to replicate findings of mirror-symmetry in high-level processing stages, as well as view-tuning in early processing stages. Data analysis decisions-pattern dissimilarity measure and data recentering-accounted for the variable observation of mirror-symmetry in late processing stages. The model provides a unifying explanation of MVPA studies of viewpoint selectivity. We also show how common analysis choices can lead to erroneous conclusions.

Quantitative MRI maps of human neocortex explored using cell type-specific gene expression analysis

Quantitative magnetic resonance imaging (qMRI) allows extraction of reproducible and robust parameter maps. However, the connection to underlying biological substrates remains murky, especially in the complex, densely packed cortex. We investigated associations in human neocortex between qMRI parameters and neocortical cell types by comparing the spatial distribution of the qMRI parameters longitudinal relaxation rate (${R_{1}}$), effective transverse relaxation rate (${R_{2}}^{\ast }$), and magnetization transfer saturation (MTsat) to gene expression from the Allen Human Brain Atlas, then combining this with lists of genes enriched in specific cell types found in the human brain. As qMRI parameters are magnetic field strength-dependent, the analysis was performed on MRI data at 3T and 7T. All qMRI parameters significantly covaried with genes enriched in GABA- and glutamatergic neurons, i.e. they were associated with cytoarchitecture. The qMRI parameters also significantly covaried with the distribution of genes enriched in astrocytes (${R_{2}}^{\ast }$ at 3T, ${R_{1}}$ at 7T), endothelial cells (${R_{1}}$ and MTsat at 3T), microglia (${R_{1}}$ and MTsat at 3T, ${R_{1}}$ at 7T), and oligodendrocytes and oligodendrocyte precursor cells (${R_{1}}$ at 7T). These results advance the potential use of qMRI parameters as biomarkers for specific cell types.

Ultraliser: a framework for creating multiscale, high-fidelity and geometrically realistic 3D models for in silico neuroscience

Ultraliser is a neuroscience-specific software framework capable of creating accurate and biologically realistic 3D models of complex neuroscientific structures at intracellular (e.g. mitochondria and endoplasmic reticula), cellular (e.g. neurons and glia) and even multicellular scales of resolution (e.g. cerebral vasculature and minicolumns). Resulting models are exported as triangulated surface meshes and annotated volumes for multiple applications in in silico neuroscience, allowing scalable supercomputer simulations that can unravel intricate cellular structure-function relationships. Ultraliser implements a high-performance and unconditionally robust voxelization engine adapted to create optimized watertight surface meshes and annotated voxel grids from arbitrary non-watertight triangular soups, digitized morphological skeletons or binary volumetric masks. The framework represents a major leap forward in simulation-based neuroscience, making it possible to employ high-resolution 3D structural models for quantification of surface areas and volumes, which are of the utmost importance for cellular and system simulations. The power of Ultraliser is demonstrated with several use cases in which hundreds of models are created for potential application in diverse types of simulations. Ultraliser is publicly released under the GNU GPL3 license on GitHub (BlueBrain/Ultraliser).

SIGNIFICANCE

There is crystal clear evidence on the impact of cell shape on its signaling mechanisms. Structural models can therefore be insightful to realize the function; the more realistic the structure can be, the further we get insights into the function. Creating realistic structural models from existing ones is challenging, particularly when needed for detailed subcellular simulations. We present Ultraliser, a neuroscience-dedicated framework capable of building these structural models with realistic and detailed cellular geometries that can be used for simulations.

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.

Laminar-specific functional connectivity mapping with multi-slice line-scanning fMRI

Despite extensive studies detecting laminar functional magnetic resonance imaging (fMRI) signals to illustrate the canonical microcircuit, the spatiotemporal characteristics of laminar-specific information flow across cortical regions remain to be fully investigated in both evoked and resting conditions at different brain states. Here, we developed a multislice line-scanning fMRI (MS-LS) method to detect laminar fMRI signals in adjacent cortical regions with high spatial (50 μm) and temporal resolution (100 ms) in anesthetized rats. Across different trials, we detected either laminar-specific positive or negative blood-oxygen-level-dependent (BOLD) responses in the surrounding cortical region adjacent to the most activated cortex under the evoked condition. Specifically, in contrast to typical Layer (L) 4 correlation across different regions due to the thalamocortical projections for trials with positive BOLD, a strong correlation pattern specific in L2/3 was detected for trials with negative BOLD in adjacent regions, which indicated brain state-dependent laminar-fMRI responses based on corticocortical interaction. Also, in resting-state (rs-) fMRI study, robust lag time differences in L2/3, 4, and 5 across multiple cortices represented the low-frequency rs-fMRI signal propagation from caudal to rostral slices. In summary, our study provided a unique laminar fMRI mapping scheme to better characterize trial-specific intra- and inter-laminar functional connectivity in evoked and resting-state MS-LS.

Modelling the depth-dependent VASO and BOLD responses in human primary visual cortex

Functional magnetic resonance imaging (fMRI) using a blood-oxygenation-level-dependent (BOLD) contrast is a common method for studying human brain function noninvasively. Gradient-echo (GRE) BOLD is highly sensitive to the blood oxygenation change in blood vessels; however, the spatial signal specificity can be degraded due to signal leakage from activated lower layers to superficial layers in depth-dependent (also called laminar or layer-specific) fMRI. Alternatively, physiological variables such as cerebral blood volume using the VAscular-Space-Occupancy (VASO) contrast have shown higher spatial specificity compared to BOLD. To better understand the physiological mechanisms such as blood volume and oxygenation changes and to interpret the measured depth-dependent responses, models are needed which reflect vascular properties at this scale. For this purpose, we extended and modified the "cortical vascular model" previously developed to predict layer-specific BOLD signal changes in human primary visual cortex to also predict a layer-specific VASO response. To evaluate the model, we compared the predictions with experimental results of simultaneous VASO and BOLD measurements in a group of healthy participants. Fitting the model to our experimental data provided an estimate of CBV change in different vascular compartments upon neural activity. We found that stimulus-evoked CBV change mainly occurs in small arterioles, capillaries, and intracortical arteries and that the contribution from venules and ICVs is smaller. Our results confirm that VASO is less susceptible to large vessel effects compared to BOLD, as blood volume changes in intracortical arteries did not substantially affect the resulting depth-dependent VASO profiles, whereas depth-dependent BOLD profiles showed a bias towards signal contributions from intracortical veins.

Human-in-the-Loop Optimization of Transcranial Electrical Stimulation at the Point of Care: A Computational Perspective

Individual differences in the responsiveness of the brain to transcranial electrical stimulation (tES) are increasingly demonstrated by the large variability in the effects of tES. Anatomically detailed computational brain models have been developed to address this variability; however, static brain models are not “realistic” in accounting for the dynamic state of the brain. Therefore, human-in-the-loop optimization at the point of care is proposed in this perspective article based on systems analysis of the neurovascular effects of tES. First, modal analysis was conducted using a physiologically detailed neurovascular model that found stable modes in the 0 Hz to 0.05 Hz range for the pathway for vessel response through the smooth muscle cells, measured with functional near-infrared spectroscopy (fNIRS). During tES, the transient sensations can have arousal effects on the hemodynamics, so we present a healthy case series for black-box modeling of fNIRS−pupillometry of short-duration tDCS effects. The block exogeneity test rejected the claim that tDCS is not a one-step Granger cause of the fNIRS total hemoglobin changes (HbT) and pupil dilation changes (p < 0.05). Moreover, grey-box modeling using fNIRS of the tDCS effects in chronic stroke showed the HbT response to be significantly different (paired-samples t-test, p < 0.05) between the ipsilesional and contralesional hemispheres for primary motor cortex tDCS and cerebellar tDCS, which was subserved by the smooth muscle cells. Here, our opinion is that various physiological pathways subserving the effects of tES can lead to state−trait variability, which can be challenging for clinical translation. Therefore, we conducted a case study on human-in-the-loop optimization using our reduced-dimensions model and a stochastic, derivative-free covariance matrix adaptation evolution strategy. We conclude from our computational analysis that human-in-the-loop optimization of the effects of tES at the point of care merits investigation in future studies for reducing inter-subject and intra-subject variability in neuromodulation.

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic-haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

Fiber-based lactate recordings with fluorescence resonance energy transfer sensors by applying an magnetic resonance-informed correction of hemodynamic artifacts

Significance: Fluorescence resonance energy transfer (FRET) sensors offer enormous benefits when studying neurophysiology through confocal microscopy. Yet, their use for fiber-based in vivo recordings is hampered by massive confounding effects and has therefore been scarcely reported. Aim: We aim to investigate whether in vivo fiber-based lactate recordings in the rodent brain are feasible with FRET sensors and implement a correction algorithm for the predominant hemodynamic artifact. Approach: We performed fiber-based FRET recordings of lactate (Laconic) and calcium (Twitch-2B) simultaneously with functional MRI and pharmacological MRI. MR-derived parameters were applied to correct hemodynamic artifacts. Results of FRET measurements were validated by local field potential, magnetic resonance spectroscopy, and blood analysis. Results: Hemodynamic artifacts dominated fiber-based in vivo FRET measurements with both Laconic and Twitch-2B. Our MR-based correction algorithm enabled to remove the artifacts and detect lactate and calcium changes during sensory stimulation or intravenous lactate injections. Conclusions: In vivo fiber-based lactate recordings are feasible using FRET-based sensors. However, signal corrections are required. MR-derived hemodynamic parameters can successfully be applied for artifact correction.

Toward an integrative neurovascular framework for studying brain networks

Brain functional connectivity based on the measure of blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signals has become one of the most widely used measurements in human neuroimaging. However, the nature of the functional networks revealed by BOLD fMRI can be ambiguous, as highlighted by a recent series of experiments that have suggested that typical resting-state networks can be replicated from purely vascular or physiologically driven BOLD signals. After going through a brief review of the key concepts of brain network analysis, we explore how the vascular and neuronal systems interact to give rise to the brain functional networks measured with BOLD fMRI. This leads us to emphasize a view of the vascular network not only as a confounding element in fMRI but also as a functionally relevant system that is entangled with the neuronal network. To study the vascular and neuronal underpinnings of BOLD functional connectivity, we consider a combination of methodological avenues based on multiscale and multimodal optical imaging in mice, used in combination with computational models that allow the integration of vascular information to explain functional connectivity.

Mathematical modelling of haemorrhagic transformation within a multi-scale microvasculature network

Objective Haemorrhagic transformation (HT) is one of the most common complications after ischaemic stroke caused by damage to the blood-brain barrier (BBB) that could be the result of stroke progression or a complication of stroke treatment with reperfusion therapy. The aim of this study is to develop further a previous simple HT mathematical model into an enlarged multi-scale microvasculature model in order to investigate the effects of HT on the surrounding tissue and vasculature. In addition, this study investigates the relationship between tissue displacement and vascular geometry. Approach By modelling tissue displacement, capillary compression, hydraulic conductivity in tissue and vascular permeability, we establish a mathematical model to describe the change of intracranial pressure (ICP) surrounding the damaged vascular bed after HT onset applied to a 3D multi-scale microvasculature. The use of a voxel-scale model then enables us to compare our HT simulation with available clinical imaging data for perfusion and cerebral blood volume (CBV) in the multi-scale microvasculature network. Main results We showed that the haematoma diameter and the maximum tissue displacement are approximately proportional to the diameter of the breakdown vessel. Based on the voxel-scale model, we found that perfusion reduces by approximately 13-17 % and CBV reduces by around 20-25 % after HT onset due to the effect of capillary compression caused by increased interstitial pressure. The results are in good agreement with the limited experimental data. Significance This model, by enabling us to bridge the gap between the microvascular scale and clinically measurable parameters, thus provides a foundation for more detailed validation and understanding of HT in patients.

Imaging of the pial arterial vasculature of the human brain in vivo using high-resolution 7T time-of-flight angiography

The pial arterial vasculature of the human brain is the only blood supply to the neocortex, but quantitative data on the morphology and topology of these mesoscopic arteries (diameter 50–300 µm) remains scarce. Because it is commonly assumed that blood flow velocities in these vessels are prohibitively slow, non-invasive time-of-flight magnetic resonance angiography (TOF-MRA)—which is well suited to high 3D imaging resolutions—has not been applied to imaging the pial arteries. Here, we provide a theoretical framework that outlines how TOF-MRA can visualize small pial arteries in vivo, by employing extremely small voxels at the size of individual vessels. We then provide evidence for this theory by imaging the pial arteries at 140 µm isotropic resolution using a 7 Tesla (T) magnetic resonance imaging (MRI) scanner and prospective motion correction, and show that pial arteries one voxel width in diameter can be detected. We conclude that imaging pial arteries is not limited by slow blood flow, but instead by achievable image resolution. This study represents the first targeted, comprehensive account of imaging pial arteries in vivo in the human brain. This ultra-high-resolution angiography will enable the characterization of pial vascular anatomy across the brain to investigate patterns of blood supply and relationships between vascular and functional architecture.

A unified 3D map of microscopic architecture and MRI of the human brain

We present the first three-dimensional (3D) concordance maps of cyto- and fiber architecture of the human brain, combining histology, immunohistochemistry, and 7-T quantitative magnetic resonance imaging (MRI), in two individual specimens. These 3D maps each integrate data from approximately 800 microscopy sections per brain, showing neuronal and glial cell bodies, nerve fibers, and interneuronal populations, as well as ultrahigh-field quantitative MRI, all coaligned at the 200-μm scale to the stacked blockface images obtained during sectioning. These unprecedented 3D multimodal datasets are shared without any restrictions and provide a unique resource for the joint study of cell and fiber architecture of the brain, detailed anatomical atlasing, or modeling of the microscopic underpinnings of MRI contrasts.

Vascularizing the brain in vitro

The brain is arguably the most fascinating and complex organ in the human body. Recreating the brain in vitro is an ambition restricted by our limited understanding of its structure and interacting elements. One of these interacting parts, the brain microvasculature, is distinguished by a highly selective barrier known as the blood-brain barrier (BBB), limiting the transport of substances between the blood and the nervous system. Numerous in vitro models have been used to mimic the BBB and constructed by implementing a variety of microfabrication and microfluidic techniques. However, currently available models still cannot accurately imitate the in vivo characteristics of BBB. In this article, we review recent BBB models by analyzing each parameter affecting the accuracy of these models. Furthermore, we propose an investigation of the synergy between BBB models and neuronal tissue biofabrication, which results in more advanced models, including neurovascular unit microfluidic models and vascularized brain organoid-based models.

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.

Seeing the Forest and Its Trees Together: Implementing 3D Light Microscopy Pipelines for Cell Type Mapping in the Mouse Brain

The brain is composed of diverse neuronal and non-neuronal cell types with complex regional connectivity patterns that create the anatomical infrastructure underlying cognition. Remarkable advances in neuroscience techniques enable labeling and imaging of these individual cell types and their interactions throughout intact mammalian brains at a cellular resolution allowing neuroscientists to examine microscopic details in macroscopic brain circuits. Nevertheless, implementing these tools is fraught with many technical and analytical challenges with a need for high-level data analysis. Here we review key technical considerations for implementing a brain mapping pipeline using the mouse brain as a primary model system. Specifically, we provide practical details for choosing methods including cell type specific labeling, sample preparation (e.g., tissue clearing), microscopy modalities, image processing, and data analysis (e.g., image registration to standard atlases). We also highlight the need to develop better 3D atlases with standardized anatomical labels and nomenclature across species and developmental time points to extend the mapping to other species including humans and to facilitate data sharing, confederation, and integrative analysis. In summary, this review provides key elements and currently available resources to consider while developing and implementing high-resolution mapping methods.

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.

Digital Reconstruction of the Neuro-Glia-Vascular Architecture

Astrocytes connect the vasculature to neurons mediating the supply of nutrients and biochemicals. They are involved in a growing number of physiological and pathophysiological processes that result from biophysical, physiological, and molecular interactions in this neuro-glia-vascular ensemble (NGV). The lack of a detailed cytoarchitecture severely restricts the understanding of how they support brain function. To address this problem, we used data from multiple sources to create a data-driven digital reconstruction of the NGV at micrometer anatomical resolution. We reconstructed 0.2 mm3 of the rat somatosensory cortex with 16 000 morphologically detailed neurons, 2500 protoplasmic astrocytes, and its microvasculature. The consistency of the reconstruction with a wide array of experimental measurements allows novel predictions of the NGV organization, allowing the anatomical reconstruction of overlapping astrocytic microdomains and the quantification of endfeet connecting each astrocyte to the vasculature, as well as the extent to which they cover the latter. Structural analysis showed that astrocytes optimize their positions to provide uniform vascular coverage for trophic support and signaling. However, this optimal organization rapidly declines as their density increases. The NGV digital reconstruction is a resource that will enable a better understanding of the anatomical principles and geometric constraints, which govern how astrocytes support brain function.

Grey-box modeling and hypothesis testing of functional near-infrared spectroscopy-based cerebrovascular reactivity to anodal high-definition tDCS in healthy humans

Transcranial direct current stimulation (tDCS) has been shown to evoke hemodynamics response; however, the mechanisms have not been investigated systematically using systems biology approaches. Our study presents a grey-box linear model that was developed from a physiologically detailed multi-compartmental neurovascular unit model consisting of the vascular smooth muscle, perivascular space, synaptic space, and astrocyte glial cell. Then, model linearization was performed on the physiologically detailed nonlinear model to find appropriate complexity (Akaike information criterion) to fit functional near-infrared spectroscopy (fNIRS) based measure of blood volume changes, called cerebrovascular reactivity (CVR), to high-definition (HD) tDCS. The grey-box linear model was applied on the fNIRS-based CVR during the first 150 seconds of anodal HD-tDCS in eleven healthy humans. The grey-box linear models for each of the four nested pathways starting from tDCS scalp current density that perturbed synaptic potassium released from active neurons for Pathway 1, astrocytic transmembrane current for Pathway 2, perivascular potassium concentration for Pathway 3, and voltage-gated ion channel current on the smooth muscle cell for Pathway 4 were fitted to the total hemoglobin concentration (tHb) changes from optodes in the vicinity of 4x1 HD-tDCS electrodes as well as on the contralateral sensorimotor cortex. We found that the tDCS perturbation Pathway 3 presented the least mean square error (MSE, median <2.5%) and the lowest Akaike information criterion (AIC, median -1.726) from the individual grey-box linear model fitting at the targeted-region. Then, minimal realization transfer function with reduced-order approximations of the grey-box model pathways was fitted to the ensemble average tHb time series. Again, Pathway 3 with nine poles and two zeros (all free parameters), provided the best Goodness of Fit of 0.0078 for Chi-Square difference test of nested pathways. Therefore, our study provided a systems biology approach to investigate the initial transient hemodynamic response to tDCS based on fNIRS tHb data. Future studies need to investigate the steady-state responses, including steady-state oscillations found to be driven by calcium dynamics, where transcranial alternating current stimulation may provide frequency-dependent physiological entrainment for system identification. We postulate that such a mechanistic understanding from system identification of the hemodynamics response to transcranial electrical stimulation can facilitate adequate delivery of the current density to the neurovascular tissue under simultaneous portable imaging in various cerebrovascular diseases.

Large-scale characterization of the microvascular geometry in development and disease by tissue clearing and quantitative ultramicroscopy

Three-dimensional assessment of optically cleared, entire organs and organisms has recently become possible by tissue clearing and selective plane illumination microscopy ("ultramicroscopy"). Resulting datasets can be highly complex, encompass over a thousand images with millions of objects and data of several gigabytes per acquisition. This constitutes a major challenge for quantitative analysis. We have developed post-processing tools to quantify millions of microvessels and their distribution in three-dimensional datasets from ultramicroscopy and demonstrate the capabilities of our pipeline within entire mouse brains and embryos. Using our developed acquisition, segmentation, and analysis platform, we quantify physiological vascular networks in development and the healthy brain. We compare various geometric vessel parameters (e.g. vessel density, radius, tortuosity) in the embryonic spinal cord and brain as well as in different brain regions (basal ganglia, corpus callosum, cortex). White matter tract structures (corpus callosum, spinal cord) showed lower microvascular branch densities and longer vessel branch length compared to grey matter (cortex, basal ganglia). Furthermore, we assess tumor neoangiogenesis in a mouse glioma model to compare tumor core and tumor border. The developed methodology allows rapid quantification of three-dimensional datasets by semi-automated segmentation of fluorescently labeled objects with conventional computer hardware. Our approach can aid preclinical investigations and paves the way towards "quantitative ultramicroscopy".

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.

Quantitative relations between BOLD responses, cortical energetics and impulse firing across cortical depth

The blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal arises as a consequence of changes in cerebral blood flow (CBF) and cerebral metabolic rate of oxygen ( CMR O 2 ) that in turn are modulated by changes in neural activity. Recent advances in imaging have achieved sub-millimetre resolution and allowed investigation of the BOLD response as a function of cortical depth. Here, we adapt our previous theory relating the BOLD signal to neural activity to produce a quantitative model that incorporates venous blood draining between cortical layers. The adjustable inputs to the model are the neural activity and a parameter governing this blood draining. A three-layer version for transient neural inputs and a multi-layer version for constant or tonic neural inputs are able to account for a variety of experimental results, including negative BOLD signals.

Brain microvasculature has a common topology with local differences in geometry that match metabolic load

The microvasculature underlies the supply networks that support neuronal activity within heterogeneous brain regions. What are common versus heterogeneous aspects of the connectivity, density, and orientation of capillary networks? To address this, we imaged, reconstructed, and analyzed the microvasculature connectome in whole adult mice brains with sub-micrometer resolution. Graph analysis revealed common network topology across the brain that leads to a shared structural robustness against the rarefaction of vessels. Geometrical analysis, based on anatomically accurate reconstructions, uncovered a scaling law that links length density, i.e., the length of vessel per volume, with tissue-to-vessel distances. We then derive a formula that connects regional differences in metabolism to differences in length density and, further, predicts a common value of maximum tissue oxygen tension across the brain. Last, the orientation of capillaries is weakly anisotropic with the exception of a few strongly anisotropic regions; this variation can impact the interpretation of fMRI data.

Modelling the leptomeningeal collateral circulation during acute ischaemic stroke

A novel model of the leptomeningeal collateral circulation is created by combining data from multiple sources with statistical scaling laws. The extent of the collateral circulation is varied by defining a collateral vessel probability. Blood flow and pressure are simulated using a one-dimensional steady state blood flow model. The leptomeningeal collateral vessels provide significant flow during a stroke. The pressure drop over an occlusion predicted by the model ranges between 60 and 85 mmHg depending on the extent of the collateral circulation. The linear transport of contrast material was simulated in the circulatory network. The time delay of peak contrast over an occlusion is 3.3 s in the model, and 2.1 s (IQR 0.8-4.0 s) when measured in dynamic CTA data of acute ischaemic stroke patients. Modelling the leptomeningeal collateral circulation could lead to better estimates of infarct volume and patient outcome.