A Rapid Capillary-Pressure Driven Micro-Channel to Demonstrate Newtonian Fluid Behavior of Zebrafish Blood at High Shear Rates

Role of tissue biomechanics in the formation and function of myocardial trabeculae in zebrafish embryos

Cardiac trabeculae are uneven ventricular muscular structures that develop during early embryonic heart development at the outer curvature of the ventricle. Their biomechanical function is not completely understood, and while their formation is known to be mechanosensitive, it is unclear whether ventricular tissue internal stresses play an important role in their formation. Here, we performed imaging and image-based cardiac biomechanics simulations on zebrafish embryonic ventricles to investigate these issues. Microscopy-based ventricular strain measurements show that the appearance of trabeculae coincided with enhanced deformability of the ventricular wall. Image-based biomechanical simulations reveal that the presence of trabeculae reduces ventricular tissue internal stresses, likely acting as structural support in response to the geometry of the ventricle. Passive ventricular pressure-loading experiments further reveal that the formation of trabeculae is associated with a spatial homogenization of ventricular tissue stiffnesses in healthy hearts, but gata1 morphants with a disrupted trabeculation process retain a spatial stiffness heterogeneity. Our findings thus suggest that modulating ventricular wall deformability, stresses, and stiffness are among the biomechanical functions of trabeculae. Further, experiments with gata1 morphants reveal that a reduction in fluid pressures and consequently ventricular tissue internal stresses can disrupt trabeculation, but a subsequent restoration of ventricular tissue internal stresses via vasopressin rescues trabeculation, demonstrating that tissue stresses are important to trabeculae formation. Overall, we find that tissue biomechanics is important to the formation and function of embryonic heart trabeculation. KEY POINTS: Trabeculations are fascinating and important cardiac structures and their abnormalities are linked to embryonic demise. However, their function in the heart and their mechanobiological formation processes are not completely understood. Our imaging and modelling show that tissue biomechanics is the key here. We find that trabeculations enhance cardiac wall deformability, reduce fluid pressure stresses, homogenize wall stiffness, and have alignments that are optimal for providing load-bearing structural support for the heart. We further discover that high ventricular tissue internal stresses consequent to high fluid pressures are needed for trabeculation formation through a rescue experiment, demonstrating that myocardial tissue stresses are as important as fluid flow wall shear stresses for trabeculation formation.

Microfluidic techniques for mechanical measurements of biological samples

The use of microfluidics to make mechanical property measurements is increasingly common. Fabrication of microfluidic devices has enabled various types of flow control and sensor integration at micrometer length scales to interrogate biological materials. For rheological measurements of biofluids, the small length scales are well suited to reach high rates, and measurements can be made on droplet-sized samples. The control of flow fields, constrictions, and external fields can be used in microfluidics to make mechanical measurements of individual bioparticle properties, often at high sampling rates for high-throughput measurements. Microfluidics also enables the measurement of bio-surfaces, such as the elasticity and permeability properties of layers of cells cultured in microfluidic devices. Recent progress on these topics is reviewed, and future directions are discussed.

Biomechanical Assessment of Red Blood Cells in Pulsatile Blood Flows

As rheological properties are substantially influenced by red blood cells (RBCs) and plasma, the separation of their individual contributions in blood is essential. The estimation of multiple rheological factors is a critical issue for effective early detection of diseases. In this study, three rheological properties (i.e., viscoelasticity, RBC aggregation, and blood junction pressure) are measured by analyzing the blood velocity and image intensity in a microfluidic device. Using a single syringe pump, the blood flow rate sets to a pulsatile flow pattern (Q_b[t] = 1 + 0.5 sin(2πt/240) mL/h). Based on the discrete fluidic circuit model, the analytical formula of the time constant (λ_b) as viscoelasticity is derived and obtained at specific time intervals by analyzing the pulsatile blood velocity. To obtain RBC aggregation by reducing blood velocity substantially, an air compliance unit (ACU) is used to connect polyethylene tubing (i.d. = 250 µm, length = 150 mm) to the blood channel in parallel. The RBC aggregation index (AI) is obtained by analyzing the microscopic image intensity. The blood junction pressure (β) is obtained by integrating the blood velocity within the ACU. As a demonstration, the present method is then applied to detect either RBC-aggregated blood with different concentrations of dextran solution or hardened blood with thermally shocked RBCs. Thus, it can be concluded that the present method has the ability to consistently detect differences in diluent or RBCs in terms of three rheological properties.

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.

Fluid mechanics of the zebrafish embryonic heart trabeculation

