In Silico Design of Heterogeneous Microvascular Trees Using Generative Adversarial Networks and Constrained Constructive Optimization.

OBJECTIVE: Designing physiologically adequate microvascular trees is of crucial relevance for bioengineering functional tissues and organs. Yet, currently available methods are poorly suited to replicate the morphological and topological heterogeneity of real microvascular trees because the parameters used to control tree generation are too simplistic to mimic results of the complex angiogenetic and structural adaptation processes in vivo. METHODS: We propose a method to overcome this limitation by integrating a conditional deep convolutional generative adversarial network (cDCGAN) with a local fractal dimension-oriented constrained constructive optimization (LFDO-CCO) strategy. The cDCGAN learns the patterns of real microvascular bifurcations allowing for their artificial replication. The LFDO-CCO strategy connects the generated bifurcations hierarchically to form microvascular trees with a vessel density corresponding to that observed in healthy tissues. RESULTS: The generated artificial microvascular trees are consistent with real microvascular trees regarding characteristics such as fractal dimension, vascular density, and coefficient of variation of diameter, length, and tortuosity. CONCLUSIONS: These results support the adoption of the proposed strategy for the generation of artificial microvascular trees in tissue engineering as well as for computational modeling and simulations of microcirculatory physiology.

Coalescent angiogenesis-evidence for a novel concept of vascular network maturation.

Angiogenesis describes the formation of new blood vessels from pre-existing vascular structures. While the most studied mode of angiogenesis is vascular sprouting, specific conditions or organs favor intussusception, i.e., the division or splitting of an existing vessel, as preferential mode of new vessel formation. In the present study, sustained (33-h) intravital microscopy of the vasculature in the chick chorioallantoic membrane (CAM) led to the hypothesis of a novel non-sprouting mode for vessel generation, which we termed "coalescent angiogenesis." In this process, preferential flow pathways evolve from isotropic capillary meshes enclosing tissue islands. These preferential flow pathways progressively enlarge by coalescence of capillaries and elimination of internal tissue pillars, in a process that is the reverse of intussusception. Concomitantly, less perfused segments regress. In this way, an initially mesh-like capillary network is remodeled into a tree structure, while conserving vascular wall components and maintaining blood flow. Coalescent angiogenesis, thus, describes the remodeling of an initial, hemodynamically inefficient mesh structure, into a hierarchical tree structure that provides efficient convective transport, allowing for the rapid expansion of the vasculature with maintained blood supply and function during development.

Visualizing the spatiotemporal pattern of yolk sac membrane vascular network by enhanced local fractal analysis.

OBJECTIVE: To establish methods for providing a comprehensive and detailed description of the spatial distribution of the vascular networks, and to reveal the spatiotemporal pattern of the yolk sac membrane vascular network during the angiogenic procedure. METHODS: Addressing the limitations in the conventional local fractal analysis, an improved approach, named scanning average local fractal dimension, was proposed. This method was conducted on 6 high-resolution vascular images of the yolk sac membrane for 3 eggs at two stages (E3 and E4) to characterize the spatial distribution of the complexity of the vascular network. RESULTS: With the proposed method, the spatial distribution of the complexity of the yolk sac membrane vascular network was visualized. From E3 to E4, the local fractal dimension increased in 3 eggs, 1.80 ± 0.02 vs. 1.85 ± 0.02, 1.72 ± 0.03 vs. 1.83 ± 0.02, and 1.77 ± 0.03 vs. 1.82 ± 0.02, respectively. The mean local fractal dimension in the most distal area from the embryo proper was the lowest at E3 while the highest at E4. At E3, the most peaks of the local fractal dimension were located in the vein territories and shifted to artery territories at E4. CONCLUSIONS: The spatial distribution of the complexity of the yolk sac membrane vascular network exhibited diverse patterns at different stages. In addition from E3 to E4, the increment of complexity at the intersection areas between arteries and sinus terminalis was with the most advance. This is consistent with the physiologic evidence. The present work provides a potential approach for investigating the spatiotemporal pattern of the angiogenic process.

