Computational simulations of the 4D micro-circulatory network in zebrafish tail amputation and regeneration

Hemodynamic regulation allows stable growth of microvascular networks

How do vessels find optimal radii? Capillaries are known to adapt their radii to maintain the shear stress of blood flow at the vessel wall at a set point, yet models of adaptation purely based on average shear stress have not been able to produce complex loopy networks that resemble real microvascular systems. For narrow vessels where red blood cells travel in a single file, the shear stress on vessel endothelium peaks sharply when a red blood cell passes through. We show that stable shear-stress-based adaptation is possible if vessel shear stress set points are cued to the stress peaks. Model networks that respond to peak stresses alone can quantitatively reproduce the observed zebrafish trunk microcirculation, including its adaptive trajectory when hematocrit changes or parts of the network are amputated. Our work reveals the potential for mechanotransduction alone to generate stable hydraulically tuned microvascular networks.

A cell-and-plasma numerical model reveals hemodynamic stress and flow adaptation in zebrafish microvessels after morphological alteration

The development of a functional cardiovascular system ensures a sustainable oxygen, nutrient and hormone delivery system for successful embryonic development and homeostasis in adulthood. While early vessels are formed by biochemical signaling and genetic programming, the onset of blood flow provides mechanical cues that participate in vascular remodeling of the embryonic vascular system. The zebrafish is a prolific animal model for studying the quantitative relationship between blood flow and vascular morphogenesis due to a combination of favorable factors including blood flow visualization in optically transparent larvae. In this study, we have developed a cell-and-plasma blood transport model using computational fluid dynamics (CFD) to understand how red blood cell (RBC) partitioning affect lumen wall shear stress (WSS) and blood pressure in zebrafish trunk blood vascular networks with altered rheology and morphology. By performing live imaging of embryos with reduced hematocrit, we discovered that cardiac output and caudal artery flow rates were maintained. These adaptation trends were recapitulated in our CFD models, which showed reduction in network WSS via viscosity reduction in the caudal artery/vein and via pressure gradient weakening in the intersegmental vessels (ISVs). Embryos with experimentally reduced lumen diameter showed reduced cardiac output and caudal artery flow rate. Factoring in this trend into our CFD models, simulations highlighted that lumen diameter reduction increased vessel WSS but this increase was mitigated by flow reduction due to the adaptive network pressure gradient weakening. Additionally, hypothetical network CFD models with different vessel lumen diameter distribution characteristics indicated the significance of axial variation in lumen diameter and cross-sectional shape for establishing physiological WSS gradients along ISVs. In summary, our work demonstrates how both experiment-driven and hypothetical CFD modeling can be employed for the study of blood flow physiology during vascular remodeling.

Integrating 4-D light-sheet fluorescence microscopy and genetic zebrafish system to investigate ambient pollutants-mediated toxicity

Ambient air pollutants, including PM2.5 (aerodynamic diameter d ~2.5 μm), PM10 (d ~10 μm), and ultrafine particles (UFP: d < 0.1 μm) impart both short- and long-term toxicity to various organs, including cardiopulmonary, central nervous, and gastrointestinal systems. While rodents have been the principal animal model to elucidate air pollution-mediated organ dysfunction, zebrafish (Danio rerio) is genetically tractable for its short husbandry and life cycle to study ambient pollutants. Its electrocardiogram (ECG) resembles that of humans, and the fluorescent reporter-labeled tissues in the zebrafish system allow for screening a host of ambient pollutants that impair cardiovascular development, organ regeneration, and gut-vascular barriers. In parallel, the high spatiotemporal resolution of light-sheet fluorescence microscopy (LSFM) enables investigators to take advantage of the transparent zebrafish embryos and genetically labeled fluorescent reporters for imaging the dynamic cardiac structure and function at a single-cell resolution. In this context, our review highlights the integrated strengths of the genetic zebrafish system and LSFM for high-resolution and high-throughput investigation of ambient pollutants-mediated cardiac and intestinal toxicity.

Combined computational analysis and cytology show limited depth osteogenic effect on bone defects in negative pressure wound therapy

Background: The treatment of bone defects remains a clinical challenge. The effect of negative pressure wound therapy (NPWT) on osteogenesis in bone defects has been recognized; however, bone marrow fluid dynamics under negative pressure (NP) remain unknown. In this study, we aimed to examine the marrow fluid mechanics within trabeculae by computational fluid dynamics (CFD), and to verify osteogenic gene expression, osteogenic differentiation to investigate the osteogenic depth under NP. Methods: The human femoral head is scanned using micro-CT to segment the volume of interest (VOI) trabeculae. The VOI trabeculae CFD model simulating the bone marrow cavity is developed by combining the Hypermesh and ANSYS software. The effect of trabecular anisotropy is investigated, and bone regeneration effects are simulated under NP scales of -80, -120, -160, and -200 mmHg. The working distance (WD) is proposed to describe the suction depth of the NP. Finally, gene sequence analysis, cytological experiments including bone mesenchymal stem cells (BMSCs) proliferation and osteogenic differentiation are conducted after the BMSCs are cultured under the same NP scale. Results: The pressure, shear stress on trabeculae, and marrow fluid velocity decrease exponentially with an increase in WD. The hydromechanics of fluid at any WD inside the marrow cavity can be theoretically quantified. The NP scale significantly affects the fluid properties, especially those fluid close to the NP source; however, the effect of the NP scale become marginal as WD deepens. Anisotropy of trabecular structure coupled with the anisotropic hydrodynamic behavior of bone marrow; An NP of -120 mmHg demonstrates the majority of bone formation-related genes, as well as the most effective proliferation and osteogenic differentiation of BMSCs compared to the other NP scales. Conclusion: An NP of -120 mmHg may have the optimal activated ability to promote osteogenesis, but the effective WD may be limited to a certain depth. These findings help improve the understanding of fluid mechanisms behind NPWT in treating bone defects.

High-Throughput Imaging of Blood Flow Reveals Developmental Changes in Distribution Patterns of Hemodynamic Quantities in Developing Zebrafish

Mechanical forces from blood flow and pressure (hemodynamic forces) contribute to the formation and shaping of the blood vascular network during embryonic development. Previous studies have demonstrated that hemodynamic forces regulate signaling and gene expression in endothelial cells that line the inner surface of vascular tubes, thereby modifying their cellular state and behavior. Given its important role in vascular development, we still know very little about the quantitative aspects of hemodynamics that endothelial cells experience due to the difficulty in measuring forces in vivo. In this study, we sought to determine the magnitude of wall shear stress (WSS) exerted on ECs by blood flow in different vessel types and how it evolves during development. Utilizing the zebrafish as a vertebrate model system, we have established a semi-automated high-throughput fluorescent imaging system to capture the flow of red blood cells in an entire zebrafish between 2- and 6-day post-fertilization (dpf). This system is capable of imaging up to 50 zebrafish at a time. A semi-automated analysis method was developed to calculate WSS in zebrafish trunk vessels. This was achieved by measuring red blood cell flow using particle tracking velocimetry analysis, generating a custom-made script to measure lumen diameter, and measuring local tube hematocrit levels to calculate the effective blood viscosity at each developmental stage. With this methodology, we were able to determine WSS magnitude in different vessels at different stages of embryonic and larvae growth and identified developmental changes in WSS, with absolute levels of peak WSS in all vessel types falling to levels below 0.3 Pa at 6 dpf. Additionally, we discovered that zebrafish display an anterior-to-posterior trend in WSS at each developmental stage.