Embryonic heart development is a mechanosensitive process, where specific fluid forces are needed for the correct development, and abnormal mechanical stimuli can lead to malformations. It is thus important to understand the nature of embryonic heart fluid forces. However, the fluid dynamical behaviour close to the embryonic endocardial surface is very sensitive to the geometry and motion dynamics of fine-scale cardiac trabecular surface structures. Here, we conducted image-based computational fluid dynamics (CFD) simulations to quantify the fluid mechanics associated with the zebrafish embryonic heart trabeculae. To capture trabecular geometric and motion details, we used a fish line that expresses fluorescence at the endocardial cell membrane, and high resolution 3D confocal microscopy. Our endocardial wall shear stress (WSS) results were found to exceed those reported in existing literature, which were estimated using myocardial rather than endocardial boundaries. By conducting simulations of single intra-trabecular spaces under varied scenarios, where the translational or deformational motions (caused by contraction) were removed, we found that a squeeze flow effect was responsible for most of the WSS magnitude in the intra-trabecular spaces, rather than the shear interaction with the flow in the main ventricular chamber. We found that trabecular structures were responsible for the high spatial variability of the magnitude and oscillatory nature of WSS, and for reducing the endocardial deformational burden. We further found cells attached to the endocardium within the intra-trabecular spaces, which were likely embryonic hemogenic cells, whose presence increased endocardial WSS. Overall, our results suggested that a complex multi-component consideration of both anatomic features and motion dynamics were needed to quantify the trabeculated embryonic heart fluid mechanics.

In the embryonic heart, the mechanical forces that blood fluid imposes on the cardiac tissues are known to be important biological stimuli that affect the proper heart development. We thus perform careful quantification of these forces, using the zebrafish embryo as a model. To do this, we perform high resolution imaging of zebrafish embryonic hearts and image-based flow simulations. We find that the use of a particular fish line that expresses fluorescence at the exact boundary between heart tissue and blood, that is the endocardial cell membrane boundary, is important to give high quality results. The heart’s inner surface has uneven trabeculation structures. We find that they cause fluid forces to have spatial variability and an oscillatory nature. We also find that there is a squeezing motion of cardiac tissues on the trabeculation fluid spaces, which is the main mechanism that generated fluid forces. Fluid forces are also affected by a number of cardiac cells that were developing into blood cells, lodged in the trabeculation fluid spaces. Our investigations provide an understanding of the complexity of the fluid forces on the inner surface of the embryonic heart, and our quantifications will be useful to future studies on the biology elicited by these fluid forces.

In vivo acoustic manipulation of microparticles in zebrafish embryos

In vivo micromanipulation using ultrasound is an exciting technology with promises for cancer research, brain research, vasculature biology, diseases, and treatment development. In the present work, we demonstrate in vivo manipulation of gas-filled microparticles using zebrafish embryos as a vertebrate model system. Micromanipulation methods often are conducted in vitro, and they do not fully reflect the complex environment associated in vivo. Four piezoelectric actuators were positioned orthogonally to each other around an off-centered fluidic channel that allowed for two-dimensional manipulation of intravenously injected microbubbles. Selective manipulation of microbubbles inside a blood vessel with micrometer precision was achieved without interfering with circulating blood cells. Last, we studied the viability of zebrafish embryos subjected to the acoustic field. This successful high-precision, in vivo acoustic manipulation of intravenously injected microbubbles offers potentially promising therapeutic options.

Sequential quantification of blood and diluent using red cell sedimentation-based separation and pressure-induced work in a microfluidic channel

The erythrocyte sedimentation method has been widely used to detect inflammatory diseases. However, this conventional method still has several drawbacks, such as a large blood volume (∼1 mL) and difficulty in continuous monitoring. Most importantly, image-based methods cannot quantify RBC-rich blood (blood) and RBC-free blood (diluent) simultaneously. In this study, instead of visualizing interface movement in the blood syringe, a simple method is proposed to quantify blood and diluent in microfluidic channels sequentially. The hematocrit was set to 25% to enhance RBC sedimentation and form two layers (blood and diluent) in the blood syringe. An air cavity (∼300 μL) inside the blood syringe was secured to completely remove dead volumes (∼200 μL) in fluidic paths (syringe needle and tubing). Thus, a small blood volume (V_b = 50 μL) suctioned into the blood syringe is sufficient for supplying blood and diluent in the blood channel sequentially. The relative ratio of blood resident time (RBC-to-diluent separation) was quantified using λ_b, which was obtained by quantifying the image intensity of blood flow. After the junction pressure (P_j) and blood volume (V) were obtained by analyzing the interface in the coflowing channel, the averaged work (W_p [Pa mm³]) was calculated and adopted to detect blood and diluent, respectively. The proposed method was then applied with various concentrations of dextran solution to detect aggregation-elevated blood. The W_p of blood and diluent exhibited substantial differences with respect to dextran solutions ranging from C_dex = 10 to C_dex = 40 mg mL^-1. Moreover, λ_b did not exhibit substantial differences in blood with C_dex > 10 mg mL^-1. The variations in λ_b were comparable to those of the previous method based on interface movement in the blood syringe. In conclusion, the W_P could detect blood as well as diluents more effectively than λ_b. Furthermore, the proposed method substantially reduced the blood volume from 1 mL to 50 μL.

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