Pulsatility damping in the microcirculation: Basic pattern and modulating factors.

Blood flow pulsatility is an important determinant of macro- and microvascular physiology. Pulsatility is damped largely in the microcirculation, but the characteristics of this damping and the factors that regulate it have not been fully elucidated yet. Applying computational approaches to real microvascular network geometry, we examined the pattern of pulsatility damping and the role of potential damping factors, including pulse frequency, vascular viscous resistance, vascular compliance, viscoelastic behavior of the vessel wall, and wave propagation and reflection. To this end, three full rat mesenteric vascular networks were reconstructed from intravital microscopic recordings, a one-dimensional (1D) model was used to reproduce pulsatile properties within the network, and potential damping factors were examined by sensitivity analysis. Results demonstrate that blood flow pulsatility is predominantly damped at the arteriolar side and remains at a low level at the venular side. Damping was sensitive to pulse frequency, vascular viscous resistance and vascular compliance, whereas viscoelasticity of the vessel wall or wave propagation and reflection contributed little to pulsatility damping. The present results contribute to our understanding of mechanical forces and their regulation in the microcirculation.

Gap junctions regulate vessel diameter in chick chorioallantoic membrane vasculature by both tone-dependent and structural mechanisms.

OBJECTIVE: In this study, we examined the impact of gap junction blockade on chick chorioallantoic membrane microvessels. METHODS: Expression of Cx37, Cx40/42, and Cx43 in chick chorioallantoic membrane tissue was studied by PCR, Western blot, and confocal immunofluorescence microscopy. Vessel diameter changes occurring under gap junction blockade with carbenoxolone (175 µmol/L), palmitoleic acid (100 µmol/L), 43 GAP27 (1 mmol/L) were analyzed by intravital microscopy. To analyze vascular tone, chick chorioallantoic membrane vessels were exposed to a vasodilator cocktail consisting of acetylcholine (10 µmol/L), adenosine (100 µmol/L), papaverine (200 µmol/L), and sodium nitroprusside (10 µmol/L). RESULTS: In chick chorioallantoic membrane lysates, Western blot analysis revealed the expression of Cx40 and Cx43. Immunofluorescence in intact chick chorioallantoic membrane vasculature showed only Cx43, limited to arterial vessel walls. Upon gap junction blockade (3 hours) arterial and venous diameters decreased to 0.50 ± 0.03 and 0.36 ± 0.06 (carbenoxolone), 0.72 ± 0.08 and 0.63 ± 0.15 (palmitoleic acid) and 0.77 ± 0.004 and 0.58 ± 0.05 (GAP27), relative to initial values. Initially, diameter decrease was dominated by increasing vascular tone. After 6 hours, however, vessel tone was reduced, suggesting structural network remodeling. CONCLUSIONS: Our findings suggest a major role for connexins in mediating acute and chronic diameter changes in developing vascular networks.

Investigation into the diversity in the fractal dimensions of arterioles and venules in a microvascular network - A quantitative analysis.

Fractal dimension is a robust fractal parameter for estimating the morphology of vascular networks. It reflects the property of vascular networks that may vary and thus, differentiate between individual networks and/or identify physiological and pathological conditions. As such, fractal dimension differs also between arteriolar and venular compartments, yet the underlying reason is so far unclear. In order to understand the mechanisms behind these differences, we quantitatively analyzed the impacts of vessel attributes on the fractal dimension. Fractal dimension and vessel attributes given by vessel density (VD), vessel length density (VL), and diameter index (DI=VD/VL) were analyzed in three microvascular networks of the rat mesentery, which were reconstructed from experimental data. The results show that differences in diameter between arterioles and venules are primarily responsible for arterio-venous differences in fractal dimension. Moreover, multiple linear regression analysis demonstrates that the sensitivity of the variation of fractal dimension to vessel length and diameter varies with the type of the vessels. While the change of vessel length contributes 57.8 ±â€¯3.4% to the variation of arteriolar dimension, vessel diameter contributes 63.9 ±â€¯4.8% to the variation of venular dimension. The present study provides an explanation for the different fractal dimension and dimension variation in arteriolar and venular compartments. It highlights the importance of estimating the fractal dimensions of arterioles and venules separately, which will enhance the ability of feature extraction by fractal analysis in physiological and clinical application.

