Vectorization of optically sectioned brain microvasculature: learning aids completion of vascular graphs by connecting gaps and deleting open-ended segments

Mesh-free high-resolution simulation of cerebrocortical oxygen supply with fast Fourier preconditioning

Oxygen transfer from blood vessels to cortical brain tissue is representative of a class of problems with mixed-domain character. Large-scale efficient computation of tissue oxygen concentration is dependent on the manner in which the tubular network of blood vessels is coupled to the tissue mesh. Models which explicitly resolve the interface between the tissue and vasculature with a contiguous mesh are prohibitively expensive for very dense cerebral microvasculature. We propose a mixed-domain mesh-free technique whereby a vascular anatomical network (VAN) represented as a thin directed graph serves for convection of blood oxygen, and the surrounding extravascular tissue is represented as a Cartesian grid of 3D voxels throughout which oxygen is transported by diffusion. We split the network and tissue meshes by the Schur complement method of domain decomposition to obtain a reduced set of system equations for the tissue oxygen concentration at steady state. The use of a Cartesian grid allows the corresponding matrix equation to be solved approximately with a fast Fourier transform-based Poisson solver, which serves as an effective preconditioner for Krylov subspace iteration. The performance of this method enables the steady-state simulation of cortical oxygen perfusion for anatomically accurate vascular networks down to single micron resolution without the need for supercomputers.

MESH-FREE HIGH-RESOLUTION SIMULATION OF CEREBROCORTICAL OXYGEN SUPPLY WITH FAST FOURIER PRECONDITIONING

Oxygen transfer from blood vessels to cortical brain tissue is representative of a class of problems with mixed-domain character. Large-scale efficient computation of tissue oxygen concentration is dependent on the manner in which the tubular network of blood vessels is coupled to the tissue mesh. Models which explicitly resolve the interface between the tissue and vasculature with a contiguous mesh are prohibitively expensive for very dense cerebral microvasculature. We propose a mixed-domain mesh-free technique whereby a vascular anatomical network (VAN) represented as a thin directed graph serves for convection of blood oxygen, and the surrounding extravascular tissue is represented as a Cartesian grid of 3D voxels throughout which oxygen is transported by diffusion. We split the network and tissue meshes by the Schur complement method of domain decomposition to obtain a reduced set of system equations for the tissue oxygen concentration. The use of a Cartesian grid allows the corresponding matrix equation to be solved approximately with a fast Fourier transform based Poisson solver, which serves as an effective preconditioner for Krylov subspace iteration. The performance of this method enables the steady state simulation of cortical oxygen perfusion for anatomically accurate vascular networks down to single micron resolution without the need for supercomputers.

Practitioner Points

We present a novel mixed-domain framework for efficiently modeling O ₂ extraction kinetics in the brain. Model equations are generated by graph-theoretic methods for mixed domains.Dual mesh domain decomposition with FFT preconditioning yields very fast simulation times for extremely high spatial resolution.

Mathematical synthesis of the cortical circulation for the whole mouse brain-part II: Microcirculatory closure

Recent advancements in multiphoton imaging and vascular reconstruction algorithms have increased the amount of data on cerebrovascular circulation for statistical analysis and hemodynamic simulations. Experimental observations offer fundamental insights into capillary network topology but mainly within a narrow field of view typically spanning a small fraction of the cortical surface (less than 2%). In contrast, larger-resolution imaging modalities, such as computed tomography (CT) or magnetic resonance imaging (MRI), have whole-brain coverage but capture only larger blood vessels, overlooking the microscopic capillary bed. To integrate data acquired at multiple length scales with different neuroimaging modalities and to reconcile brain-wide macroscale information with microscale multiphoton data, we developed a method for synthesizing hemodynamically equivalent vascular networks for the entire cerebral circulation. This computational approach is intended to aid in the quantification of patterns of cerebral blood flow and metabolism for the entire brain. In part I, we described the mathematical framework for image-guided generation of synthetic vascular networks covering the large cerebral arteries from the circle of Willis through the pial surface network leading back to the venous sinuses. Here in part II, we introduce novel procedures for creating microcirculatory closure that mimics a realistic capillary bed. We demonstrate our capability to synthesize synthetic vascular networks whose morphometrics match empirical network graphs from three independent state-of-the-art imaging laboratories using different image acquisition and reconstruction protocols. We also successfully synthesized twelve vascular networks of a complete mouse brain hemisphere suitable for performing whole-brain blood flow simulations. Synthetic arterial and venous networks with microvascular closure allow whole-brain hemodynamic predictions. Simulations across all length scales will potentially illuminate organ-wide supply and metabolic functions that are inaccessible to models reconstructed from image data with limited spatial coverage.

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.