Wall shear stress (WSS) contributes to the mechanotransduction underlying microvascular development and regeneration. Using computational fluid dynamics, we elucidated the interplay between WSS and vascular remodelling in a zebrafish model of tail amputation and regeneration. The transgenic Tg (fli1:eGFP; Gata1:ds-red) zebrafish line was used to track the three-dimensional fluorescently labelled vascular endothelium for post-image segmentation and reconstruction of the fluid domain. Particle image velocimetry was used to validate the blood flow. Following amputation to the dorsal aorta and posterior cardinal vein (PCV), vasoconstriction developed in the dorsal longitudinal anastomotic vessel (DLAV) along with increased WSS in the proximal segmental vessels (SVs) from amputation. Angiogenesis ensued at the tips of the amputated DLAV and PCV where WSS was minimal. At 2 days post amputation (dpa), vasodilation occurred in a pair of SVs proximal to amputation, followed by increased blood flow and WSS; however, in the SVs distal to amputation, WSS normalized to the baseline. At 3 dpa, the blood flow increased in the arterial SV proximal to amputation and through anastomosis with DLAV formed a loop with PCV. Thus, our in silico modelling revealed the interplay between WSS and microvascular adaptation to changes in WSS and blood flow to restore microcirculation following tail amputation.

Cardiac forces regulate zebrafish heart valve delamination by modulating Nfat signaling

In the clinic, most cases of congenital heart valve defects are thought to arise through errors that occur after the endothelial–mesenchymal transition (EndoMT) stage of valve development. Although mechanical forces caused by heartbeat are essential modulators of cardiovascular development, their role in these later developmental events is poorly understood. To address this question, we used the zebrafish superior atrioventricular valve (AV) as a model. We found that cellularized cushions of the superior atrioventricular canal (AVC) morph into valve leaflets via mesenchymal–endothelial transition (MEndoT) and tissue sheet delamination. Defects in delamination result in thickened, hyperplastic valves, and reduced heart function. Mechanical, chemical, and genetic perturbation of cardiac forces showed that mechanical stimuli are important regulators of valve delamination. Mechanistically, we show that forces modulate Nfatc activity to control delamination. Together, our results establish the cellular and molecular signature of cardiac valve delamination in vivo and demonstrate the continuous regulatory role of mechanical forces and blood flow during valve formation.

Why do developing zebrafish atrioventricular heart valves become hyperplastic under certain hemodynamic conditions? This study suggests that part of the answer lies in how the mechanosensitive Nfat pathway regulates the valve mesenchymal-to-endothelial transition.

Wasp controls oriented migration of endothelial cells to achieve functional vascular patterning

Endothelial cell migration and proliferation are essential for the establishment of a hierarchical organization of blood vessels and optimal distribution of blood. However, how these cellular processes are quantitatively coordinated to drive vascular network morphogenesis remains unknown. Here, using the zebrafish vasculature as a model system, we demonstrate that the balanced distribution of endothelial cells, as well as the resulting regularity of vessel calibre, is a result of cell migration from veins towards arteries and cell proliferation in veins. We identify the Wiskott-Aldrich Syndrome protein (WASp) as an important molecular regulator of this process and show that loss of coordinated migration from veins to arteries upon wasb depletion results in aberrant vessel morphology and the formation of persistent arteriovenous shunts. We demonstrate that WASp achieves its function through the coordination of junctional actin assembly and PECAM1 recruitment and provide evidence that this is conserved in humans. Overall, we demonstrate that functional vascular patterning in the zebrafish trunk is established through differential cell migration regulated by junctional actin, and that interruption of differential migration may represent a pathomechanism in vascular malformations.

Summary: Regular diameter of developing veins and arteries in the zebrafish trunk is controlled by differential endothelial cell proliferation and WASp-driven directed cell migration.

Haemodynamic dependence of mechano-genetic evolution of the cardiovascular system in Japanese medaka

