In vivo wall shear measurements within the developing zebrafish heart

Evolution of cardiac tissue and flow mechanics in developing Japanese Medaka

The effects of pressure drop across cardiac valve cushion regions and endocardial wall strain in the early developmental stages of a teleost species heart are poorly understood. In the presented work, we utilize microscale particle image velocimetry (μPIV) flow measurements of developing medaka hearts from 3 to 14 dpf (n = 5 at each dpf) to quantify the pressure field and endocardial wall strain. Peak pressure drop at the atrioventricular canal (ΔPAVC) and outflow tract (ΔPOFT) show a steady increase with fish age progression. Pressure drops when non-dimensionalized with blood viscosity and heart rate at each dpf are comparable with measurements in zebrafish hearts. Retrograde flows captured at these regions display a negative pressure drop. A novel metric, Endocardial Work (EW), is introduced by analyzing the ΔPAVC-strain curves, which is a non-invasive measure of work required for ventricle filling. EW is a metric that can differentiate between the linear heart stage (< 100 Pa-%), cardiac looped chamber stage (< 300 Pa-%), and the fully formed chamber stage (> 300 Pa-%).

Chronotropic and vasoactive properties of the gut bacterial short-chain fatty acids in larval zebrafish

Short-chain fatty acids (SCFAs) produced by the gut bacteria have been associated with cardiovascular dysfunction in humans and rodents. However, studies exploring effects of SCFAs on cardiovascular parameters in the zebrafish, an increasingly popular model in cardiovascular research, remain limited. Here, we performed fecal bacterial 16S sequencing and gas chromatography/mass spectrometry (GC-MS) to determine the composition and abundance of gut microbiota and SCFAs in adult zebrafish. Following this, the acute effects of major SCFAs on heart rate and vascular tone were measured in anesthetized zebrafish larvae using fecal concentrations of butyrate, acetate, and propionate. Finally, we investigated if coincubation with butyrate may lessen the effects of angiotensin II (ANG II) and phenylephrine (PE) on vascular tone in anesthetized zebrafish larvae. We found that the abundance in Proteobacteria, Firmicutes, and Fusobacteria phyla in the adult zebrafish resembled those reported in rodents and humans. SCFA levels with highest concentration of acetate (27.43 µM), followed by butyrate (2.19 µM) and propionate (1.65 µM) were observed in the fecal samples of adult zebrafish. Immersion in butyrate and acetate produced a ∼20% decrease in heart rate (HR), respectively, with no observed effects of propionate. Butyrate alone also produced an ∼25% decrease in the cross-sectional width of the dorsal aorta (DA) at 60 min (*P < 0.05), suggesting compensatory vasoconstriction, with no effects of either acetate or propionate. In addition, butyrate significantly alleviated the decrease in DA cross-sectional width produced by both ANG II and PE. We demonstrate the potential for zebrafish in investigation of host-microbiota interactions in cardiovascular health.NEW & NOTEWORTHY We highlight the presence of a core gut microbiota and demonstrate in vivo short-chain fatty acid production in adult zebrafish. In addition, we show cardio-beneficial vasoactive and chronotropic properties of butyrate, and chronotropic properties of acetate in anesthetized zebrafish larvae.

Small fish, big discoveries: zebrafish shed light on microbial biomarkers for neuro-immune-cardiovascular health

Zebrafish (Danio rerio) have emerged as a powerful model to study the gut microbiome in the context of human conditions, including hypertension, cardiovascular disease, neurological disorders, and immune dysfunction. Here, we highlight zebrafish as a tool to bridge the gap in knowledge in linking the gut microbiome and physiological homeostasis of cardiovascular, neural, and immune systems, both independently and as an integrated axis. Drawing on zebrafish studies to date, we discuss challenges in microbiota transplant techniques and gnotobiotic husbandry practices. We present advantages and current limitations in zebrafish microbiome research and discuss the use of zebrafish in identification of microbial enterotypes in health and disease. We also highlight the versatility of zebrafish studies to further explore the function of human conditions relevant to gut dysbiosis and reveal novel therapeutic targets.

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.

Ultrasound deep learning for monitoring of flow-vessel dynamics in murine carotid artery

