3D Reconstruction of Chick Embryo Vascular Geometries Using Non-invasive High-Frequency Ultrasound for Computational Fluid Dynamics Studies

Myocardial Biomechanics and the Consequent Differentially Expressed Genes of the Left Atrial Ligation Chick Embryonic Model of Hypoplastic Left Heart Syndrome

Left atrial ligation (LAL) of the chick embryonic heart is a model of the hypoplastic left heart syndrome (HLHS) where a purely mechanical intervention without genetic or pharmacological manipulation is employed to initiate cardiac malformation. It is thus a key model for understanding the biomechanical origins of HLHS. However, its myocardial mechanics and subsequent gene expressions are not well-understood. We performed finite element (FE) modeling and single-cell RNA sequencing to address this. 4D high-frequency ultrasound imaging of chick embryonic hearts at HH25 (ED 4.5) were obtained for both LAL and control. Motion tracking was performed to quantify strains. Image-based FE modeling was conducted, using the direction of the smallest strain eigenvector as the orientations of contractions, the Guccione active tension model and a Fung-type transversely isotropic passive stiffness model that was determined via micro-pipette aspiration. Single-cell RNA sequencing of left ventricle (LV) heart tissues was performed for normal and LAL embryos at HH30 (ED 6.5) and differentially expressed genes (DEG) were identified.After LAL, LV thickness increased by 33%, strains in the myofiber direction increased by 42%, while stresses in the myofiber direction decreased by 50%. These were likely related to the reduction in ventricular preload and underloading of the LV due to LAL. RNA-seq data revealed potentially related DEG in myocytes, including mechano-sensing genes (Cadherins, NOTCH1, etc.), myosin contractility genes (MLCK, MLCP, etc.), calcium signaling genes (PI3K, PMCA, etc.), and genes related to fibrosis and fibroelastosis (TGF-β, BMP, etc.). We elucidated the changes to the myocardial biomechanics brought by LAL and the corresponding changes to myocyte gene expressions. These data may be useful in identifying the mechanobiological pathways of HLHS.

Ultrasound-Based Three-Dimensional Microangiography for Repeated Noninvasive Imaging of Neovascularization

Neovascularization plays an essential part in the progression of several diseases as well as recovery processes such as angiogenesis-guided neurogenesis after a spinal cord injury. To visualize and quantify neovascularization in an experimental setup, it is desirable to apply a noninvasive, repeatable, and harmless technique. Here, we apply a newly developed ultrasound-based technique to construct angiographies and quantify neovascularization in the regenerating spinal cord of a regeneration competent animal model, the Mexican axolotl. We measured vessel volume fraction prior to the induction of contusion and transection spinal cord injury and repeated this measurement directly after injury and at 3, 6, 12, and 63 days after injury. Although neither of the injury types resulted in statistically significant differences in vessel volume fraction relative to sham-operated animals, there was a statistically significant increase in neovascularization over time in all groups. Additionally, vessel volume fraction at the final time point (63 days after injury) was quantified with micro-CT imaging after vascular perfusion with a contrast agent, confirming no statistically significant difference in neovascularization between injury types. Ex vivo vessel volume fraction measured by micro-CT was significantly different from the in vivo ultrasound-based measurement at the same time point. This is likely a result of incomplete vascular perfusion with the contrast agent before micro-CT imaging, which was supported by subsequent histological evaluation. In summary, the results suggest that the ultrasound-based angiographic procedure, we demonstrate here, is applicable to visualize and quantify neovascularization in a noninvasive and harmless fashion in longitudinal experiments circumventing the limitations of contrast agent-dependent techniques.

Fluid mechanics of the left atrial ligation chick embryonic model of hypoplastic left heart syndrome

Left atrial ligation (LAL) of the chick embryonic heart at HH21 is a model of the hypoplastic left heart syndrome (HLHS) disease, demonstrating morphological and hemodynamic features similar to human HLHS cases. Since it relies on mechanical intervention without genetic or pharmacological manipulations, it is a good model for understanding the biomechanics origins of such HLHS malformations. To date, however, the fluid mechanical environment of this model is poorly understood. In the current study, we performed 4D ultrasound imaging of LAL and normal chick embryonic hearts and 4D cardiac flow simulations to help shed light on the mechanical environment that may lead to the HLHS morphology. Results showed that the HH25 LAL atrial function was compromised, and velocities in the ventricle were reduced. The HH25 LAL ventricles developed a more triangular shape with a sharper apex, and in some cases, the atrioventricular junction shifted medially. These changes led to more sluggish flow near the ventricular free wall and apex, where more fluid particles moved in an oscillatory manner with the motion of the ventricular wall, while slowly being washed out, resulting in lower wall shear stresses and higher oscillatory indices. Consequent to these flow conditions, at HH28, even before septation is complete, the left ventricle was found to be hypoplastic while the right ventricle was found to be larger in compensation. Our results suggest that the low and oscillatory flow near the left side of the heart may play a role in causing the HLHS morphology in the LAL model.