The progression of cardiac gene expression-wall shear stress (WSS) interplay is critical to identifying developmental defects during cardiovascular morphogenesis. However, mechano-genetics from the embryonic to larval stages are poorly understood in vertebrates. We quantified peak WSS in the heart and tail vessels of Japanese medaka from 3 days post fertilization (dpf) to 14 dpf using in vivo micro-particle image velocimetry flow measurements, and in parallel analysed the expression of five cardiac genes (fgf8, hoxb6b, bmp4, nkx2.5, smyd1). Here, we report that WSS in the atrioventricular canal (AVC), ventricular outflow tract (OFT), and the caudal vessels in medaka peak with inflection points at 6 dpf and 10-11 dpf instead of a monotonic trend. Retrograde flows are captured at the AVC and OFT of the medaka heart for the first time. In addition, all genes were upregulated at 3 dpf and 7 dpf, indicating a possible correlation between the two, with the cardiac gene upregulation preceding WSS increase in order to facilitate cardiac wall remodelling.

Microfluidic Viscometer Using a Suspending Micromembrane for Measurement of Biosamples

The viscosity of biofluids such as blood and saliva can reflect an individual’s health conditions, and viscosity measurements are therefore considered in health monitoring and disease diagnosis. However, conventional viscometers can only handle a larger liquid volume beyond the quantity that can be extracted from a person. Though very effective, micro-sensors based on electrokinetic, ultrasonic, or other principles often have strict requirements for the supporting equipment and complicated procedures and signal processing. Sample contamination is always an important issue. In this paper, we report a microfluidic viscometer requiring a small volume of biosamples (<50 µL) and straightforward operation procedures. It is fabricated with low-cost and biocompatible polymeric materials as one-time-use devices, such that contamination is no longer the concern. It contains a suspending micromembrane located along a microchannel. Under a steady driving pressure, the membrane displacement is a function of viscosity of the liquid sample being tested. We derived a simple analytical relation and perform a simulation for converting the membrane displacement to the sample viscosity. We conducted experiments with liquids (water and mineral oil) with defined properties to verify such a relation. We further applied the micro-viscometer to measure bovine blood samples with different hematocrit levels. It can be concluded that the microfluidic viscometer has a high compatibility with a broad range of biomedical applications.

Capillary-Driven Flow Microfluidics Combined with Smartphone Detection: An Emerging Tool for Point-of-Care Diagnostics

Point-of-care (POC) or near-patient testing allows clinicians to accurately achieve real-time diagnostic results performed at or near to the patient site. The outlook of POC devices is to provide quicker analyses that can lead to well-informed clinical decisions and hence improve the health of patients at the point-of-need. Microfluidics plays an important role in the development of POC devices. However, requirements of handling expertise, pumping systems and complex fluidic controls make the technology unaffordable to the current healthcare systems in the world. In recent years, capillary-driven flow microfluidics has emerged as an attractive microfluidic-based technology to overcome these limitations by offering robust, cost-effective and simple-to-operate devices. The internal wall of the microchannels can be pre-coated with reagents, and by merely dipping the device into the patient sample, the sample can be loaded into the microchannel driven by capillary forces and can be detected via handheld or smartphone-based detectors. The capabilities of capillary-driven flow devices have not been fully exploited in developing POC diagnostics, especially for antimicrobial resistance studies in clinical settings. The purpose of this review is to open up this field of microfluidics to the ever-expanding microfluidic-based scientific community.

Design and Fabrication of Capillary-Driven Flow Device for Point-Of-Care Diagnostics

Point-of-care (POC) diagnostics enables the diagnosis and monitoring of patients from the clinic or their home. Ideally, POC devices should be compact, portable and operatable without the requirement of expertise or complex fluid mechanical controls. This paper showcases a chip-and-dip device, which works on the principle of capillary-driven flow microfluidics and allows analytes’ detection by multiple microchannels in a single microchip via smartphone imaging. The chip-and-dip device, fabricated with inexpensive materials, works by simply dipping the reagents-coated microchip consisting of microchannels into a fluidic sample. The sample is loaded into the microchannels via capillary action and reacts with the reagents to produce a colourimetric signal. Unlike dipstick tests, this device allows the loading of bacterial/pathogenic samples for antimicrobial testing. A single device can be coated with multiple reagents, and more analytes can be detected in one sample. This platform could be used for a wide variety of assays. Here, we show the design, fabrication and working principle of the chip-and-dip flow device along with a specific application consisting in the determination of β-lactamase activity and cortisol. The simplicity, robustness and multiplexing capability of the chip-and-dip device will allow it to be used for POC diagnostics.

4D modelling of fluid mechanics in the zebrafish embryonic heart