Voxelized simulation of cerebral oxygen perfusion elucidates hypoxia in aged mouse cortex

Departures of normal blood flow and metabolite distribution from the cerebral microvasculature into neuronal tissue have been implicated with age-related neurodegeneration. Mathematical models informed by spatially and temporally distributed neuroimage data are becoming instrumental for reconstructing a coherent picture of normal and pathological oxygen delivery throughout the brain. Unfortunately, current mathematical models of cerebral blood flow and oxygen exchange become excessively large in size. They further suffer from boundary effects due to incomplete or physiologically inaccurate computational domains, numerical instabilities due to enormous length scale differences, and convergence problems associated with condition number deterioration at fine mesh resolutions. Our proposed simple finite volume discretization scheme for blood and oxygen microperfusion simulations does not require expensive mesh generation leading to the critical benefit that it drastically reduces matrix size and bandwidth of the coupled oxygen transfer problem. The compact problem formulation yields rapid and stable convergence. Moreover, boundary effects can effectively be suppressed by generating very large replica of the cortical microcirculation in silico using an image-based cerebrovascular network synthesis algorithm, so that boundaries of the perfusion simulations are far removed from the regions of interest. Massive simulations over sizeable portions of the cortex with feature resolution down to the micron scale become tractable with even modest computer resources. The feasibility and accuracy of the novel method is demonstrated and validated with in vivo oxygen perfusion data in cohorts of young and aged mice. Our oxygen exchange simulations quantify steep gradients near penetrating blood vessels and point towards pathological changes that might cause neurodegeneration in aged brains. This research aims to explain mechanistic interactions between anatomical structures and how they might change in diseases or with age. Rigorous quantification of age-related changes is of significant interest because it might aide in the search for imaging biomarkers for dementia and Alzheimer’s disease.

Brain function critically depends on the maintenance of physiological blood supply and metabolism in the cortex. Disturbances to adequate perfusion have been linked to age-related neurodegeneration. However, the precise correlation between age-related hemodynamic changes and the resulting decline in oxygen delivery is not well understood and has not been quantified. Therefore, we introduce a new compact, and therefore highly scalable, computational method for predicting the physiological relationship between hemodynamics and cortical oxygen perfusion for large sections of the cortical microcirculation. We demonstrate the novel mesh generation-free (MGF), multi-scale simulation approach through realistic in vivo case studies of cortical microperfusion in the mouse brain. We further validate mechanistic correlations and a quantitative relationship between blood flow and brain oxygenation using experimental data from cohorts of young, middle aged and old mouse brains. Our computational approach overcomes size and performance limitations of previous unstructured meshing techniques to enable the prediction of oxygen tension with a spatial resolution of least two orders of magnitude higher than previously possible. Our simulation results support the hypothesis that structural changes in the microvasculature induce hypoxic pockets in the aged brain that are absent in the healthy, young mouse.

Deep learning toolbox for automated enhancement, segmentation, and graphing of cortical optical coherence tomography microangiograms