Computational Modeling of Blood Flow Hemodynamics for Biomechanical Investigation of Cardiac Development and Disease

The heart is the first functional organ in a developing embryo. Cardiac development continues throughout developmental stages while the heart goes through a serious of drastic morphological changes. Previous animal experiments as well as clinical observations showed that disturbed hemodynamics interfere with the development of the heart and leads to the formation of a variety of defects in heart valves, heart chambers, and blood vessels, suggesting that hemodynamics is a governing factor for cardiogenesis, and disturbed hemodynamics is an important source of congenital heart defects. Therefore, there is an interest to image and quantify the flowing blood through a developing heart. Flow measurement in embryonic fetal heart can be performed using advanced techniques such as magnetic resonance imaging (MRI) or echocardiography. Computational fluid dynamics (CFD) modeling is another approach especially useful when the other imaging modalities are not available and in-depth flow assessment is needed. The approach is based on numerically solving relevant physical equations to approximate the flow hemodynamics and tissue behavior. This approach is becoming widely adapted to simulate cardiac flows during the embryonic development. While there are few studies for human fetal cardiac flows, many groups used zebrafish and chicken embryos as useful models for elucidating normal and diseased cardiogenesis. In this paper, we explain the major steps to generate CFD models for simulating cardiac hemodynamics in vivo and summarize the latest findings on chicken and zebrafish embryos as well as human fetal hearts.

High resolution, label-free photoacoustic imaging of live chicken embryo developing in bioengineered eggshell

Chicken embryos have been proven to be an attractive vertebrate model for biomedical research. They have helped in making significant contributions for advancements in various fields like developmental biology, cancer research and cardiovascular studies. However, a non-invasive, label-free method of imaging live chicken embryo at high resolution still needs to be developed and optimized. In this work, we have shown the potential of photoacoustic tomography (PAT) for imaging live chicken embryos cultured in bioengineered eggshells. Laser pulses at wavelengths of 532 and 740 nm were used for attaining cross-sectional images of chicken embryos at different developmental stages. Cross-sections along different depths were imaged to gain knowledge of the relative depth of different vessels and organs. Due to high optical absorption of vasculature and embryonic eye, images with good optical contrast could be acquired using this method. We have thus reported a label-free method of performing cross-sectional imaging of chicken embryos at high resolution demonstrating the capacity of PAT as a promising tool for avian embryo imaging.

Retinal oxygen supply shaped the functional evolution of the vertebrate eye

The retina has a very high energy demand but lacks an internal blood supply in most vertebrates. Here we explore the hypothesis that oxygen diffusion limited the evolution of retinal morphology by reconstructing the evolution of retinal thickness and the various mechanisms for retinal oxygen supply, including capillarization and acid-induced haemoglobin oxygen unloading. We show that a common ancestor of bony fishes likely had a thin retina without additional retinal oxygen supply mechanisms and that three different types of retinal capillaries were gained and lost independently multiple times during the radiation of vertebrates, and that these were invariably associated with parallel changes in retinal thickness. Since retinal thickness confers multiple advantages to vision, we propose that insufficient retinal oxygen supply constrained the functional evolution of the eye in early vertebrates, and that recurrent origins of additional retinal oxygen supply mechanisms facilitated the phenotypic evolution of improved functional eye morphology.

Organ Dynamics and Hemodynamic of the Whole HH25 Avian Embryonic Heart, Revealed by Ultrasound Biomicroscopy, Boundary Tracking, and Flow Simulations

Congenital heart malformations occur to substantial number of pregnancies. Studies showed that abnormal flow biomechanical environments could lead to malformations, making it important to understand the biomechanical environment of the developing heart. We performed 4D high-frequency ultrasound scans of chick embryonic hearts at HH25 to study the biomechanics of the whole heart (atria and ventricle). A novel and high-fidelity motion estimation technique, based on temporal motion model and non-rigid image registration algorithm, allowed automatic tracking of fluid-structure boundaries from scan images, and supported flow simulations. Results demonstrated that atrial appendages were the most contractile portion of the atria, having disproportionately high contribution to atrial blood pumping for its volume in the atria. However, the atria played a small role in blood pumping compared to the ventricle, as it had much lower ejection energy expenditure, and as the ventricle appeared to be able to draw inflow from the veins directly during late diastole. Spatially and temporally averaged wall shear stresses (WSS) for various cardiac structures were 0.062–0.068 Pa, but spatial-averaged WSS could be as high as 0.54 Pa in the RV. WSS was especially elevated at the atrial inlet, atrioventricular junction, regions near to the outflow tract, and at dividing lines between the left and right atrium and left and right side of the ventricle, where septation had begun and the lumen had narrowed. Elevated WSS could serve as biomechanics stimulation for proper growth and development.

Fluid dynamics and forces in the HH25 avian embryonic outflow tract