Several arterial diseases are closely related with mechanical properties of the blood vessel and interactions of flow-vessel dynamics such as mean flow velocity, wall shear stress (WSS) and vascular strain. However, there is an opportunity to improve the measurement accuracy of vascular properties and hemodynamics by adopting deep learning-based ultrasound imaging for flow-vessel dynamics (DL-UFV). In this study, the DL-UFV is proposed by devising an integrated neural network for super-resolved localization and vessel wall segmentation, and it is also combined with tissue motion estimation and flow measurement techniques such as speckle image velocimetry and speckle tracking velocimetry for measuring velocity field information of blood flow. Performance of the DL-UFV is verified by comparing with other conventional techniques in tissue-mimicking phantoms. After the performance verification, in vivo feasibility is demonstrated in the murine carotid artery with different pathologies: aging and diabetes mellitus (DM). The mutual comparison of flow-vessel dynamics and histological analyses shows correlations between the immunoreactive region and abnormal flow-vessel dynamics interactions. The DL-UFV improves biases in measurements of velocity, WSS, and strain with up to 4.6-fold, 15.1-fold, and 22.2-fold in the tissue-mimicking phantom, respectively. Mean flow velocities and WSS values of the DM group decrease by 30% and 20% of those of the control group, respectively. Mean flow velocities and WSS values of the aging group (34.11 cm/s and 13.17 dyne/cm²) are slightly smaller than those of the control group (36.22 cm/s and 14.25 dyne/cm²). However, the strain values of the aging and DM groups are much smaller than those of the control group (p < 0.05). This study shows that the DL-UFV performs better than the conventional ultrasound-based flow and strain measurement techniques for measuring vascular stiffness and complicated flow-vessel dynamics. Furthermore, the DL-UFV demonstrates its excellent performance in the analysis of the hemodynamic and hemorheological effects of DM and aging on the flow and vascular characteristics. This work provides useful hemodynamic information, including mean flow velocity, WSS and strain with high-resolution for diagnosing the pathogenesis of arterial diseases. This information can be used for monitoring progression and regression of atherosclerotic diseases in clinical practice.

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.

A hybrid of light-field and light-sheet imaging to study myocardial function and intracardiac blood flow during zebrafish development

Biomechanical forces intimately contribute to cardiac morphogenesis. However, volumetric imaging to investigate the cardiac mechanics with high temporal and spatial resolution remains an imaging challenge. We hereby integrated light-field microscopy (LFM) with light-sheet fluorescence microscopy (LSFM), coupled with a retrospective gating method, to simultaneously access myocardial contraction and intracardiac blood flow at 200 volumes per second. While LSFM allows for the reconstruction of the myocardial function, LFM enables instantaneous acquisition of the intracardiac blood cells traversing across the valves. We further adopted deformable image registration to quantify the ventricular wall displacement and particle tracking velocimetry to monitor intracardiac blood flow. The integration of LFM and LSFM enabled the time-dependent tracking of the individual blood cells and the differential rates of segmental wall displacement during a cardiac cycle. Taken together, we demonstrated a hybrid system, coupled with our image analysis pipeline, to simultaneously capture the myocardial wall motion with intracardiac blood flow during cardiac development.

During the conception of the heart, cardiac muscular contraction and blood flow generate biomechanical forces to influence the functional and structural development. To elucidate the underlying biomechanical mechanisms, we have embraced the zebrafish system for the ease of genetic and pharmacological manipulations and its rapidity for organ development. However, acquiring the dynamic processes (space + time domain) in the small beating zebrafish heart remains a challenge. In the presence of a rapid heartbeat, microscopy is confined by temporal resolution to image the cardiac contraction and blood flow. In this context, we demonstrated an integrated light-sheet and light-field imaging system to visualize cardiac contraction along with the flowing blood cells inside the cardiac chambers. Assuming the periodicity of the cardiac cycle, we synchronized the image data in post-processing for 3-D reconstruction. We further quantified the velocity of the various regions of cardiac muscular contraction, and tracked the individual blood cells during the cardiac cycles. The time-dependent velocity maps allow for uncovering differential segments of cardiac contraction and relaxation, and for revealing the patterns of blood flow. Thus, our integrated light-sheet and light-field imaging system provides an experimental basis to further investigate cardiac function and development.

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.

