Pharmacological Modulation of Hemodynamics in Adult Zebrafish In Vivo

Introduction

Hemodynamic parameters in zebrafish receive increasing attention because of their important role in cardiovascular processes such as atherosclerosis, hematopoiesis, sprouting and intussusceptive angiogenesis. To study underlying mechanisms, the precise modulation of parameters like blood flow velocity or shear stress is centrally important. Questions related to blood flow have been addressed in the past in either embryonic or ex vivo-zebrafish models but little information is available for adult animals. Here we describe a pharmacological approach to modulate cardiac and hemodynamic parameters in adult zebrafish in vivo.

Materials and Methods

Adult zebrafish were paralyzed and orally perfused with salt water. The drugs isoprenaline and sodium nitroprusside were directly applied with the perfusate, thus closely resembling the preferred method for drug delivery in zebrafish, namely within the water. Drug effects on the heart and on blood flow in the submental vein were studied using electrocardiograms, in vivo-microscopy and mathematical flow simulations.

Results

Under control conditions, heart rate, blood flow velocity and shear stress varied less than ± 5%. Maximal chronotropic effects of isoprenaline were achieved at a concentration of 50 μmol/L, where it increased the heart rate by 22.6 ± 1.3% (n = 4; p < 0.0001). Blood flow velocity and shear stress in the submental vein were not significantly increased. Sodium nitroprusside at 1 mmol/L did not alter the heart rate but increased blood flow velocity by 110.46 ± 19.64% (p = 0.01) and shear stress by 117.96 ± 23.65% (n = 9; p = 0.03).

Discussion

In this study, we demonstrate that cardiac and hemodynamic parameters in adult zebrafish can be efficiently modulated by isoprenaline and sodium nitroprusside. Together with the suitability of the zebrafish for in vivo-microscopy and genetic modifications, the methodology described permits studying biological processes that are dependent on hemodynamic alterations.     

Analysis of early embryonic great-vessel microcirculation in zebrafish using high-speed confocal μPIV

In the developing cardiovascular system, hemodynamic vascular loading is critical for angiogenesis and cardiovascular adaptation. Normal zebrafish embryos with transgenically-labeled endothelial and red blood cells provide an excellent in vivo model for studying the fluid-flow induced vascular loading. To characterize the developmental hemodynamics of early embryonic great-vessel microcirculation in the zebrafish embryo, two complementary studies (experimental and numerical) are presented. Quantitative comparison of the wall shear stress (WSS) at the first aortic arch (AA1) of wild-type zebrafish embryos during two consecutive developmental stages is presented, using time-resolved confocal micro-particle image velocimetry (μPIV). Analysis showed that there was significant WSS difference between 32 and 48 h post-fertilization (hpf) wild-type embryos, which correlates with normal arch morphogenesis. The vascular distensibility of the arch wall at systole and the acceleration/deceleration rates of time-lapse phase-averaged streamwise blood flow curves were also analyzed. To estimate the influence of a novel intermittent red-blood cell (RBC) loading on the endothelium, a numerical two-phase, volume of fluid (VOF) flow model was further developed with realistic in vivo conditions. These studies showed that near-wall effects and cell clustering increased WSS augmentation at a minimum of 15% when the distance of RBC from arch vessel wall was less than 3 μm or when RBC cell-to-cell distance was less than 3 μm. When compared to a smooth wall, the WSS augmentation increased by a factor of ~1.4 due to the roughness of the wall created by the endothelial cell profile. These results quantitatively highlight the contribution of individual RBC flow patterns on endothelial WSS in great-vessel microcirculation and will benefit the quantitative understanding of mechanotransduction in embryonic great vessel biology, including arteriovenous malformations (AVM).     

Flow-pressure drop measurement and calculation in a tapered femoral artery of a dog

This study describes the in vivo measurement of pressure drop and flow during the cardiac cycle in the femoral artery of a dog, and the computer simulation of the experiment based on the use of the measured flow, vessel dimensions and blood viscosity. In view of the experimental uncertainty in obtaining the accurate velocity profile at the wall region, the velocity pulse at the center was measured and numerical calculations were performed for the center Une instantaneous velocity and within the two limits of spatial distribution of inlet flow conditions: uniform and parabolic. Temporal and spatial variations of flow parameters, i.e., velocity profile, shear rate, non-Newtonian viscosity, wall shear stress, and pressure drop were calculated. There existed both positive and negative shear rates during a pulse cycle, i.e., the arterial wall experiences zero shear three times during a cardiac cycle. For the parabolic inlet condition, the taper of the artery not only increased the magnitude of the positive and negative shear rates, but caused a steep gradient in shear rate, a phenomenon which in turn affects wall shear stress and pressure. In contrast, for the uniform inlet condition, the flow through the tapered artery was predominantly the developing type, which resulted in reduction in magnitude of wall shear rate along the axial direction.     

