Vascular graph model to simulate the cerebral blood flow in realistic vascular networks

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.

Cerebral perfusion simulation using realistically generated synthetic trees for healthy and stroke patients

BACKGROUND AND OBJECTIVE

Cerebral vascular diseases are among the most burdensome diseases faced by society. However, investigating the pathophysiology of diseases as well as developing future treatments still relies heavily on expensive in-vivo and in-vitro studies. The generation of realistic, patient-specific models of the cerebrovascular system capable of simulating hemodynamics and perfusion promises the ability to simulate diseased states, therefore accelerating development cycles using in silico studies and opening opportunities for the individual assessment of diseased states, treatment planning, and the prediction of outcomes. By providing a patient-specific, anatomically detailed and validated model of the human cerebral vascular system, we aim to provide the basis for future in silico investigations of the cerebral physiology and pathology.

METHODS

In this retrospective study, a processing pipeline for patient-specific quantification of cerebral perfusion was developed and applied to healthy individuals and a stroke patient. Major arteries are segmented from 3T MR angiography data. A synthetic tree generation algorithm titled tissue-growth based optimization (GBO)1 is used to extend vascular trees beyond the imaging resolution. To investigate the anatomical accuracy of the generated trees, morphological parameters are compared against those of 7 T MRI, 9.4 T MRI, and dissection data. Using the generated vessel model, hemodynamics and perfusion are simulated by solving one-dimensional blood flow equations combined with Darcy flow equations.

RESULTS

Morphological data of three healthy individuals (mean age 47 years ± 15.9 [SD], 2 female) was analyzed. Bifurcation and physiological characteristics of the synthetically generated vessels are comparable to those of dissection data. The inability of MRI based segmentation to resolve small branches and the small volume investigated cause a mismatch in the comparison to MRI data. Cerebral perfusion was estimated for healthy individuals and a stroke patient. The simulated perfusion is compared against Arterial-Spin-Labeling MRI perfusion data. Good qualitative agreement is found between simulated and measured cerebral blood flow (CBF)2. Ischemic regions are predicted well, however ischemia severity is overestimated.

CONCLUSIONS

GBO successfully generates detailed cerebral vascular models with realistic morphological parameters. Simulations based on the resulting networks predict perfusion territories and ischemic regions successfully.

Modeling of blood flow in cerebral arterial circulation and its dynamic impact on electrical conductivity in a realistic multi-compartment head model

BACKGROUND AND OBJECTIVE

This study aims to assess the dynamic impact of non-Newtonian cerebral arterial circulation on electrical conductivity within a realistic multi-compartment head model. Evaluating this research question is crucial and challenging due to its relevance to electrophysiological modalities like transcranial electrical stimulation (tES), electro-/magnetoencephalography (EEG/MEG), and electrical impedance tomography (EIT). In these modalities, accurate forward modeling depends on the electrical conductivity, which is affected by complex tortuous vessel networks, limited data acquisition in Magnetic Resonance Imaging (MRI), and non-linear blood flow phenomena, including shear rate and viscosity in non-Newtonian fluid.

METHODS

To obtain an approximation for the blood concentration, we first use Navier-Stokes equations (NSEs) to solve for the pressure and velocity of the blood in the major vessels. Then Fick's law is used to solve for the blood concentration in the tissues. Finally, Archie's law is used to estimate the electrical conductivity distribution based on the blood concentration.

RESULTS

The results, obtained with an open 7 Tesla MRI dataset, suggest that a dynamic model of cerebral blood flow (CBF) for both arterial and microcirculation can be established; we find blood pressure and electrical conductivity distributions given a numerically simulated pulse sequence and approximate the blood concentration and electrical conductivity inside the brain based on those.

CONCLUSIONS

Our model provides an approximation of the dynamical blood flow and the corresponding electrical conductivity distribution in the different parts of the brain. The advantage of our approach is that it is applicable with limited a priori information about the blood flow and with an arbitrary head model distinguishing the arteries.

Pressure-Poisson equation in numerical simulation of cerebral arterial circulation and its effect on the electrical conductivity of the brain

BACKGROUND AND OBJECTIVE

This study considers dynamic modeling of the cerebral arterial circulation and reconstructing an atlas for the electrical conductivity of the brain. Electrical conductivity is a governing parameter in several electrophysiological modalities applied in neuroscience, such as electroencephalography (EEG), transcranial electrical stimulation (tES), and electrical impedance tomography (EIT). While high-resolution 7-Tesla (T) Magnetic Resonance Imaging (MRI) data allow for reconstructing the cerebral arteries with a cross-sectional diameter larger than the voxel size, electrical conductivity cannot be directly inferred from MRI data. Brain models of electrophysiology typically associate each brain tissue compartment with a constant electrical conductivity, omitting any dynamic effects of cerebral blood circulation. Incorporating those effects poses the challenge of solving a system of incompressible Navier-Stokes equations (NSEs) in a realistic multi-compartment head model. However, using a simplified circulation model is well-motivated since, on the one hand, the complete system does not always have a numerically stable solution and, on the other hand, the full set of arteries cannot be perfectly reconstructed from the MRI data, meaning that any solution will be approximative.

METHODS

We postulate that circulation in the distinguishable arteries can be estimated via the pressure-Poisson equation (PPE), which is coupled with Fick's law of diffusion for microcirculation. To establish a fluid exchange model between arteries and microarteries, a boundary condition derived from the Hagen-Poisseuille model is applied. The relationship between the estimated volumetric blood concentration and the electrical conductivity of the brain tissue is approximated through Archie's law for fluid flow in porous media.

RESULTS

Through the formulation of the PPE and a set of boundary conditions (BCs) based on the Hagen-Poisseuille model, we obtained an equivalent formulation of the incompressible Stokes equation (SE). Thus, allowing effective blood pressure estimation in cerebral arteries segmented from open 7T MRI data.

CONCLUSIONS

As a result of this research, we developed and built a useful modeling framework that accounts for the effects of dynamic blood flow on a novel MRI-based electrical conductivity atlas. The electrical conductivity perturbation obtained in numerical experiments has an appropriate overall match with previous studies on this subject. Further research to validate these results will be necessary.

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.

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.

Graph theory applications in congenital heart disease

Graph theory can be used to address problems with complex network structures. Congenital heart diseases (CHDs) involve complex abnormal connections between chambers, vessels, and organs. We proposed a new method to represent CHDs based on graph theory, wherein vertices were defined as the spaces through which blood flows and edges were defined by the blood flow between the spaces and direction of the blood flow. The CHDs of tetralogy of Fallot (TOF) and transposition of the great arteries (TGA) were selected as examples for constructing directed graphs and binary adjacency matrices. Patients with totally repaired TOF, surgically corrected d-TGA, and Fontan circulation undergoing four-dimensional (4D) flow magnetic resonance imaging (MRI) were included as examples for constructing the weighted adjacency matrices. The directed graphs and binary adjacency matrices of the normal heart, extreme TOF undergoing a right modified Blalock-Taussig shunt, and d-TGA with a ventricular septal defect were constructed. The weighted adjacency matrix of totally repaired TOF was constructed using the peak velocities obtained from 4D flow MRI. The developed method is promising for representing CHDs and may be helpful in developing artificial intelligence and conducting future research on CHD.