Effects of crude oil vapors on the cardiovascular flow of embryonic Gulf killifish

Direct contact with toxicants in crude oil during embryogenesis causes cardiovascular defects, but the effects of exposure to airborne volatile organic compounds released from spilled oil are not well understood. The effects of crude oil-derived airborne toxicants on peripheral blood flow were examined in Gulf killifish (Fundulus grandis) since this model completes embryogenesis in the air. Particle image velocimetry was used to measure in vivo blood flow in intersegmental arteries of control and oil-exposed embryos. Significant effects in oil-exposed embryos included increased pulse rate, reduced mean blood flow speed and volumetric flow rate, and decreased pulsatility, demonstrating that normal-appearing oil-exposed embryos retain underlying cardiovascular defects. Further, hematocrit moderately increased in oil-exposed embryos. This study highlights the potential for fine-scale physiological measurement techniques to better understand the sub-lethal effects of oil exposure and demonstrates the efficacy of Gulf killifish as a unique teleost model for aerial toxicant exposure studies.

Ultrasound Deep Learning for Wall Segmentation and Near-Wall Blood Flow Measurement

Studies of medical flow imaging have technical limitations for accurate analysis of blood flow dynamics and vessel wall interaction at arteries. We propose a new deep learning-based boundary detection and compensation (DL-BDC) technique in ultrasound (US) imaging. It can segment vessel boundaries by harnessing the convolutional neural network and wall motion compensation in the analysis of near-wall flow dynamics. The network enables training from real and synthetic US images together. The performance of the technique is validated through synthetic US images and tissue-mimicking phantom experiments. The neural network performs well with high Dice coefficients of over 0.94 and 0.9 for lumens and walls, outperforming previous segmentation techniques. Then, the performance of the wall motion compensation is examined for compliant phantoms. When DL-BDC is applied to flow influenced by wall motion, root-mean-square errors are less than 0.07%. The technique is utilized to analyze flow dynamics and wall interaction with varying elastic moduli of the phantoms. The results show that the flow dynamics and wall shear stress values are consistent with the expected values of the compliant phantoms, and their wall motion behavior is observed with pulse wave propagation. This strategy makes US imaging capable of simultaneous measurement of blood flow and vessel dynamics in human arteries for their accurate interaction analysis. DL-BDC can segment vessel walls fast, accurately, and robustly. It enables to measure the near-wall flow precisely by determining the vessel boundary dynamics. This approach can be beneficial in flow dynamics and wall interaction analyses in various biomedical applications.

Validating the Paradigm That Biomechanical Forces Regulate Embryonic Cardiovascular Morphogenesis and Are Fundamental in the Etiology of Congenital Heart Disease

The goal of this review is to provide a broad overview of the biomechanical maturation and regulation of vertebrate cardiovascular (CV) morphogenesis and the evidence for mechanistic relationships between function and form relevant to the origins of congenital heart disease (CHD). The embryonic heart has been investigated for over a century, initially focusing on the chick embryo due to the opportunity to isolate and investigate myocardial electromechanical maturation, the ability to directly instrument and measure normal cardiac function, intervene to alter ventricular loading conditions, and then investigate changes in functional and structural maturation to deduce mechanism. The paradigm of "Develop and validate quantitative techniques, describe normal, perturb the system, describe abnormal, then deduce mechanisms" was taught to many young investigators by Dr. Edward B. Clark and then validated by a rapidly expanding number of teams dedicated to investigate CV morphogenesis, structure-function relationships, and pathogenic mechanisms of CHD. Pioneering studies using the chick embryo model rapidly expanded into a broad range of model systems, particularly the mouse and zebrafish, to investigate the interdependent genetic and biomechanical regulation of CV morphogenesis. Several central morphogenic themes have emerged. First, CV morphogenesis is inherently dependent upon the biomechanical forces that influence cell and tissue growth and remodeling. Second, embryonic CV systems dynamically adapt to changes in biomechanical loading conditions similar to mature systems. Third, biomechanical loading conditions dynamically impact and are regulated by genetic morphogenic systems. Fourth, advanced imaging techniques coupled with computational modeling provide novel insights to validate regulatory mechanisms. Finally, insights regarding the genetic and biomechanical regulation of CV morphogenesis and adaptation are relevant to current regenerative strategies for patients with CHD.

