Modeling of Cerebral Oxygen Transport Based on In vivo Microscopic Imaging of Microvascular Network Structure, Blood Flow, and Oxygenation
- 1Optics Division, Department of Radiology, MHG/MIT/HMS Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital and Harvard Medical School, Charlestown, MA, USA
- 2Institut de Mécanique des Fluides de Toulouse, Toulouse, France
- 3Department of Physiology, University of Arizona, Tucson, AZ, USA
- 4Departments of Neurosciences and Radiology, University of California, San Diego, La Jolla, CA, USA
Oxygen is delivered to brain tissue by a dense network of microvessels, which actively control cerebral blood flow (CBF) through vasodilation and contraction in response to changing levels of neural activity. Understanding these network-level processes is immediately relevant for (1) interpretation of functional Magnetic Resonance Imaging (fMRI) signals, and (2) investigation of neurological diseases in which a deterioration of neurovascular and neuro-metabolic physiology contributes to motor and cognitive decline. Experimental data on the structure, flow and oxygen levels of microvascular networks are needed, together with theoretical methods to integrate this information and predict physiologically relevant properties that are not directly measurable. Recent progress in optical imaging technologies for high-resolution in vivo measurement of the cerebral microvascular architecture, blood flow, and oxygenation enables construction of detailed computational models of cerebral hemodynamics and oxygen transport based on realistic three-dimensional microvascular networks. In this article, we review state-of-the-art optical microscopy technologies for quantitative in vivo imaging of cerebral microvascular structure, blood flow and oxygenation, and theoretical methods that utilize such data to generate spatially resolved models for blood flow and oxygen transport. These “bottom-up” models are essential for the understanding of the processes governing brain oxygenation in normal and disease states and for eventual translation of the lessons learned from animal studies to humans.
The energy requirements of the cerebral cortex are met almost exclusively by oxidative metabolism of glucose, necessitating continuous delivery of oxygen by diffusion into the tissue from flowing blood (Buxton and Frank, 1997; Attwell and Laughlin, 2001; Raichle and Gusnard, 2002; Lin et al., 2010; Hall et al., 2012). The distance that oxygen can diffuse from blood into oxygen-consuming tissue is limited. An approximate upper bound is the one-dimensional diffusion distance, LD = (2DαP0/M)1/2, where D and α are the diffusivity and solubility of oxygen in tissue, P0 is the partial pressure of oxygen (pO2) in blood, and M is the rate of oxygen consumption (Secomb et al., 1993). Typical values for rodent cerebral cortex are Dα = 6 × 10−10 cm3O2/cm/s/mmHg, P0 = 50 mmHg and M = 6 cm3O2/100 g/min, resulting in LD ≈ 77 μm (Secomb et al., 2000). The tissue must therefore be supplied by a dense network of microvessels (Hirsch et al., 2012; Blinder et al., 2013), so that each point in the tissue is within a distance LD of the nearest vessel. A further implication of this result is that significant gradients in PO2, over a scale of tens of μm, are present in the tissue surrounding each microvessel that is contributing to the diffusive delivery of oxygen between blood and tissue (Popel, 1989; Pittman, 2005; Goldman, 2008). This does not necessarily imply that significant gradients in PO2 are present around all microvessels. For instance, if the oxygen demand of the tissue is largely met by diffusion from arterioles (Sakadžić et al., 2014), then PO2 gradients around neighboring capillaries are correspondingly reduced. This is more likely to occur at lower oxygen consumption rates (Secomb and Hsu, 1994). As a further example, heterogeneous levels of oxygen within capillaries can lead to situations where some capillaries do not participate in diffusive oxygen delivery and may even take up oxygen (Secomb and Hsu, 1994).
The oxygen demand of cerebral cortex depends on the level of neuronal activity, which varies in space and time across cortical regions. However, in spite of significant research efforts, our current understanding of the neurovascular coupling mechanisms by which blood flow and oxygen delivery are coordinated with spatially and temporally varying metabolic demands in the brain remains incomplete in both healthy and diseased brain (Iadecola, 2004; Girouard and Iadecola, 2006; Raichle and Mintun, 2006; Iadecola and Nedergaard, 2007; Gordon et al., 2008; Cauli and Hamel, 2010; Lecrux et al., 2011; Devor et al., 2012). In the context of task-induced hemodynamic activity in a healthy subject, neurovascular coupling is commonly associated with the causal chain of events, where changes in neural activity drive changes in energy metabolism which then drive changes in blood flow (Raichle and Mintun, 2006). However, recent experimental data suggests an alternative hypothesis that much of the acute vasodilation and constriction under healthy conditions appears to be driven by molecules related to neural signaling itself, i.e., release of vasoactive signaling agents such as neuropeptides, prostaglandins, nitric oxide and, possibly, K+ from active neurons (Attwell and Iadecola, 2002; Cauli and Hamel, 2010; Kleinfeld et al., 2011). This implies that CBF and CMRO2 are driven in parallel by neural activity, and potentially by different aspects of neural activity (Devor et al., 2012; Buxton et al., 2014).
The mechanisms by which blood flow and oxygen delivery are coordinated with spatially and temporally varying metabolic demand in the brain critically depend on vascular network-level interactions. In the present review, we focus on the microvascular networks (Figure 1) where the blood flow rate in any given segment depends not only on the flow resistance of that segment but also on the distributions of resistance and blood flow in the vessel segments to which it is connected. An increase in flow in a given segment requires dilation of the upstream arterioles supplying that segment, in addition to possible dilation of the segment itself and downstream segments. Furthermore, the oxygen content of blood arriving at a given segment depends on the extraction that has already occurred, and is therefore affected by oxygen transport phenomena occurring along the upstream part of the flow pathway. Conducted responses, in which vasodilator or vasoconstrictor signals travel upstream from capillaries along arterioles, play an important role by coordinating upstream dilation (Jensen and Holstein-Rathlou, 2013).
Figure 1. Schematic illustration of cortical vasculature, showing network structures involved in neurovascular coupling. The cerebral cortex receives its blood supply from a mesh of pial arteries and veins lying on the cortical surface. Penetrating arterioles and venules branch off the pial vessels and traverse the thickness of the cortex, supplying dense arrays of capillaries. Transient activation of neuronal signaling leads to an increased metabolic demand in a tissue region, necessitating increased blood flow, which is achieved by adjustment of vascular diameters in the upstream vessels feeding that region and possibly also in the capillaries and downstream vessels. Open arrows indicate direction of blood flow.
A further challenge in understanding the dynamics of tissue oxygenation is the heterogeneity of structural and hemodynamic parameters in the microcirculation (Pries et al., 1995). The significance of heterogeneous blood flow, resulting in high capillary transit time heterogeneity (CTH), has been highlighted by Jespersen and Østergaard (2012). By introducing variation in capillary velocities in an array of capillaries, they showed that high CTH can lead to reduced oxygen extraction. A thorough analysis of the consequences of heterogeneity requires consideration of microvascular network structure.
The value of theoretical models for gaining quantitative understanding of oxygen transport to tissue has long been recognized (Krogh, 1919). A variety of approaches have been used. In the Krogh cylinder model (Krogh, 1919), the exchange of oxygen between a single capillary and a surrounding cylinder of tissue is analyzed, taking into account the radial gradients of oxygen concentration in the tissue. The application of this approach to cerebral oxygen transport is discussed by Mintun et al. (2001). Other approaches have been based on the assumption that oxygen efflux from capillaries is proportional to plasma oxygen concentration (Buxton and Frank, 1997) or proportional to the difference between intravascular and tissue oxygen level, where the tissue space is represented as a well-mixed compartment, characterized by a specific tissue oxygen level (Jespersen and Østergaard, 2012). In order to represent the contributions of different classes of vessels to oxygen transport, a compartmental model was developed by Sharan et al. (1989).
While these modeling approaches provide useful insights into aspects of cerebral oxygen transport, they do not explicitly represent the effects of network structure and flow distribution. Detailed network-level models are valuable because they permit analysis of the effects of microvascular network architecture and flow distribution on the distribution of oxygen in the tissue (Popel, 1989; Secomb et al., 2000, 2004; Goldman, 2008). Such models require extensive experimental data as input, including information on the three-dimensional structure of microvascular networks containing a large number (typically hundreds or thousands) of segments, together with at least partial information on vessel blood flow rates and oxygen content, and estimates of key oxygen transport parameters. Until recently, however, experimental data available for modeling were mostly from planar or restricted three-dimensional structures, flow rates were generally not measured in the same structures, and data on tissue oxygenation were rarely available.
Recently, high-spatial resolution optical imaging technologies capable of penetrating up to 1 mm into brain tissue have been developed and applied to in vivo brain imaging, providing enhanced capabilities to obtain the data needed for such models. The advantages of optical imaging include high temporal and spatial resolution within reasonable thicknesses of tissue, numerous methods for enhancing contrast, and modest cost of instrumentation. To provide the data needed for network-level modeling of microvascular blood flow and oxygen transport, the imaging modality should meet the following requirements. (i) It should allow imaging of depths of several hundred μm into the tissue, to enable three-dimensional mapping of microvascular structures spanning a substantial fraction of the cortical thickness. (ii) It should have high spatial resolution, sufficient to resolve individual capillaries. (iii) It should allow estimation of blood flow rates in individual microvessels. (iv) It should allow estimation of spatially resolved oxygen levels in microvessels and/or the surrounding tissue. (v) For studies of neurovascular coupling, it should have high temporal resolution (~1 s), sufficient to resolve the time course of hemodynamic changes. In this article, we review state-of-the-art optical microscopy technologies for quantitative in vivo imaging of cerebral microvascular structure, blood flow, and oxygenation. We then review strategies for network-level modeling of cerebral microcirculation using data obtained from such in vivo imaging. These methods will contribute to future investigations of neurovascular and neurometabolic coupling and provide the essential bridge for translation from the microscopic to the network/systems level needed for human translation.
Optical Imaging of Cerebral Microvascular Structure, Hemodynamics and Oxygen Transport
Optical Imaging of Cerebral Microvascular Structure
Optical microscopy techniques that have been applied in vivo to image the cerebral microvasculature in rodents include multi-photon laser scanning microscopy (MPM), optical coherence tomography (OCT), and photoacoustic imaging (PAI). Other “methods to watch” are light-sheet (Keller and Ahrens, 2015) microscopy and super-resolution ultrasound (US) (Errico et al., 2015). The former has been initially targeted to relatively small and transparent organisms and tissue samples (Huisken et al., 2004) but was recently adapted for imaging of neurovascular activity in the mouse cerebral cortex in vivo albeit with limited depth penetration (Bouchard et al., 2015). The latter started as a slow imaging modality (Viessmann et al., 2013; Christensen-Jeffries et al., 2015), but the most recent report achieved mapping of the vasculature in the whole rat brain with ~8 μm resolution within tens of seconds (Errico et al., 2015).
Among these techniques, MPM is the most mature and widely applied for in vivo acquisition of high- resolution cortical angiograms (Helmchen and Denk, 2005; Svoboda and Yasuda, 2006). MPM detects photoluminescence (fluorescence or phosphorescence) from endogenous or exogenous chromophores, and images are typically created from the sequences of sequential point measurements. Brain imaging with MPM can be performed through the thinned skull in mice (Drew et al., 2010), although more commonly a glass-covered cranial window is used allowing increased penetration depth into the brain. Images with a high signal-to-noise ratio are routinely obtained down to a 600–700 μm cortical depth by commercially available MPM setups, and recent technological advances have allowed imaging at a depth of >1 mm through the entire depth of the mouse cortex (Kobat et al., 2011; Horton et al., 2013). Spatial resolution of MPM in vivo imaging is typically 1 μm or better, which enables accurate imaging of all microvessels, including capillaries. However, the method is relatively slow, typically needing tens of minutes to acquire a three-dimensional angiogram over a ~1 mm field of view. Shadows beneath the large cortical surface vessels can cause discontinuities in the acquired data, and optical wave front distortions can impair imaging at larger depths. An example of a mouse cortical microangiogram obtained in vivo by MPM is shown in Figure 2 (Dorr et al., 2012). This study examined changes in cortical microvascular structure and function in a transgenic mouse model of Alzheimer's disease. Structural microvascular impairment (increased tortuosity and decreased diameter of the penetrating arterioles) was associated with accumulation of the amyloid-β and rescued by the administration of scyllo-inositol.
Figure 2. MPM microangiography. In vivo multi-photon fluorescence microscopy images of the cortical microvasculature in a 3-month old non-transgenic mouse. Acquired in the primary somatosensory cortex, this maximum intensity projection image (A) shows penetrating arterioles running vertically down from the cortical surface (top of the bottom row image) as well as ascending venules interspersed with a dense capillary network. (B) A tubular model obtained by segmentation of the branching vessel network shown in (A). The color bar indicates the coding of average vessel diameters. Top row images are parallel with the cortical surface, while bottom row images are perpendicular to the cortical surface. Adapted with permission from Dorr et al. (2012).
In comparison with MPM, other modalities for high-resolution in vivo imaging of cortical microvasculature have been explored less. Among the other modalities, OCT has significant potential (Wang et al., 2007; Vakoc et al., 2009; Srinivasan et al., 2010). An optical analog of ultrasonic pulse-echo imaging (Drexler and Fujimoto, 2008), OCT detects scattering within the tissue to form cross-sectional images, and does not require exogenous contrast agents. Spatial resolution generally ranges from ~1 μm to a few tens of μm. Signals reflected from multiple locations along the optical axis (A-scans) can be acquired simultaneously, which gives OCT a significant speed advantage over MPM in volumetric imaging, with OCT microangiography taking several seconds compared with several minutes for MPM in comparable volumes of tissue. The typical penetration depth in rodent brain is ~1 mm (Srinivasan et al., 2011). The high-speed volumetric imaging is achieved by utilizing low numerical aperture (NA) imaging objectives, providing long depths of focus but reduced lateral spatial resolution. While OCT is currently subject to this tradeoff between lateral resolution and imaging depth, advances in hardware and processing algorithms may alleviate this limitation in the future (Leitgeb et al., 2006; Grulkowski et al., 2014; Mo et al., 2015). Shadows beneath large cortical surface vessels and optical wave front distortions at larger imaging depths are also present in OCT angiograms. In a study conducted by Vakoc et al. (2009), OCT microangiography was utilized to assess angiogenesis induced by a brain tumor. OCT provided microangiograms comparable to those obtained by MPM, including microvascular diameters (Figure 3).
Figure 3. (A) The depth-projected vasculature within the first 2 mm of mouse brain bearing a xenotransplanted U87 human glioblastoma multiforme tumor imaged with OFDI. Depth is denoted by color: yellow (superficial) to red (deep). Scale bar, 500 μm. Right panels: Validation of the morphological measurements obtained from OFDI and MPM. (B,C) Normal brain vasculature acquired by OFDI (B) and MPM (C). The automated vascular tracing was applied to registered vascular data sets to quantify the resolution of OFDI angiography and (D) validation of the morphological measurements obtained from OFDI. Scale bars, 250 μm. Adapted with permission from Vakoc et al. (2009).
PAI's contrast mechanism comes from the absorption of the excitation light by endogenous or exogenous chromophores, thereby creating heat. Images are formed by detecting ultrasonic waves generated by thermoelastic tissue expansion (Wang and Gao, 2014; Yao and Wang, 2014) The high-resolution application of PAI (optical-resolution photoacoustic microscopy, OR-PAM) relies on tight focusing of excitation light (Maslov et al., 2008) and can achieve a spatial resolution of several μm with ~1 mm penetration depth in a rodent brain (Hu et al., 2009; Wang et al., 2013a). The detection of ultrasonic waves in OR-PAM allows rapid acquisition of A-scans and, analogous to OCT, generation of three-dimensional images from two-dimensional scans.
Relative to in vivo microangiography, ex vivo imaging uses more invasive sample preparation and imaging approaches, and permits imaging of larger microvascular networks at high resolution (Yuan et al., 2015). Such methods include physical sample sectioning combined with confocal microscopy (Cassot et al., 2006; Lee et al., 2009) or two-photon microscopy (Blinder et al., 2013), tissue optical clearing (Chung and Deisseroth, 2013), and vascular casting with subsequent chemical tissue clearing and imaging with scanning electron microscopy (Konerding et al., 2001) or ionizing radiation (Plouraboue et al., 2004). Disadvantages of ex vivo microangiography include possible changes in vascular diameters compared to the in vivo state, difficulty in achieving complete filing of the vasculature when labeling the intravascular lumens, and spatial distortions due to sample preparation and/or physical sectioning.
Segmentation of Cerebral Microvascular Angiograms
Regardless of the imaging method used, the resulting data set typically consists of a volumetric intensity map. To provide a basis for analyzing functional properties such as blood flow or oxygen transport, this map must be converted into a defined vascular structure (Plouraboue et al., 2004; Wang et al., 2007; Reichold et al., 2009; Blinder et al., 2013; Yuan et al., 2015). Typically, this requires two main steps. First, the intravascular space must be distinguished from the extravascular space by segmentation, usually achieved setting an intensity threshold. From the results of this step, basic parameters such as vascular volume and vascular surface area can be deduced. Second, the intravascular space must be represented as a set of connected vessel segments, with identified lengths, diameters and nodal connection points. This information is required for the analysis of blood flow in the network, and provides a basis for modeling oxygen transport to tissue. Numerous vascular segmentation algorithms have been developed in recent years (Fridman et al., 2004; Kirbas and Quek, 2004; McIntosh and Hamarneh, 2006; Tyrrell et al., 2007; Shikata et al., 2009; Tsai et al., 2009; Rennie et al., 2011; Fraz et al., 2012; Hong et al., 2014). Detailed discussion of these algorithms is beyond the scope of this review.
Although the process of segmentation of the volumetric microvascular intensity images is simple in principle, it is subject to a number of practical difficulties. The spatial resolution of the imaging system is an important factor. If vessel boundaries are not sharply resolved, estimates of vessel diameters may be subject to substantial uncertainty. Furthermore, the choice of intensity threshold affects the region identified as being intravascular. According to Poiseuille's law, the flow rate in a vessel segment resulting from a given pressure drop varies as the fourth power of diameter. Predicted blood flow rates in microvascular networks are therefore very sensitive to estimated diameters. Variations in image contrast or intensity across the region of interest can lead to systematic deviations. In regions of low image intensity, some segments may fall below the chosen threshold, leading to incomplete mapping of the network and the appearance of multiple blind-ended vessels. The type of fluorescent labeling used for red blood cells, blood plasma or endothelial cells must also be considered when estimating vessel diameters. If a macromolecular plasma marker is used, the apparent diameter may be less than the anatomical diameter because large molecules may be excluded from the endothelial surface layer or glycocalyx, a layer up to about 1 μm thick of glycoproteins, polysaccharides and proteins that covers the inner surface of the vasculature (Pries et al., 2000; Vink and Duling, 2000).
The process of describing the vascular network in terms of segments and nodes, based on a three-dimensional map of the intravascular space, presents further challenges, particularly if the map is not of high quality. Available segmentation algorithms are unable to achieve high accuracy in an automated mode, and manual intervention by the user is often needed to achieve satisfactory results. This is highly time-consuming for cortical volumes such as those obtained in vivo by MPM (~700 × 700 × 700 μm3). Development of improved methods for analysis of three-dimensional in vivo microvascular images is a pressing need (Wang et al., 2013a; Leahy et al., 2015).
Imaging of Cerebral Microvessel Blood Flows
Information on microvessel blood flow rates is essential for the analysis of oxygen delivery to cerebral tissue. Measurement of flow in all segments of extensive microvascular networks is currently not feasible. However, flow measurements in a subset of segments can be combined with theoretical models of blood flow in microvascular networks, permitting estimation of flow in all segments, as discussed below.
MPM can be used for high-resolution depth-resolved measurements of vessel diameters, RBC fluxes and flow velocities in small animal models (Kleinfeld et al., 1998; Kamoun et al., 2010; Shih et al., 2012; Figure 4). MPM utilizes fluorescent labeling of either plasma or RBCs and tracks RBC speed and flux in individual microvessels, including capillaries. There is a tradeoff between the number of simultaneously monitored vessels and temporal resolution. Fast blood flow transients that are characteristic of functional hyperemia or spreading depressions and require ~1 s imaging temporal resolution can be continuously monitored in only a few vessel segments at a time, while all vessels in the field of view can be sampled with ~12 s temporal resolution (Kamoun et al., 2010). However, multiplexing strategies for MPM are being explored (Ducros et al., 2013; Lecoq et al., 2014; Yang et al., 2016).
Figure 4. TPM flow. Simultaneous measurement of lumen diameter and red blood cell velocity in multiple vessels using spatially optimized line scans. (A) Image of fluorescein-dextran-labeled vessels in the rat somatosensory cortex taken with a 4X objective. (B) High-magnification image of a surface arteriole and venule in the forelimb region collected with a 40X objective. The pattern of the two-photon scanning laser is superimposed. Portions of the scan path along the length of the vessel are used to calculate red blood cell (RBC) velocity, whereas portions moving across the width of the vessels are used to calculate lumen diameter (Driscoll et al., 2011). Scans were acquired at a rate of 733 scan cycles per second. (C) The scan path is colored to show the error between the desired scan path and the actual path the mirrors traversed. The error along linear portions of the image is typically < 1 μm and increases when the mirrors undergo rapid acceleration. The error between successive scans of the same path is < 0.15 μm, several times lower than the point-spread function of a two-photon laser scanning microscopy. (D) The upper traces show the scan path and mirror speed as a function of time. Note that portions used to acquire diameter and velocity data are constant speed. The line scans generated from the path can be stacked sequentially as a space–time plot as shown in the lower image panels. Each image panel shows ~100 ms of data collected before, during, and after an electrical stimulation of the contralateral forelimb of the anesthetized rat. The stimulus was a 1 mA current, delivered for 3 s at 3 Hz with a 100-ms pulse width (Devor et al., 2007). (E) Vessel diameter is calculated as the full-width at half-maximum of a time-average of several scans across the width of a vessel (left). Red blood cell velocity is calculated from the angle of the RBC streaks (right; Drew et al., 2010). (F) Data traces of lumen diameter, RBC velocity, and RBC flux for the arteriole and venule after processing to remove heartbeat and smoothing with a running window. Figure adapted from Driscoll et al. (2011, Book chapter).
OCT enables depth-resolved imaging of the absolute blood flow in individual cortical arterioles and venules (Wang et al., 2007; Wang and An, 2009; Srinivasan et al., 2010; Bouwens et al., 2013; Lee et al., 2013, 2014) as well as measurements of RBC flux in capillaries (Ren et al., 2012a; Srinivasan et al., 2012; Lee et al., 2013; Tokayer et al., 2013; Weiss et al., 2013; Figure 5). Compared to MPM measurements of blood flow, the advantages of OCT include increased penetration depth through the thinned skull in mice and >1 mm penetration depth through a cranial window, reliance on endogenous contrast (i.e., optical scattering) instead of exogenous contrast agents, and improved acquisition speed. OCT measurements of cortical blood flow have been validated (Srinivasan et al., 2011) and full volumetric imaging of blood flow over a cortical surface area of 1 mm2 is possible in ~1 min (Srinivasan et al., 2010). In addition to CBF measurements in non-capillary vessels, significant progress was made in the quantitative OCT measurement of capillary blood flow (Srinivasan et al., 2012; Lee et al., 2013, 2014). Volumetric images of capillary RBC flux can be obtained in a few minutes, opening the possibility of studying capillary perfusion in the brain and providing extensive empirical data for theoretical models of RBC flux distribution within capillaries. In combination with MPM, this technology will enable studies of the role of pericytes in CBF control. Commercial systems are available, facilitating widespread adoption of OCT. However, the data analysis for the OCT blood flow measurements is significantly more complex than for the MPM flow measurements. Therefore, widespread adoption of OCT flow measurements may also require the availability of sophisticated data processing tools. In short, while OCT has only recently emerged as a tool for studying brain CBF (Jia et al., 2009; Srinivasan et al., 2011; Ren et al., 2012b; Shen et al., 2014), its importance is expected to grow rapidly.
Figure 5. OCT RBC flow. (A) The enface MIP of the 3D angiogram with color indicating the depth from the cortical surface. Bar = 100 μm. (B–D) Estimated RBC speed, flux, and density are presented as color spots on the MIP angiogram. We used relatively comparable ranges (median ± 40%) for all quantities. (E) Histograms (top) and depth profiles (bottom) of the RBC speed, flux, and density measured from three animals (n = 2; 259 measures in total). Comparable histogram ranges (median ± 40%) are used. In the depth profiles, data are presented as mean ± s.d. (F) Correlations between the flux vs. the speed and density (n = 2259). Dashed lines indicate median ± 40% of the density (left) and speed (right). Adapted with permission from Lee et al. (2013).
PAI can be also applied for measuring blood flow. By utilizing this technology, the characteristics of cerebral blood flow (CBF) can be measured using several different approaches, including blood velocity measurement in arterioles and venules based on Doppler broadening of the ultrasonic bandwidth (Yao et al., 2010), combining ultrasonic thermal tagging with PAI (Wang et al., 2013b), and even single RBC velocity measurement (Wang et al., 2013a). A combination of PAI measurements of cortical blood flow and oxygen saturation may allow estimation of local oxygen consumption rates.
Imaging of Cerebral Oxygenation
Theoretical models can be used to predict the distribution of tissue oxygen levels surrounding a microvascular network, as discussed later. However, direct measurement of tissue and/or blood oxygen levels is needed to provide boundary conditions for such models, to allow estimation of unknown parameters, and to test the resulting predictions. A large body of experimental studies have established the basic properties of intra- and extravascular oxygen changes during increases in neuronal activity and associated functional hyperemia (vasodilation and an increase in CBF and volume) (Vanzetta and Grinvald, 1999; Vovenko, 1999; Ances et al., 2001; Erecińska and Silver, 2001; Masamoto et al., 2003; Thompson et al., 2003; Offenhauser et al., 2005; Viswanathan and Freeman, 2007; Sharan et al., 2008; Yaseen et al., 2009; Vazquez et al., 2010). However, only recently several optical microscopy imaging technologies were developed opening the door for microscopic in vivo imaging of intravascular and tissue oxygenation with unprecedented spatial and temporal resolution within cortex.
Two-photon microscopy-based phosphorescence lifetime imaging (PLIO2) allows mapping intravascular and tissue partial pressure of O2 (PO2) with micrometer spatial resolution (Finikova et al., 2008; Sakadžić et al., 2010; Lecoq et al., 2011; Figure 6). The technique measures oxygen-dependent phosphorescence lifetimes of an exogenous contrast agent (Vanderkooi et al., 1987; Rumsey et al., 1988). The phosphorescence lifetime of a probe depends on the PO2 in the immediate vicinity of the probe, providing a spatially localized measurement of dissolved oxygen. Probe molecules for MPM were specially designed for two-photon excitation, with a high degree of encapsulation that ensures stability of the lifetime calibration in a complex biological environment (Finikova et al., 2008; Lebedev et al., 2009; Roussakis et al., 2014). Unlike spectroscopy-based hemoglobin saturation measurements, PLIO2 lifetime imaging is insensitive to changes in tissue optical properties during imaging. The acquisition speed is currently limited to 0.2–0.5 s per measurement point by relatively long phosphorescence lifetimes and the number of decay averages required at each point. Further development of oxygen sensitive dyes with significantly higher quantum yield, two-photon absorption cross section, and dynamic range will enable significant improvements in acquisition speed and precision of PO2 imaging (Esipova and Vinogradov, 2014). This technology has been used to obtain high-resolution maps of oxygen concentration distribution in both the microvasculature and tissue under various conditions (Devor et al., 2011; Kazmi et al., 2013; Parpaleix et al., 2013; Sakadžić et al., 2014; Spencer et al., 2014) and the awake mouse brain (Lyons et al., 2016).
Figure 6. Two-photon microscopy imaging of PO2. (A) Maximum intensity projection (MIP) of the 200−μm-thick cortical microvascular stack obtained by TPM. Blood plasma was labeled by FITC. (B) Top-view projection of the segmented microvasculature with the intravascular PO2 measurements obtained by TPM. Mean vascular segment PO2 measurements were color-coded and overlaid on the segmented microvascular structure. Scale bars, 200 μm. Adapted with permission from (Sakadžić et al., 2014). (C) Simultaneous measurement of PO2 in cortical vasculature and tissue. Measured PO2 values overlaid with the gray scale phosphorescence intensity image at 60 μm depth. Measurements were performed at the location of an ascending venule. Measurement location is marked with the white rectangle in the inset (bottom right), showing MIP of 80 μm-thick FITC-labeled microvasculature stack. Ascending venule (blue) and branches of the descending arteriole (red) are color-coded for easier identification. Scale bars, 50 μm. Adapted with permission from Sakadžić et al. (2010).
Intravascular oxygenation can be also measured with PAI because hemoglobin is the dominant optical absorber of visible and near-infrared radiation in the brain. PAI provides quantification of both oxygen saturation and total hemoglobin concentration (Figure 7). In particular, high spatial resolution variants of PAI can achieve ~1 mm penetration depth in the brain and were successfully applied to measure arteriolar, venular and capillary oxygen saturation (SO2) as well as single RBC oxygenation (Wang et al., 2013a; Yao and Wang, 2014). Recently, high acquisition rate PAI was used to image the hemodynamic response to hind limb stimulation in the mouse somatosensory cortex (Yao et al., 2015). PAI has also been applied to image brain pathologies including edema after cold injury (Xu et al., 2011), and dynamics of the microvascular oxygenation after middle cerebral artery occlusion (Hu et al., 2011) and during epileptic seizure (Tsytsarev et al., 2013). The ability of PAI to measure RBC oxygen saturation and flux in capillaries has potential applications for analyzing capillary flow and oxygenation under various conditions.
Figure 7. Fast functional photoacoustic microscopy (PAM) of the mouse brain. (A) Representative xy projected brain vasculature image through an intact skull. (B) Representative enhanced xz projected brain vasculature image acquired over a 0.6 × 0.6 mm2 region with depth scanning, where the signal amplitude was normalized depthwise. (C) PAM of oxygen saturation of hemoglobin (SO2) in the same mouse brain as in (A), acquired by using the single-wavelength pulse-width-based method (PW-SO2) with two lasers. SV, skull vessel. Adapted with permission from Yao et al. (2015).
Another imaging technique, visible light OCT (vis-OCT), is emerging as an alternative that may provide rapid non-invasive measurements of microvascular oxygenation (Robles et al., 2011; Pan et al., 2014; Chong et al., 2015a). This will permit measurements of tissue and vascular structure, blood flow and oxygen transport, all obtained using OCT (Figure 8).
Figure 8. OCT measurement of SO2. Quantification of chromophores in the mouse brain in an en face view. (A) Maximum intensity projection of R2 values from the fit (Equation 19) shows the highest values near the centers of vessels, with a decrease at the edges. (B) The parameter Φ accounts for RBC scattering effects. (C) Saturation map, showing clear distinctions between arteries and veins. (D) Map of the maximum of the product of oxygenated hemoglobin concentration and distance shows that veins and arteries contain oxyhemoglobin. (E) By comparison, under the given experimental conditions, most of deoxyhemoglobin is contained in the veins. (F) The map of the maximum of the product of total hemoglobin concentration and distance shows larger values in larger vessels, with localized increases at vessel crossings. It should be noted that quantitative measurements of chromophores can be achieved by integrating the maps (D–F) in the transverse plane (x and y dimensions). All maps were displayed with transparency based on the local R2 values at each transverse location, averaged over depth. An artery (a) and vein (v) are labeled. Adapted with permission from Chong et al. (2015b).
Theoretical Modeling of Cerebral Microvascular Hemodynamics and Oxygen Transport
Blood Flow Simulation
Analysis of oxygen transport in a microvascular network requires knowledge of blood flow and red blood cell (RBC) fluxes in each vessel segment, since they govern the convective transport of oxygen. To model microvascular hemodynamics, the network is generally represented as a set of interconnected cylindrical segments, in which each segment is characterized by its conductance, defined as the ratio of volume flow rate (Q) to pressure drop (ΔP) (Lipowsky and Zweifach, 1974). According to Poiseuille's law, the conductance is given by πd4/(128Lμ), where d is the diameter, L is the length and μ is the viscosity of blood. This formula shows that conductance is strongly dependent on vessel diameter: a small increase in diameter leads to a significant flow increase. In network flow simulations, small inaccuracies in diameters derived from image processing can therefore significantly affect the predicted flow distributions. The effective viscosity of blood, μ, has been determined experimentally, and varies strongly with vessel diameter as a result of the particulate nature of blood (the Fåhraeus-Lindqvist effect). Empirical equations derived from experimental observations give apparent blood viscosity as a function of vessel diameter and hematocrit (Pries et al., 1992, 1994). These equations take into account the inhomogeneous characteristics of blood flow in microvessels, including the formation of a cell-free or cell-depleted layer near vessel walls, resulting from the fact that it is a concentrated suspension of RBCs. The effects of the endothelial surface layer adjacent to the vessel walls are also included.
As a result of the particulate nature of blood, the partition of RBC flux at diverging microvessel bifurcations generally differs from the partition of plasma flow. This implies that the hematocrit in each daughter vessel may differ from that in the parent vessel. Empirical equations have been developed to describe this “phase separation” effect, including its dependence on the flow split in the bifurcation, the vessel diameters, and the discharge hematocrit in the parent vessel (Pries et al., 1989; Pries and Secomb, 2005). The discharge hematocrit (HD) is defined as the volume flux of RBCs as a fraction of the total volume flow rate, and should be distinguished from the tube hematocrit (HT), which is defined as the fraction of the vessel volume occupied by RBCs. Generally HD is larger than HT, because RBCs preferentially flow near the vessel center-line where the velocity is higher, a phenomenon known as the Fåhraeus effect. According to the above definitions, the volume flux of RBCs is given by QHD, where Q is the volume flow rate of the segment. The empirical equations referred to above for phase separation refer to situations where a single vessel splits into exactly two branches. An alternative approach, which has the advantage of applicability to branch points with more than two outflowing segments, has been proposed (Gould and Linninger, 2015b). This approach excludes the possibility of zero hematocrit in low-flow branches, whereas such behavior can be observed in vivo, and fits experimental data less closely than the earlier empirical equations (Pries et al., 1989; Gould and Linninger, 2015a). A further alternative approach has been proposed based on the assumption that the distribution of RBCs in diverging bifurcations satisfies an optimality principle (Sriram et al., 2014).
When network models of this type are applied to observed microvascular network structures, the problem arises that many vessels typically cross the boundaries of the sampled region. To obtain a unique flow solution, it is necessary to impose flow or pressure boundary conditions (BCs) on these inflowing and outflowing segments, and hematocrit boundary conditions on inflowing segments. Not all these flow or pressure BCs are generally available from experimental observations, although a subset may be available. Several approaches have been used to address this difficulty. Lorthois et al. (2011a,b) considered two regimes: (i) uniform pressure BCs applied to boundary capillaries, with a pressure value chosen to give zero net flux into the capillaries; (ii) zero flow conditions imposed at boundary capillaries. These two cases gave upper and lower bounds on the total network perfusion. The latter approach was also used by Gagnon et al. (2015a,b). Fry et al. (2012) developed an approach to address this problem by solving for the flows throughout the network while minimizing the deviation from an expected “target” pressure and shear stress in each vessel segment. These target values must be estimated based on observations of typical values in the tissue under consideration. In an analysis of flow distribution in cortical networks (Gagnon et al., 2015a), pressure boundary conditions were imposed using literature values based on vessel diameters. The computation was constrained using measurements of flow in a subset of vessels (Figure 9). This was shown to reduce the impact of pressure boundary conditions and estimated vessel resistances. The error in flow computation was reduced by 50% when flow measurements in major outflowing vessels were included in the set of constraints.
Figure 9. Reconstruction of microvascular cerebral blood flow in the mouse cortex using Doppler OCT flow measurements in outflowing venules to constrain the computation. (A) MPM angiogram. (B) Doppler OCT velocity projection map. (C) Reconstructed flow. Adapted with permission from Gagnon et al. (2015a).
Models for Oxygen Transport
Oxygen transport takes place over a range of spatial scales and involves multiple biophysical processes. The challenge of developing theoretical models for oxygen transport to tissue has stimulated a variety of approaches. Here we briefly review the development of the field, focusing on cerebral microcirculation. For a more detailed review, see Goldman (2008).
A key early contribution was a model inspired by the somewhat regular structure of skeletal muscle (Krogh, 1919). The Krogh cylinder model assumes an array of parallel, evenly spaced capillaries, each of which supplies oxygen to a tissue cylinder surrounding it, with a fixed rate of consumption. This leads to a simple steady-state reaction-diffusion equation, through which the radial variation in tissue PO2 can be computed as a function of distance from the vessel. The Krogh model has been extended to include other effects including time-dependent transport, PO2-dependent oxygen consumption according to Michaelis-Menten kinetics, myoglobin-facilitated tissue transport and intravascular resistance to radial oxygen diffusion. Hudetz et al. (1999) adapted the Krogh model to study oxygenation in the cerebral cortex, assuming the capillary supplied a conical tissue region that tapered downstream, such that venous capillaries supply smaller regions. This model exhibited a non-linear relation between consumption and blood flow, i.e., for a significant rise in consumption rate, the blood flow increase required to avoid hypoxia was disproportionately larger.
Using a finite-element model of oxygen transport by discrete moving RBCs, Lücker et al. (2015) analyzed the transient variation in tissue oxygenation due to the passage of RBCs and the difference in PO2 between RBCs and plasma. These variations have been termed erythrocyte-associated transients (EAT). The authors suggested that a higher capillary density on the venular side of the network must be present to ensure tissue oxygenation. Network-scale models of oxygenation including transient effects associated with each RBC would be very expensive computationally. An approximate approach based on an effective blood PO2 is generally used, as described in the next section.
In general, Krogh-type approaches are appropriate for muscle, particularly under conditions of high oxygen demand, but do not provide a good approximation in situations where the geometry is not approximated by parallel capillaries. The Krogh model overestimates minimum tissue PO2 and thus fails to predict the onset of hypoxia in such geometries (Secomb et al., 1993, 2000). Therefore, approaches have been developed for modeling oxygen delivery by networks of vessels. One relatively simple approximate approach is to treat the extravascular space as a well-mixed compartment with respect to oxygen transport, with the resistance to oxygen transport residing in the vessel wall (Jespersen and Østergaard, 2012; Angleys et al., 2015). Finite element or finite difference methods have been employed to simulate oxygenation in arrays of parallel capillaries (Hoofd, 1995; Goldman and Popel, 2000; Beard and Bassingthwaighte, 2001), networks with interconnected parallel capillaries (Wieringa et al., 1993; Goldman and Popel, 2000; Beard and Bassingthwaighte, 2001), synthetic networks generated based on measured distributions of anatomical properties (Lauwers et al., 2008; Linninger et al., 2013; Safaeian and David, 2013; Park and Payne, 2016) and networks extracted from three-dimensional angiograms (Fang et al., 2008; Gagnon et al., 2015b). An efficient computational method applicable to arbitrary vascular network geometries has been developed using a Green's function approach (Hsu and Secomb, 1989; Secomb et al., 2004), as discussed later.
Intravascular Oxygen Transport
The rate of oxygen transport along a microvessel by convection is given by:
where Pb is the PO2 in the blood, S is the fractional oxyhemoglobin saturation, HD is the discharge hematocrit, C0 is the concentration of hemoglobin-bound oxygen in a fully saturated RBC and αpl and αrbc are the solubilities of oxygen in plasma and RBCs respectively. This may be simplified to give:
where αeff, the effective solubility of oxygen in blood, depends on the discharge hematocrit HD and on the solubilities in plasma and red blood cells. Under normal conditions, the majority of transported oxygen is bound to hemoglobin in the RBCs. The contribution of dissolved oxygen to convective oxygen transport is normally small, but may become significant if Pb is very high or HD is very low. The binding of oxygen to hemoglobin is a complex process which depends on several factors including blood pH, temperature and CO2 concentration. However, the saturation can be fairly accurately represented by the Hill equation:
where P50 is the PO2 at 50% saturation, which is in the range 29–40.5 mmHg depending on the species, and n is the Hill exponent, which is typically in the range 2–3 (Ellsworth et al., 1988; Uchida et al., 1998; Ellis et al., 2002). The value of P50 is important in determining the critical range of PO2 values for oxygen unbinding, due to the sharp gradient in the Hill equation curve around this value. Conservation of mass implies that:
where s is distance along the segment and qv is the rate of diffusive oxygen efflux per unit vessel length. In simulations of oxygen delivery by microvascular networks, the above equations are applied to each segment in the network.
The above equation for convective oxygen flux J assumes that the blood PO2 at each point along the vessel is described by a single value, Pb. In reality, significant gradients in PO2 occur within the vessel cross-section, associated with the radial diffusion of oxygen from blood to tissue (Hellums, 1977). Therefore, the quantity Pb represents an effective average blood PO2, such that the convective oxygen flux is correctly specified.
Oxygen Diffusion and Consumption in Tissue
If the diffusivity D and solubility α of oxygen in the tissue are uniform, the tissue PO2, P(x, y, z) satisfies:
under steady-state conditions, where ∇2 is the Laplacian operator and M(P) is the oxygen consumption rate. Experimentally, oxygen consumption rates are approximately constant above a certain tissue PO2 value, but drop sharply below this. This behavior is commonly represented by Michaelis-Menten kinetics:
where M0 is the rate when oxygen is not rate-limiting (maximal consumption rate), and P0 is the PO2 value at half-maximal oxygen consumption. In previous work, this has usually been assumed to be about 1 mmHg, but recent experimental results put it closer to 5–10 mmHg (Golub and Pittman, 2012).
At the blood-tissue interface, the tissue oxygen field P(x, y, z) must satisfy boundary conditions representing continuity of diffusive oxygen flux and continuity of partial pressure. Continuity of flux implies that:
where r is radial distance from the vessel centerline, a is the vessel radius and the integral is around the vessel circumference. The matching of PO2 values at the interface requires consideration of the radial intracapillary gradients, mentioned earlier. Since these gradients drive the radial flux of oxygen, they can be assumed to be proportional to the rate of oxygen loss from the vessel. This relationship can be expressed as
where Pv(s) is the tissue PO2 averaged around the circumference of the vessel and K represents intravascular resistance to radial oxygen transport. This resistance parameter shows a strong dependence on vessel diameter (Hellums et al., 1996) and depends on hematocrit (Roy and Secomb, 2014). As a result of these matching conditions, the equations for oxygen transport by convection in the vessels and by diffusion in the tissue are strongly coupled and must be solved simultaneously. This is a challenging problem when large networks containing many segments with a complex geometry are considered.
Green's Function Methods
One approach to address this challenge is to represent the tissue oxygen field as a superposition of fields resulting from a distribution of oxygen sources (the vessel segments) and oxygen sinks (consumption in the tissue). In an infinite domain, under steady-state conditions, the PO2 field resulting from a unit point source of oxygen is given by:
where r is the distance from the source. This is the Green's function for the three-dimensional Laplacian operator. This approach permits computation of spatially resolved oxygen fields for complicated vascular network structures with a reasonable computational cost (Hsu and Secomb, 1989; Secomb et al., 2004) and has been applied to a various tissues including brain (Secomb et al., 2000). The Green's function approach uses an iterative solution method and can be applied when the reaction rate kinetics are non-linear, as in the case of Michaelis-Menten kinetics, or the solute binding characteristics in blood are non-linear, as in the case of the Hill equation for oxygen binding to hemoglobin. An analogous approach has also been developed for time-dependent diffusion problems (Secomb, 2015).
Applications and Future Directions
Estimation of Microvascular Blood Flow
As already described, current imaging approaches can provide three-dimensional maps of extensive network structures, including vessel diameters, together with flow rate measurements in a subset of segments. Theoretical models can then be used to reconstruct the entire flow distribution. The limitation of incomplete flow information for boundary segments can be addressed by constraining flows in selected segments (Figure 9) to match measured values (Gagnon et al., 2015a) or by using an optimization approach based on typical distributions of pressures and shear rates (Fry et al., 2012). Combining these two methods represents a logical future development. The emergence of new imaging modalities for measuring flow in capillaries, such as OCT (Srinivasan et al., 2012; Lee et al., 2013, 2014), MPM (Kleinfeld et al., 1998; Desjardins et al., 2014) and PAM (Wang et al., 2013a) will provide flow data for more vessel segments, and therefore allow more accurate flow reconstructions. Improvements in acquisition speed will allow for dynamic measurements during brain activation, which will permit studies aiming to improve interpretation of fMRI signals that are routinely used in human neuroscience studies (Bouchard et al., 2015).
In a study of blood flow in the microvessel networks of the rat mesentery, Pries et al. (1994) showed that the apparent viscosity of blood differed strongly from that in glass tubes of corresponding diameters. The main cause of this discrepancy was found to be the presence of a relatively thick (~1 μm) glycocalyx or endothelial surface layer (Pries et al., 2000). However, it is not known to what extent this effect occurs in other tissues, including brain. The methods discussed here provide a potential approach for investigating the hemodynamic effects of the ESL in brain microcirculation.
Assessment of Tissue Oxygenation
Overall oxygen transport rates to cortical tissue can be estimated from tissue-level parameters. For instance, the rate of oxygen consumption per unit volume is given by:
where CaO2 is the arterial oxygen concentration and OEF is the oxygen extraction fraction given by:
and SaO2 and SvO2 are the oxygen saturations in arteries and veins, respectively. With the development of cortical measurements of PO2, SO2, and CaO2 with MPM (Sakadžić et al., 2010; Lecoq et al., 2011), PAM (Yao et al., 2013), and OCT (Chong et al., 2015b), SaO2 and SvO2 can be measured in vivo and OEF can be computed from these parameters, allowing assessment of local variations in these quantities.
However, these parameters are not sufficient to allow estimation of tissue oxygen levels. Due to the heterogeneous structure of the microcirculation and the significant gradients in PO2, the spatial frequency distribution of tissue PO2 is much wider than would be the case in an ideal structure such as a Krogh cylinder (Secomb et al., 2000). Theoretical models for oxygen transport are needed to compute the spatially resolved tissue oxygen field (Secomb et al., 2000, 2004; Fang et al., 2008; Linninger et al., 2013; Gagnon et al., 2015b). However, the predictions of such models depend on oxygen consumption rate, oxygen levels in inflowing vessels, vessel flow rates, and other parameters that are not accurately known. A combination of modeling and measurements is needed to constrain these parameters (Gagnon et al., 2015b; Figure 10). The increasing availability of micro-scale measurements of tissue and vessel PO2, combined with the computationally efficient Green's function method for simulating oxygen transport, will make possible increasingly detailed and robust predictions of network-scale tissue PO2 distributions under a range of normal and pathological conditions.
Figure 10. (A) Distribution of the partial pressure of oxygen (PO2) simulated across the vascular network using the FEM model. (B) MPM experimental measurement of PO2 in vivo. (C) Quantitative comparison of simulated and experimental PO2 and SO2 distributions across the vascular network. Adapted with permission from Gagnon et al. (2015b).
Modeling of Neurovascular Coupling and fMRI Signals
On the macroscopic level of vascular dynamics, a number of theoretical models have been developed to describe the active regulation of blood flow to the brain. Among them, Ursino and Lodi (1998) used a compartmental model to analyze the regulation of CBF in response to changes in arterial pressure, intracranial pressure and carbon dioxide level. Other models such as the Balloon model (Buxton et al., 1998) or the double-compliance model (Mandeville et al., 1999) have also been developed.
The mechanisms of blood flow regulation and neurovascular coupling at the microscopic level are still not fully understood. Recently, a number of groups produced controversial evidence regarding the role of pericytes in controlling dilation and constriction of the capillary bed (Hall et al., 2014; Hill et al., 2015). The population of pericytes is heterogeneous, and not all of them are contractile as documented by many groups in various tissues including brain (Armulik et al., 2011; Attwell et al., 2015; Hartmann et al., 2015; Trost et al., 2016). There is a general agreement that high-branching-order capillaries cannot actively dilate or constrict (Hall et al., 2014; Hill et al., 2015). The debate is focused on low order branches off diving arterioles in a mouse cortex and whether these branches should be considered capillaries or precapillary arterioles. Hall at al. classified all branches off diving cerebral arterioles as capillaries based on the presence of pericytes. However, pericytes are present also on top of the smooth muscle wall in arterioles including larger pial arteries. Therefore, capillary definition based on the presence of pericytes may be ambiguous. Two-photon imaging indicates that the first couple of branching orders off diving cerebral arterioles are noticeably thicker that the average high order capillary supporting the existence of precapillary arterioles in between diving arterioles and the capillary bed. Regardless of nomenclature issues, there is a consensus across studies that these low-order branches (which may include true capillaries) can dilate and constrict (Tian et al., 2010; Hall et al., 2014; Uhlirova et al., 2016).
The distinction between precapillary arterioles and capillaries at sites of control is functionally relevant because terminal arterioles typically supply several capillaries, and so do not provide control at the individual capillary level. Also, an arteriole is spatially separated from the capillaries that it feeds, implying the need for mechanisms to coordinate arteriolar contraction with metabolic needs at the capillary level. The modeling approaches described above can be used to investigate the implications of different proposed mechanisms of neurovascular coupling. For example, network models can be used to investigate the effects of local vascular dilation on flow distribution (Reichold et al., 2009). Development of spatially resolved network-level models of neurovascular coupling is at an early stage and represents an important challenge for future work. Complexity arises from the fact that multiple biological mechanisms are involved. The roles of pressure-dependent (myogenic), shear-dependent and metabolic responses, and of conduction of signals along the vessel wall, have been explored in an integrated model for flow regulation in a simplified vascular network in skeletal muscle (Arciero et al., 2008; Carlson et al., 2008). In the brain, the roles that potassium ion fluxes in synaptic regions, astrocytes and arterioles play in vasodilation have been investigated theoretically (Witthoft et al., 2013). Progress in developing spatially resolved models for neurovascular coupling, network blood flow and oxygen transport has potential to lead to better understanding of the mechanisms determining CBF and oxygen transport (Boas et al., 2008; Fang et al., 2008) and to provide a basis for bottom-up modeling of fMRI signals (Gagnon et al., 2015b).
As a step in this direction, the experimental and theoretical methods reviewed here provide a basis for modeling the effects of prescribed changes in vascular diameter on blood flow and tissue oxygenation (Figure 11). For instance, Gagnon et al. (2015b) used experimental time courses of arterial dilation obtained with MPM (Tian et al., 2010) as a basis for simulating the resulting changes in blood flow and tissue oxygenation. Given the dynamics of CBF in the vasculature and the time course of CMRO2 in the cortical tissue, the dynamic distribution of SO2 was computed (Figure 11). Good agreement was obtained between the simulations and the experimental confocal microscopic measurements of SO2 during functional activation. The results of such simulations can be extended to model fMRI signals from first principles with Monte Carlo simulations of proton diffusion. This approach allowed quantification of the compartmental microvascular origin of BOLD-fMRI for different magnet strengths and MR pulse sequences with high accuracy. This approach also predicted that the BOLD response measured in human studies is influenced by the local folding of the cortex and its orientation in the magnet, a prediction that was confirmed by measuring the amplitude of the BOLD response to hypercapnia for different cortical orientations (Gagnon et al., 2015b).
Figure 11. (A) Spatiotemporal evolution of SO2 changes following a 2 s forepaw stimulus. (B) Comparison of simulated SO2 changes (n = 6 animals) with experimental SO2 changes (n = 10 animals) measured in pial vessels during a forepaw stimulus with confocal microscopy. Adapted with permission from Gagnon et al. (2015b).
Progress in understanding neurovascular function (or its dysfunction in disease) requires experimental technologies and computational tools to accurately account for spatially distributed and dynamic processes of vasodilation/constriction and O2 transport. These tools—on their own—will not be sufficient for reconstruction of the underlying neuronal activity. Eventual progress toward this goal, however, will not take place without understanding the normal vascular network dynamics and how it is altered in disease. In the future, coupling these realistic 3D vascular network models to large-scale models of cortical circuits (Markram et al., 2015) would enable testing specific mechanistic hypotheses of functional hyperemia.
All authors listed, have made substantial, direct and intellectual contribution to the work, and approved it for publication.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer DMT declared a shared affiliation, though no other collaboration, with one of the authors AD to the handling Editor, who ensured that the process nevertheless met the standards of a fair and objective review.
Supported by NIH grants NS091230, EB018464, NS092986, NS055104, EB021018 and HL070657.
Ances, B. M., Greenberg, J. H., and Detre, J. A. (2001). The effects of graded hypercapnia on the activation flow coupling response due to forepaw stimulation in alpha-chloralose anesthetized rats. Brain Res. 911, 82–88. doi: 10.1016/S0006-8993(01)02721-4
Angleys, H., Østergaard, L., and Jespersen, S. N. (2015). The effects of capillary transit time heterogeneity (CTH) on brain oxygenation. J. Cereb. Blood Flow Metab. 35, 806–817. doi: 10.1038/jcbfm.2014.254
Arciero, J. C., Carlson, B. E., and Secomb, T. W. (2008). Theoretical model of metabolic blood flow regulation: roles of ATP release by red blood cells and conducted responses. Am. J. Physiol. Heart Circ. Physiol. 295, H1562–H1571. doi: 10.1152/ajpheart.00261.2008
Armulik, A., Genové, G., and Betsholtz, C. (2011). Pericytes: developmental, physiological, and pathological perspectives, problems, and promises. Dev. Cell 21, 193–215. doi: 10.1016/j.devcel.2011.07.001
Blinder, P., Tsai, P. S., Kaufhold, J. P., Knutsen, P. M., Suhl, H., and Kleinfeld, D. (2013). The cortical angiome: an interconnected vascular network with noncolumnar patterns of blood flow. Nat. Neurosci. 16, 889–897. doi: 10.1038/nn.3426
Boas, D. A., Jones, S. R., Devor, A., Huppert, T. J., and Dale, A. M. (2008). A vascular anatomical network model of the spatio-temporal response to brain activation. Neuroimage 40, 1116–1129. doi: 10.1016/j.neuroimage.2007.12.061
Bouchard, M. B., Voleti, V., Mendes, C. S., Lacefield, C., Grueber, W. B., Mann, R. S., et al. (2015). Swept confocally-aligned planar excitation (scape) microscopy for high speed volumetric imaging of behaving organisms. Nat. Photon. 9, 113–119. doi: 10.1038/nphoton.2014.323
Bouwens, A., Szlag, D., Szkulmowski, M., Bolmont, T., Wojtkowski, M., and Lasser, T. (2013). Quantitative lateral and axial flow imaging with optical coherence microscopy and tomography. Opt. Express 21, 17711–17729. doi: 10.1364/OE.21.017711
Buxton, R. B., and Frank, L. R. (1997). A model for the coupling between cerebral blood flow and oxygen metabolism during neural stimulation. J. Cereb. Blood Flow Metab. 17, 64–72. doi: 10.1097/00004647-199701000-00009
Buxton, R. B., Griffeth, V. E. M., Simon, A. B., Moradi, F. A., and Shmuel, A. (2014). Variability of the coupling of blood flow and oxygen metabolism responses in the brain: a problem for interpreting BOLD studies but potentially a new window on the underlying neural activity. Front. Neurosci. 8:139. doi: 10.3389/fnins.2014.00139
Buxton, R. B., Wong, E. C., and Frank, L. R. (1998). Dynamics of blood flow and oxygenation changes during brain activation: the balloon model. Magn. Reson. Med. 39, 855–864. doi: 10.1002/mrm.1910390602
Carlson, B. E., Arciero, J. C., and Secomb, T. W. (2008). Theoretical model of blood flow autoregulation: roles of myogenic, shear-dependent, and metabolic responses. Am. J. Physiol. Heart Circ. Physiol. 295, H1572–H1579. doi: 10.1152/ajpheart.00262.2008
Cassot, F., Lauwers, F., Fouard, C., Prohaska, S., and Lauwers-Cances, V. (2006). A novel three-dimensional computer-assisted method for a quantitative study of microvascular networks of the human cerebral cortex. Microcirculation 13, 1–18. doi: 10.1080/10739680500383407
Chong, S. P., Merkle, C. W., Leahy, C., Radhakrishnan, H., and Srinivasan, V. J. (2015b). Quantitative microvascular hemoglobin mapping using visible light spectroscopic Optical Coherence Tomography. Biomed. Opt. Express 6, 1429–1422. doi: 10.1364/BOE.6.001429
Chong, S. P., Merkle, C. W., Leahy, C., and Srinivasan, V. J. (2015a). Cerebral metabolic rate of oxygen (CMRO_2) assessed by combined Doppler and spectroscopic OCT. Biomed. Opt. Express 6, 3941–3911. doi: 10.1364/BOE.6.003941
Christensen-Jeffries, K., Browning, R. J., Tang, M. X., Dunsby, C., and Eckersley, R. J. (2015). In vivo acoustic super-resolution and super-resolved velocity mapping using microbubbles. IEEE Trans. Med. Imaging 34, 433–440. doi: 10.1109/TMI.2014.2359650
Desjardins, M., Berti, R., Lefebvre, J., Dubeau, S., and Lesage, F. (2014). Aging-related differences in cerebral capillary blood flow in anesthetized rats. Neurobiol. Aging 35, 1947–1955. doi: 10.1016/j.neurobiolaging.2014.01.136
Devor, A. A., Sakadzic, S. S., Saisan, P. A. P., Yaseen, M. A. M., Roussakis, E. E., Srinivasan, V. J. V., et al. (2011). “Overshoot” of O2 is required to maintain baseline tissue oxygenation at locations distal to blood vessels. J. Neurosci. 31, 13676–13681. doi: 10.1523/JNEUROSCI.1968-11.2011
Devor, A., Boas, D. A., Einevoll, G. T., Buxton, R. B., and Dale, A. M. (2012). “Neuronal basis of non-invasive functional imaging: from microscopic neurovascular dynamics to BOLD fMRI,” in Neural Metabolism In vivo Advances in Neurobiology, eds I.-Y. Choi and R. Gruetter (Boston, MA: Springer), 433–500.
Devor, A., Tian, P., Nishimura, N., Teng, I. C., Hillman, E. M. C., Narayanan, S. N., et al. (2007). Suppressed neuronal activity and concurrent arteriolar vasoconstriction may explain negative blood oxygenation level-dependent signal[s]. J. Neurosci. 27, 4452–4459. doi: 10.1523/JNEUROSCI.0134-07.2007
Dorr, A., Sahota, B., Chinta, L. V., Brown, M. E., Lai, A. Y., Ma, K., et al. (2012). Amyloid-β-dependent compromise of microvascular structure and function in a model of Alzheimer's disease. Brain 135, 3039–3050. doi: 10.1093/brain/aws243
Drew, P. J., Shih, A. Y., Driscoll, J. D., Knutsen, P. M., Blinder, P., Davalos, D., et al. (2010). Chronic optical access through a polished and reinforced thinned skull. Nat. Meth. 7, 981–984. doi: 10.1038/nmeth.1530
Driscoll, J. D., Shih, A. Y., Drew, P. J., Valmianski, I., Cauwenberghs, G., and Kleinfeld, D. (2011). “Two-photon imaging of blood flow in rat cortex,” in Imaging in Neuroscience: A Laboratory Manual (Book 2), Vol. 84, eds F. Helmchen A. Konnerth, and R. Yuste (Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press), 927–938.
Ellis, C. G., Bateman, R. M., Sharpe, M. D., Sibbald, W. J., and Gill, R. (2002). Effect of a maldistribution of microvascular blood flow on capillary O(2) extraction in sepsis. Am. J. Physiol. Heart Circ. Physiol. 282, H156–H164. doi: 10.1186/cc1589
Ellsworth, M. L., Popel, A. S., and Pittman, R. N. (1988). Assessment and impact of heterogeneities of convective oxygen transport parameters in capillaries of striated muscle: experimental and theoretical. Microvasc. Res. 35, 341–362. doi: 10.1016/0026-2862(88)90089-1
Errico, C., Pierre, J., Pezet, S., Desailly, Y., Lenkei, Z., Couture, O., et al. (2015). Ultrafast ultrasound localization microscopy for deep super-resolution vascular imaging. Nature 527, 499–502. doi: 10.1038/nature16066
Fang, Q., Sakadžić, S., Ruvinskaya, L., Devor, A., Dale, A. M., and Boas, D. A. (2008). Oxygen advection and diffusion in a three- dimensional vascular anatomical network. Opt. Express 16, 17530–17541. doi: 10.1364/OE.16.017530
Finikova, O. S., Lebedev, A. Y., Aprelev, A., Troxler, T., Gao, F., Garnacho, C., et al. (2008). Oxygen microscopy by two-photon-excited phosphorescence. Chemphyschem 9, 1673–1679. doi: 10.1002/cphc.200800296
Fraz, M. M., Remagnino, P., Hoppe, A., Uyyanonvara, B., Rudnicka, A. R., Owen, C. G., et al. (2012). Blood vessel segmentation methodologies in retinal images–a survey. Comput. Methods Prog. Biomed. 108, 407–433. doi: 10.1016/j.cmpb.2012.03.009
Gagnon, L., Sakadžić, S., Lesage, F., Mandeville, E. T., Fang, Q., Yaseen, M. A., et al. (2015a). Multimodal reconstruction of microvascular-flow distributions using combined two-photon microscopy and Doppler optical coherence tomography. Neurophoton 2:015008. doi: 10.1117/1.NPh.2.1.015008
Gagnon, L., Sakadžić, S., Lesage, F., Musacchia, J. J., Lefebvre, J., Fang, Q., et al. (2015b). Quantifying the microvascular origin of BOLD-fMRI from first principles with two-photon microscopy and an oxygen-sensitive nanoprobe. J. Neurosci. 35, 3663–3675. doi: 10.1523/JNEUROSCI.3555-14.2015
Girouard, H., and Iadecola, C. (2006). Neurovascular coupling in the normal brain and in hypertension, stroke, and Alzheimer disease. J. Appl. Physiol. 100, 328–335. doi: 10.1152/japplphysiol.00966.2005
Goldman, D., and Popel, A. S. (2000). A computational study of the effect of capillary network anastomoses and tortuosity on oxygen transport. J. Theor. Biol. 206, 181–194. doi: 10.1006/jtbi.2000.2113
Gordon, G. R. J., Choi, H. B., Rungta, R. L., Ellis-Davies, G. C. R., and MacVicar, B. A. (2008). Brain metabolism dictates the polarity of astrocyte control over arterioles. Nature 456, 745–749. doi: 10.1038/nature07525
Grulkowski, I., Szulzycki, K., and Wojtkowski, M. (2014). Microscopic OCT imaging with focus extension by ultrahigh-speed acousto-optic tunable lens and stroboscopic illumination. Opt. Express 22, 31746–31760. doi: 10.1364/OE.22.031746
Hall, C. N., Klein-Flügge, M. C., Howarth, C., and Attwell, D. (2012). Oxidative phosphorylation, not glycolysis, powers presynaptic and postsynaptic mechanisms underlying brain information processing. J. Neurosci. 32, 8940–8951. doi: 10.1523/JNEUROSCI.0026-12.2012
Hall, C. N., Reynell, C., Gesslein, B., Hamilton, N. B., Mishra, A., Sutherland, B. A., et al. (2014). Capillary pericytes regulate cerebral blood flow in health and disease. Nature 508, 55–60. doi: 10.1038/nature13165
Hartmann, D. A., Underly, R. G., Grant, R. I., Watson, A. N., Lindner, V., and Shih, A. Y. (2015). Pericyte structure and distribution in the cerebral cortex revealed by high-resolution imaging of transgenic mice. Neurophotonics 2:041402. doi: 10.1117/1.NPh.2.4.041402
Hill, R. A., Tong, L., Yuan, P., Murikinati, S., Gupta, S., and Grutzendler, J. (2015). Regional blood flow in the normal and ischemic brain is controlled by arteriolar smooth muscle cell contractility and not by capillary pericytes. Neuron 87, 95–110. doi: 10.1016/j.neuron.2015.06.001
Hirsch, S., Reichold, J., Schneider, M., Székely, G., and Weber, B. (2012). Topology and hemodynamics of the cortical cerebrovascular system. J. Cereb. Blood Flow Metab. 32, 952–967. doi: 10.1038/jcbfm.2012.39
Hoofd, L. (1995). Calculation of oxygen pressures in tissue with anisotropic capillary orientation. I. Two-dimensional analytical solution for arbitrary capillary characteristics. Math. Biosci. 129, 1–23. doi: 10.1016/0025-5564(94)00049-6
Horton, N. G., Wang, K., Kobat, D., Clark, C. G., Wise, F. W., Schaffer, C. B., et al. (2013). In vivo three-photon microscopy of subcortical structures within an intact mouse brain. Nature Photon. 7, 205–209. doi: 10.1038/nphoton.2012.336
Hudetz, A. G., Wood, J. D., Biswal, B. B., Krolo, I., and Kampine, J. P. (1999). Effect of hemodilution on RBC velocity, supply rate, and hematocrit in the cerebral capillary network. J. Appl. Physiol. 87, 505–509.
Huisken, J., Swoger, J., Del Bene, F., Wittbrodt, J., and Stelzer, E. H. K. (2004). Optical sectioning deep inside live embryos by selective plane illumination microscopy. Science 305, 1007–1009. doi: 10.1126/science.1100035
Jespersen, S. N., and Østergaard, L. (2012). The roles of cerebral blood flow, capillary transit time heterogeneity, and oxygen tension in brain oxygenation and metabolism. J. Cereb. Blood Flow Metab. 32, 264–277. doi: 10.1038/jcbfm.2011.153
Jia, Y., Alkayed, N., and Wang, R. K. (2009). Potential of optical microangiography to monitor cerebral blood perfusion and vascular plasticity following traumatic brain injury in mice in vivo. J. Biomed. Opt. 14, 040505. doi: 10.1117/1.3207121
Kamoun, W. S., Chae, S.-S., Lacorre, D. A., Tyrrell, J. A., Mitre, M., Gillissen, M. A., et al. (2010). Simultaneous measurement of RBC velocity, flux, hematocrit and shear rate in vascular networks. Nat. Meth. 7, 655–660. doi: 10.1038/nmeth.1475
Kazmi, S. M. S., Parthasarthy, A. B., Song, N. E., Jones, T. A., and Dunn, A. K. (2013). Chronic imaging of cortical blood flow using multi-exposure speckle imaging. J. Cereb. Blood Flow Metab. 33, 798–808. doi: 10.1038/jcbfm.2013.57
Kleinfeld, D., Blinder, P., Drew, P. J., Driscoll, J. D., Muller, A., Tsai, P. S., et al. (2011). A guide to delineate the logic of neurovascular signaling in the brain. Front. Neuroenerg. 3:1. doi: 10.3389/fnene.2011.00001
Kleinfeld, D., Mitra, P. P., Helmchen, F., and Denk, W. (1998). Fluctuations and stimulus-induced changes in blood flow observed in individual capillaries in layers 2 through 4 of rat neocortex. Proc. Natl. Acad. Sci. U.S.A. 95, 15741–15746. doi: 10.1073/pnas.95.26.15741
Konerding, M. A., Fait, E., and Gaumann, A. (2001). 3D microvascular architecture of pre-cancerous lesions and invasive carcinomas of the colon. Br. J. Cancer 84, 1354–1362. doi: 10.1054/bjoc.2001.1809
Krogh, A. (1919). The number and distribution of capillaries in muscles with calculations of the oxygen pressure head necessary for supplying the tissue. J. Physiol. 52, 409–415. doi: 10.1113/jphysiol.1919.sp001839
Lauwers, F., Cassot, F., Lauwers-Cances, V., Puwanarajah, P., and Duvernoy, H. (2008). Morphometry of the human cerebral cortex microcirculation: general characteristics and space-related profiles. Neuroimage 39, 936–948. doi: 10.1016/j.neuroimage.2007.09.024
Leahy, C., Radhakrishnan, H., Weiner, G., Goldberg, J. L., and Srinivasan, V. J. (2015). Mapping the 3D connectivity of the rat inner retinal vascular network using OCT angiography. Invest. Ophthalmol. Vis. Sci. 56, 5785–5793. doi: 10.1167/iovs.15-17210
Lebedev, A. Y., Cheprakov, A. V., Sakadžić, S., Boas, D. A., Wilson, D. F., and Vinogradov, S. A. (2009). Dendritic phosphorescent probes for oxygen imaging in biological systems. ACS Appl. Mater. Interfaces 1, 1292–1304. doi: 10.1021/am9001698
Lecoq, J., Parpaleix, A., Roussakis, E., Ducros, M., Houssen, Y. G., Vinogradov, S. A., et al. (2011). Simultaneous two-photon imaging of oxygen and blood flow in deep cerebral vessels. Nat. Med. 17, 893–898. doi: 10.1038/nm.2394
Lecoq, J., Savall, J., Vučinić, D., Grewe, B. F., Kim, H., Li, J. Z., et al. (2014). Visualizing mammalian brain area interactions by dual-axis two-photon calcium imaging. Nat. Neurosci. 17, 1825–1829. doi: 10.1038/nn.3867
Lecrux, C., Toussay, X., Kocharyan, A., Fernandes, P., Neupane, S., Lévesque, M., et al. (2011). Pyramidal neurons are “neurogenic hubs” in the neurovascular coupling response to whisker stimulation. J. Neurosci. 31, 9836–9847. doi: 10.1523/JNEUROSCI.4943-10.2011
Lee, J., Jiang, J. Y., Wu, W., Lesage, F., and Boas, D. A. (2014). Statistical intensity variation analysis for rapid volumetric imaging of capillary network flux. Biomed. Opt. Express 5:1160. doi: 10.1364/BOE.5.001160
Lee, J., Niederer, S., Nordsletten, D., Le Grice, I., Smaill, B., Kay, D., Smith, N., et al. (2009). Coupling contraction, excitation, ventricular and coronary blood flow across scale and physics in the heart. Philos. Trans. Math. Phys. Eng. Sci. 367, 2311–2331. doi: 10.1098/rsta.2008.0311
Lee, J., Wu, W., Lesage, F., and Boas, D. A. (2013). Multiple-capillary measurement of RBC speed, flux, and density with optical coherence tomography. J. Cereb. Blood Flow Metab. 33, 1707–1710. doi: 10.1038/jcbfm.2013.158
Leitgeb, R. A., Villiger, M., Bachmann, A. H., Steinmann, L., and Lasser, T. (2006). Extended focus depth for Fourier domain optical coherence microscopy. Opt. Lett. 31, 2450–2452. doi: 10.1364/OL.31.002450
Lin, A.-L., Fox, P. T., Hardies, J., Duong, T. Q., and Gao, J.-H. (2010). Nonlinear coupling between cerebral blood flow, oxygen consumption, and ATP production in human visual cortex. Proc. Natl. Acad. Sci. U.S.A. 107, 8446–8451. doi: 10.1073/pnas.0909711107
Linninger, A. A., Gould, I. G., Marinnan, T., Hsu, C.-Y., Chojecki, M., and Alaraj, A. (2013). Cerebral microcirculation and oxygen tension in the human secondary cortex. Ann. Biomed. Eng. 41, 2264–2284. doi: 10.1007/s10439-013-0828-0
Lorthois, S., Cassot, F., and Lauwers, F. (2011a). Simulation study of brain blood flow regulation by intra-cortical arterioles in an anatomically accurate large human vascular network. Part II: flow variations induced by global or localized modifications of arteriolar diameters. Neuroimage 54, 2840–2853. doi: 10.1016/j.neuroimage.2010.10.040
Lorthois, S., Cassot, F., and Lauwers, F. (2011b). Simulation study of brain blood flow regulation by intra-cortical arterioles in an anatomically accurate large human vascular network: Part I: Methodology and baseline flow. Neuroimage 54, 12–12. doi: 10.1016/j.neuroimage.2010.09.032
Lücker, A., Weber, B., and Jenny, P. (2015). A dynamic model of oxygen transport from capillaries to tissue with moving red blood cells. Am. J. Physiol. Heart Circ. Physiol. 308, H206–H216. doi: 10.1152/ajpheart.00447.2014
Mandeville, J. B., Marota, J. J., Ayata, C., Zaharchuk, G., Moskowitz, M. A., Rosen, B. R., et al. (1999). Evidence of a cerebrovascular postarteriole windkessel with delayed compliance. J. Cereb. Blood Flow Metab. 19, 679–689. doi: 10.1097/00004647-199906000-00012
Markram, H., Muller, E., Ramaswamy, S., Reimann, M. W., Abdellah, M., Sanchez, C. A., et al. (2015). Reconstruction and simulation of neocortical microcircuitry. Cell 163, 456–492. doi: 10.1016/j.cell.2015.09.029
Masamoto, K., Omura, T., Takizawa, N., Kobayashi, H., Katura, T., Maki, A., et al. (2003). Biphasic changes in tissue partial pressure of oxygen closely related to localized neural activity in guinea pig auditory cortex. J. Cereb. Blood Flow Metab. 23, 1075–1084. doi: 10.1097/01.WCB.0000084248.20114.B3
McIntosh, C., and Hamarneh, G. (2006). “Vessel crawlers: 3D physically-based deformable organisms for vasculature segmentation and analysis,” in 2006 IEEE Computer Society Conference on Computer Vision and Pattern Recognition, Vol. 1 (New York, NY), 1084–1091.
Mintun, M. A., Lundstrom, B. N., Snyder, A. Z., Vlassenko, A. G., Shulman, G. L., and Raichle, M. E. (2001). Blood flow and oxygen delivery to human brain during functional activity: theoretical modeling and experimental data. Proc. Natl. Acad. Sci. U.S.A. 98, 6859–6864. doi: 10.1073/pnas.111164398
Mo, J., de Groot, M., and de Boer, J. F. (2015). Depth-encoded synthetic aperture optical coherence tomography of biological tissues with extended focal depth. Opt. Express 23, 4935–4945. doi: 10.1364/OE.23.004935
Offenhauser, N., Thomsen, K., Caesar, K., and Lauritzen, M. (2005). Activity-induced tissue oxygenation changes in rat cerebellar cortex: interplay of postsynaptic activation and blood flow. J. Physiol. 565, 279–294. doi: 10.1113/jphysiol.2005.082776
Pan, Y., You, J., Volkow, N. D., Park, K., and Du, C. (2014). Ultrasensitive detection of 3D cerebral microvascular network dynamics in vivo. Neuroimage 103, 492–501. doi: 10.1016/j.neuroimage.2014.08.051
Plouraboue, F., Cloetens, P., Fonta, C., Steyer, A., Lauwers, F., and Marc-Vergnes, J. P. (2004). X-ray high-resolution vascular network imaging. J. Microsc. 215, 139–148. doi: 10.1111/j.0022-2720.2004.01362.x
Pries, A. R., and Secomb, T. W. (2005). Microvascular blood viscosity in vivo and the endothelial surface layer. Am. J. Physiol. Heart Circ. Physiol. 289, H2657–H2664. doi: 10.1152/ajpheart.00297.2005
Pries, A. R., Secomb, T. W., Gessner, T., Sperandio, M. B., Gross, J. F., and Gaehtgens, P. (1994). Resistance to blood flow in microvessels in vivo. Circ. Res. 75, 904–915. doi: 10.1161/01.RES.75.5.904
Reichold, J., Stampanoni, M., Lena Keller, A., Buck, A., Jenny, P., and Weber, B. (2009). Vascular graph model to simulate the cerebral blood flow in realistic vascular networks. J. Cereb. Blood Flow Metab. 29, 1429–1443. doi: 10.1038/jcbfm.2009.58
Ren, H., Du, C., Park, K., Volkow, N. D., and Pan, Y. (2012a). Quantitative imaging of red blood cell velocity in vivo using optical coherence Doppler tomography. Appl. Phys. Lett. 100, 233702–2337024. doi: 10.1063/1.4726115
Ren, H., Du, C., Yuan, Z., Park, K., Volkow, N. D., and Pan, Y. (2012b). Cocaine-induced cortical microischemia in the rodent brain: clinical implications. Mol. Psychiatry 17, 1017–1025. doi: 10.1038/mp.2011.160
Rennie, M. Y., Detmar, J., Whiteley, K. J., Yang, J., Jurisicova, A., Adamson, S. L., et al. (2011). Vessel tortuousity and reduced vascularization in the fetoplacental arterial tree after maternal exposure to polycyclic aromatic hydrocarbons. Am. J. Physiol. Heart Circ. Physiol. 300, H675–H684. doi: 10.1152/ajpheart.00510.2010
Rumsey, W. L., Vanderkooi, J. M., and Wilson, D. F. (1988). Imaging of phosphorescence: a novel method for measuring oxygen distribution in perfused tissue. Science 241, 1649–1651. doi: 10.1126/science.3420417
Safaeian, N., and David, T. (2013). A computational model of oxygen transport in the cerebrocapillary levels for normal and pathologic brain function. J. Cereb. Blood Flow Metab. 33, 1633–1641. doi: 10.1038/jcbfm.2013.119
Sakadžić, S., Mandeville, E. T., Gagnon, L., Musacchia, J. J., Yaseen, M. A., Yücel, M. A., et al. (2014). Large arteriolar component of oxygen delivery implies a safe margin of oxygen supply to cerebral tissue. Nat. Comms 5, 5734. doi: 10.1038/ncomms6734
Sakadžić, S., Roussakis, E., Yaseen, M. A., Mandeville, E. T., Srinivasan, V. J., Arai, K., et al. (2010). Two-photon high-resolution measurement of partial pressure of oxygen in cerebral vasculature and tissue. Nat. Meth. 7, 755–759. doi: 10.1038/nmeth.1490
Secomb, T. W. (2015). A Green's function method for simulation of time-dependent solute transport and reaction in realistic microvascular geometries. Math. Med. Biol. doi: 10.1093/imammb/dqv031. [Epub ahead of print].
Secomb, T. W., Hsu, R., Beamer, N. B., and Coull, B. M. (2000). Theoretical simulation of oxygen transport to brain by networks of microvessels: effects of oxygen supply and demand on tissue hypoxia. Microcirculation 7, 237–247. doi: 10.1111/j.1549-8719.2000.tb00124.x
Secomb, T. W., Hsu, R., Dewhirst, M. W., Klitzman, B., and Gross, J. F. (1993). Analysis of oxygen transport to tumor tissue by microvascular networks. Int. J. Radiat. Oncol. Biol. Phys. 25, 481–489. doi: 10.1016/0360-3016(93)90070-C
Secomb, T. W., Hsu, R., Park, E. Y. H., and Dewhirst, M. W. (2004). Green's function methods for analysis of oxygen delivery to tissue by microvascular networks. Ann. Biomed. Eng. 32, 1519–1529. doi: 10.1114/B:ABME.0000049036.08817.44
Sharan, M., Jones, M. D., Koehler, R. C., Traystman, R. J., and Popel, A. S. (1989). A compartmental model for oxygen transport in brain microcirculation. Ann. Biomed. Eng. 17, 13–38. doi: 10.1007/BF02364271
Sharan, M., Vovenko, E. P., Vadapalli, A., Popel, A. S., and Pittman, R. N. (2008). Experimental and theoretical studies of oxygen gradients in rat pial microvessels. J. Cereb. Blood Flow Metab. 28, 1597–1604. doi: 10.1038/jcbfm.2008.51
Shen, H.-Y., Sun, H., Hanthorn, M. M., Zhi, Z., Lan, J.-Q., Poulsen, D. J., et al. (2014). Overexpression of adenosine kinase in cortical astrocytes and focal neocortical epilepsy in mice. J. Neurosurg. 120, 628–638. doi: 10.3171/2013.10.JNS13918
Shih, A. Y., Driscoll, J. D., Drew, P. J., Nishimura, N., Schaffer, C. B., and Kleinfeld, D. (2012). Two-photon microscopy as a tool to study blood flow and neurovascular coupling in the rodent brain. J. Cereb. Blood Flow Metab. 32, 1277–1309. doi: 10.1038/jcbfm.2011.196
Shikata, H., McLennan, G., Hoffman, E. A., and Sonka, M. (2009). Segmentation of pulmonary vascular trees from thoracic 3D CT images. Int. J. Biomed. Imaging 2009, 636240–636211. doi: 10.1155/2009/636240
Spencer, J. A., Ferraro, F., Roussakis, E., Klein, A., Wu, J., Runnels, J. M., et al. (2014). Direct measurement of local oxygen concentration in the bone marrow of live animals. Nature 508, 269–273. doi: 10.1038/nature13034
Srinivasan, V. J., Atochin, D. N., Radhakrishnan, H., Jiang, J. Y., Ruvinskaya, S., Wu, W., et al. (2011). Optical coherence tomography for the quantitative study of cerebrovascular physiology. J. Cereb. Blood Flow Metab. 31, 1339–1345. doi: 10.1038/jcbfm.2011.19
Srinivasan, V. J., Radhakrishnan, H., Lo, E. H., Mandeville, E. T., Jiang, J. Y., Barry, S., et al. (2012). OCT methods for capillary velocimetry. Biomed. Opt. Express 3, 612–629. doi: 10.1364/BOE.3.000612
Srinivasan, V. J., SakadŽiæ, S., Gorczynska, I., Ruvinskaya, S., Wu, W., Fujimoto, J. G., et al. (2010). Quantitative cerebral blood flow with optical coherence tomography. Opt. Express 18, 2477–2494. doi: 10.1364/OE.18.002477
Sriram, K., Intaglietta, M., and Tartakovsky, D. M. (2014). Hematocrit dispersion in asymmetrically bifurcating vascular networks. Am. J. Physiol. Heart Circ. Physiol. 307, H1576–H1586. doi: 10.1152/ajpheart.00283.2014
Tian, P., Teng, I. C., May, L. D., Kurz, R., Lu, K., Scadeng, M., et al. (2010). Cortical depth-specific microvascular dilation underlies laminar differences in blood oxygenation level-dependent functional MRI signal. Proc. Natl. Acad. Sci. U.S.A. 107, 15246–15251. doi: 10.1073/pnas.1006735107
Tokayer, J., Jia, Y., Dhalla, A.-H., and Huang, D. (2013). Blood flow velocity quantification using split-spectrum amplitude-decorrelation angiography with optical coherence tomography. Biomed. Opt. Express 4, 1909–1924. doi: 10.1364/BOE.4.001909
Trost, A., Lange, S., Schroedl, F., Bruckner, D., Motloch, K. A., Bogner, B., et al. (2016). Brain and retinal pericytes: origin, function and role. Front. Cell. Neurosci. 10:20. doi: 10.3389/fncel.2016.00020
Tsai, P. S., Kaufhold, J. P., Blinder, P., Friedman, B., Drew, P. J., Karten, H. J., et al. (2009). Correlations of neuronal and microvascular densities in murine cortex revealed by direct counting and colocalization of nuclei and vessels. J. Neurosci. 29, 14553–14570. doi: 10.1523/JNEUROSCI.3287-09.2009
Tsytsarev, V., Rao, B., Maslov, K. I., Li, L., and Wang, L. V. (2013). Photoacoustic and optical coherence tomography of epilepsy with high temporal and spatial resolution and dual optical contrasts. J. Neurosci. Methods 216, 142–145. doi: 10.1016/j.jneumeth.2013.04.001
Tyrrell, J. A., di Tomaso, E., Fuja, D., Tong, R., Kozak, K., Jain, R. K., et al. (2007). Robust 3-D modeling of vasculature imagery using superellipsoids. IEEE Trans. Med. Imaging 26, 223–237. doi: 10.1109/TMI.2006.889722
Uhlirova, H., Kiliç, K., Tian, P., Thunemann, M., Desjardins, M., Saisan, P. A., et al. (2016). Cell type specificity of neurovascular coupling in cerebral cortex. Elife 5:e14315. doi: 10.7554/eLife.14315
Vakoc, B. J., Lanning, R. M., Tyrrell, J. A., Padera, T. P., Bartlett, L. A., Stylianopoulos, T., et al. (2009). Three-dimensional microscopy of the tumor microenvironment in vivo using optical frequency domain imaging. Nat. Med. 15, 1219–1223. doi: 10.1038/nm.1971
Vanderkooi, J. M., Maniara, G., Green, T. J., and Wilson, D. F. (1987). An optical method for measurement of dioxygen concentration based upon quenching of phosphorescence. J. Biol. Chem. 262, 5476–5482.
Vanzetta, I., and Grinvald, A. (1999). Increased cortical oxidative metabolism due to sensory stimulation: implications for functional brain imaging. Science 286, 1555–1558. doi: 10.1126/science.286.5444.1555
Viessmann, O. M., Eckersley, R. J., Christensen-Jeffries, K., Tang, M. X., and Dunsby, C. (2013). Acoustic super-resolution with ultrasound and microbubbles. Phys. Med. Biol. 58, 6447–6458. doi: 10.1088/0031-9155/58/18/6447
Vovenko, E. (1999). Distribution of oxygen tension on the surface of arterioles, capillaries and venules of brain cortex and in tissue in normoxia: an experimental study on rats. Pflugers Arch. 437, 617–623. doi: 10.1007/s004240050825
Wang, L., Xia, J., Yao, J., Maslov, K. I., and Wang, L. V. (2013b). Ultrasonically encoded photoacoustic flowgraphy in biological tissue. Phys. Rev. Lett. 111:204301. doi: 10.1103/PhysRevLett.111.204301
Weiss, N., van Leeuwen, T. G., and Kalkman, J. (2013). Localized measurement of longitudinal and transverse flow velocities in colloidal suspensions using optical coherence tomography. Phys. Rev. E. 88:042312. doi: 10.1103/physreve.88.042312
Yang, W., Miller, J.-E. K., Carrillo-Reid, L., Pnevmatikakis, E., Paninski, L., Yuste, R., et al. (2016). Simultaneous multi-plane imaging of neural circuits. Neuron 89, 269–284. doi: 10.1016/j.neuron.2015.12.012
Yao, J., Maslov, K. I., Shi, Y., Taber, L. A., and Wang, L. V. (2010). In vivo photoacoustic imaging of transverse blood flow by using Doppler broadening of bandwidth. Opt. Lett. 35, 1419–1421. doi: 10.1364/OL.35.001419
Yao, J., Wang, L., Yang, J.-M., Maslov, K. I., Wong, T. T. W., Li, L., et al. (2015). High-speed label-free functional photoacoustic microscopy of mouse brain in action. Nat. Meth. 12, 407–410. doi: 10.1038/nmeth.3336
Yao, J., Xia, J., Maslov, K. I., Nasiriavanaki, M., Tsytsarev, V., Demchenko, A. V., et al. (2013). Noninvasive photoacoustic computed tomography of mouse brain metabolism in vivo. Neuroimage 64, 257–266. doi: 10.1016/j.neuroimage.2012.08.054
Yaseen, M. A., Srinivasan, V. J., Sakadžić, S., Wu, W., Ruvinskaya, S., Vinogradov, S. A., et al. (2009). Optical monitoring of oxygen tension in cortical microvessels with confocal microscopy. Opt. Express 17, 22341–22350. doi: 10.1364/OE.17.022341
Keywords: cerebral blood flow (CBF), cerebral blood flow measurement, cerebrovascular circulation, brain imaging methods, modeling and simulations
Citation: Gagnon L, Smith AF, Boas DA, Devor A, Secomb TW and Sakadžić S (2016) Modeling of Cerebral Oxygen Transport Based on In vivo Microscopic Imaging of Microvascular Network Structure, Blood Flow, and Oxygenation. Front. Comput. Neurosci. 10:82. doi: 10.3389/fncom.2016.00082
Received: 03 March 2016; Accepted: 25 July 2016;
Published: 31 August 2016.
Edited by:Franck Plouraboué, Centre National de la Recherche Scientifique, France
Reviewed by:Daniel M. Tartakovsky, University of California, San Diego, USA
Andreas A. Linninger, University of Illinois at Chicago, USA
Copyright © 2016 Gagnon, Smith, Boas, Devor, Secomb and Sakadžić. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Louis Gagnon, firstname.lastname@example.org