In vivo imaging of embryonic vascular development using transgenic zebrafish

In this study we describe a model system that allows continuous in vivo observation of the vertebrate embryonic vasculature. We find that the zebrafish fli1 promoter is able to drive expression of enhanced green fluorescent protein (EGFP) in all blood vessels throughout embryogenesis. We demonstrate the utility of vascular-specific transgenic zebrafish in conjunction with time-lapse multiphoton laser scanning microscopy by directly observing angiogenesis within the brain of developing embryos. Our images reveal that blood vessels undergoing active angiogenic growth display extensive filopodial activity and pathfinding behavior similar to that of neuronal growth cones. We further show, using the zebrafish mindbomb mutant as an example, that the expression of EGFP within developing blood vessels permits detailed analysis of vascular defects associated with genetic mutations. Thus, these transgenic lines allow detailed analysis of both wild type and mutant embryonic vasculature and, together with the ability to perform large scale forward-genetic screens in zebrafish, will facilitate identification of new mutants affecting vascular development.     

Joint dynamic imaging of morphogenesis and function in the developing heart

In the developing heart, time-lapse imaging is particularly challenging. Changes in heart morphology due to tissue growth or long-term reorganization are difficult to follow because they are much subtler than the rapid shape changes induced by the heartbeat. Therefore, imaging heart development usually requires slowing or stopping the heart. This, however, leads to information loss about the unperturbed heart shape and the dynamics of heart function. To overcome this limitation, we have developed a non-invasive heart imaging technique to jointly document heart function (at fixed stages of development) as well as its morphogenesis (at any fixed phase in the heartbeat) that does not require stopping or slowing the heart. We review the challenges for imaging heart development and our methodology, which is based on computationally combining and analyzing multiple high-speed image sequences acquired throughout the course of development. We present results obtained in the developing zebrafish heart. Image analysis of the acquired data yielded blood flow velocity maps and made it possible to follow the relative movement of individual cells over several hours.Key words: cardiac imaging, zebrafish, fluorescence imaging, heart development, registration, fast imaging     

Joint dynamic imaging of morphogenesis and function in the developing heart

In the developing heart, time-lapse imaging is particularly challenging. Changes in heart morphology due to tissue growth or long-term reorganization are difficult to follow because they are much subtler than the rapid shape changes induced by the heartbeat. Therefore, imaging heart development usually requires slowing or stopping the heart. This, however, leads to information loss about the unperturbed heart shape and the dynamics of heart function. To overcome this limitation, we have developed a non-invasive heart imaging technique to jointly document heart function (at fixed stages of development) as well as its morphogenesis (at any fixed phase in the heartbeat) that does not require stopping or slowing the heart. We review the challenges for imaging heart development and our methodology, which is based on computationally combining and analyzing multiple high-speed image sequences acquired throughout the course of development. We present results obtained in the developing zebrafish heart. Image analysis of the acquired data yielded blood flow velocity maps and made it possible to follow the relative movement of individual cells over several hours.     

Radial artery wall shear stress evaluation in patients with arteriovenous fistula for hemodialysis access

It has been extensively documented that changes in blood flow induce vascular remodeling and this phenomenon seems to be correlated to the shear forces imposed on the vessel wall by motion of blood. Wall shear stress, the tractive force that acts on the endothelium, has been shown to influence endothelial cell function. To study changes in wall shear stress that develop on the vessel wall upon changes of blood flow, we set up a technique that allows estimation of shear stress in the radial artery of patients on chronic hemodialysis therapy. The technique is based on color-flow Doppler examination of the radial artery before and after surgical creation of radiocephalic fistula for hemodialysis. Calculation of time function wall shear stress and blood flow rate in the radial artery is performed on the basis of arterial diameter, center-line velocity waveform and blood viscosity, using a numerical method developed according to Womersley's theory for pulsatile flow in tubes. The results presented confirm that the model developed is suitable for calculation of the wall shear stress that develops in the radial artery of patients before and after surgical creation of an arteriovenous fistula for hemodialysis. This methodology was developed for characterization of wall shear stress in the radial artery but may be well applied to other vessels that can be examined by echo-Doppler technique.     

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.     