Contractile and hemodynamic forces coordinate Notch1b-mediated outflow tract valve formation

Biomechanical forces and endothelial-to-mesenchymal transition (EndoMT) are known to mediate valvulogenesis. However, the relative contributions of myocardial contractile and hemodynamic shear forces remain poorly understood. We integrated 4-D light-sheet imaging of transgenic zebrafish models with moving-domain computational fluid dynamics to determine effects of changes in contractile forces and fluid wall shear stress (WSS) on ventriculobulbar (VB) valve development. Augmentation of myocardial contractility with isoproterenol increased both WSS and Notch1b activity in the developing outflow tract (OFT) and resulted in VB valve hyperplasia. Increasing WSS in the OFT, achieved by increasing blood viscosity through EPO mRNA injection, also resulted in VB valve hyperplasia. Conversely, decreasing myocardial contractility by Tnnt2a morpholino oligonucleotide (MO) administration, 2,3-butanedione monoxime treatment, or Plcγ1 inhibition completely blocked VB valve formation, which could not be rescued by increasing WSS or activating Notch. Decreasing WSS in the OFT, achieved by slowing heart rate with metoprolol or reducing viscosity with Gata1a MO, did not affect VB valve formation. Immunofluorescent staining with the mesenchymal marker, DM-GRASP, revealed that biomechanical force-mediated Notch1b activity is implicated in EndoMT to modulate valve morphology. Altogether, increases in WSS result in Notch1b- EndoMT-mediated VB valve hyperplasia, whereas decreases in contractility result in reduced Notch1b activity, absence of EndoMT, and VB valve underdevelopment. Thus, we provide developmental mechanotransduction mechanisms underlying Notch1b-mediated EndoMT in the OFT.

Fluid forces shape the embryonic heart: Insights from zebrafish

Heart formation involves a complex series of tissue rearrangements, during which regions of the developing organ expand, bend, converge, and protrude in order to create the specific shapes of important cardiac components. Much of this morphogenesis takes place while cardiac function is underway, with blood flowing through the rapidly contracting chambers. Fluid forces are therefore likely to influence the regulation of cardiac morphogenesis, but it is not yet clear how these biomechanical cues direct specific cellular behaviors. In recent years, the optical accessibility and genetic amenability of zebrafish embryos have facilitated unique opportunities to integrate the analysis of flow parameters with the molecular and cellular dynamics underlying cardiogenesis. Consequently, we are making progress toward a comprehensive view of the biomechanical regulation of cardiac chamber emergence, atrioventricular canal differentiation, and ventricular trabeculation. In this review, we highlight a series of studies in zebrafish that have provided new insight into how cardiac function can shape cardiac morphology, with a particular focus on how hemodynamics can impact cardiac cell behavior. Over the long-term, this knowledge will undoubtedly guide our consideration of the potential causes of congenital heart disease.

Fluid dynamics in heart development: effects of hematocrit and trabeculation

Recent in vivo experiments have illustrated the importance of understanding the haemodynamics of heart morphogenesis. In particular, ventricular trabeculation is governed by a delicate interaction between haemodynamic forces, myocardial activity, and morphogen gradients, all of which are coupled to genetic regulatory networks. The underlying haemodynamics at the stage of development in which the trabeculae form is particularly complex, given the balance between inertial and viscous forces. Small perturbations in the geometry, scale, and steadiness of the flow can lead to changes in the overall flow structures and chemical morphogen gradients, including the local direction of flow, the transport of morphogens, and the formation of vortices. The immersed boundary method was used to solve the two-dimensional fluid-structure interaction problem of fluid flow moving through a two chambered heart of a zebrafish (Danio rerio), with a trabeculated ventricle, at 96 hours post fertilization (hpf). Trabeculae heights and hematocrit were varied, and simulations were conducted for two orders of magnitude of Womersley number, extending beyond the biologically relevant range (0.2-12.0). Both intracardial and intertrabecular vortices formed in the ventricle for biologically relevant parameter values. The bifurcation from smooth streaming flow to vortical flow depends upon the trabeculae geometry, hematocrit, and Womersley number, $Wo$. This work shows the importance of hematocrit and geometry in determining the bulk flow patterns in the heart at this stage of development.