Optical coherence tomography angiography (OCTA) is becoming increasingly popular for neuroscientific study, but it remains challenging to objectively quantify angioarchitectural properties from 3D OCTA images. This is mainly due to projection artifacts or "tails" underneath vessels caused by multiple-scattering, as well as the relatively low signal-to-noise ratio compared to fluorescence-based imaging modalities. Here, we propose a set of deep learning approaches based on convolutional neural networks (CNNs) to automated enhancement, segmentation and gap-correction of OCTA images, especially of those obtained from the rodent cortex. Additionally, we present a strategy for skeletonizing the segmented OCTA and extracting the underlying vascular graph, which enables the quantitative assessment of various angioarchitectural properties, including individual vessel lengths and tortuosity. These tools, including the trained CNNs, are made publicly available as a user-friendly toolbox for researchers to input their OCTA images and subsequently receive the underlying vascular network graph with the associated angioarchitectural properties.

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.

Simulations of blood as a suspension predicts a depth dependent hematocrit in the circulation throughout the cerebral cortex

Recent advances in modeling oxygen supply to cortical brain tissue have begun to elucidate the functional mechanisms of neurovascular coupling. While the principal mechanisms of blood flow regulation after neuronal firing are generally known, mechanistic hemodynamic simulations cannot yet pinpoint the exact spatial and temporal coordination between the network of arteries, arterioles, capillaries and veins for the entire brain. Because of the potential significance of blood flow and oxygen supply simulations for illuminating spatiotemporal regulation inside the cortical microanatomy, there is a need to create mathematical models of the entire cerebral circulation with realistic anatomical detail. Our hypothesis is that an anatomically accurate reconstruction of the cerebrocirculatory architecture will inform about possible regulatory mechanisms of the neurovascular interface. In this article, we introduce large-scale networks of the murine cerebral circulation spanning the Circle of Willis, main cerebral arteries connected to the pial network down to the microcirculation in the capillary bed. Several multiscale models were generated from state-of-the-art neuroimaging data. Using a vascular network construction algorithm, the entire circulation of the middle cerebral artery was synthesized. Blood flow simulations indicate a consistent trend of higher hematocrit in deeper cortical layers, while surface layers with shorter vascular path lengths seem to carry comparatively lower red blood cell (RBC) concentrations. Moreover, the variability of RBC flux decreases with cortical depth. These results support the notion that plasma skimming serves a self-regulating function for maintaining uniform oxygen perfusion to neurons irrespective of their location in the blood supply hierarchy. Our computations also demonstrate the practicality of simulating blood flow for large portions of the mouse brain with existing computer resources. The efficient simulation of blood flow throughout the entire middle cerebral artery (MCA) territory is a promising milestone towards the final aim of predicting blood flow patterns for the entire brain.

The brain’s astonishing cognitive capacity depends on the coordination between neurons and the cerebral circulation, a system known as the neurovascular unit. The spatial and temporal coupling between these two networks is the object of intense research. However, the concise anatomical description of the cerebral circulation has so far been intractable. This paper introduces a methodology for the in silico creation of realistic models for the entire cerebral circulation. This innovation incorporates topological data from several neuroimaging modalities covering three lengths scales as input into a computer algorithm, which assembles anatomically accurate circulatory networks. When simulating blood flow as red blood cells suspended in plasma for experimental and synthetic cortical network models, we discovered that red blood cells tend to be more concentrated in deeper layers of the cortex compared to the surface. RBC fluxes are more homogenous in deeper layers. The phenomenon of depth dependent red blood cell supply supports the notion that the intricate architecture of the cortical microcirculation serves a self-regulating function to maintain uniform oxygen perfusion to neurons. We also demonstrate the practicality of predicting blood flow patterns for the entire brain with existing computer power.

Whole-Brain Vasculature Reconstruction at the Single Capillary Level

The distinct organization of the brain’s vascular network ensures that it is adequately supplied with oxygen and nutrients. However, despite this fundamental role, a detailed reconstruction of the brain-wide vasculature at the capillary level remains elusive, due to insufficient image quality using the best available techniques. Here, we demonstrate a novel approach that improves vascular demarcation by combining CLARITY with a vascular staining approach that can fill the entire blood vessel lumen and imaging with light-sheet fluorescence microscopy. This method significantly improves image contrast, particularly in depth, thereby allowing reliable application of automatic segmentation algorithms, which play an increasingly important role in high-throughput imaging of the terabyte-sized datasets now routinely produced. Furthermore, our novel method is compatible with endogenous fluorescence, thus allowing simultaneous investigations of vasculature and genetically targeted neurons. We believe our new method will be valuable for future brain-wide investigations of the capillary network.