An anatomically detailed arterial-venous network model. Cerebral and coronary circulation

In recent years, several works have addressed the problem of modeling blood flow phenomena in veins, as a response to increasing interest in modeling pathological conditions occurring in the venous network and their connection with the rest of the circulatory system. In this context, one-dimensional models have proven to be extremely efficient in delivering predictions in agreement with in-vivo observations. Pursuing the increase of anatomical accuracy and its connection to physiological principles in haemodynamics simulations, the main aim of this work is to describe a novel closed-loop Anatomically-Detailed Arterial-Venous Network (ADAVN) model. An extremely refined description of the arterial network consisting of 2,185 arterial vessels is coupled to a novel venous network featuring high level of anatomical detail in cerebral and coronary vascular territories. The entire venous network comprises 189 venous vessels, 79 of which drain the brain and 14 are coronary veins. Fundamental physiological mechanisms accounting for the interaction of brain blood flow with the cerebro-spinal fluid and of the coronary circulation with the cardiac mechanics are considered. Several issues related to the coupling of arterial and venous vessels at the microcirculation level are discussed in detail. Numerical simulations are compared to patient records published in the literature to show the descriptive capabilities of the model. Furthermore, a local sensitivity analysis is performed, evidencing the high impact of the venous circulation on main cardiovascular variables.

A network-based model of dynamic cerebral autoregulation

Cerebrovascular diseases continue to be one of the leading causes of morbidity and mortality in humans. Abnormalities in dynamic cerebral autoregulation (dCA) have been implicated in many of these disease conditions. Accurate models are therefore needed to better understand the complex pathophysiology behind impaired dCA. We thus present here a simple framework for modelling a vessel-driven network model of dCA in the microvasculature, as opposed to the conventional compartmental modelling approach. Network models incorporate the actual connectivity and anatomy of the vasculature, thereby allowing us to include and trace changes in the calibre and morphology of individual vessels, investigate the spatial specificity and heterogeneity of the various control mechanisms to help disentangle their contributions, and link the model parameters to the actual network physiology. The proposed control feedback mechanisms are incorporated at the level of the individual vessel, and the dynamic pressure and flow fields are solved for here within a simple vessel network. In response to an upstream pressure drop, the network is found to be able to recover cerebral blood flow (CBF) while exhibiting the characteristic autoregulatory behaviour in terms of changes in vessel calibre and the biphasic flow response. We assess the feasibility of our formulation in larger networks by comparing the simulation results to those obtained using a one-dimensional (1D) model of CBF applied to the same microvasculature network and find that our model results are in very good agreement with the 1D solution, while significantly reducing the computational cost, thus enabling more detailed models of network behaviour to be adopted in the future. Accurate and computationally feasible models of dCA that are more representative of the vasculature can help increase the translatability of haemodynamic models into the clinical environment, which would help develop more informed treatment guidelines for patients with cerebrovascular diseases.

Tumor Evolution Models of Phase-Field Type with Nonlocal Effects and Angiogenesis

In this survey article, a variety of systems modeling tumor growth are discussed. In accordance with the hallmarks of cancer, the described models incorporate the primary characteristics of cancer evolution. Specifically, we focus on diffusive interface models and follow the phase-field approach that describes the tumor as a collection of cells. Such systems are based on a multiphase approach that employs constitutive laws and balance laws for individual constituents. In mathematical oncology, numerous biological phenomena are involved, including temporal and spatial nonlocal effects, complex nonlinearities, stochasticity, and mixed-dimensional couplings. Using the models, for instance, we can express angiogenesis and cell-to-matrix adhesion effects. Finally, we offer some methods for numerically approximating the models and show simulations of the tumor's evolution in response to various biological effects.

A skeleton context-aware 3D fully convolutional network for abdominal artery segmentation

PURPOSE

This paper aims to propose a deep learning-based method for abdominal artery segmentation. Blood vessel structure information is essential to diagnosis and treatment. Accurate blood vessel segmentation is critical to preoperative planning. Although deep learning-based methods perform well on large organs, segmenting small organs such as blood vessels is challenging due to complicated branching structures and positions. We propose a 3D deep learning network from a skeleton context-aware perspective to improve segmentation accuracy. In addition, we propose a novel 3D patch generation method which could strengthen the structural diversity of a training data set.

METHOD

The proposed method segments abdominal arteries from an abdominal computed tomography (CT) volume using a 3D fully convolutional network (FCN). We add two auxiliary tasks to the network to extract the skeleton context of abdominal arteries. In addition, our skeleton-based patch generation (SBPG) method further enables the FCN to segment small arteries. SBPG generates a 3D patch from a CT volume by leveraging artery skeleton information. These methods improve the segmentation accuracies of small arteries.

RESULTS

We used 20 cases of abdominal CT volumes to evaluate the proposed method. The experimental results showed that our method outperformed previous segmentation accuracies. The averaged precision rate, recall rate, and F-measure were 95.5%, 91.0%, and 93.2%, respectively. Compared to a baseline method, our method improved 1.5% the averaged recall rate and 0.7% the averaged F-measure.

CONCLUSIONS

We present a skeleton context-aware 3D FCN to segment abdominal arteries from an abdominal CT volume. In addition, we propose a 3D patch generation method. Our fully automated method segmented most of the abdominal artery regions. The method produced competitive segmentation performance compared to previous methods.

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.

VascuViz: a multimodality and multiscale imaging and visualization pipeline for vascular systems biology

Despite advances in imaging, image-based vascular systems biology has remained challenging because blood vessel data are often available only from a single modality or at a given spatial scale, and cross-modality data are difficult to integrate. Therefore, there is an exigent need for a multimodality pipeline that enables ex vivo vascular imaging with magnetic resonance imaging, computed tomography and optical microscopy of the same sample, while permitting imaging with complementary contrast mechanisms from the whole-organ to endothelial cell spatial scales. To achieve this, we developed 'VascuViz'-an easy-to-use method for simultaneous three-dimensional imaging and visualization of the vascular microenvironment using magnetic resonance imaging, computed tomography and optical microscopy in the same intact, unsectioned tissue. The VascuViz workflow permits multimodal imaging with a single labeling step using commercial reagents and is compatible with diverse tissue types and protocols. VascuViz's interdisciplinary utility in conjunction with new data visualization approaches opens up new vistas in image-based vascular systems biology.