Modeling the hematocrit distribution in microcirculatory networks: A quantitative evaluation of a phase separation model.

OBJECTIVE: Theoretical models are essential tools for studying microcirculatory function. Recently, the validity of a well-established phase separation model was questioned and it was claimed that it produces problematically low hematocrit predictions and lack of red cells in small diameter vessels. We conducted a quantitative evaluation of this phase separation model to establish common ground for future research. METHODS: Model predictions were validated against a comprehensive database with measurements from 4 mesenteric networks. A Bayesian data analysis framework was used to integrate measurements and network model simulations into a combined analysis and to model uncertainties related to network boundary conditions as well as phase separation model parameters. The model evaluation was conducted within a cross-validation scheme. RESULTS: Unlike the recently reported results, our analysis demonstrated good correspondence in global characteristics between measurements and predictions. In particular, predicted hematocrits for vessels with small diameters were consistent with measurements. Incorporating phase separation model parameter uncertainties further reduced the hematocrit validation error by 17% and led to the absence of red-cell-free segments. Corresponding model parameters are presented as alternatives to standard parameters. CONCLUSIONS: Consistent with earlier studies, our quantitative model evaluation supports the continued use of the established phase separation model.

Pulse wave velocity in the microcirculation reflects both vascular compliance and resistance: Insights from computational approaches.

OBJECTIVE: PWV is the speed of pulse wave propagation through the circulatory system. mPWV emerges as a novel indicator of hypertension, yet it remains unclear how different vascular properties affect mPWV. We aim to identify the biomechanical determinants of mPWV. METHODS: A 1D model was used to simulate PWV in a rat mesenteric microvascular network and, for comparison, in a human macrovascular arterial network. Sensitivity analysis was performed to assess the relationship between PWV and vascular compliance and resistance. RESULTS: The 1D model enabled adequate simulation of PWV in both micro- and macrovascular networks. Simulated arterial PWV changed as a function of vascular compliance but not resistance, in that arterial PWV varied at a rate of 0.30 m/s and -6.18 × 10-3  m/s per 10% increase in vascular compliance and resistance, respectively. In contrast, mPWV depended on both vascular compliance and resistance, as it varied at a rate of 2.79 and -2.64 cm/s per 10% increase in the respective parameters. CONCLUSIONS: The present study identifies vascular compliance and resistance in microvascular networks as critical determinants of mPWV. We anticipate that mPWV can be utilized as an effective indicator for the assessment of microvascular biomechanical properties.

Coronary microcirculatory pathophysiology: can we afford it to remain a black box?

Coronary microvascular networks play the key role in determining blood flow distribution in the heart. Matching local blood supply to tissue metabolic demand entails continuous adaptation of coronary vessels via regulation of smooth muscle tone and structural dilated vessel diameter. The importance of coronary microcirculation for relevant pathological conditions including angina in patients with normal or near-normal coronary angiograms [microvascular angina (MVA)] and heart failure with preserved ejection fraction (HFpEF) is increasingly recognized. For MVA, clinical studies have shown a prevalence of up to 40% in patients with suspected coronary artery disease and a relevant impact on adverse cardiovascular events including cardiac death, stroke, and heart failure. Despite a continuously increasing number of corresponding clinical studies, the knowledge on pathophysiological cause-effect relations involving coronary microcirculation is, however, still very limited. A number of pathophysiological hypotheses for MVA and HFpEF have been suggested but are not established to a degree, which would allow definition of nosological entities, stratification of affected patients, or development of effective therapeutic strategies. This may be related to a steep decline in experimental (animal) pathophysiological studies in this area during the last 15 years. Since technology to experimentally investigate microvascular pathophysiology in the beating heart is increasingly, in principle, available, a concerted effort to build 'coronary microcirculatory observatories' to close this gap and to accelerate clinical progress in this area is suggested.