The impact of vessel size, orientation and intravascular contribution on the neurovascular fingerprint of BOLD bSSFP fMRI

Monte Carlo simulations have been used to analyze oxygenation-related signal changes in pass-band balanced steady state free precession (bSSFP) as well as in gradient echo (GE) and spin echo (SE) sequences. Signal changes were calculated for artificial cylinders and neurovascular networks acquired from the mouse parietal cortex by two-photon laser scanning microscopy at 1 μm isotropic resolution. Signal changes as a function of vessel size, blood volume, vessel orientation to the main magnetic field B₀ as well as relations of intra- and extravascular and of micro- and macrovascular contributions have been analyzed. The results show that bSSFP is highly sensitive to extravascular and microvascular components. Furthermore, GE and bSSFP, and to a lesser extent SE, exhibit a strong dependence of their signal change on the orientation of the vessel network to B₀.

In vivo multiphoton imaging of a diverse array of fluorophores to investigate deep neurovascular structure

We perform high-resolution, non-invasive, in vivo deep-tissue imaging of the mouse neocortex using multiphoton microscopy with a high repetition rate optical parametric amplifier laser source tunable between λ=1,100 and 1,400 nm. By combining the high repetition rate (511 kHz) and high pulse energy (400 nJ) of our amplifier laser system, we demonstrate imaging of vasculature labeled with Texas Red and Indocyanine Green, and neurons expressing tdTomato and yellow fluorescent protein. We measure the blood flow speed of a single capillary at a depth of 1.2 mm, and image vasculature to a depth of 1.53 mm with fine axial steps (5 μm) and reasonable acquisition times. The high image quality enabled analysis of vascular morphology at depths to 1.45 mm.

NET: a new framework for the vectorization and examination of network data

Background

The analysis of complex networks both in general and in particular as pertaining to real biological systems has been the focus of intense scientific attention in the past and present. In this paper we introduce two tools that provide fast and efficient means for the processing and quantification of biological networks like Drosophila tracheoles or leaf venation patterns: the Network Extraction Tool (NET) to extract data and the Graph-edit-GUI (GeGUI) to visualize and modify networks.

Results

NET is especially designed for high-throughput semi-automated analysis of biological datasets containing digital images of networks. The framework starts with the segmentation of the image and then proceeds to vectorization using methodologies from optical character recognition. After a series of steps to clean and improve the quality of the extracted data the framework produces a graph in which the network is represented only by its nodes and neighborhood-relations. The final output contains information about the adjacency matrix of the graph, the width of the edges and the positions of the nodes in space. NET also provides tools for statistical analysis of the network properties, such as the number of nodes or total network length. Other, more complex metrics can be calculated by importing the vectorized network to specialized network analysis packages.

GeGUI is designed to facilitate manual correction of non-planar networks as these may contain artifacts or spurious junctions due to branches crossing each other. It is tailored for but not limited to the processing of networks from microscopy images of Drosophila tracheoles.

Conclusion

The networks extracted by NET closely approximate the network depicted in the original image. NET is fast, yields reproducible results and is able to capture the full geometry of the network, including curved branches. Additionally GeGUI allows easy handling and visualization of the networks.

Electronic supplementary material

The online version of this article (doi:10.1186/s13029-017-0064-3) contains supplementary material, which is available to authorized users.

Depth-dependent flow and pressure characteristics in cortical microvascular networks