Notch1b and neuregulin are required for specification of central cardiac conduction tissue

Normal heart function is critically dependent on the timing and coordination provided by a complex network of specialized cells: the cardiac conduction system. We have employed functional assays in zebrafish to explore early steps in the patterning of the conduction system that previously have been inaccessible. We demonstrate that a ring of atrioventricular conduction tissue develops at 40 hours post-fertilization in the zebrafish heart. Analysis of the mutant cloche reveals a requirement for endocardial signals in the formation of this tissue. The differentiation of these specialized cells, unlike that of adjacent endocardial cushions and valves, is not dependent on blood flow or cardiac contraction. Finally, both neuregulin and notch1b are necessary for the development of atrioventricular conduction tissue. These results are the first demonstration of the endocardial signals required for patterning central ;slow' conduction tissue, and they reveal the operation of distinct local endocardial-myocardial interactions within the developing heart tube.     

4D subject-specific inverse modeling of the chick embryonic heart outflow tract hemodynamics

Blood flow plays a critical role in regulating embryonic cardiac growth and development, with altered flow leading to congenital heart disease. Progress in the field, however, is hindered by a lack of quantification of hemodynamic conditions in the developing heart. In this study, we present a methodology to quantify blood flow dynamics in the embryonic heart using subject-specific computational fluid dynamics (CFD) models. While the methodology is general, we focused on a model of the chick embryonic heart outflow tract (OFT), which distally connects the heart to the arterial system, and is the region of origin of many congenital cardiac defects. Using structural and Doppler velocity data collected from optical coherence tomography, we generated 4D (\(\hbox {3D}\,+\,\hbox {time}\)) embryo-specific CFD models of the heart OFT. To replicate the blood flow dynamics over time during the cardiac cycle, we developed an iterative inverse-method optimization algorithm, which determines the CFD model boundary conditions such that differences between computed velocities and measured velocities at one point within the OFT lumen are minimized. Results from our developed CFD model agree with previously measured hemodynamics in the OFT. Further, computed velocities and measured velocities differ by \(<\)15 % at locations that were not used in the optimization, validating the model. The presented methodology can be used in quantifications of embryonic cardiac hemodynamics under normal and altered blood flow conditions, enabling an in-depth quantitative study of how blood flow influences cardiac development.     

Fluid Dynamics of Heart Development

The morphology, muscle mechanics, fluid dynamics, conduction properties, and molecular biology of the developing embryonic heart have received much attention in recent years due to the importance of both fluid and elastic forces in shaping the heart as well as the striking relationship between the heart’s evolution and development. Although few studies have directly addressed the connection between fluid dynamics and heart development, a number of studies suggest that fluids may play a key role in morphogenic signaling. For example, fluid shear stress may trigger biochemical cascades within the endothelial cells of the developing heart that regulate chamber and valve morphogenesis. Myocardial activity generates forces on the intracardiac blood, creating pressure gradients across the cardiac wall. These pressures may also serve as epigenetic signals. In this article, the fluid dynamics of the early stages of heart development is reviewed. The relevant work in cardiac morphology, muscle mechanics, regulatory networks, and electrophysiology is also reviewed in the context of intracardial fluid dynamics.     

In vivo micro particle image velocimetry measurements of blood-plasma in the embryonic avian heart

The measurement of blood-plasma velocity distributions with spatial and temporal resolution in vivo is inevitable for the determination of shear stress distributions in complex geometries at unsteady flow conditions like in the beating heart. A non-intrusive, whole-field velocity measurement technique is required that is capable of measuring instantaneous flow fields at sub-millimeter scales in highly unsteady flows. Micro particle image velocimetry (muPIV) meets these demands, but requires special consideration and methodologies in order to be utilized for in vivo studies in medical and biological research. We adapt muPIV to measure the blood-plasma velocity in the beating heart of a chicken embryo. In the current work, bio-inert, fluorescent liposomes with a nominal diameter of 400 nm are added to the flow as a tracer. Because of their small dimension and neutral buoyancy the liposomes closely follow the movement of the blood-plasma and allow the determination of the velocity gradient close to the wall. The measurements quantitatively resolve the velocity distribution in the developing ventricle and atrium of the embryo at nine different stages within the cardiac cycle. Up to 400 velocity vectors per measurement give detailed insight into the fluid dynamics of the primitive beating heart. A rapid peristaltic contraction accelerates the flow to peak velocities of 26 mm/s, with the velocity distribution showing a distinct asymmetrical profile in the highly curved section of the outflow tract. In relation to earlier published gene-expression experiments, the results underline the significance of fluid forces for embryonic cardiogenesis. In general, the measurements demonstrate that muPIV has the potential to develop into a general tool for instationary flow conditions in complex flow geometries encountered in cardiovascular research.     

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.     

