Development of an image-based model for capillary vasculature of retina

Simulation of murine retinal hemodynamics in response to tail suspension

The etiology of spaceflight-associated neuro-ocular syndrome (SANS) remains unclear. Recent murine studies indicate there may be a link between the space environment and retinal endothelial dysfunction. Post-fixed control (N = 4) and 14-day tail-suspended (TS) (N = 4) mice eye samples were stained and imaged for the vessel plexus and co-located regions of endothelial cell death. A custom workflow combined whole-mounted and tear reconstructed three-dimensional (3D) spherical retinal plexus models with computational fluid dynamics (CFD) simulation that accounted for the Fåhræus-Lindqvist effect and boundary conditions that accommodated TS fluid pressure measurements and deeper capillary layer blood flow distribution. TS samples exhibited reduced surface area (4.6 ± 0.5 mm2 vs. 3.5 ± 0.3 mm2, P = 0.010) and shorter lengths between branches in small vessels (<10 μm, 69.5 ± 0.6 μm vs. 60.4 ± 1.1 μm, P < 0.001). Wall shear stress (WSS) and pressure were higher in TS mice compared to controls, particularly in smaller vessels (<10 μm, WSS: 6.57 ± 1.08 Pa vs. 4.72 ± 0.67 Pa, P = 0.034, Pressure: 72.04 ± 3.14 mmHg vs. 50.64 ± 6.74 mmHg, P = 0.004). Rates of retinal endothelial cell death were variable in TS mice compared to controls. WSS and pressure were generally higher in cell death regions, both within and between cohorts, but significance was variable and limited to small to medium-sized vessels (<20 μm). These findings suggest a link may exist between emulated microgravity and retinal endothelial dysfunction that may have implications for SANS development. Future work with increased sample sizes of larger species or spaceflight cohorts should be considered.

Methods to measure blood flow and vascular reactivity in the retina

Disturbances of retinal perfusion are involved in the onset and maintenance of several ocular diseases, including diabetic retinopathy, glaucoma, and retinal vascular occlusion. Hence, knowledge on ocular vascular anatomy and function is highly relevant for basic research studies and for clinical judgment and treatment. The retinal vasculature is composed of the superficial, intermediate, and deep vascular layer. Detection of changes in blood flow and vascular diameter especially in smaller vessels is essential to understand and to analyze vascular diseases. Several methods to evaluate blood flow regulation in the retina have been described so far, but no gold standard has been established. For highly reliable assessment of retinal blood flow, exact determination of vessel diameter is necessary. Several measurement methods have already been reported in humans. But for further analysis of retinal vascular diseases, studies in laboratory animals, including genetically modified mice, are important. As for mice, the small vessel size is challenging requiring devices with high optic resolution. In this review, we recapitulate different methods for retinal blood flow and vessel diameter measurement. Moreover, studies in humans and in experimental animals are described.

Ocular Fluid Mechanics and Drug Delivery: A Review of Mathematical and Computational Models

The human eye is a complex biomechanical structure with a range of biomechanical processes involved in various physiological as well as pathological conditions. Fluid flow inside different domains of the eye is one of the most significant biomechanical processes that tend to perform a wide variety of functions and when combined with other biophysical processes play a crucial role in ocular drug delivery. However, it is quite difficult to comprehend the effect of these processes on drug transport and associated treatment experimentally because of ethical constraints and economic feasibility. Computational modeling on the other hand is an excellent means to understand the associated complexity between these aforementioned processes and drug delivery. A wide range of computational models specific to different types of fluids present in different domains of the eye as well as varying drug delivery modes has been established to understand the fluid flow behavior and drug transport phenomenon in an insilico manner. These computational models have been used as a non-invasive tool to aid ophthalmologists in identifying the challenges associated with a particular drug delivery mode while treating particular eye diseases and to advance the understanding of the biomechanical behavior of the eye. In this regard, the author attempts to summarize the existing computational and mathematical approaches proposed in the last two decades for understanding the fluid mechanics and drug transport associated with different domains of the eye, together with their application to modify the existing treatment processes.