A nonlinear multi-scale model for blood circulation in a realistic vascular system

In the last decade, numerical models have become an increasingly important tool in biological and medical science. Numerical simulations contribute to a deeper understanding of physiology and are a powerful tool for better diagnostics and treatment. In this paper, a nonlinear multi-scale model framework is developed for blood flow distribution in the full vascular system of an organ. We couple a quasi one-dimensional vascular graph model to represent blood flow in larger vessels and a porous media model to describe flow in smaller vessels and capillary bed. The vascular model is based on Poiseuille's Law, with pressure correction by elasticity and pressure drop estimation at vessels' junctions. The porous capillary bed is modelled as a two-compartment domain (artery and venous) using Darcy's Law. The fluid exchange between the artery and venous capillary bed compartments is defined as blood perfusion. The numerical experiments show that the proposed model for blood circulation: (i) is closely dependent on the structure and parameters of both the larger vessels and of the capillary bed, and (ii) provides a realistic blood circulation in the organ. The advantage of the proposed model is that it is complex enough to reliably capture the main underlying physiological function, yet highly flexible as it offers the possibility of incorporating various local effects. Furthermore, the numerical implementation of the model is straightforward and allows for simulations on a regular desktop computer.

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.

Parallel generation of extensive vascular networks with application to an archetypal human kidney model

Given the relevance of the inextricable coupling between microcirculation and physiology, and the relation to organ function and disease progression, the construction of synthetic vascular networks for mathematical modelling and computer simulation is becoming an increasingly broad field of research. Building vascular networks that mimic in vivo morphometry is feasible through algorithms such as constrained constructive optimization (CCO) and variations. Nevertheless, these methods are limited by the maximum number of vessels to be generated due to the whole network update required at each vessel addition. In this work, we propose a CCO-based approach endowed with a domain decomposition strategy to concurrently create vascular networks. The performance of this approach is evaluated by analysing the agreement with the sequentially generated networks and studying the scalability when building vascular networks up to 200 000 vascular segments. Finally, we apply our method to vascularize a highly complex geometry corresponding to the cortex of a prototypical human kidney. The technique presented in this work enables the automatic generation of extensive vascular networks, removing the limitation from previous works. Thus, we can extend vascular networks (e.g. obtained from medical images) to pre-arteriolar level, yielding patient-specific whole-organ vascular models with an unprecedented level of detail.

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.

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.

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.

Adaptive constrained constructive optimisation for complex vascularisation processes

Mimicking angiogenetic processes in vascular territories acquires importance in the analysis of the multi-scale circulatory cascade and the coupling between blood flow and cell function. The present work extends, in several aspects, the Constrained Constructive Optimisation (CCO) algorithm to tackle complex automatic vascularisation tasks. The main extensions are based on the integration of adaptive optimisation criteria and multi-staged space-filling strategies which enhance the modelling capabilities of CCO for specific vascular architectures. Moreover, this vascular outgrowth can be performed either from scratch or from an existing network of vessels. Hence, the vascular territory is defined as a partition of vascular, avascular and carriage domains (the last one contains vessels but not terminals) allowing one to model complex vascular domains. In turn, the multi-staged space-filling approach allows one to delineate a sequence of biologically-inspired stages during the vascularisation process by exploiting different constraints, optimisation strategies and domain partitions stage by stage, improving the consistency with the architectural hierarchy observed in anatomical structures. With these features, the aDaptive CCO (DCCO) algorithm proposed here aims at improving the modelled network anatomy. The capabilities of the DCCO algorithm are assessed with a number of anatomically realistic scenarios.

Validating Linear Systems Analysis for Laminar fMRI: Temporal Additivity for Stimulus Duration Manipulations

Advancements in ultra-high field (7 T and higher) magnetic resonance imaging (MRI) scanners have made it possible to investigate both the structure and function of the human brain at a sub-millimeter scale. As neuronal feedforward and feedback information arrives in different layers, sub-millimeter functional MRI has the potential to uncover information processing between cortical micro-circuits across cortical depth, i.e. laminar fMRI. For nearly all conventional fMRI analyses, the main assumption is that the relationship between local neuronal activity and the blood oxygenation level dependent (BOLD) signal adheres to the principles of linear systems theory. For laminar fMRI, however, directional blood pooling across cortical depth stemming from the anatomy of the cortical vasculature, potentially violates these linear system assumptions, thereby complicating analysis and interpretation. Here we assess whether the temporal additivity requirement of linear systems theory holds for laminar fMRI. We measured responses elicited by viewing stimuli presented for different durations and evaluated how well the responses to shorter durations predicted those elicited by longer durations. We find that BOLD response predictions are consistently good predictors for observed responses, across all cortical depths, and in all measured visual field maps (V1, V2, and V3). Our results suggest that the temporal additivity assumption for linear systems theory holds for laminar fMRI. We thus show that the temporal additivity assumption holds across cortical depth for sub-millimeter gradient-echo BOLD fMRI in early visual cortex.

A Graph-Based Framework for Multiscale Modeling of Physiological Transport

Trillions of chemical reactions occur in the human body every second, where the generated products are not only consumed locally but also transported to various locations in a systematic manner to sustain homeostasis. Current solutions to model these biological phenomena are restricted in computability and scalability due to the use of continuum approaches in which it is practically impossible to encapsulate the complexity of the physiological processes occurring at diverse scales. Here, we present a discrete modeling framework defined on an interacting graph that offers the flexibility to model multiscale systems by translating the physical space into a metamodel. We discretize the graph-based metamodel into functional units composed of well-mixed volumes with vascular and cellular subdomains; the operators defined over these volumes define the transport dynamics. We predict glucose drift governed by advective-dispersive transport in the vascular subdomains of an islet vasculature and cross-validate the flow and concentration fields with finite-element-based COMSOL simulations. Vascular and cellular subdomains are coupled to model the nutrient exchange occurring in response to the gradient arising out of reaction and perfusion dynamics. The application of our framework for modeling biologically relevant test systems shows how our approach can assimilate both multi-omics data from in vitro-in vivo studies and vascular topology from imaging studies for examining the structure-function relationship of complex vasculatures. The framework can advance simulation of whole-body networks at user-defined levels and is expected to find major use in personalized medicine and drug discovery.

Pre-operative Assessment of Ablation Margins for Variable Blood Perfusion Metrics in a Magnetic Resonance Imaging Based Complex Breast Tumour Anatomy: Simulation Paradigms in Thermal Therapies

BACKGROUND AND OBJECTIVES