Model-based inference from microvascular measurements: Combining experimental measurements and model predictions using a Bayesian probabilistic approach.

OBJECTIVE: In vivo imaging of the microcirculation and network-oriented modeling have emerged as powerful means of studying microvascular function and understanding its physiological significance. Network-oriented modeling may provide the means of summarizing vast amounts of data produced by high-throughput imaging techniques in terms of key, physiological indices. To estimate such indices with sufficient certainty, however, network-oriented analysis must be robust to the inevitable presence of uncertainty due to measurement errors as well as model errors. METHODS: We propose the Bayesian probabilistic data analysis framework as a means of integrating experimental measurements and network model simulations into a combined and statistically coherent analysis. The framework naturally handles noisy measurements and provides posterior distributions of model parameters as well as physiological indices associated with uncertainty. RESULTS: We applied the analysis framework to experimental data from three rat mesentery networks and one mouse brain cortex network. We inferred distributions for more than 500 unknown pressure and hematocrit boundary conditions. Model predictions were consistent with previous analyses, and remained robust when measurements were omitted from model calibration. CONCLUSION: Our Bayesian probabilistic approach may be suitable for optimizing data acquisition and for analyzing and reporting large data sets acquired as part of microvascular imaging studies.

Dynamic remodeling of arteriolar collaterals after acute occlusion in chick chorioallantoic membrane.

OBJECTIVE: After arteriolar occlusion, collaterals enlarge and initially elevated WSS normalizes. While most previous studies focused on endpoints of such adaptive changes in larger collaterals, the present investigation aimed to continuously determine the relation between WSS and diameter in microvascular collaterals during adaptive reactions. METHODS: In Hamburger-Hamilton stage 40 CAMs, junction points between arteriolar segments were identified and the third upstream segment on one side was occluded. Intravital microscopy recordings were taken for 24 hours post-occlusion. Segment diameter and blood velocity were measured: WSS and capillary density were calculated. RESULTS: After occlusion, vascular diameters exhibited an immediate decrease, then increased with a time constant of 2.5 ± 0.8 hours and reached a plateau of up to 60% above baseline after about 7 hours. Vascular tone showed no significant change. WSS exhibited an immediate increase post-occlusion and linearly returned to baseline after about 12 hours. Local WSS change and diameter change rate showed similar patterns during the initial but not the later phase of post-occlusive adaptation. CONCLUSIONS: CAM collaterals undergo fast structural remodeling within 24 hours post-occlusion. This remodeling might be driven by local WSS and by other regulators within the vascular network.

Modeling of pulsatile flow-dependent nitric oxide regulation in a realistic microvascular network.

Hemodynamic pulsatility has been reported to regulate microcirculatory function. To quantitatively assess the impact of flow pulsatility on the microvasculature, a mathematical model was first developed to simulate the regulation of NO production by pulsatile flow in the microcirculation. Shear stress and pressure pulsatility were selected as regulators of endothelial NO production and NO-dependent vessel dilation as feedback to control microvascular hemodynamics. The model was then applied to a real microvascular network of the rat mesentery consisting of 546 microvessels. As compared to steady flow conditions, pulsatile flow increased the average NO concentration in arterioles from 256.8±93.1nM to 274.8±101.1nM (P<0.001), with a corresponding increase in vessel dilation by approximately 7% from 27.5±10.6% to 29.4±11.4% (P<0.001). In contrast, NO concentration and vessel size showed a far lesser increase (about 1.7%) in venules under pulsatile flow as compared to steady flow conditions. Network perfusion and flow heterogeneity were improved under pulsatile flow conditions, and vasodilation within the network was more sensitive to heart rate changes than pulse pressure amplitude. The proposed model simulates the role of flow pulsatility in the regulation of a complex microvascular network in terms of NO concentration and hemodynamics under varied physiological conditions.