The embryonic outflow tract (OFT) eventually undergoes aorticopulmonary septation to form the aorta and pulmonary artery, and it is hypothesized that blood flow mechanical forces guide this process. We performed detailed studies of the geometry, wall motions, and fluid dynamics of the HH25 chick embryonic OFT just before septation, using noninvasive 4D high-frequency ultrasound and computational flow simulations. The OFT exhibited expansion and contraction waves propagating from proximal to distal end, with periods of luminal collapse at locations of the two endocardial cushions. This, combined with periods of reversed flow, resulted in the OFT cushions experiencing wall shear stresses (WSS or flow drag forces) with elevated oscillatory characteristics, which could be important to signal for further development of cushions into valves and septum. Furthermore, the OFT exhibits interesting double-helical flow during systole, where a pair of helical flow structures twisted about each other from the proximal to distal end. This coincided with the location of the future aorticopulmonary septum, which also twisted from the proximal to distal end, suggesting that this flow pattern may be guiding OFT septation.

A method to quantify mechanobiologic forces during zebrafish cardiac development using 4-D light sheet imaging and computational modeling

Blood flow and mechanical forces in the ventricle are implicated in cardiac development and trabeculation. However, the mechanisms of mechanotransduction remain elusive. This is due in part to the challenges associated with accurately quantifying mechanical forces in the developing heart. We present a novel computational framework to simulate cardiac hemodynamics in developing zebrafish embryos by coupling 4-D light sheet imaging with a stabilized finite element flow solver, and extract time-dependent mechanical stimuli data. We employ deformable image registration methods to segment the motion of the ventricle from high resolution 4-D light sheet image data. This results in a robust and efficient workflow, as segmentation need only be performed at one cardiac phase, while wall position in the other cardiac phases is found by image registration. Ventricular hemodynamics are then quantified by numerically solving the Navier-Stokes equations in the moving wall domain with our validated flow solver. We demonstrate the applicability of the workflow in wild type zebrafish and three treated fish types that disrupt trabeculation: (a) chemical treatment using AG1478, an ErbB2 signaling inhibitor that inhibits proliferation and differentiation of cardiac trabeculation; (b) injection of gata1a morpholino oligomer (gata1aMO) suppressing hematopoiesis and resulting in attenuated trabeculation; (c) weak-atriumm58 mutant (wea) with inhibited atrial contraction leading to a highly undeveloped ventricle and poor cardiac function. Our simulations reveal elevated wall shear stress (WSS) in wild type and AG1478 compared to gata1aMO and wea. High oscillatory shear index (OSI) in the grooves between trabeculae, compared to lower values on the ridges, in the wild type suggest oscillatory forces as a possible regulatory mechanism of cardiac trabeculation development. The framework has broad applicability for future cardiac developmental studies focused on quantitatively investigating the role of hemodynamic forces and mechanotransduction during morphogenesis.

Organ Dynamics and Fluid Dynamics of the HH25 Chick Embryonic Cardiac Ventricle as Revealed by a Novel 4D High-Frequency Ultrasound Imaging Technique and Computational Flow Simulations

Past literature has provided evidence that a normal mechanical force environment of blood flow may guide normal development while an abnormal environment can lead to congenital malformations, thus warranting further studies on embryonic cardiovascular flow dynamics. In the current study, we developed a non-invasive 4D high-frequency ultrasound technique, and use it to analyze cardiovascular organ dynamics and flow dynamics. Three chick embryos at stage HH25 were scanned with high frequency ultrasound in cine-B-mode at multiple planes spaced at 0.05 mm. 4D images of the heart and nearby arteries were generated via temporal and spatial correlation coupled with quadratic mean ensemble averaging. Dynamic mesh CFD was performed to understand the flow dynamics in the ventricle of the 2 hearts. Our imaging technique has sufficiently high resolution to enable organ dynamics quantification and CFD. Fine structures such as the aortic arches and details such as the cyclic distension of the carotid arteries were captured. The outflow tract completely collapsed during ventricular diastole, possible serving the function of a valve to prevent regurgitation. CFD showed that ventricular wall shear stress (WSS) were in the range of 0.1-0.5 Pa, and that the left side of the common ventricle experienced lower WSS than the right side. The pressure gradient from the inlet to the outlet of the ventricle was positive over most of the cardiac cycle, and minimal regurgitation flow was observed, despite the absence of heart valves. We developed a new image-based CFD method to elucidate cardiac organ dynamics and flow dynamics of embryonic hearts. The embryonic heart appeared to be optimized to generate net forward flow despite the absence of valves, and the WSS environment appeared to be side-specific.

A method to study the hemodynamics of chicken embryo's aortic arches using optical coherence tomography

Congenital cardiovascular defects are the leading cause of birth defect related death. It has been hypothesized that fluid mechanical forces of embryonic blood flow affect cardiovascular development and play a role in congenital malformations. Studies in small animal embryos can improve our understanding of congenital malformations and can lead to better treatment. We present a feasibility study in which high-resolution optical coherence tomography (OCT) and computational fluid dynamics (CFD) are combined to provide quantitative analysis of the embryonic flow mechanics and the associated anatomy in a small animal model.