Image-guided medical interventions facilitates precise visualization at treatment site. The conformal prediction for sparing healthy tissue fringes precisely in the vicinity of irregular tumour anatomy remains clinically challenging. Pre-clinical image-based computational modelling is imperative as it helps in enhancement of treatment quality, augmenting clinical-decision making, while planning, targeting, controlling, monitoring and assessing treatment response with an effective risk assessment before the onset of treatment in clinical settings. In this study, the influence of heat deposition rate (SAR), exposure duration, and variable blood perfusion metrics for a patient-specific breast tumour is quantified considering the tumour margins thereby suggesting need of geometrically accurate models.

METHODS

A three-dimensional realistic model mimicking dimensions of a female breast, comprising ~1.7 cm irregular tumour, was generated from patient specific two-dimensional DICOM format MRI images through image segmentation tools MIMICS 19.0® and 3-Matic 11.0® which is finally exported to COMSOL Multiphysics 5.2® as a volumetric mesh for finite element analysis. The Pennes bioheat transfer model and Arrhenius thermal damage model of cell-death are integrated to simulate a coupled biophysics problem. A comparative blood perfusion analysis is done to evaluate the response of tumour during heating considering thermal damage extent, including the tumour margins while sparing critical adjoining healthy tissues.

RESULTS

The evaluated thermal damage zones for 1 mm, 2 mm and 3 mm fringe heating region (beyond tumour boundary) reveals 0.09%, 0.21% and 0.34% thermal damage to the healthy tissue (which is <1%) and thus successful necrosis of the tumour. The iterative computational experiments suggests treatment margins < 5 mm are sufficient enough as heating beyond 3 mm fringe layer leads to higher damage surrounding the tumour approximately 1.5 times the tumour volume. Further, the heat-dosage requirements are 22% more for highly perfused tumour as compared to moderately perfused tumour with an approximate double time to ablate the whole tumour volume.

CONCLUSIONS

Depending on the blood perfusion characteristics of a tumour, it is a trade-off between heat-dosage (SAR) and exposure/treatment duration to get desired thermal damage including the irregular tumour boundaries while taking into account, the margin of healthy tissue. The suggested patient-specific integrated multiphysics-model based on MRI-Images may be implemented for pre-treatment planning based on the tumour blood perfusion to evaluate the thermal ablation zone dimensions clinically and thereby avoiding the damage of off-target tissues. Thus, risks involving underestimation or overestimation of thermal coagulation zones may be minimised while preserving the surrounding normal breast parenchyma.

DeepVesselNet: Vessel Segmentation, Centerline Prediction, and Bifurcation Detection in 3-D Angiographic Volumes

We present DeepVesselNet, an architecture tailored to the challenges faced when extracting vessel trees and networks and corresponding features in 3-D angiographic volumes using deep learning. We discuss the problems of low execution speed and high memory requirements associated with full 3-D networks, high-class imbalance arising from the low percentage (<3%) of vessel voxels, and unavailability of accurately annotated 3-D training data-and offer solutions as the building blocks of DeepVesselNet. First, we formulate 2-D orthogonal cross-hair filters which make use of 3-D context information at a reduced computational burden. Second, we introduce a class balancing cross-entropy loss function with false-positive rate correction to handle the high-class imbalance and high false positive rate problems associated with existing loss functions. Finally, we generate a synthetic dataset using a computational angiogenesis model capable of simulating vascular tree growth under physiological constraints on local network structure and topology and use these data for transfer learning. We demonstrate the performance on a range of angiographic volumes at different spatial scales including clinical MRA data of the human brain, as well as CTA microscopy scans of the rat brain. Our results show that cross-hair filters achieve over 23% improvement in speed, lower memory footprint, lower network complexity which prevents overfitting and comparable accuracy that does not differ from full 3-D filters. Our class balancing metric is crucial for training the network, and transfer learning with synthetic data is an efficient, robust, and very generalizable approach leading to a network that excels in a variety of angiography segmentation tasks. We observe that sub-sampling and max pooling layers may lead to a drop in performance in tasks that involve voxel-sized structures. To this end, the DeepVesselNet architecture does not use any form of sub-sampling layer and works well for vessel segmentation, centerline prediction, and bifurcation detection. We make our synthetic training data publicly available, fostering future research, and serving as one of the first public datasets for brain vessel tree segmentation and analysis.

Predicting Vessel Diameter Changes to Up-Regulate Biphasic Blood Flow During Activation in Realistic Microvascular Networks

A dense network of blood vessels distributes blood to different regions of the brain. To meet the temporarily and spatially varying energy demand resulting from changes in neuronal activity, the vasculature is able to locally up-regulate the blood supply. However, to which extent diameter changes of different vessel types contribute to the up-regulation, as well as the spatial and temporal characteristics of their changes, are currently unknown. Here, we present a new simulation method, which solves an inverse problem to calculate diameter changes of individual blood vessels needed to achieve predefined blood flow distributions in microvascular networks. This allows us to systematically compare the impact of different vessel types in various regulation scenarios. Moreover, the method offers the advantage that it handles the stochastic nature of blood flow originating from tracking the movement of individual red blood cells. Since the inverse problem is formulated for time-averaged pressures and flow rates, a deterministic approach for calculating the diameter changes is used, which allows us to apply the method for large realistic microvascular networks with high-dimensional parameter spaces. Our results obtained in both artificial and realistic microvascular networks reveal that diameter changes at the level of capillaries enable a very localized regulation of blood flow. In scenarios where only larger vessels, i.e., arterioles, are allowed to adapt, the flow increase cannot be confined to a specific activated region and flow changes spread into neighboring regions. Furthermore, relatively small dilations and constrictions of all vessel types can lead to substantial changes of capillary blood flow distributions. This suggests that small scale regulation is necessary to obtain a localized increase in blood flow.

A 3D-1D coupled blood flow and oxygen transport model to generate microvascular networks

In this work, we introduce an algorithmic approach to generate microvascular networks starting from larger vessels that can be reconstructed without noticeable segmentation errors. Contrary to larger vessels, the reconstruction of fine-scale components of microvascular networks shows significant segmentation errors, and an accurate mapping is time and cost intense. Thus there is a need for fast and reliable reconstruction algorithms yielding surrogate networks having similar stochastic properties as the original ones. The microvascular networks are constructed in a marching way by adding vessels to the outlets of the vascular tree from the previous step. To optimise the structure of the vascular trees, we use Murray's law to determine the radii of the vessels and bifurcation angles. In each step, we compute the local gradient of the partial pressure of oxygen and adapt the orientation of the new vessels to this gradient. At the same time, we use the partial pressure of oxygen to check whether the considered tissue block is supplied sufficiently with oxygen. Computing the partial pressure of oxygen, we use a 3D-1D coupled model for blood flow and oxygen transport. To decrease the complexity of a fully coupled 3D model, we reduce the blood vessel network to a 1D graph structure and use a bi-directional coupling with the tissue which is described by a 3D homogeneous porous medium. The resulting surrogate networks are analysed with respect to morphological and physiological aspects.