Acute Lung Injury Causes Asynchronous Alveolar Ventilation That Can Be Corrected by Individual Sighs.

RATIONALE: Improved ventilation strategies have been the mainstay for reducing mortality in acute respiratory distress syndrome. Their unique clinical effectiveness is, however, unmatched by our understanding of the underlying mechanobiology, and their impact on alveolar dynamics and gas exchange remains largely speculative. OBJECTIVES: To assess changes in alveolar dynamics and associated effects on local gas exchange in experimental models of acute lung injury (ALI) and their responsiveness to sighs. METHODS: Alveolar dynamics and local gas exchange were studied in vivo by darkfield microscopy and multispectral oximetry in experimental murine models of ALI induced by hydrochloric acid, Tween instillation, or in antibody-mediated transfusion-related ALI. MEASUREMENTS AND MAIN RESULTS: Independent of injury mode, ALI resulted in asynchronous alveolar ventilation characteristic of alveolar pendelluft, which either spontaneously resolved or progressed to a complete cessation or even inversion of alveolar ventilation. The functional relevance of the latter phenomena was evident as impaired blood oxygenation in juxtaposed lung capillaries. Individual sighs (2 × 10 s at inspiratory plateau pressure of 30 cm H2O) largely restored normal alveolar dynamics and gas exchange in acid-induced ALI, yet not in Tween-induced surfactant depletion. CONCLUSIONS: We describe for the first time in detail the different forms and temporal sequence of impaired alveolar dynamics in the acutely injured lung and report the first direct visualization of alveolar pendelluft. Moreover, we identify individual sighs as an effective strategy to restore intact alveolar ventilation by a mechanism independent of alveolar collapse and reopening.

Structure and hemodynamics of vascular networks in the chorioallantoic membrane of the chicken.

The chick chorioallantoic membrane (CAM) is extensively used as an in vivo model. Here, structure and hemodynamics of CAM vessel trees were analyzed and compared with predictions of Murray's law. CAM microvascular networks of Hamburger-Hamilton stage 40 chick embryos were scanned by videomicroscopy. Three networks with ∼3,800, 580, and 480 segments were digitally reconstructed, neglecting the capillary mesh. Vessel diameters (D) and segment lengths were measured, and generation numbers and junctional exponents at bifurcations were derived. In selected vessels, flow velocities (v) and hematocrit were measured. Hemodynamic simulations, incorporating the branching of capillaries from preterminal vessels, were used to estimate v, volume flow, shear stress (τ), and pressure for all segments of the largest network. For individual arteriovenous flow pathways, terminal arterial and venous generation numbers are negatively correlated, leading to low variability of total topological and morphological pathway lengths. Arteriolar velocity is proportional to diameter (v∝D1.03 measured, v∝D0.93 modeling), giving nearly uniform τ levels (τ∝D0.05). Venular trees exhibit slightly higher exponents (v∝D1.3, τ∝D0.38). Junctional exponents at divergent and convergent bifurcations were 2.05 ± 1.13 and 1.97 ± 0.95 (mean ± SD) in contrast to the value 3 predicted by Murray's law. In accordance with Murray's law, τ levels are (nearly) maintained in CAM arterial (venular) trees, suggesting vascular adaptation to shear stress. Arterial and venous trees show an interdigitating arrangement providing homogeneous flow pathway properties and have preterminal capillary branches. These properties may facilitate efficient oxygen exchange in the CAM during rapid embryonic growth.

Microvascular Plasticity: Angiogenesis in Health and Disease--Preface.

This Special Topic Issue is concerned with the mechanisms that determine the structure of microvascular networks. The vast number of vessels and the highly plastic character of the microcirculation give evidence that microvascular network structures emerge as a result of responses of individual vessels and cells to the local stimuli that they experience, through a combination of angiogenesis, remodeling and pruning. The articles in this issue of Microcirculation address a range of cellular and molecular mechanisms involved in these processes.

Microvascular hemodynamics in the chick chorioallantoic membrane.