Fluid Dynamics of Ventricular Filling in the Embryonic Heart

The vertebrate embryonic heart first forms as a valveless tube that pumps blood using waves of contraction. As the heart develops, the atrium and ventricle bulge out from the heart tube, and valves begin to form through the expansion of the endocardial cushions. As a result of changes in geometry, conduction velocities, and material properties of the heart wall, the fluid dynamics and resulting spatial patterns of shear stress and transmural pressure change dramatically. Recent work suggests that these transitions are significant because fluid forces acting on the cardiac walls, as well as the activity of myocardial cells that drive the flow, are necessary for correct chamber and valve morphogenesis. In this article, computational fluid dynamics was used to explore how spatial distributions of the normal forces acting on the heart wall change as the endocardial cushions grow and as the cardiac wall increases in stiffness. The immersed boundary method was used to simulate the fluid-moving boundary problem of the cardiac wall driving the motion of the blood in a simplified model of a two-dimensional heart. The normal forces acting on the heart walls increased during the period of one atrial contraction because inertial forces are negligible and the ventricular walls must be stretched during filling. Furthermore, the force required to fill the ventricle increased as the stiffness of the ventricular wall was increased. Increased endocardial cushion height also drastically increased the force necessary to contract the ventricle. Finally, flow in the moving boundary model was compared to flow through immobile rigid chambers, and the forces acting normal to the walls were substantially different.     

Comparison between a generalized Newtonian model and a network-type multiscale model for hemodynamic behavior in the aortic arch: Validation with 4D MRI data for a case study

Blood is a complex fluid in which the presence of the various constituents leads to significant changes in its rheological properties. Thus, an appropriate non-Newtonian model is advisable; and we choose a Modified version of the rheological model of Phan-Thien and Tanner (MPTT). The different parameters of this model, derived from the rheology of polymers, allow characterization of the non-Newtonian nature of blood, taking into account the behavior of red blood cells in plasma. Using the MPTT model that we implemented in the open access software OpenFOAM, numerical simulations have been performed on blood flow in the thoracic aorta for a healthy patient. We started from a patient-specific model which was constructed from medical images. Exiting flow boundary conditions have been developped, based on a 3-element Windkessel model to approximate physiological conditions. The parameters of the Windkessel model were calibrated with in vivo measurements of flow rate and pressure. The influence of the selected viscosity of red blood cells on the flow and wall shear stress (WSS) was investigated. Results obtained from this model were compared to those of the Newtonian model, and to those of a generalized Newtonian model, as well as to in vivo dynamic data from 4D MRI during a cardiac cycle. Upon evaluating the results, the MPTT model shows better agreement with the MRI data during the systolic and diastolic phases than the Newtonian or generalized Newtonian model, which confirms our interest in using a complex viscoelastic model.     

Vascular remodeling of the mouse yolk sac requires hemodynamic force

The embryonic heart and vessels are dynamic and form and remodel while functional. Much has been learned about the genetic mechanisms underlying the development of the cardiovascular system, but we are just beginning to understand how changes in heart and vessel structure are influenced by hemodynamic forces such as shear stress. Recent work has shown that vessel remodeling in the mouse yolk sac is secondarily effected when cardiac function is reduced or absent. These findings indicate that proper circulation is required for vessel remodeling, but have not defined whether the role of circulation is to provide mechanical cues, to deliver oxygen or to circulate signaling molecules. Here, we used time-lapse confocal microscopy to determine the role of fluid-derived forces in vessel remodeling in the developing murine yolk sac. Novel methods were used to characterize flows in normal embryos and in embryos with impaired contractility (Mlc2a(-/-)). We found abnormal plasma and erythroblast circulation in these embryos, which led us to hypothesize that the entry of erythroblasts into circulation is a key event in triggering vessel remodeling. We tested this by sequestering erythroblasts in the blood islands, thereby lowering the hematocrit and reducing shear stress, and found that vessel remodeling and the expression of eNOS (Nos3) depends on erythroblast flow. Further, we rescued remodeling defects and eNOS expression in low-hematocrit embryos by restoring the viscosity of the blood. These data show that hemodynamic force is necessary and sufficient to induce vessel remodeling in the mammalian yolk sac.     