A better knowledge of the flow and pressure distribution in realistic microvascular networks is needed for improving our understanding of neurovascular coupling mechanisms and the related measurement techniques. Here, numerical simulations with discrete tracking of red blood cells (RBCs) are performed in three realistic microvascular networks from the mouse cerebral cortex. Our analysis is based on trajectories of individual RBCs and focuses on layer-specific flow phenomena until a cortical depth of 1 mm. The individual RBC trajectories reveal that in the capillary bed RBCs preferentially move in plane. Hence, the capillary flow field shows laminar patterns and a layer-specific analysis is valid. We demonstrate that for RBCs entering the capillary bed close to the cortical surface (< 400 μm) the largest pressure drop takes place in the capillaries (37%), while for deeper regions arterioles are responsible for 61% of the total pressure drop. Further flow characteristics, such as capillary transit time or RBC velocity, also vary significantly over cortical depth. Comparison of purely topological characteristics with flow-based ones shows that a combined interpretation of topology and flow is indispensable. Our results provide evidence that it is crucial to consider layer-specific differences for all investigations related to the flow and pressure distribution in the cortical vasculature. These findings support the hypothesis that for an efficient oxygen up-regulation at least two regulation mechanisms must be playing hand in hand, namely cerebral blood flow increase and microvascular flow homogenization. However, the contribution of both regulation mechanisms to oxygen up-regulation likely varies over depth.

Gap-free segmentation of vascular networks with automatic image processing pipeline

Current image processing techniques capture large vessels reliably but often fail to preserve connectivity in bifurcations and small vessels. Imaging artifacts and noise can create gaps and discontinuity of intensity that hinders segmentation of vascular trees. However, topological analysis of vascular trees require proper connectivity without gaps, loops or dangling segments. Proper tree connectivity is also important for high quality rendering of surface meshes for scientific visualization or 3D printing. We present a fully automated vessel enhancement pipeline with automated parameter settings for vessel enhancement of tree-like structures from customary imaging sources, including 3D rotational angiography, magnetic resonance angiography, magnetic resonance venography, and computed tomography angiography. The output of the filter pipeline is a vessel-enhanced image which is ideal for generating anatomical consistent network representations of the cerebral angioarchitecture for further topological or statistical analysis. The filter pipeline combined with computational modeling can potentially improve computer-aided diagnosis of cerebrovascular diseases by delivering biometrics and anatomy of the vasculature. It may serve as the first step in fully automatic epidemiological analysis of large clinical datasets. The automatic analysis would enable rigorous statistical comparison of biometrics in subject-specific vascular trees. The robust and accurate image segmentation using a validated filter pipeline would also eliminate operator dependency that has been observed in manual segmentation. Moreover, manual segmentation is time prohibitive given that vascular trees have more than thousands of segments and bifurcations so that interactive segmentation consumes excessive human resources. Subject-specific trees are a first step toward patient-specific hemodynamic simulations for assessing treatment outcomes.

The capillary bed offers the largest hemodynamic resistance to the cortical blood supply

The cortical angioarchitecture is a key factor in controlling cerebral blood flow and oxygen metabolism. Difficulties in imaging the complex microanatomy of the cortex have so far restricted insight about blood flow distribution in the microcirculation. A new methodology combining advanced microscopy data with large scale hemodynamic simulations enabled us to quantify the effect of the angioarchitecture on the cerebral microcirculation. High-resolution images of the mouse primary somatosensory cortex were input into with a comprehensive computational model of cerebral perfusion and oxygen supply ranging from the pial vessels to individual brain cells. Simulations of blood flow, hematocrit and oxygen tension show that the wide variation of hemodynamic states in the tortuous, randomly organized capillary bed is responsible for relatively uniform cortical tissue perfusion and oxygenation. Computational analysis of microcirculatory blood flow and pressure drops further indicates that the capillary bed, including capillaries adjacent to feeding arterioles (d < 10 µm), are the largest contributors to hydraulic resistance.

Imaging and graphing of cortical vasculature using dynamically focused optical coherence microscopy angiography

Recently, optical coherence tomography (OCT) angiography has enabled label-free imaging of vasculature based on dynamic scattering in vessels. However, quantitative volumetric analysis of the vascular networks depicted in OCT angiography data has remained challenging. Multiple-scattering tails (artifacts specific to the imaging geometry) make automated assessment of vascular morphology problematic. We demonstrate that dynamically focused optical coherence microscopy (OCM) angiography with a high numerical aperture, chosen so the scattering length greatly exceeds the depth-of-field, significantly reduces the deleterious effect of multiple-scattering tails in synthesized angiograms. Capitalizing on the improved vascular image quality, we devised and tailored a self-correcting automated graphing approach that achieves a reconstruction of cortical microvasculature from OCM angiography data sets with accuracy approaching that attained by trained operators. The automated techniques described here will facilitate more widespread study of vascular network topology in health and disease.