High-speed extended-volume blood flow measurement using engineered point-spread function

Experimental characterization of blood flow in living organisms is crucial for understanding the development and function of cardiovascular systems, but there has been no technique reported for snapshot imaging of thick samples in large volumes with high precision. We have combined computational microscopy and the diffraction-free, self-bending property of Airy-beams to track fluorescent beads with sub-micron precision through an extended axial range (up to 600 μm) within the flowing blood of 3 days post-fertilization (dpf) zebrafish embryos. The spatial trajectories of the tracer beads within flowing blood were recorded during transit through both cardinal and intersegmental vessels, and the trajectories were found to be consistent with the segmentation of the vasculature recorded using selective-plane illumination microscopy (SPIM). This method provides sufficiently precise spatial and temporal measurement of 3D blood flow that has the potential for directly probing key biomechanical quantities such as wall shear stress, as well as exploring the fluidic repercussions of cardiovascular diseases. Although we demonstrate the technique for blood flow, the ten-fold better enhancement in the depth range offers improvements in a wide range of applications of high-speed precision measurement of fluid flow, from microfluidics through measurement of cell dynamics to macroscopic aerosol characterizations.

Hemodynamics in Cardiac Development

The beating heart is subject to intrinsic mechanical factors, exerted by contraction of the myocardium (stretch and strain) and fluid forces of the enclosed blood (wall shear stress). The earliest contractions of the heart occur already in the 10-somite stage in the tubular as yet unsegmented heart. With development, the looping heart becomes asymmetric providing varying diameters and curvatures resulting in unequal flow profiles. These flow profiles exert various wall shear stresses and as a consequence different expression patterns of shear responsive genes. In this paper we investigate the morphological alterations of the heart after changing the blood flow by ligation of the right vitelline vein in a model chicken embryo and analyze the extended expression in the endocardial cushions of the shear responsive gene Tgfbeta receptor III. A major phenomenon is the diminished endocardial-mesenchymal transition resulting in hypoplastic (even absence of) atrioventricular and outflow tract endocardial cushions, which might be lethal in early phases. The surviving embryos exhibit several cardiac malformations including ventricular septal defects and malformed semilunar valves related to abnormal development of the aortopulmonary septal complex and the enclosed neural crest cells. We discuss the results in the light of the interactions between several shear stress responsive signaling pathways including an extended review of the involved Vegf, Notch, Pdgf, Klf2, eNos, Endothelin and Tgfβ/Bmp/Smad networks.

Using Zebrafish for Investigating the Molecular Mechanisms of Drug-Induced Cardiotoxicity

Over the last decade, the zebrafish (Danio rerio) has emerged as a model organism for cardiovascular research. Zebrafish have several advantages over mammalian models. For instance, the experimental cost of using zebrafish is comparatively low; the embryos are transparent, develop externally, and have high fecundity making them suitable for large-scale genetic screening. More recently, zebrafish embryos have been used for the screening of a variety of toxic agents, particularly for cardiotoxicity testing. Zebrafish has been shown to exhibit physiological responses that are similar to mammals after exposure to medicinal drugs including xenobiotics, hormones, cancer drugs, and also environmental pollutants, including pesticides and heavy metals. In this review, we provided a summary for recent studies that have used zebrafish to investigate the molecular mechanisms of drug-induced cardiotoxicity. More specifically, we focused on the techniques that were exploited by us and others for cardiovascular toxicity assessment and described several microscopic imaging and analysis protocols that are being used for the estimation of a variety of cardiac hemodynamic parameters.

Influence of blood flow on cardiac development

The role of hemodynamics in cardiovascular development is not well understood. Indeed, it would be remarkable if it were, given the dauntingly complex array of intricately synchronized genetic, molecular, mechanical, and environmental factors at play. However, with congenital heart defects affecting around 1 in 100 human births, and numerous studies pointing to hemodynamics as a factor in cardiovascular morphogenesis, this is not an area in which we can afford to remain in the dark. This review seeks to present the case for the importance of research into the biomechanics of the developing cardiovascular system. This is accomplished by i) illustrating the basics of some of the highly complex processes involved in heart development, and discussing the known influence of hemodynamics on those processes; ii) demonstrating how altered hemodynamic environments have the potential to bring about morphological anomalies, citing studies in multiple animal models with a variety of perturbation methods; iii) providing examples of widely used technological innovations which allow for accurate measurement of hemodynamic parameters in embryos; iv) detailing the results of studies in avian embryos which point to exciting correlations between various hemodynamic manipulations in early development and phenotypic defect incidence in mature hearts; and finally, v) stressing the relevance of uncovering specific biomechanical pathways involved in cardiovascular formation and remodeling under adverse conditions, to the potential treatment of human patients. The time is ripe to unravel the contributions of hemodynamics to cardiac development, and to recognize their frequently neglected role in the occurrence of heart malformation phenotypes.