Age dependence of retinal vascular plexus attenuation in the triple transgenic mouse model of Alzheimer's disease

The influence of Alzheimer's disease (AD) progression and severity on the structural and functional integrity of the cerebral vasculature is well recognized. The retina is an extension of the brain; thus, changes in retinal vascular features may serve as markers of AD cerebrovascular pathologies. However, differentiating normal aging-versus AD-induced retinal vascular changes is unresolved. Therefore, we compared and quantified changes in superficial (SVP), intermediate (IVP), and deep (DVP) retinal vascular plexuses in young, middle-age, and old triple transgenic mouse model of AD (3xT-AD) to the changes that occur in age-matched controls (C57BL/6j). We used immunostaining combined with a novel tissue optical clearing approach along with a computational tool for quantitative analysis of vascular network alterations (vessel length and density) in SVP, IVP, and DVP. All three layers had comparable structural features and densities in young 3xTg-AD and control animals. In controls, IVP and DVP densities decreased with aging (-14% to -32% change from young to old, p < 0.05), while no changes were observed in SVP. In contrast, vascular parameters in the transgenic group decreased in all three layers with aging (-12% to -49% change from young to old, p < 0.05). Furthermore, in the old group, SVP and DVP vascular parameters were lower in the transgenics compared to age-matched controls (p < 0.05). Our analysis demonstrates that normal aging and progression of AD lead to various degrees of vascular alterations in the retina. Specifically, compared to normal aging, changes in vascular features of SVP and DVP regions of the retina are accelerated during AD progression. Considering recent advances in the field of depth-resolved imaging of retinal capillary network and microangiography, noninvasive quantitative monitoring of changes in retinal vascular network parameters of SVP and DVP may serve as markers for diagnosis and staging of Alzheimer's disease and discriminating AD-induced vascular attenuation from age-related vasculopathy.

Metabolic Signaling in a Theoretical Model of the Human Retinal Microcirculation

Impaired blood flow and oxygenation contribute to many ocular pathologies, including glaucoma. Here, a mathematical model is presented that combines an image-based heterogeneous representation of retinal arterioles with a compartmental description of capillaries and venules. The arteriolar model of the human retina is extrapolated from a previous mouse model based on confocal microscopy images. Every terminal arteriole is connected in series to compartments for capillaries and venules, yielding a hybrid model for predicting blood flow and oxygenation throughout the retinal microcirculation. A metabolic wall signal is calculated in each vessel according to blood and tissue oxygen levels. As expected, a higher average metabolic signal is generated in pathways with a lower average oxygen level. The model also predicts a wide range of metabolic signals dependent on oxygen levels and specific network location. For example, for high oxygen demand, a threefold range in metabolic signal is predicted despite nearly identical PO2 levels. This whole-network approach, including a spatially nonuniform structure, is needed to describe the metabolic status of the retina. This model provides the geometric and hemodynamic framework necessary to predict ocular blood flow regulation and will ultimately facilitate early detection and treatment of ischemic and metabolic disorders of the eye.

Concurrent OCT and OCT angiography of retinal neurovascular degeneration in the 5XFAD Alzheimer's disease mice

Significance: As one part of the central nervous system, the retina manifests neurovascular defects in Alzheimer's disease (AD). Quantitative imaging of retinal neurovascular abnormalities may promise a new method for early diagnosis and treatment assessment of AD. Previous imaging studies of transgenic AD mouse models have been limited to the central part of the retina. Given that the pathological hallmarks of AD frequently appear in different peripheral quadrants, a comprehensive regional investigation is needed for a better understanding of the retinal degeneration associated with AD-like pathology. Aim: We aim to demonstrate concurrent optical coherence tomography (OCT) and OCT angiography (OCTA) of retinal neuronal and vascular abnormalities in the 5XFAD mouse model and to investigate region-specific retinal degeneration. Approach: A custom-built OCT system was used for retinal imaging. Retinal thickness, vessel width, and vessel density were quantitatively measured. The artery and vein (AV) were classified for differential AV analysis, and trilaminar vascular plexuses were segmented for depth-resolved density measurement. Results: It was observed that inner and outer retinal thicknesses were explicitly reduced in the dorsal and temporal quadrants, respectively, in 5XFAD mice. A significant arterial narrowing in 5XFAD mice was also observed. Moreover, overall capillary density consistently showed a decreasing trend in 5XFAD mice, but regional specificity was not identified. Conclusions: Quadrant- and layer-specific neurovascular degeneration was observed in 5XFAD mice. Concurrent OCT and OCTA promise a noninvasive method for quantitative monitoring of AD progression and treatment assessment.