Mapping the 3D Connectivity of the Rat Inner Retinal Vascular Network Using OCT Angiography

PURPOSE

The purpose of this study is to demonstrate three-dimensional (3D) graphing based on optical coherence tomography (OCT) angiography for characterization of the inner retinal vascular architecture and determination of its topologic principles.

METHODS

Rat eyes (N = 3) were imaged with a 1300-nm spectral/Fourier domain OCT microscope. A topologic model of the inner retinal vascular network was obtained from OCT angiography data using a combination of automated and manually-guided image processing techniques. Using a resistive network model, with experimentally-quantified flow in major retinal vessels near the optic nerve head as boundary conditions, theoretical changes in the distribution of flow induced by vessel dilations were inferred.

RESULTS

A topologically-representative 3D vectorized graph of the inner retinal vasculature, derived from OCT angiography data, is presented. The laminar and compartmental connectivity of the vasculature are characterized. In contrast to sparse connectivity between the superficial vitreal vasculature and capillary plexuses of the inner retina, connectivity between the two capillary plexus layers is dense. Simulated dilation of single arterioles is shown to produce both localized and lamina-specific changes in blood flow, while dilation of capillaries in a given retinal vascular layer is shown to lead to increased total flow in that layer.

CONCLUSIONS

Our graphing and modeling data suggest that vascular architecture enables both local and lamina-specific control of blood flow in the inner retina. The imaging, graph analysis, and modeling approach presented here will help provide a detailed characterization of vascular changes in a variety of retinal diseases, both in experimental preclinical models and human subjects.

Sequential Monte Carlo tracking of the marginal artery by multiple cue fusion and random forest regression

Given the potential importance of marginal artery localization in automated registration in computed tomography colonography (CTC), we have devised a semi-automated method of marginal vessel detection employing sequential Monte Carlo tracking (also known as particle filtering tracking) by multiple cue fusion based on intensity, vesselness, organ detection, and minimum spanning tree information for poorly enhanced vessel segments. We then employed a random forest algorithm for intelligent cue fusion and decision making which achieved high sensitivity and robustness. After applying a vessel pruning procedure to the tracking results, we achieved statistically significantly improved precision compared to a baseline Hessian detection method (2.7% versus 75.2%, p<0.001). This method also showed statistically significantly improved recall rate compared to a 2-cue baseline method using fewer vessel cues (30.7% versus 67.7%, p<0.001). These results demonstrate that marginal artery localization on CTC is feasible by combining a discriminative classifier (i.e., random forest) with a sequential Monte Carlo tracking mechanism. In so doing, we present the effective application of an anatomical probability map to vessel pruning as well as a supplementary spatial coordinate system for colonic segmentation and registration when this task has been confounded by colon lumen collapse.

The cortical angiome: an interconnected vascular network with noncolumnar patterns of blood flow

What is the nature of the vascular architecture in the cortex that allows the brain to meet the energy demands of neuronal computations? We used high-throughput histology to reconstruct the complete angioarchitecture and the positions of all neuronal somata of multiple cubic millimeter regions of vibrissa primary sensory cortex in mouse. Vascular networks were derived from the reconstruction. In contrast with the standard model of cortical columns that are tightly linked with the vascular network, graph-theoretical analyses revealed that the subsurface microvasculature formed interconnected loops with a topology that was invariant to the position and boundary of columns. Furthermore, the calculated patterns of blood flow in the networks were unrelated to location of columns. Rather, blood sourced by penetrating arterioles was effectively drained by the penetrating venules to limit lateral perfusion. This analysis provides the underpinning to understand functional imaging and the effect of penetrating vessels strokes on brain viability.