Computational feasibility of calculating the steady-state blood flow rate through the vasculature of the entire human body

The human body contains approximately 20 billion individual blood vessels that deliver nutrients and oxygen to tissues. While blood flow is a well-developed field of research, no previous studies have calculated the blood flow rates through more than 5 million connected vessels. The goal of this study was to test if it is computationally feasible to calculate the blood flow rates through a vasculature equal in size to that of the human body. We designed and implemented a two-step algorithm to calculate the blood flow rates using principles of steady-state fluid dynamics. Steady-state fluid dynamics is an accurate approximation for the microvascular and venous structures in the human body. To determine the computational feasibility, we measured and evaluated the execution time, scalability, and memory usage to quantify the computational requirements. We demonstrated that it is computationally feasible to calculate the blood flow rate through 17 billion vessels in 6.5 hours using 256 compute nodes. The computational modeling of blood flow rate in entire organisms may find application in research on drug delivery, treatment of cancer metastases, and modulation of physiological performance.

3D Digital Anatomic Angioarchitecture of the Rat Spinal Cord: A Synchrotron Radiation Micro-CT Study

Comprehensive analysis of 3D angioarchitecture within the intact rat spinal cord remains technically challenging due to its sophisticated anatomical properties. In this study, we aim to present a framework for ultrahigh-resolution digitalized mapping of the normal rat spinal cord angioarchitecture and to determine the physiological parameters using synchrotron radiation micro-CT (SRμCT). Male SD rats were used in this ex vivo study. After a proportional mixture of contrast agents perfusion, the intact spinal cord covered the cervical spinal from the upper of the 1st cervical vertebra to the 5th lumbar vertebra was harvested and cut into proper lengths within three distinct regions: Cervical 3–5 levels, Thoracic 10–12 levels, Lumbar 3–5 levels spinal cord and examined using SRμCT. This method enabled the replication of the complicated microvasculature network of the normal rat spinal cord at the ultrahigh-resolution level, allowing for the precise quantitative analysis of the vascular morphological difference among cervical, thoracic and lumbar spinal cord in a 3D manner. Apart from a series of delicate 3D digital anatomical maps of the rat spinal cord angioarchitecture ranging from the cervical and thoracic to the lumbar spinal cord were presented, the 3D reconstruction data of SRμCT made the 3D printing of the spinal cord targeted selected microvasculature reality, that possibly provided deep insight into the nature and role of spinal cord intricate angioarchitecture. Our data proposed a new approach to outline systematic visual and quantitative evaluations on the 3D arrangement of the entire hierarchical microvasculature of the normal rat spinal cord at ultrahigh resolution. The technique may have great potential and become useful for future research on the poorly understood nature and function of the neurovascular interaction, particularly to investigate their pathology changes in various models of neurovascular disease.

Interactive visualization and analysis of morphological skeletons of brain vasculature networks with VessMorphoVis

MOTIVATION

Accurate morphological models of brain vasculature are key to modeling and simulating cerebral blood flow in realistic vascular networks. This in silico approach is fundamental to revealing the principles of neurovascular coupling. Validating those vascular morphologies entails performing certain visual analysis tasks that cannot be accomplished with generic visualization frameworks. This limitation has a substantial impact on the accuracy of the vascular models employed in the simulation.

RESULTS

We present VessMorphoVis, an integrated suite of toolboxes for interactive visualization and analysis of vast brain vascular networks represented by morphological graphs segmented originally from imaging or microscopy stacks. Our workflow leverages the outstanding potentials of Blender, aiming to establish an integrated, extensible and domain-specific framework capable of interactive visualization, analysis, repair, high-fidelity meshing and high-quality rendering of vascular morphologies. Based on the initial feedback of the users, we anticipate that our framework will be an essential component in vascular modeling and simulation in the future, filling a gap that is at present largely unfulfilled.

AVAILABILITY AND IMPLEMENTATION

VessMorphoVis is freely available under the GNU public license on Github at https://github.com/BlueBrain/VessMorphoVis. The morphology analysis, visualization, meshing and rendering modules are implemented as an add-on for Blender 2.8 based on its Python API (application programming interface). The add-on functionality is made available to users through an intuitive graphical user interface, as well as through exhaustive configuration files calling the API via a feature-rich command line interface running Blender in background mode.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Linear systems analysis for laminar fMRI: Evaluating BOLD amplitude scaling for luminance contrast manipulations

A fundamental assumption of nearly all functional magnetic resonance imaging (fMRI) analyses is that the relationship between local neuronal activity and the blood oxygenation level dependent (BOLD) signal can be described as following linear systems theory. With the advent of ultra-high field (7T and higher) MRI scanners, it has become possible to perform sub-millimeter resolution fMRI in humans. A novel and promising application of sub-millimeter fMRI is measuring responses across cortical depth, i.e. laminar imaging. However, the cortical vasculature and associated directional blood pooling towards the pial surface strongly influence the cortical depth-dependent BOLD signal, particularly for gradient-echo BOLD. This directional pooling may potentially affect BOLD linearity across cortical depth. Here we assess whether the amplitude scaling assumption for linear systems theory holds across cortical depth. For this, we use stimuli with different luminance contrasts to elicit different BOLD response amplitudes. We find that BOLD amplitude across cortical depth scales with luminance contrast, and that this scaling is identical across cortical depth. Although nonlinearities may be present for different stimulus configurations and acquisition protocols, our results suggest that the amplitude scaling assumption for linear systems theory across cortical depth holds for luminance contrast manipulations in sub-millimeter laminar BOLD fMRI.

From whole-organ imaging to in-silico blood flow modeling: A new multi-scale network analysis for revisiting tissue functional anatomy

We present a multi-disciplinary image-based blood flow perfusion modeling of a whole organ vascular network for analyzing both its structural and functional properties. We show how the use of Light-Sheet Fluorescence Microscopy (LSFM) permits whole-organ micro-vascular imaging, analysis and modelling. By using adapted image post-treatment workflow, we could segment, vectorize and reconstruct the entire micro-vascular network composed of 1.7 million vessels, from the tissue-scale, inside a ∼ 25 × 5 × 1 = 125mm3 volume of the mouse fat pad, hundreds of times larger than previous studies, down to the cellular scale at micron resolution, with the entire blood perfusion modeled. Adapted network analysis revealed the structural and functional organization of meso-scale tissue as strongly connected communities of vessels. These communities share a distinct heterogeneous core region and a more homogeneous peripheral region, consistently with known biological functions of fat tissue. Graph clustering analysis also revealed two distinct robust meso-scale typical sizes (from 10 to several hundred times the cellular size), revealing, for the first time, strongly connected functional vascular communities. These community networks support heterogeneous micro-environments. This work provides the proof of concept that in-silico all-tissue perfusion modeling can reveal new structural and functional exchanges between micro-regions in tissues, found from community clusters in the vascular graph.