Microvascular Tissue Engineering-A Review

Tissue engineering and regenerative medicine have come a long way in recent decades, but the lack of functioning vasculature is still a major obstacle preventing the development of thicker, physiologically relevant tissue constructs. A large part of this obstacle lies in the development of the vessels on a microscale-the microvasculature-that are crucial for oxygen and nutrient delivery. In this review, we present the state of the art in the field of microvascular tissue engineering and demonstrate the challenges for future research in various sections of the field. Finally, we illustrate the potential strategies for addressing some of those challenges.

Blood flow regulation and oxygen transport in a heterogeneous model of the mouse retina

Elevated intraocular pressure is the primary risk factor for glaucoma, yet vascular health and ocular hemodynamics have also been established as important risk factors for the disease. The precise physiological mechanisms and processes by which flow impairment and reduced tissue oxygenation relate to retinal ganglion cell death are not fully known. Mathematical modeling has emerged as a useful tool to help decipher the role of hemodynamic alterations in glaucoma. Several previous models of the retinal microvasculature and tissue have investigated the individual impact of spatial heterogeneity, flow regulation, and oxygen transport on the system. This study combines all three of these components into a heterogeneous mathematical model of retinal arterioles that includes oxygen transport and acute flow regulation in response to changes in pressure, shear stress, and oxygen demand. The metabolic signal (S_i) is implemented as a wall-derived signal that reflects the oxygen deficit along the network, and three cases of conduction are considered: no conduction, a constant signal, and a flow-weighted signal. The model shows that the heterogeneity of the downstream signal serves to regulate flow better than a constant conducted response. In fact, the increases in average tissue PO₂ due to a flow-weighted signal are often more significant than if the entire level of signal is increased. Such theoretical work supports the importance of the non-uniform structure of the retinal vasculature when assessing the capability and/or dysfunction of blood flow regulation in the retinal microcirculation.

Regulation of blood flow in the retinal trilaminar vascular network

Light stimulation evokes neuronal activity in the retina, resulting in the dilation of retinal blood vessels and increased blood flow. This response, named functional hyperemia, brings oxygen and nutrients to active neurons. However, it remains unclear which vessels mediate functional hyperemia. We have characterized blood flow regulation in the rat retina in vivo by measuring changes in retinal vessel diameter and red blood cell (RBC) flux evoked by a flickering light stimulus. We found that, in first- and second-order arterioles, flicker evoked large (7.5 and 5.0%), rapid (0.73 and 0.70 s), and consistent dilations. Flicker-evoked dilations in capillaries were smaller (2.0%) and tended to have a slower onset (0.97 s), whereas dilations in venules were smaller (1.0%) and slower (1.06 s) still. The proximity of pericyte somata did not predict capillary dilation amplitude. Expression of the contractile protein α-smooth muscle actin was high in arterioles and low in capillaries. Unexpectedly, we found that blood flow in the three vascular layers was differentially regulated. Flicker stimulation evoked far larger dilations and RBC flux increases in the intermediate layer capillaries than in the superficial and deep layer capillaries (2.6 vs 0.9 and 0.7% dilation; 25.7 vs 0.8 and 11.3% RBC flux increase). These results indicate that functional hyperemia in the retina is driven primarily by active dilation of arterioles. The dilation of intermediate layer capillaries is likely mediated by active mechanisms as well. The physiological consequences of differential regulation in the three vascular layers are discussed.