3D + time blood flow mapping using SPIM-microPIV in the developing zebrafish heart

We present SPIM-μPIV as a flow imaging system, capable of measuring in vivo flow information with 3D micron-scale resolution. Our system was validated using a phantom experiment consisting of a flow of beads in a 50 μm diameter FEP tube. Then, with the help of optical gating techniques, we obtained 3D + time flow fields throughout the full heartbeat in a ∼3 day old zebrafish larva using fluorescent red blood cells as tracer particles. From this we were able to recover 3D flow fields at 31 separate phases in the heartbeat. From our measurements of this specimen, we found the net pumped blood volume through the atrium to be 0.239 nL per beat. SPIM-μPIV enables high quality in vivo measurements of flow fields that will be valuable for studies of heart function and fluid-structure interaction in a range of small-animal models.

Hemodynamic Studies for Analyzing the Teratogenic Effects of Drugs in the Zebrafish Embryo

Investigations of teratogenic effects of drugs generally involve testing the drug on animals and zebrafish embryo is a commonly used animal model for that purpose. In these studies, cardiovascular function of the animals needs to be evaluated to reveal the influence of exposure on the development of the cardiovascular system as well as on the growth of the whole animal. Here, relevant microscopy imaging and analysis protocols are described to calculate a variety of hemodynamic parameters for zebrafish embryos exposed to clinical drugs.

Anisotropic shear stress patterns predict the orientation of convergent tissue movements in the embryonic heart

Myocardial contractility and blood flow provide essential mechanical cues for the morphogenesis of the heart. In general, endothelial cells change their migratory behavior in response to shear stress patterns, according to flow directionality. Here, we assessed the impact of shear stress patterns and flow directionality on the behavior of endocardial cells, the specialized endothelial cells of the heart. At the early stages of zebrafish heart valve formation, we show that endocardial cells are converging to the valve-forming area and that this behavior depends upon mechanical forces. Quantitative live imaging and mathematical modeling allow us to correlate this tissue convergence with the underlying flow forces. We predict that tissue convergence is associated with the direction of the mean wall shear stress and of the gradient of harmonic phase-averaged shear stresses, which surprisingly do not match the overall direction of the flow. This contrasts with the usual role of flow directionality in vascular development and suggests that the full spatial and temporal complexity of the wall shear stress should be taken into account when studying endothelial cell responses to flow in vivo.

Summary: Blood flow modeling shows that dynamic shear stress patterns, rather than mean flow direction, predict the stereotypical behavior of endocardial cells during the early steps of heart valve formation.

Variation in wall shear stress in channel networks of zebrafish models

Physiological functions of vascular endothelial cells (ECs) vary depending on wall shear stress (WSS) magnitude, and the functional change affects the pathologies of various cardiovascular systems. Several in vitro and in vivo models have been used to investigate the functions of ECs under different WSS conditions. However, these models have technical limitations in precisely mimicking the physiological environments of ECs and monitoring temporal variations of ECs in detail. Although zebrafish (Danio rerio) has several strategies to overcome these technical limitations, zebrafish cannot be used as a perfect animal model because applying various WSS conditions on blood vessels of zebrafish is difficult. This study proposes a new zebrafish model in which various WSS can be applied to the caudal vein. The WSS magnitude is controlled by blocking some parts of blood-vessel networks. The accuracy and reproducibility of the proposed method are validated using an equivalent circuit model of blood vessels in zebrafish. The proposed method is applied to lipopolysaccharide (LPS)-stimulated zebrafish as a typical application. The proposed zebrafish model can be used as an in vivo animal model to investigate the relationship between WSS and EC physiology or WSS-induced cardiovascular diseases.

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.