UTX/KDM6A Deletion Promotes Recovery of Spinal Cord Injury by Epigenetically Regulating Vascular Regeneration

The regeneration of the blood vessel system post spinal cord injury (SCI) is essential for the repair of neurological function. As a significant means to regulate gene expression, epigenetic regulation of angiogenesis in SCI is still largely unknown. Here, we found that Ubiquitously Transcribed tetratricopeptide repeat on chromosome X (UTX), the histone H3K27 demethylase, increased significantly in endothelial cells post SCI. Knockdown of UTX can promote the migration and tube formation of endothelial cells. The specific knockout of UTX in endothelial cells enhanced angiogenesis post SCI accompanied with improved neurological function. In addition, we found regulation of UTX expression can change the level of microRNA 24 (miR-24) in vitro. The physical binding of UTX to the promotor of miR-24 was indicated by chromatin immunoprecipitation (ChIP) assay. Meanwhile, methylation sequencing of endothelial cells demonstrated that UTX could significantly decrease the level of methylation in the miR-24 promotor. Furthermore, miR-24 significantly abolished the promoting effect of UTX deletion on angiogenesis in vitro and in vivo. Finally, we predicted the potential target mRNAs of miR-24 related to angiogenesis. We indicate that UTX deletion can epigenetically promote the vascular regeneration and functional recovery post SCI by forming a regulatory network with miR-24.

Emergent dynamics of neuromorphic nanowire networks

Neuromorphic networks are formed by random self-assembly of silver nanowires. Silver nanowires are coated with a polymer layer after synthesis in which junctions between two nanowires act as resistive switches, often compared with neurosynapses. We analyze the role of single junction switching in the dynamical properties of the neuromorphic network. Network transitions to a high-conductance state under the application of a voltage bias higher than a threshold value. The stability and permanence of this state is studied by shifting the voltage bias in order to activate or deactivate the network. A model of the electrical network with atomic switches reproduces the relation between individual nanowire junctions switching events with current pathway formation or destruction. This relation is further manifested in changes in 1/f power-law scaling of the spectral distribution of current. The current fluctuations involved in this scaling shift are considered to arise from an essential equilibrium between formation, stochastic-mediated breakdown of individual nanowire-nanowire junctions and the onset of different current pathways that optimize power dissipation. This emergent dynamics shown by polymer-coated Ag nanowire networks places this system in the class of optimal transport networks, from which new fundamental parallels with neural dynamics and natural computing problem-solving can be drawn.

Red blood cells stabilize flow in brain microvascular networks

Capillaries are the prime location for oxygen and nutrient exchange in all tissues. Despite their fundamental role, our knowledge of perfusion and flow regulation in cortical capillary beds is still limited. Here, we use in vivo measurements and blood flow simulations in anatomically accurate microvascular network to investigate the impact of red blood cells (RBCs) on microvascular flow. Based on these in vivo and in silico experiments, we show that the impact of RBCs leads to a bias toward equating the values of the outflow velocities at divergent capillary bifurcations, for which we coin the term “well-balanced bifurcations”. Our simulation results further reveal that hematocrit heterogeneity is directly caused by the RBC dynamics, i.e. by their unequal partitioning at bifurcations and their effect on vessel resistance. These results provide the first in vivo evidence of the impact of RBC dynamics on the flow field in the cortical microvasculature. By structural and functional analyses of our blood flow simulations we show that capillary diameter changes locally alter flow and RBC distribution. A dilation of 10% along a vessel length of 100 μm increases the flow on average by 21% in the dilated vessel downstream a well-balanced bifurcation. The number of RBCs rises on average by 27%. Importantly, RBC up-regulation proves to be more effective the more balanced the outflow velocities at the upstream bifurcation are. Taken together, we conclude that diameter changes at capillary level bear potential to locally change the flow field and the RBC distribution. Moreover, our results suggest that the balancing of outflow velocities contributes to the robustness of perfusion. Based on our in silico results, we anticipate that the bi-phasic nature of blood and small-scale regulations are essential for a well-adjusted oxygen and energy substrate supply.

Glucose and oxygen are key energy sources of the brain. As energy storage capabilities are limited in the brain, a continuous supply of oxygen and glucose via the bloodstream is crucial for the brain’s functioning. The bulk of discharge occurs at the level of capillaries, which are the smallest and most frequent vessels of the cortical vasculature. Nonetheless, our understanding of perfusion and topology of the capillary bed is still limited. Here, we use in vivo two-photon based blood flow measurements and numerical simulations in large realistic microvascular networks to study the flow in the cortical microvasculature. Our results reveal that the impact of red blood cells enhances the robustness of microvascular perfusion and increases the heterogeneity in red blood cell distribution. It is well established that higher neuronal activity leads to an increase in blood flow. However, the precise regulation mechanisms and their spatial extent remain largely unknown. We show that small-scale regulations locally alter flow and red blood cell distribution. We suggest that these mechanisms are key for an efficient and flexible circulatory system. Moreover, our results reveal a novel role of the bi-phasic nature of blood.

Glioblastoma multiforme restructures the topological connectivity of cerebrovascular networks

Glioblastoma multiforme alters healthy tissue vasculature by inducing angiogenesis and vascular remodeling. To fully comprehend the structural and functional properties of the resulting vascular network, it needs to be studied collectively by considering both geometric and topological properties. Utilizing Single Plane Illumination Microscopy (SPIM), the detailed capillary structure in entire healthy and tumor-bearing mouse brains could be resolved in three dimensions. At the scale of the smallest capillaries, the entire vascular systems of bulk U87- and GL261-glioblastoma xenografts, their respective cores, and healthy brain hemispheres were modeled as complex networks and quantified with fundamental topological measures. All individual vessel segments were further quantified geometrically and modular clusters were uncovered and characterized as meta-networks, facilitating an analysis of large-scale connectivity. An inclusive comparison of large tissue sections revealed that geometric properties of individual vessels were altered in glioblastoma in a relatively subtle way, with high intra- and inter-tumor heterogeneity, compared to the impact on the vessel connectivity. A network topology analysis revealed a clear decomposition of large modular structures and hierarchical network organization, while preserving most fundamental topological classifications, in both tumor models with distinct growth patterns. These results augment our understanding of cerebrovascular networks and offer a topological assessment of glioma-induced vascular remodeling. The findings may help understand the emergence of hypoxia and necrosis, and prove valuable for therapeutic interventions such as radiation or antiangiogenic therapy.