Effect of intraocular pressure on the hemodynamics of the central retinal artery: a mathematical model

Retinal hemodynamics plays a crucial role in the pathophysiology of several ocular diseases. There are clear evidences that the hemodynamics of the central retinal artery (CRA) is strongly affected by the level of intraocular pressure (IOP), which is the pressure inside the eye globe. However, the mechanisms through which this occurs are still elusive. The main goal of this paper is to develop a mathematical model that combines the mechanical action of IOP and the blood flow in the CRA to elucidate the mechanisms through which IOP elevation affects the CRA hemodynamics. Our model suggests that the development of radial compressive regions in the lamina cribrosa (a collagen structure in the optic nerve pierced by the CRA approximately in its center) might be responsible for the clinically-observed blood velocity reduction in the CRA following IOP elevation. The predictions of the mathematical model are in very good agreement with experimental and clinical data. Our model also identifies radius and thickness of the lamina cribrosa as major factors affecting the IOP-CRA relationship, suggesting that anatomical differences among individuals might lead to different hemodynamic responses to IOP elevation.

Intraocular pressure, blood pressure, and retinal blood flow autoregulation: a mathematical model to clarify their relationship and clinical relevance

PURPOSE

This study investigates the relationship between intraocular pressure (IOP) and retinal hemodynamics and predicts how arterial blood pressure (BP) and blood flow autoregulation (AR) influence this relationship.

METHODS

A mathematical model is developed to simulate blood flow in the central retinal vessels and retinal microvasculature as current flowing through a network of resistances and capacitances. Variable resistances describe active and passive diameter changes due to AR and IOP. The model is validated by using clinically measured values of retinal blood flow and velocity. The model simulations for six theoretical patients with high, normal, and low BP (HBP-, NBP-, LBP-) and functional or absent AR (-wAR, -woAR) are compared with clinical data.

RESULTS

The model predicts that NBPwAR and HBPwAR patients can regulate retinal blood flow (RBF) as IOP varies between 15 and 23 mm Hg and between 23 and 29 mm Hg, respectively, whereas LBPwAR patients do not adequately regulate blood flow if IOP is 15 mm Hg or higher. Hemodynamic alterations would be noticeable only if IOP changes occur outside of the regulating range, which, most importantly, depend on BP. The model predictions are consistent with clinical data for IOP reduction via surgery and medications and for cases of induced IOP elevation.

CONCLUSIONS

The theoretical model results suggest that the ability of IOP to induce noticeable changes in retinal hemodynamics depends on the levels of BP and AR of the individual. These predictions might help to explain the inconsistencies found in the clinical literature concerning the relationship between IOP and retinal hemodynamics.

An imaging-based stochastic model for simulation of tumour vasculature

A mathematical model which reconstructs the structure of existing vasculature using patient-specific anatomical, functional and molecular imaging as input was developed. The vessel structure is modelled according to empirical vascular parameters, such as the mean vessel branching angle. The model is calibrated such that the resultant oxygen map modelled from the simulated microvasculature stochastically matches the input oxygen map to a high degree of accuracy (R(2) ≈ 1). The calibrated model was successfully applied to preclinical imaging data. Starting from the anatomical vasculature image (obtained from contrast-enhanced computed tomography), a representative map of the complete vasculature was stochastically simulated as determined by the oxygen map (obtained from hypoxia [(64)Cu]Cu-ATSM positron emission tomography). The simulated microscopic vasculature and the calculated oxygenation map successfully represent the imaged hypoxia distribution (R(2) = 0.94). The model elicits the parameters required to simulate vasculature consistent with imaging and provides a key mathematical relationship relating the vessel volume to the tissue oxygen tension. Apart from providing an excellent framework for visualizing the imaging gap between the microscopic and macroscopic imagings, the model has the potential to be extended as a tool to study the dynamics between the tumour and the vasculature in a patient-specific manner and has an application in the simulation of anti-angiogenic therapies.