Abnormal blood flow mechanics can result in pathological heart malformation, underlining the importance of understanding embryonic cardiac fluid mechanics. In the current study, we performed image-based computational fluid dynamics simulation of the zebrafish embryonic heart ventricles and characterized flow mechanics, organ dynamics, and energy dynamics in detail. 4D scans of 5 days post-fertilization embryonic hearts with GFP-labelled myocardium were acquired using line-scan focal modulation microscopy. This revealed that the zebrafish hearts exhibited a wave-like contractile/relaxation motion from the inlet to the outlet during both systole and diastole, which we showed to be an energy efficient configuration. No impedance pumping effects of pressure and velocity waves were observed. Due to its tube-like configuration, inflow velocities were higher near the inlet and smaller at the outlet and vice versa for outflow velocities. This resulted in an interesting spatial wall shear stress (WSS) pattern where WSS waveforms near the inlet and those near the outlet were out of phase. There was large spatial variability in WSS magnitudes. Peak WSS was in the range of 47.5-130 dyne/cm² at the inflow and outflow tracts, but were much smaller, in the range of 4-11 dyne/cm², in the mid-ventricular segment. Due to very low Reynolds number and the highly viscous environment, intraventricular pressure gradients were high, suggesting substantial energy losses of flow through the heart.

Blood flow-induced Notch activation and endothelial migration enable vascular remodeling in zebrafish embryos

Arteries and veins are formed independently by different types of endothelial cells (ECs). In vascular remodeling, arteries and veins become connected and some arteries become veins. It is unclear how ECs in transforming vessels change their type and how fates of individual vessels are determined. In embryonic zebrafish trunk, vascular remodeling transforms arterial intersegmental vessels (ISVs) into a functional network of arteries and veins. Here we find that, once an ISV is connected to venous circulation, venous blood flow promotes upstream migration of ECs that results in displacement of arterial ECs by venous ECs, completing the transformation of this ISV into a vein without trans-differentiation of ECs. Arterial blood flow initiated in two neighboring ISVs prevents their transformation into veins by activating Notch signaling in ECs. Together, different responses of ECs to arterial and venous blood flow lead to formation of a balanced network with equal numbers of arteries and veins.

Multiscale light-sheet for rapid imaging of cardiopulmonary system

The ability to image tissue morphogenesis in real-time and in 3-dimensions (3-D) remains an optical challenge. The advent of light-sheet fluorescence microscopy (LSFM) has advanced developmental biology and tissue regeneration research. In this review, we introduce a LSFM system in which the illumination lens reshapes a thin light-sheet to rapidly scan across a sample of interest while the detection lens orthogonally collects the imaging data. This multiscale strategy provides deep-tissue penetration, high-spatiotemporal resolution, and minimal photobleaching and phototoxicity, allowing in vivo visualization of a variety of tissues and processes, ranging from developing hearts in live zebrafish embryos to ex vivo interrogation of the microarchitecture of optically cleared neonatal hearts. Here, we highlight multiple applications of LSFM and discuss several studies that have allowed better characterization of developmental and pathological processes in multiple models and tissues. These findings demonstrate the capacity of multiscale light-sheet imaging to uncover cardiovascular developmental and regenerative phenomena.

Spatial and temporal variations in hemodynamic forces initiate cardiac trabeculation

Hemodynamic shear force has been implicated as modulating Notch signaling-mediated cardiac trabeculation. Whether the spatiotemporal variations in wall shear stress (WSS) coordinate the initiation of trabeculation to influence ventricular contractile function remains unknown. Using light-sheet fluorescent microscopy, we reconstructed the 4D moving domain and applied computational fluid dynamics to quantify 4D WSS along the trabecular ridges and in the groves. In WT zebrafish, pulsatile shear stress developed along the trabecular ridges, with prominent endocardial Notch activity at 3 days after fertilization (dpf), and oscillatory shear stress developed in the trabecular grooves, with epicardial Notch activity at 4 dpf. Genetic manipulations were performed to reduce hematopoiesis and inhibit atrial contraction to lower WSS in synchrony with attenuation of oscillatory shear index (OSI) during ventricular development. γ-Secretase inhibitor of Notch intracellular domain (NICD) abrogated endocardial and epicardial Notch activity. Rescue with NICD mRNA restored Notch activity sequentially from the endocardium to trabecular grooves, which was corroborated by observed Notch-mediated cardiomyocyte proliferations on WT zebrafish trabeculae. We also demonstrated in vitro that a high OSI value correlated with upregulated endothelial Notch-related mRNA expression. In silico computation of energy dissipation further supports the role of trabeculation to preserve ventricular structure and contractile function. Thus, spatiotemporal variations in WSS coordinate trabecular organization for ventricular contractile function.

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.