A new framework for assessing subject-specific whole brain circulation and perfusion using MRI-based measurements and a multi-scale continuous flow model

A large variety of severe medical conditions involve alterations in microvascular circulation. Hence, measurements or simulation of circulation and perfusion has considerable clinical value and can be used for diagnostics, evaluation of treatment efficacy, and for surgical planning. However, the accuracy of traditional tracer kinetic one-compartment models is limited due to scale dependency. As a remedy, we propose a scale invariant mathematical framework for simulating whole brain perfusion. The suggested framework is based on a segmentation of anatomical geometry down to imaging voxel resolution. Large vessels in the arterial and venous network are identified from time-of-flight (ToF) and quantitative susceptibility mapping (QSM). Macro-scale flow in the large-vessel-network is accurately modelled using the Hagen-Poiseuille equation, whereas capillary flow is treated as two-compartment porous media flow. Macro-scale flow is coupled with micro-scale flow by a spatially distributing support function in the terminal endings. Perfusion is defined as the transition of fluid from the arterial to the venous compartment. We demonstrate a whole brain simulation of tracer propagation on a realistic geometric model of the human brain, where the model comprises distinct areas of grey and white matter, as well as large vessels in the arterial and venous vascular network. Our proposed framework is an accurate and viable alternative to traditional compartment models, with high relevance for simulation of brain perfusion and also for restoration of field parameters in clinical brain perfusion applications.

Dynamic Changes in Microvascular Flow Conductivity and Perfusion After Myocardial Infarction Shown by Image-Based Modeling

Background Microcirculation is a decisive factor in tissue reperfusion inadequacy following myocardial infarction ( MI ). Nonetheless, experimental assessment of blood flow in microcirculation remains a bottleneck. We sought to model blood flow properties in coronary microcirculation at different time points after MI and to compare them with healthy conditions to obtain insights into alterations in cardiac tissue perfusion. Methods and Results We developed an image-based modeling framework that permitted feeding a continuum flow model with anatomical data previously obtained from the pig coronary microvasculature to calculate physiologically meaningful permeability tensors. The tensors encompassed the microvascular conductivity and were also used to estimate the arteriole-venule drop in pressure and myocardial blood flow. Our results indicate that the tensors increased in a bimodal pattern at infarcted areas on days 1 and 7 after MI while a nonphysiological decrease in arteriole-venule drop in pressure was observed; contrary, the tensors and the arteriole-venule drop in pressure on day 3 after MI , and in remote areas, were closer to values for healthy tissue. Myocardial blood flow calculated using the condition-dependent arteriole-venule drop in pressure decreased in infarcted areas. Last, we simulated specific modes of vascular remodeling, such as vasodilation, vasoconstriction, or pruning, and quantified their distinct impact on microvascular conductivity. Conclusions Our study unravels time- and region-dependent alterations of tissue perfusion related to the structural changes occurring in the coronary microvasculature due to MI . It also paves the way for conducting simulations in new therapeutic interventions in MI and for image-based microvascular modeling by applying continuum flow models in other biomedical scenarios.

Brain Capillary Networks Across Species: A few Simple Organizational Requirements Are Sufficient to Reproduce Both Structure and Function

Despite the key role of the capillaries in neurovascular function, a thorough characterization of cerebral capillary network properties is currently lacking. Here, we define a range of metrics (geometrical, topological, flow, mass transfer, and robustness) for quantification of structural differences between brain areas, organs, species, or patient populations and, in parallel, digitally generate synthetic networks that replicate the key organizational features of anatomical networks (isotropy, connectedness, space-filling nature, convexity of tissue domains, characteristic size). To reach these objectives, we first construct a database of the defined metrics for healthy capillary networks obtained from imaging of mouse and human brains. Results show that anatomical networks are topologically equivalent between the two species and that geometrical metrics only differ in scaling. Based on these results, we then devise a method which employs constrained Voronoi diagrams to generate 3D model synthetic cerebral capillary networks that are locally randomized but homogeneous at the network-scale. With appropriate choice of scaling, these networks have equivalent properties to the anatomical data, demonstrated by comparison of the defined metrics. The ability to synthetically replicate cerebral capillary networks opens a broad range of applications, ranging from systematic computational studies of structure-function relationships in healthy capillary networks to detailed analysis of pathological structural degeneration, or even to the development of templates for fabrication of 3D biomimetic vascular networks embedded in tissue-engineered constructs.

A computational model for microcirculation including Fahraeus-Lindqvist effect, plasma skimming and fluid exchange with the tissue interstitium

We present a two-phase model for microcirculation that describes the interaction of plasma with red blood cells. The model takes into account of typical effects characterizing the microcirculation, such as the Fahraeus-Lindqvist effect and plasma skimming. Besides these features, the model describes the interaction of capillaries with the surrounding tissue. More precisely, the model accounts for the interaction of capillary transmural flow with the surrounding interstitial pressure. Furthermore, the capillaries are represented as one-dimensional channels with arbitrary, possibly curved configuration. The latter two features rely on the unique ability of the model to account for variations of flow rate and pressure along the axis of the capillary, according to a local differential formulation of mass and momentum conservation. Indeed, the model stands on a solid mathematical foundation, which is also addressed in this work. In particular, we present the model derivation, the variational formulation, and its approximation using the finite element method. Finally, we conclude the work with a comparative computational study of the importance of the Fahraeus-Lindqvist, plasma skimming, and capillary leakage effects on the distribution of flow in a microvascular network.

Accurate flow in augmented networks (AFAN): an approach to generating three-dimensional biomimetic microfluidic networks with controlled flow

In vivo, microvasculature provides oxygen, nutrients, and soluble factors necessary for cell survival and function. The highly tortuous, densely-packed, and interconnected three-dimensional (3D) architecture of microvasculature ensures that cells receive these crucial components. The ability to duplicate microvascular architecture in tissue-engineered models could provide a means to generate large-volume constructs as well as advanced microphysiological systems. Similarly, the ability to induce realistic flow in engineered microvasculature is crucial to recapitulating in vivo-like flow and transport. Advanced biofabrication techniques are capable of generating 3D, biomimetic microfluidic networks in hydrogels, however, these models can exhibit systemic aberrations in flow due to incorrect boundary conditions. To overcome this problem, we developed an automated method for generating synthetic augmented channels that induce the desired flow properties within three-dimensional microfluidic networks. These augmented inlets and outlets enforce the appropriate boundary conditions for achieving specified flow properties and create a three-dimensional output useful for image-guided fabrication techniques to create biomimetic microvascular networks.

Branching Pattern of the Cerebral Arterial Tree

Quantitative data on branching patterns of the human cerebral arterial tree are lacking in the 1.0–0.1 mm radius range. We aimed to collect quantitative data in this range, and to study if the cerebral artery tree complies with the principle of minimal work (Law of Murray).