Endoglin controls blood vessel diameter through endothelial cell shape changes in response to haemodynamic cues

The hierarchical organization of properly sized blood vessels ensures the correct distribution of blood to all organs of the body, and is controlled via haemodynamic cues. In current concepts, an endothelium-dependent shear stress set point causes blood vessel enlargement in response to higher flow rates, while lower flow would lead to blood vessel narrowing, thereby establishing homeostasis. We show that during zebrafish embryonic development increases in flow, after an initial expansion of blood vessel diameters, eventually lead to vessel contraction. This is mediated via endothelial cell shape changes. We identify the transforming growth factor beta co-receptor endoglin as an important player in this process. Endoglin mutant cells and blood vessels continue to enlarge in response to flow increases, thus exacerbating pre-existing embryonic arterial-venous shunts. Together, our data suggest that cell shape changes in response to biophysical cues act as an underlying principle allowing for the ordered patterning of tubular organs.

Heart function and hemodynamic analysis for zebrafish embryos

The Zebrafish has emerged to become a powerful vertebrate animal model for cardiovascular research in recent years. Its advantages include easy genetic manipulation, transparency, small size, low cost, and the ability to survive without active circulation at early stages of development. Sequencing the whole genome and identifying ortholog genes with human genome made it possible to induce clinically relevant cardiovascular defects via genetic approaches. Heart function and disturbed hemodynamics need to be assessed in a reliable manner for these disease models in order to reveal the mechanobiology of induced defects. This effort requires precise determination of blood flow patterns as well as hemodynamic stress (i.e., wall shear stress and pressure) levels within the developing heart. While traditional approach involves time-lapse brightfield microscopy to track cell and tissue movements, in more recent studies fast light-sheet fluorescent microscopes are utilized for that purpose. Integration of more complicated techniques like particle image velocimetry and computational fluid dynamics modeling for hemodynamic analysis holds a great promise to the advancement of the Zebrafish studies. Here, we discuss the latest developments in heart function and hemodynamic analysis for Zebrafish embryos and conclude with our future perspective on dynamic analysis of the Zebrafish cardiovascular system. Developmental Dynamics 246:868-880, 2017. © 2017 Wiley Periodicals, Inc.

The Impact of Fluid Inertia on In Vivo Estimation of Mitral Valve Leaflet Constitutive Properties and Mechanics

We examine the influence of the added mass effect (fluid inertia) on mitral valve leaflet stress during isovolumetric phases. To study this effect, oscillating flow is applied to a flexible membrane at various frequencies to control inertia. Resulting membrane strain is calculated through a three-dimensional reconstruction of markers from stereo images. To investigate the effect in vivo, the analysis is repeated on a published dataset for an ovine mitral valve (Journal of Biomechanics 42(16): 2697-2701). The membrane experiment demonstrates that the relationship between pressure and strain must be corrected with a fluid inertia term if the ratio of inertia to pressure differential approaches 1. In the mitral valve, this ratio reaches 0.7 during isovolumetric contraction for an acceleration of 6 m/s(2). Acceleration is reduced by 72% during isovolumetric relaxation. Fluid acceleration also varies along the leaflet during isovolumetric phases, resulting in spatial variations in stress. These results demonstrate that fluid inertia may be the source of the temporally and spatially varying stiffness measurements previously seen through inverse finite element analysis of in vivo data during isovolumetric phases. This study demonstrates that there is a need to account for added mass effects when analyzing in vivo constitutive relationships of heart valves.

Interplay between cardiac function and heart development

Mechanotransduction refers to the conversion of mechanical forces into biochemical or electrical signals that initiate structural and functional remodeling in cells and tissues. The heart is a kinetic organ whose form changes considerably during development and disease. This requires cardiomyocytes to be mechanically durable and able to mount coordinated responses to a variety of environmental signals on different time scales, including cardiac pressure loading and electrical and hemodynamic forces. During physiological growth, myocytes, endocardial and epicardial cells have to adaptively remodel to these mechanical forces. Here we review some of the recent advances in the understanding of how mechanical forces influence cardiac development, with a focus on fluid flow forces. This article is part of a Special Issue entitled: Cardiomyocyte Biology: Integration of Developmental and Environmental Cues in the Heart edited by Marcus Schaub and Hughes Abriel.