An Investigation on Platelet Transport during Thrombus Formation at Micro-Scale Stenosis

This paper reports on an investigation of mass transport of blood cells at micro-scale stenosis where local strain-rate micro-gradients trigger platelet aggregation. Using a microfluidic flow focusing platform we investigate the blood flow streams that principally contribute to platelet aggregation under shear micro-gradient conditions. We demonstrate that relatively thin surface streams located at the channel wall are the primary contributor of platelets to the developing aggregate under shear gradient conditions. Furthermore we delineate a role for red blood cell hydrodynamic lift forces in driving enhanced advection of platelets to the stenosis wall and surface of developing aggregates. We show that this novel microfluidic platform can be effectively used to study the role of mass transport phenomena driving platelet recruitment and aggregate formation and believe that this approach will lead to a greater understanding of the mechanisms underlying shear-gradient dependent discoid platelet aggregation in the context of cardiovascular diseases such as acute coronary syndromes and ischemic stroke.     

Increased Inlet Blood Flow Velocity Predicts Low Wall Shear Stress in the Cephalic Arch of Patients with Brachiocephalic Fistula Access

BackgroundAn autogenous arteriovenous fistula is the optimal vascular access for hemodialysis. In the case of brachiocephalic fistula, cephalic arch stenosis commonly develops leading to access failure. We have hypothesized that a contribution to fistula failure is low wall shear stress resulting from post-fistula creation hemodynamic changes that occur in the cephalic arch.MethodsTwenty-two subjects with advanced renal failure had brachiocephalic fistulae placed. The following procedures were performed at mapping (pre-operative) and at fistula maturation (8–32 weeks post-operative): venogram, Doppler to measure venous blood flow velocity, and whole blood viscosity. Geometric and computational modeling was performed to determine wall shear stress and other geometric parameters. The relationship between hemodynamic parameters and clinical findings was examined using univariate analysis and linear regression.ResultsThe percent low wall shear stress was linearly related to the increase in blood flow velocity (p < 0.01). This relationship was more significant in non-diabetic patients (p < 0.01) than diabetic patients. The change in global measures of arch curvature and asymmetry also evolve with time to maturation (p < 0.05).ConclusionsThe curvature and hemodynamic changes during fistula maturation increase the percentage of low wall shear stress regions within the cephalic arch. Low wall shear stress may contribute to subsequent neointimal hyperplasia and resultant cephalic arch stenosis. If this hypothesis remains tenable with further studies, ways of protecting the arch through control of blood flow velocity may need to be developed.     

Microinjecting recombinant rainbow trout Ea4-peptide of pro-IGF-I into zebrafish embryos causes abnormal development in heart, red blood cells, and vasculature

E-peptides and mature insulin-like growth factors (IGFs) are produced from pre-pro-IGFs during post-translational processing and co-secreted into the circulation. Previously, we reported that introduction of a transgene encoding the secreted form of rainbow trout (rt) Ea4-peptide or human (h) Eb-peptide into newly fertilized eggs of medaka (Oryzias latipes) and zebrafish (Danio rerio) resulted in developmental defects in heart, red blood cells and vasculature. In addition to vasculature and red blood cell developmental defects, multiple phenocopies of heart developmental defects categorized by developmental arrest at cardiomyocyte, heart tube and heart looping stages were also observed. These results raise a question of whether rtEa4- or hEb-peptide exerts pleiotropic inhibitory effects on heart, vasculature and red blood cell development in fish embryos. To answer this question, various amounts of recombinant rtEa4-peptide were microinjected into zebrafish eggs at 1.5, 2.5 and 5.5 h post-fertilization (hpf). Although a dose-dependent developmental defect in heart, vasculature and red blood cells was observed in embryos microinjected with rtEa4-peptide at 1.5 and 2.5 hpf, the heart development in all of the microinjected embryos was arrested at the cardiomyocyte stage. Furthermore, the mRNA levels of Nkx2.5, GATA5, VEGF, GATA1 and GATA2 genes in defective embryos were significantly reduced by rtEa4-peptide. These results confirm our previous findings that rtEa4- or hEb-peptide exhibits pleiotropic effects in inhibiting heart, vasculature and red blood cell development in zebrafish embryos.