To enable easy quantification of branching patterns a semi‐automatic method was employed to measure 1,294 bifurcations and 2,031 segments on 7 T‐MRI scans of two corrosion casts embedded in a gel. Additionally, to measure segments with a radius smaller than 0.1 mm, 9.4 T‐MRI was used on a small cast section to characterize 1,147 bifurcations and 1,150 segments. Besides MRI, traditional methods were employed. Seven hundred thirty‐three bifurcations were manually measured on a corrosion cast and 1,808 bifurcations and 1,799 segment lengths were manually measured on a fresh dissected cerebral arterial tree. Data showed a large variation in branching pattern parameters (asymmetry‐ratio, area‐ratio, length‐radius‐ratio, tapering). Part of the variation may be explained by the variation in measurement techniques, number of measurements and location of measurement in the vascular tree.

This study confirms that the cerebral arterial tree complies with the principle of minimum work. These data are essential in the future development of more accurate mathematical blood flow models. Anat Rec, 302:1434–1446, 2019. © 2018 The Authors. The Anatomical Record published by Wiley Periodicals, Inc. on behalf of American Association of Anatomists.

The Relation Between Capillary Transit Times and Hemoglobin Saturation Heterogeneity. Part 2: Capillary Networks

Brain metabolism is highly dependent on continuous oxygen supply. Cortical microvascular networks exhibit heterogeneous blood flow, leading to non-uniform tissue oxygenation and capillary hemoglobin saturation. We recently proposed capillary outflow saturation heterogeneity (COSH) to represent effects of heterogeneity on oxygen supply to tissue regions most vulnerable to hypoxia, and showed that diffusive oxygen exchange among red blood cells within capillaries and among capillaries (diffusive interaction) significantly reduces COSH in simplified geometrical configurations. Here, numerical simulations of oxygen transport in capillary network geometries derived from mouse somatosensory cortex are presented. Diffusive interaction was found to reduce COSH by 41 to 62% compared to simulations where diffusive interaction was excluded. Hemoglobin saturation drop across the microvascular network is strongly correlated with red blood cell transit time, but the coefficient of variation of saturation drop is approximately one third lower. Unexpectedly, the radius of the tissue cylinder supplied by a capillary correlates weakly with the anatomical tissue cylinder radius, but strongly with hemoglobin saturation. Thus, diffusive interaction contributes greatly to the microcirculation's ability to achieve tissue oxygenation, despite heterogeneous capillary transit time and hematocrit distribution. These findings provide insight into the effects of cerebral small vessel disease on tissue oxygenation and brain function.

Toward quantitative three-dimensional microvascular networks segmentation with multiview light-sheet fluorescence microscopy

Three-dimensional (3-D) large-scale imaging of microvascular networks is of interest in various areas of biology and medicine related to structural, functional, developmental, and pathological issues. Light-sheet fluorescence microscopy (LSFM) techniques are rapidly spreading and are now on the way to offer operational solutions for large-scale tissue imaging. This contribution describes how reliable vessel segmentation can be handled from LSFM data in very large tissue volumes using a suitable image analysis workflow. Since capillaries are tubular objects of a few microns scale radius, they represent challenging structures to reliably reconstruct without distortion and artifacts. We provide a systematic analysis of multiview deconvolution image processing workflow to control and evaluate the accuracy of the reconstructed vascular network using various low to high level, metrics. We show that even if low-level structural metrics are sensitive to isotropic imaging enhancement provided by a larger number of views, functional high-level metrics, including perfusion permeability, are less sensitive. Hence, combining deconvolution and registration onto a few number of views appears sufficient for a reliable quantitative 3-D vessel segmentation for their possible use for perfusion modeling.

Scaling behavior of drug transport and absorption in in silico cerebral capillary networks

Drug delivery to the brain is challenging due to the presence of the blood-brain barrier. Mathematical modeling and simulation are essential tools for the deeper understanding of transport processes in the blood, across the blood-brain barrier and within the tissue. Here we present a mathematical model for drug delivery through capillary networks with increasingly complex topologies with the goal to understand the scaling behavior of model predictions on a coarse-to-fine sequence of grids. We apply our model to the delivery of L-Dopa, the primary drug used in the therapy of Parkinson’s Disease. Our model replicates observed blood flow rates and ratios between plasma and tissue concentrations. We propose an optimal network grain size for the simulation of tissue volumes of 1 cm³ that allows to make reliable predictions with reasonable computational costs.

How does blood regulate cerebral temperatures during hypothermia?

Macro-modeling of cerebral blood flow can help determine the impact of thermal intervention during instances of head trauma to mitigate tissue damage. This work presents a bioheat model using a 3D fluid-porous domain coupled with intersecting 1D arterial and venous vessel trees. This combined vascular porous (VaPor) model resolves both cerebral blood flow and energy equations, including heat generated by metabolism, using vasculature extracted from MRI data and is extended using a tree generation algorithm. Counter-current flows are expected to increase thermal transfer within the brain and are enforced using either the vascular structure or flow reversal, represented by a flow reversal constant, C_R. These methods exhibit larger average brain cooling (from 0.56 °C ± <0.01 °C to 0.58 °C ± <0.01 °C) compared with previous models (0.39 °C) when scalp temperature is reduced. An greater reduction in core brain temperature is observed (from 0.29 °C ± <0.01 °C to 0.45 °C ± <0.01 °C) compared to previous models (0.11 °C) due to the inclusion of counter-current cooling effects. The VaPor model also predicts that a hypothermic average temperature (<36 °C) can be reached in core regions of neonatal models using scalp cooling alone.

Numerical modelling of a peripheral arterial stenosis using dimensionally reduced models and kernel methods

In this work, we consider 2 kinds of model reduction techniques to simulate blood flow through the largest systemic arteries, where a stenosis is located in a peripheral artery, i.e., in an artery that is located far away from the heart. For our simulations, we place the stenosis in one of the tibial arteries belonging to the right lower leg (right posterior tibial artery). The model reduction techniques that are used are on the one hand dimensionally reduced models (1-D and 0-D models, the so-called mixed-dimension model) and on the other hand surrogate models produced by kernel methods. Both methods are combined in such a way that the mixed-dimension models yield training data for the surrogate model, where the surrogate model is parametrised by the degree of narrowing of the peripheral stenosis. By means of a well-trained surrogate model, we show that simulation data can be reproduced with a satisfactory accuracy and that parameter optimisation or state estimation problems can be solved in a very efficient way. Furthermore, it is demonstrated that a surrogate model enables us to present after a very short simulation time the impact of a varying degree of stenosis on blood flow, obtaining a speedup of several orders over the full model.