Association of Early Atherosclerosis with Vascular Wall Shear Stress in Hypercholesterolemic Zebrafish

Although atherosclerosis is a multifactorial disease, the role of hemodynamic information has become more important. Low and oscillating wall shear stress (WSS) that changes its direction is associated with the early stage of atherosclerosis. Several in vitro and in vivo models were proposed to reveal the relation between the WSS and the early atherosclerosis. However, these models possess technical limitations in mimicking real physiological conditions and monitoring the developmental course of the early atherosclerosis. In this study, a hypercholesterolaemic zebrafish model is proposed as a novel experimental model to resolve these limitations. Zebrafish larvae are optically transparent, which enables temporal observation of pathological variations under in vivo condition. WSS in blood vessels of 15 days post-fertilisation zebrafish was measured using a micro particle image velocimetry (PIV) technique, and spatial distribution of lipid deposition inside the model was quantitatively investigated after feeding high cholesterol diet for 10 days. Lipids were mainly deposited in blood vessel of low WSS. The oscillating WSS was not induced by the blood flows in zebrafish models. The present hypercholesterolaemic zebrafish would be used as a potentially useful model for in vivo study about the effects of low WSS in the early atherosclerosis.

Live imaging and modeling for shear stress quantification in the embryonic zebrafish heart

Hemodynamic shear stress is sensed by the endocardial cells composing the inner cell layer of the heart, and plays a major role in cardiac morphogenesis. Yet, the underlying hemodynamics and the associated mechanical stimuli experienced by endocardial cells remains poorly understood. Progress in the field has been hampered by the need for high temporal resolution imaging allowing the flow profiles generated in the beating heart to be resolved. To fill this gap, we propose a method to analyze the wall dynamics, the flow field, and the wall shear stress of the developing zebrafish heart. This method combines live confocal imaging and computational fluid dynamics to overcome difficulties related to live imaging of blood flow in the developing heart. To provide an example of the applicability of the method, we discuss the hemodynamic frequency content sensed by endocardial cells at the onset of valve formation, and how the fundamental frequency of the wall shear stress represents a unique mechanical cue to endocardial, heart-valve precursors.

Measurement of real pulsatile blood flow using X-ray PIV technique with CO2 microbubbles

Synchrotron X-ray imaging technique has been used to investigate biofluid flows in a non-destructive manner. This study aims to investigate the feasibility of the X-ray PIV technique with CO2 microbubbles as flow tracer for measurement of pulsatile blood flows under in vivo conditions. The traceability of CO2 microbubbles in a pulsatile flow was demonstrated through in vitro experiment. A rat extracorporeal bypass loop was used by connecting a tube between the abdominal aorta and jugular vein of a rat to obtain hemodynamic information of actual pulsatile blood flows without changing the hemorheological properties. The decrease in image contrast of the surrounding tissue was also investigated for in vivo applications of the proposed technique. This technique could be used to accurately measure whole velocity field information of real pulsatile blood flows and has strong potential for hemodynamic diagnosis of cardiovascular diseases.

Fluid-dynamic optimal design of helical vascular graft for stenotic disturbed flow

Although a helical configuration of a prosthetic vascular graft appears to be clinically beneficial in suppressing thrombosis and intimal hyperplasia, an optimization of a helical design has yet to be achieved because of the lack of a detailed understanding on hemodynamic features in helical grafts and their fluid dynamic influences. In the present study, the swirling flow in a helical graft was hypothesized to have beneficial influences on a disturbed flow structure such as stenotic flow. The characteristics of swirling flows generated by helical tubes with various helical pitches and curvatures were investigated to prove the hypothesis. The fluid dynamic influences of these helical tubes on stenotic flow were quantitatively analysed by using a particle image velocimetry technique. Results showed that the swirling intensity and helicity of the swirling flow have a linear relation with a modified Germano number (Gn*) of the helical pipe. In addition, the swirling flow generated a beneficial flow structure at the stenosis by reducing the size of the recirculation flow under steady and pulsatile flow conditions. Therefore, the beneficial effects of a helical graft on the flow field can be estimated by using the magnitude of Gn*. Finally, an optimized helical design with a maximum Gn* was suggested for the future design of a vascular graft.