OBJECTIVE: The microvasculature of the CAM in the developing chick embryo is characterized by interdigitating arteriolar and venular trees, connected at multiple points along their lengths to a mesh-like capillary plexus. Theoretical modeling techniques were employed to investigate the resulting hemodynamic characteristics of the CAM. METHODS: Based on previously obtained anatomical data, a model was developed in which the capillary plexus was treated as a porous medium. Supply of blood from arterioles and drainage into venules were represented by distributions of flow sources and sinks. Predicted flow velocities were compared with measurements in arterioles and venules obtained via video microscopy. RESULTS: If it was assumed that blood flowed into and out of the capillary plexus only at the ends of terminal arterioles and venules, the predicted velocities increased with decreasing diameter in vessels below 50 µm in diameter, contrary to the observations. Distributing sources/sinks along arterioles/venules led to velocities consistent with the data. CONCLUSIONS: These results imply that connections to the capillary plexus distributed along the arterioles and venules strongly affect the hemodynamic characteristics of the CAM. The theoretical model provides a basis for quantitative simulations of structural adaptation in CAM networks in response to hemodynamic stimuli.

Deformability of transfused red blood cells is a potent determinant of transfusion-induced change in recipient's blood flow.

OBJECTIVE: There is a growing concern regarding the risks in the transfusion of PRBC, as numerous studies have reported negative transfusion outcomes, including reduced blood perfusion. In search of this phenomenon's mechanism, the effect of PRBC deformability, a major determinant of blood flow, on transfusion outcome was explored. METHODS: The effect of PRBC deformability was examined by the transfusion-induced change in recipients' ∆SBF, in ß-TM patients, who are routinely treated with lifelong frequent transfusions. SBF was determined using a laser Doppler imager. RESULTS: ∆SBF was examined vs PRBC deformability, the transfusion-induced increase in ∆Hct and the recipients' SBF before transfusion (SBFB ). ∆SBF elevated with increasing PRBC deformability, with a highly significant dependence, while its elevation with ∆Hct was much less significant. ∆SBF was inversely proportional to the SBFB . CONCLUSIONS: This study provides, for the first time in humans, direct evidence that the deformability of transfused PRBC is a potent effector of transfusion outcome. Currently, PRBC are supplied primarily by the first-in-first-out criteria, while their functionality is ignored. The testing of PRBC hemodynamic quality would introduce a new paradigm into blood banking, which would contribute substantially to improving transfusion therapy.

Making microvascular networks work: angiogenesis, remodeling, and pruning.

The adequate and efficient functioning of the microcirculation requires not only numerous vessels providing a large surface area for transport but also a structure that provides short diffusion distances from capillaries to tissue and efficient distribution of convective blood flow. Theoretical models show how a combination of angiogenesis, remodeling, and pruning in response to hemodynamic and metabolic stimuli, termed "angioadaptation," generates well organized, functional networks.

Structure-based algorithms for microvessel classification.

OBJECTIVE: Recent developments in high-resolution imaging techniques have enabled digital reconstruction of three-dimensional sections of microvascular networks down to the capillary scale. To better interpret these large data sets, our goal is to distinguish branching trees of arterioles and venules from capillaries. METHODS: Two novel algorithms are presented for classifying vessels in microvascular anatomical data sets without requiring flow information. The algorithms are compared with a classification based on observed flow directions (considered the gold standard), and with an existing resistance-based method that relies only on structural data. RESULTS: The first algorithm, developed for networks with one arteriolar and one venular tree, performs well in identifying arterioles and venules and is robust to parameter changes, but incorrectly labels a significant number of capillaries as arterioles or venules. The second algorithm, developed for networks with multiple inlets and outlets, correctly identifies more arterioles and venules, but is more sensitive to parameter changes. CONCLUSIONS: The algorithms presented here can be used to classify microvessels in large microvascular data sets lacking flow information. This provides a basis for analyzing the distinct geometrical properties and modelling the functional behavior of arterioles, capillaries, and venules.