Skip to main content
  • Methodology article
  • Open access
  • Published:

Laser-scanning velocimetry: A confocal microscopy method for quantitative measurement of cardiovascular performance in zebrafish embryos and larvae



The zebrafish Danio rerio is an important model system for drug discovery and to study cardiovascular development. Using a laser-scanning confocal microscope, we have developed a non-invasive method of measuring cardiac performance in zebrafish embryos and larvae that obtains cardiovascular parameters similar to those obtained using Doppler echocardiography in mammals. A laser scan line placed parallel to the path of blood in the dorsal aorta measures blood cell velocity, from which cardiac output and indices of vascular resistance and contractility are calculated.


This technique, called laser-scanning velocimetry, was used to quantify the effects of pharmacological, developmental, and genetic modifiers of cardiac function. Laser-scanning velocimetry was applied to analyze the cardiovascular effects of morpholino knockdown of osmosensing scaffold for MEKK3 (OSM), which when mutated causes the human vascular disease cerebral cavernous malformations. OSM-deficient embryos had a constricted aortic arch and markedly increased peak cell velocity, a characteristic indicator of aortic stenosis.


These data validate laser-scanning velocimetry as a quantitative tool to measure cardiovascular performance for pharmacological and genetic analysis in zebrafish, which requires no specialized equipment other than a laser-scanning confocal microscope.


The zebrafish Danio rerio are widely used in genome-wide screens to identify genes that regulate the development and function of the vertebrate cardiovascular system [13]. Zebrafish are powerful model organisms to define genes relevant to human cardiovascular disease because of their similarity in genetic complexity to mammals, the ease of selective gene knockdown, and their suitability for large-scale mutagenesis screens [4, 5]. A beating heart and functional circulatory system forms within 26 hours of fertilization, and within 48 hours the heart develops into a two-chambered conformation, complete with functional valves. Because of their transparency, severe morphological abnormalities in cardiac development are easily observed, but the functional consequences of these defects are difficult to quantify. Methods are needed for functional analysis of cardiovascular defects that can be utilized by any investigator studying zebrafish cardiovascular physiology. Using an unmodified, commercially-available laser-scanning confocal microscope, we have developed a novel, non-invasive method of measuring cardiac performance in zebrafish embryos that provides measurements similar to those obtained with Doppler echocardiography. The technique, called laser-scanning velocimetry provides a continuous measurement of blood cell velocity that is necessary for estimating cardiac output (velocity integral) and measuring indices of vascular resistance (peak deceleration) and ventricular contractility (peak acceleration).

We first demonstrate that laser-scanning velocimetry is a sensitive and precise tool for measuring the physiological changes during normal cardiovascular development and drug-induced cardiotoxicity. We then use laser-scanning velocimetry to evaluate the cardiovascular function of zebrafish embryos with reduced expression of osmosensing scaffold for MEKK3 (OSM) – a scaffold protein, whose mutation in humans leads to the development of a vascular disorder called cerebral cavernous malformations [6]. Using microangiography we demonstrate that OSM morphant embryos have a constricted aortic arch. Analysis by laser-scanning velocimetry functionally verified the presence of aortic stenosis through a markedly increased peak cellular velocity. These experiments establish laser-scanning velocimetry as a simple yet powerful technique that allows quantitative physiological analysis of cardiac performance in zebrafish embryos.

Results and discussion

Linescan images depict cell flux and velocity within blood vessels

Linescans are images constructed with laser-scanning microscopes by repeatedly imaging a one pixel-wide line. Whereas traditional two-dimensional images provide spatial information in both the x- and y-dimensions, linescans contain spatial information in the x-dimension and temporal information in the y-dimension. When the scan line is placed across a vessel, either perpendicular or parallel to the flow of blood, the resulting linescan images contains information about the dynamics of circulation (Figure 1). Because images of the linear region of interest can be acquired rapidly, linescans provide excellent temporal resolution. In a perpendicular linescan image, the width of a blood cell in the x-dimension of the image is proportional to the physical width of the cell, whereas the length is proportional to the time required for the cell to cross the imaging line (Figure 1b and 1c). Faster moving cells cross the scan line more quickly and appear to be shorter than slower moving cells. In vessels with high cell flux, such as the dorsal aorta, periodic changes in apparent cell length are evident and reflect the periodic change in cellular velocity caused by ventricular contraction and subsequent relaxation (Figure 1c). Cell flux can easily be determined from perpendicular linescans by simply counting the number of cells that completely cross the imaging line in a given unit of time.

Figure 1
figure 1

Acquring linescan images from Danio rerio embryos. (a) In silico microangiogram of a 48 hpf wild-type zebrafish embryo. Regions of interest used for acquiring linescan images are boxed (intersomitic vessel = ISV, dorsal aorta = DA). (b and c) A single plane (XY) image (Top) showing the orientation of scan lines used to create perpendicular linescans (Bottom) of the intersomitic vessels (b) and dorsal aorta (c). The boundary between slow moving (elongated) cells and fast moving (compressed) cells reflects a ventricular contraction and is easily observed in the Dorsal Aorta (c). (d) A parallel linescan acquired from the dorsal aorta shows a series of lines whose slope is inversely proportional to cellular velocity. Distinct regions of acceleration, constant velocity, and deceleration are observed.

Many diagnostic indices of cardiac function can be calculated from a continuous measurement of the velocity of blood cells flowing through the dorsal aorta (Table 1). While cellular velocity can be calculated from the length of a cell in a perpendicular linescan image, cell lengths are difficult to measure within high-flux vessels, and the velocity calculation assumes that the cells are all spherical and have the same dimensions. Furthermore, velocity calculations based on perpendicular linescans are discontinuous: because each measurement represents the average velocity for the duration of time the cell was crossing the scan line, the instantaneous velocity of a cell is not attainable. These limitations are overcome using a scan line placed parallel to the flow of blood (Figure 1d). Parallel linescans acquired from the dorsal aorta contain a series of lines whose slopes alternate in a sigmoidal fashion (Figure 1d, bottom panel). These lines are created when a high-contrast cellular feature (such as the edge of a cell) crosses the scan line. Each scan captures the position of these features along the line. Successively adding these single line images together creates a linescan image where the positions of these features in the x-dimension are recorded as a function of time in the y-dimension. If a cell is stationary, its position will not change along the line and it will appear as a vertical stripe, whereas a cell moving faster than the scan speed will appear as a horizontal streak. Thus velocity is inversely proportional to the slope of the line traced by the cell in the linescan. Three phases of cell traces can be observed in most parallel linescans: convex, linear, and concave, which correspond to phases of accelerating, constant and decelerating velocity, respectively. The continuity of velocity information provided by laser-scanning velocimetry makes it possible to measure blood acceleration and deceleration – parameters previously only attainable using Doppler echocardiography – that are indices of vascular resistance and contractility [79].

Table 1 Measured quantities, their physiological correlates and diagnostic utility

A scan line placed in the middle of the vessel was found to be optimal for velocimetry measurements because the flux of cells crossing the line is much higher than at the extreme periphery (Figure 1c, bottom panel). While blood flow within this long artery should be laminar, significant differences in cellular velocity were not observed as the scan line was moved from the center of the lumen to the periphery. The maximum cellular velocity obtained from a scan line placed in the middle of the lumen of a 48 hpf embryo was 0.2459 ± 0.0102 cm/s, and that obtained from a scan line at the periphery of cell flux was 0.2455 ± 0.0276 cm/s. The apparently constant flow profile may be explained by the requirement of cell motion for laser-scanning velocimetry to measure velocities. While the plasma velocity may be lower very near the vessel wall, velocities in this region cannot be determined by laser-scanning velocimetry because few cells are present (Figure 1c, bottom panel).

Laser-scanning velocimetry is a sensitive technique for assessing drug-induced cardiotoxicity

To demonstrate that laser-scanning velocimetry is a sensitive technique capable of measuring subtle differences in cardiac performance caused by pharmacological modification, we acquired linescans from a two day old zebrafish embryo successively treated with increasing concentrations of the aquatic anaesthetic tricaine, which is known to suppress cardiac function at high doses. Linescans obtained from untreated embryos and those bathed in a cardiac-suppressing dose of tricaine (500 μM) yielded velocity-versus-time profiles similar to those obtained by continuous-wave Doppler echocardiography in mammals (Figure 2a and 2b). The tricaine treatment caused a 40% reduction of peak cell velocity from 0.25 cm/s to 0.15 cm/s (Figure 2a and 2b). As the concentration of tricaine increased, peak velocity decreased with a sigmoidal dose-response relationship (Figure 2c). It has long been appreciated in Doppler echocardiography that peak acceleration is closely correlated with ventricular contractility [7, 10]. Acceleration of blood cells was determined from the derivative of the velocity-versus-time profiles, and was also decreased by tricaine (Figure 2a, b, and 2d).

Figure 2
figure 2

Laser-scanning velocimetry measures the cardiac toxicity of tricaine. (a and b) Velocity and acceleration profiles acquired from the dorsal aorta of 48 hpf embryos untreated (a) or treated (b) with 500 μM tricaine. (c-g) Cardiac parameters of a single Danio rerio embryo treated with successively increasing doses of tricaine. Dose-response relationship of peak velocity (c), peak acceleration (d), stroke volume (e), heart rate (f) and cardiac output (g). Error bars represent the SEM of values obtained from four to six heartbeats in a single embryo. Lines are a non-linear least squares fit to a sigmoidal dose-response curve. Laser-Scanning Velocimetry Data Analyzer smoothing parameters used for trace, velocity, and acceleration smoothing were 5, 60, and 80, respectively.

Stroke volume, the volume of blood ejected by the heart in one beat, can be determined by integrating the continuous velocity versus time graphs obtained by laser-scanning velocimetry and multiplying by the cross-sectional area of the aorta. Since the heart is a pump, stroke volume is an excellent functional measure of cardiac performance. Estimates of stroke volume decreased significantly over a concentration range of tricaine that caused little change in heart rate (Figure 2e and 2f). Heart rates were found to remain relatively steady at low to moderate doses of tricaine (Figure 2f). Because heart rates varied little over the dose range tested, cardiac output decreased proportionately with decreases in stroke volume (Figure 2e and 2g). Our assessment of cardiac output by laser-scanning velocimetry (~ 18.6 nL/min at 48 hpf) is in good agreement with volume-approximating methods of estimating cardiac output but our method is based on actual output and does not assume normal heart morphology [11, 12].

Absolute measurements of stroke volume can be obtained when the scan line is placed along the outflow tract prior to any arterial branches. While it is not always possible to image the outflow tract prior to arterial branching, an estimate of stroke volume can be obtained by measuring cellular velocities along the dorsal aorta immediately caudal to the fusion of the lateral dorsal aortas (Figure 1a). Relative comparisons of stroke volume and cardiac output can be made between embryos of the same developmental stage when the scan line is placed in the same position along the artery. To investigate the dependence of stroke volume and cardiac output estimates on scan line position, measurements were acquired along the dorsal aorta in a rostral location (immediately caudal to the fusion of the lateral dorsal aortas) a caudal location (above the cloaca). When the scan line is placed in the caudal position, peak velocity and acceleration are underestimated by 15–20% (see Additional file 1, panels a-c). Interestingly, stroke volume and cardiac output estimates are not significantly changed when measured with either a rostral or a caudal scan line position (see Additional file 1, panels d and e).

Correlating structure and function in the developing heart

The two-chambered heart of Danio rerio matures into an adult-like morphology by five days post fertilization (dpf). We used laser-scanning velocimetry to correlate differences in cardiac function with developmental changes in heart structure (Figure 3). Because laser-scanning velocimetry is non-invasive, the development of individual embryos can be followed over time. Representative velocimetry profiles from zebrafish at 1, 2, 4, and 5 days post fertilization are shown in Figure 3a. Changes in cardiovascular parameters between days 1 and 5 are shown for both a single fish and a population (Figure 3b–h). Trends in the development of cardiovascular function observed in individuals were representative of the changes observed in the population, indicating the legitimacy and reproducibility of laser-scanning velocimetry.

Figure 3
figure 3

Functional assessment of cardiac performance in the developing zebrafish heart.(a) Representative velocimetry and acceleration profiles from Danio rerio embryos and larvae during development. (b-h) Cardiovascular parameters measured at 1, 2, 4 and 5 days post-fertilization from a single zebrafish (dotted line with triangles) and a population average (solid line with squares). (b) Peak velocity, (c) minimum velocity, (d) peak acceleration, (e) peak deceleration, (f) heart rate (no error bars are present for the single fish measurement because beat-to-beat variability is not calculated), (g) stroke volume, and (h) cardiac output. Single-fish data presented are the average obtained by analyzing four to six individual heart beats. Error bars are the SEM and represent intrafish variability. For population average data n = 5 for days 1, 2, and 4; n = 4 for day 5. Error bars are the SEM and represent interfish variability. Laser-Scanning Velocimetry Data Analyzer smoothing parameters used for trace, velocity, and acceleration smoothing were 5, 10, and 20, respectively.

At day one, the zebrafish heart exists as an elongated, primitive tube that regularly contracts but has not yet developed into a two-chambered heart [13]. Blood is ejected from the one day heart with gradual but rather sustained changes in velocity (Figure 3a). Both acceleration and deceleration are minimal at this stage (Figure 3d and 3e). By day two the heart has differentiated into a two-chambered heart, has looped into its adult conformation, and the atrio-ventricular valves have formed. Velocity-versus-time peaks for each heart beat are sharper and more frequent than at day one (Figure 3a). Notably, by day two there are significant increases in peak velocity and peak acceleration, as well as in estimates of stroke volume, and cardiac output (Figure 3b, d, g, and 3h). Peak deceleration is an index of vascular resistance and is greatest at this time. This increase in resistance provides enough opposing force to the forward motion of blood that some retrograde motion, indicated by negative minimum velocity values, is observed (Figure 3c and 3e). After day two, vascular resistance decreases and minimum velocities are greater than zero (Figure 3c and 3e). Because there is no reversal of flow direction, peak acceleration actually decreases between days two and three. Thus, it should be noted that, as in ultrasound, peak acceleration is a load-dependent measurement and should be used as an index of contractility only when there are no significant differences in vascular resistance.

Between days 2 and 5 changes in heart chamber orientation and differentiation continue, but no overt structural changes occur [13]. This is reflected by only modest changes in peak velocity, stroke volume, and cardiac output between days two and five (Figure 3b–e, and 3g).

Laser-scanning velocimetry provides functional evidence for the presence of an aortic constriction in OSM-deficient embryos

It was recently shown that valentine mutant embryos have a mutation at the ccm2/OSM locus that prevents concentric patterning of the myocardium and results in a loss of blood flow. To demonstrate that laser-scanning velocimetry is a valuable diagnostic tool to study the functional consequences of genetic modification of heart development, we studied the cardiovascular phenotype of embryos where expression of the ccm2 gene product OSM was inhibited by morpholino suppression of translation. Translation-blocking morpholino oligonucleotides targeted at or near the translation start site of the Danio rerio ccm2/osm mRNA were injected into embryos between the 1 and 4 cell stages. OSM protein expression in morpholino-injected embryos was reduced to nearly undetectable levels as demonstrated by western blotting at 48 hpf (Figure 4a and 4b). When normalized to an ERK2 loading control by densitometry, OSM expression was reduced by 98% relative to control embryos (Figure 4b). OSM-deficient embryos had pericardial edema (Figure 4c, arrow), and significantly reduced circulation of erythrocytes. Each of two OSM morpholinos generated the same phenotypes with slightly different potencies (please see Additional file 2). To show that the reduced circulation was specific for the repression of ccm2/osm translation, we rescued the morphant phenotype by co-injecting in vitro transcribed mRNA from the entire coding region of Danio rerio ccm2/osm cDNA along with the morpholino designed to target the 5' untranslated region of the endogenous ccm2/osm mRNA (OSM MO#1). Quantification of cell flux using perpendicular linescan images indicated that erythrocyte flux in OSM-deficient embryos was reduced to 55% that of control embryos, and this phenotype was almost completely restored by co-injection of ccm2/osm mRNA (Figure 4d)

Figure 4
figure 4

Laser-scanning velocimetry provides functional evidence for the presence of an aortic arch constriction in OSM-deficient embryos. (a) Western blotting for OSM in lysates from control or OSM morpholino-injected embryos harvested at 48 hpf shows nearly complete repression of OSM expression. The abundance of ERK2 is shown as a loading control. (b) Intensity of OSM bands normalized to ERK2 by densitometry for western blots shown in Panel a. (c) OSM morphant embryos at 48 hpf exhibited pericardial edema (arrowhead). (d) Co-injection of OSM mRNA along with OSM morpholino restored cell flux to levels similar to that of control embryos. Cell flux was quantified by scan lines drawn perpendicular to the dorsal aorta above the cloaca. Error bars represent the SEM, n = 3. (e) In silico microangiography of the aortic arch (arrowheads) at 28 hpf shows constricted flow through this region of the outflow tract in OSM-deficient embryos. A region of post-stenotic dilation (arrow), the result of turbulent flow, is occasionally observed. Scale bar is 25 μm. (f) Peak blood cell velocity, measured by laser-scanning velocimetry is increased in OSM-deficient embryos. Error bars represent SEM. Asterisk denotes p = 0.0017 for the two-tailed Student's T-test, n = 3. OSM morpholino #1 was used for the data presented here. Laser-Scanning Velocimetry Data Analyzer smoothing parameters used for trace, velocity, and acceleration smoothing were 5, 10, and 20, respectively.

Because homozygous null mice for ccm1/Krit1 die at E8.5 due to a constricted aortic arch [14], we investigated whether the decreased circulation in OSM morphants was due to defective vessel formation. In silico microangiography was used to outline the functional vasculature in control and OSM morphant embryos. In silico microangiography is a simple and non-invasive technique that uses digital image processing of images obtained using widefield microscopy and a CCD camera, and delivers results similar to laborious dye injection-based microangiography methods [15]. Microangiograms of the outflow tract of 26 hpf embryos revealed constricted flow through the aortic arch of morphant but not control embryos (Figure 4e). Movies of the aortic arch clearly show erythrocytes travelling side-by-side through the arch of a control embryo while there is only room for single-file travel through an OSM-deficient embryo (please see Additional file 3). Vessel patterning and lumen formation was unaffected in other vessels (please see Additional file 2), indicating a specific defect in aortic arch patterning.

Additional File 3: restricted circulation movie. Restricted blood flow through constricted aortic arch of control and OSM morpholino-injected embryo (MOV 1 MB)

In human and veterinary cardiology, Doppler echocardiography is used to diagnose aortic stenosis (constrictions). One characteristic change observed in patients with aortic stenosis is an increase in peak velocity because the constricted passage acts as a nozzle to increase the velocity of blood flowing through it [16]. Laser-scanning velocimetry measurements indicated that the peak cellular velocity in OSM morphants was 0.182 ± 0.019 cm/s versus 0.076 ± 0.011 cm/s in control embryos (Figure 4f). These results demonstrate that laser-scanning velocimetry can be used as a diagnostic tool for studying the genetics of cardiac performance as well as identifying genes responsible for congenital heart diseases.


Despite the widespread use of zebrafish for understanding cardiac development and function, few quantitative tools exist to measure functional cardiovascular phenotypes in embryos. Laser-scanning velocimetry meets this need, using a standard laser-scanning confocal microscope for non-invasive measures of cardiac function. High-resolution ultrasound has been used to measure cardiovascular performance in adult zebrafish [17], and technically could be applied to the analysis of the zebrafish embryonic cardiovascular system. Unlike confocal microscopes that are available at many research institutions, high-resolution ultrasound requires expensive and highly specialized equipment that are not generally available to many investigators. We have shown that laser-scanning velocimetry provides sensitive and precise measures of cardiac function including contractility, cardiac output, vascular resistance, heart rate, and stroke volume. Our data demonstrate that laser-scanning velocimetry easily measures drug-induced cardiotoxicity, performance changes during cardiovascular development and can diagnose anatomic lesions such as aortic stenosis in genetically-manipulated embryos.

Laser-scanning velocimetry provides quantitative and precise measures of several indices of cardiac performance (summarized in Table 1) and offers significant improvements over existing pre-clinical cardiovascular diagnostic techniques. A recent report by Pan, et. al. combines linescans with correlation spectroscopy to measure the velocity of cells travelling through the dorsal aorta of embryonic zebrafish [18]. While determining velocity vectors by correlation spectroscopy is a significant advance, this technique requires two sequential measurements that severely limit its temporal resolution. Discontinuous measurements of cellular velocities have been made in zebrafish embryos using video microscopy, and in transgenic mouse embryos expressing GFP-labeled erythrocytes using perpendicular linescans [15, 19]. Dual spot cross correlation methods of measuring blood cell velocity have existed since the mid-1970s and should be applicable to measuring blood cell velocity in zebrafish embryos with good temporal resolution [2023]. These methods use two photodiodes spaced a fixed distance apart along the path of blood flow, and cell velocities are calculated based on the time required for a cell to pass between the two detectors. Because cross-correlation methods only measure the average speed of a cell as it passes between the two detectors, they do not continuously measure velocity. This continuity is required for calculating acceleration and deceleration, which contains important information about contractility and peripheral resistance. The use of a parallel scan line allows us to continuously measure the velocity of the same object so long as it is located within the linear scanning region. This feature is what makes laser-scanning velocimetry unique in its temporal resolution for both velocity and acceleration measurements.

Video microscopy has also been used to measure cardiac output in Xenopus laevis, Gallus gallus, and Danio rerio embryos [11, 12, 2426]. This technique estimates the volume of blood ejected from the ventricle based on the observed changes in its diameter. Our measurements of cardiac output using laser-scanning velocimetry are in good agreement with these published values [12]. While the cardiac output measurements described here are estimates that assume a circular vessel cross-section, and can be underestimated if linescans are obtained from the aorta after arterial branches, they do not require normal heart morphology and valve function. A more exact determination of cardiac output may be obtained using laser-scanning velocimetry in transgenic zebrafish animals where the vessel dimensions are delineated by an endothelial expressed GFP such as fli1-GFP [27].

Recently, the genes responsible for two mutant zebrafish embryos that have similar cardiovascular phenotypes, santa and valentine, were mapped to ccm1 and ccm2, respectively [28]. The valentine mutant contains a mutation in OSM and exhibits pericardial edema, a dilated heart, and a complete lack of blood flow [28]. Morpholino knockdown of OSM resulted in significantly reduced circulation relative to control zebrafish. The fact that morpholino knockdown of OSM did not completely inhibit blood flow is probably a result of low level OSM expression in contrast to complete loss of functional OSM protein in the valentine mutant.

The observation that OSM plays a role in aortic arch lumen formation has important implications for the pathobiology of cerebral cavernous malformations. Mice homozygous for a gene-trap insertion of β-galactosidase into the ccm1/Krit1 locus also have a constricted aortic arch [14], indicating that a loss of Krit1 in mouse phenocopies a loss of OSM in zebrafish. When viewed in light of evidence that OSM and Krit1 physically associate with each other at sites of actin reorganization [29], these data suggest that Krit1 and OSM may act as a complex to organize developmental signals required for proper aortic arch lumen formation.

Zebrafish are a popular model for drug screens because of their genetic similarity to humans, the speed of screening in a physiological context and cost-effectiveness [30]. In addition to providing a diagnostic tool for molecular geneticists, laser-scanning velocimetry has obvious applications in pre-clinical drug development. For example, laser-scanning velocimetry could easily be used to identify anti-arhythmic compounds capable of reversing the fibrillating phenotype of the tremblor or reggae mutants [31]. Similarly, the potency and efficacy of potential inotropic drugs, which increase the force of ventricular contractions in heart failure, could easily be screened in wild-type zebrafish embryos to find those compounds that best increase the peak acceleration of blood cells. As demonstrated in Figure 2, in addition to supporting pre-clinical efficacy studies, laser-scanning velocimetry can also be used to readily identify molecules that possess cardiotoxicity.

Laser-scanning velocimetry quantifies the effects of genetic and pharmacological modifiers of cardiovascular function. With the increasing popularity of zebrafish in basic science and the pharmaceutical industry, the development of additional techniques to correlate physiological function with genes and their chemical modifiers promises to accelerate drug discovery and development in the post-genomic area.


Zebrafish husbandry

Zebrafish (Danio rerio) embryos used for velocimetry measurements and morpholino injections were obtained by natural spawning of a wild-type AB* line. Embryos were raised at 28.5°C in 30% Danieau solution(17 mM NaCl, 0.21 mM KCl, 0.12 mM MgSO4, 0.18 mM Ca(NO3)2, 1.5 mM HEPES pH 7.6). Zebrafish were handled in accordance with IACUC approved protocols.

Laser-scanning velocimetry

Embryos for imaging were placed in a 4-well coverglass-bottom chamber slide (Lab-Tek) in 30% Danieau solution containing the anaesthetic tricaine at a final concentration of 100 μg/ml or as indicated. After a three-minute delay to ensure a steady-state level of anaesthesia, images were captured using an Olympus Fluoview FV1000 laser-scanning microscope with an Olympus PlanApo N (60 ×/1.42 oil) lens using DIC optics. Scan lines were oriented parallel to the flow of blood within the dorsal aorta, just above the yolk ball. Approximately 15,000 lines were acquired per sample at a rate of 0.488 ms/line allowing a maximum velocity of 4.10 cm/s to be measured. Using a 20 μm line and the optical system described above, a minimum scan rate of 8 ms/line is required to record the maximum cellular velocity of 0.25 cm/s observed in these experiments. Cell paths within parallel linescans were traced using the NeuronJ plugin v.1.1.0 for ImageJ v.1.37 [32]. High-quality DIC images and proper adjustment of NeuronJ parameters were critical to obtaining quality tracings. Good results were obtained with the following NeuronJ settings: Hessian smoothing scale = 1.0, cost weight factor = 0.7, snap window size = 1 × 1, path-search window size = 800 × 800, tracing smoothing range = 10, tracing subsampling factor = 1, and line width = 1. NeuronJ trace coordinates were exported as text files then processed to calculate average velocity and acceleration versus time data with Laser-Scanning Velocimetry Data Analyzer software we developed using Visual Basic v6.0 and Access 2003 (Microsoft). Briefly, the position versus time data exported from NeuronJ was smoothed using a moving average function. The instantaneous slope of this smoothed trace was determined for each time point. Where multiple traces were measured for a single point in time, the calculated velocities were averaged. Acceleration data were obtained from the instantaneous slope of velocity-versus-time profiles. A moving average function was used throughout the calculation workflow to minimize noise. The interval of the moving average can be specified as smoothing factors in the data analyzer software. Higher smoothing factors will underestimate measured values but will increase precision. Data are only comparable when calculated using the same smoothing factors and have been specified in the figure legends.

Tricaine dose-response treatment

A single zebrafish embryo at 30 hpf was placed in a coverglass-bottom chamber slide containing 30% Danieau solution without tricaine. Parallel linescans were acquired from the free-floating, non-anesthetized embryo to obtain baseline values of an untreated embryo. After a baseline image was collected, the 30% Danieau was sequentially replaced with 30% Danieau solutions containing increasing concentrations of tricaine from 300 μM to 1 mM. After each solution change, the embryo was allowed to equilibrate with the drug for three minutes prior to imaging.

Stroke volume, heart rate, and cardiac output calculations

Borrowing principles from human ultrasonic Doppler Imaging, stroke volume and cardiac output were extracted from linescans of the dorsal aorta close to the heart, where there is little influence from peripheral vessels. Stroke volume (SV), the volume of blood ejected by the ventricle in one heartbeat, was determined by multiplying the cross sectional area of the aorta (Aaorta.) by the distance blood is pumped along the aorta in one beat (dx) (Equation 1).

SV = Aaorta × dx

Assuming a circular dorsal aorta, the cross-sectional area is given by Equation 2, where the diameter of the aorta (daorta.) was measured from the x-dimension of the perpendicular linescan.

A a o r t a = π ( d a o r t a 2 ) 2 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGbbqqdaWgaaWcbaGaemyyaeMaem4Ba8MaemOCaiNaemiDaqNaemyyaegabeaakiabg2da9GGaciab=b8aWnaabmaabaWaaSaaaeaacqWGKbazdaWgaaWcbaGaemyyaeMaem4Ba8MaemOCaiNaemiDaqNaemyyaegabeaaaOqaaiabikdaYaaaaiaawIcacaGLPaaadaahaaWcbeqaaiabikdaYaaaaaa@439E@

Because velocity is defined as distance travelled per unit time, integrating cell velocities as a function of time gives the distance travelled (Equation 3).

d x = 0 t v ( t ) d t MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGKbazcqWG4baEcqGH9aqpdaWdXaqaaiabdAha2jabcIcaOiabdsha0jabcMcaPiabdsgaKjabdsha0bWcbaGaeGimaadabaGafmiDaqNbauaaa0Gaey4kIipaaaa@3C6A@

Substitution of Equations 2 and 3 into equation 1 gives the final equation for stroke volume (Equation 4).

S V = π ( d a o r t a 2 ) 2 × 0 t v ( t ) d t MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaacqWGtbWucqWGwbGvcqGH9aqpiiGacqWFapaCdaqadaqaamaalaaabaGaemizaq2aaSbaaSqaaiabdggaHjabd+gaVjabdkhaYjabdsha0jabdggaHbqabaaakeaacqaIYaGmaaaacaGLOaGaayzkaaWaaWbaaSqabeaacqaIYaGmaaGccqGHxdaTdaWdXaqaaiabdAha2jabcIcaOiabdsha0jabcMcaPiabdsgaKjabdsha0bWcbaGaeGimaadabaGafmiDaqNbauaaa0Gaey4kIipaaaa@4BF5@

The integral of v(t)dt was approximated by calculating the area under the curve of a plot of velocity as a function of time for each heart beat using Origin v.6.0 (Microcal). Cardiac output, the volume of blood ejected by the heart in one minute, was calculated by multiplying the stroke volume by the heart rate (Equation 5).

CO = SV × HR

Heart rates were determined by dividing 60 seconds by the average duration of a heart beat, which was measured by the average peak-to-peak distance of the velocity-versus-time profile.

Morpholino knockdown of OSM

Morpholino oligonucleotides targeting the translation initiation site of the Danio rerio ccm2/osm mRNA were purchased from Gene Tools, Inc. Two OSM-specific morpholinos were used (MO#1: ATCTAATACAGCGAAAATGAAGAGC; MO#2: ATTTGTACGTAGAGATGGAGGAGGA). Morpholino sequences were designed from the +75 to -25 bp region surrounding the initiating ATG of Danio rerio ccm2/osm cDNA (IMAGE Clone 7396516) by Gene Tools, Inc. A non-targeting morpholino or 1× Danieau solution (57 mM NaCl, 0.7 mM KCl, 0.4 mM MgSO4, 0.6 mM Ca(NO3)2, 5 mM HEPES pH 7.6) were used as controls. Unless otherwise specified, morpholinos were dissolved to 2 ng/nl in 1× Danieau solution. Approximately 5 nl of morpholino solution were injected into embryos between the 1- and 4-cell stage. Repression of OSM protein was confirmed by immunoblotting using a rabbit polyclonal anti-OSM antibody produced against recombinant full-length mouse OSM [33].

Capped OSM mRNA synthesis and rescue

Capped, in vitro transcribed ccm2/osm mRNA was synthesized using Mmessage Mmachine T7 kit (Ambion) from the coding region of a Danio rerio OSM cDNA (IMAGE Clone 7396516). Embryos between the 1- and 4-cell stages were injected with 5 nl ccm2/osm mRNA (20–100 pg/nl) together with an OSM morpholino that targets the 5' UTR (MO#1) at 2–4 ng/nl).

In silicomicroangiography

In silico microangiography is an image processing approach to representing motion in a series of two-dimensional images. A DIC movie of cells flowing through the vasculature of anesthetized embryos at 28 or 48 hpf was collected as a set of 100 12-bit tiff images recorded at 15 ms intervals (X, Y).t. Using Metamorph IP software (Molecular Devices) a difference image stack was created by dividing adjacent time points; that is, for time t = n, Difference.n = (X, Y).n+1. *1000/(X, Y).n.. The resultant difference stack has intensity values only where changes occurred from one time point to the next. Maximum intensity projections through the differences stack produce an image correlating to blood cell movement over the duration of the original movie.


  1. Beis D, Bartman T, Jin SW, Scott IC, D'Amico LA, Ober EA, Verkade H, Frantsve J, Field HA, Wehman A, Baier H, Tallafuss A, Bally-Cuif L, Chen JN, Stainier DY, Jungblut B: Genetic and cellular analyses of zebrafish atrioventricular cushion and valve development. Development. 2005, 132: 4193-4204. 10.1242/dev.01970.

    Article  CAS  Google Scholar 

  2. Chen JN, Haffter P, Odenthal J, Vogelsang E, Brand M, van Eeden FJ, Furutani-Seiki M, Granato M, Hammerschmidt M, Heisenberg CP, Jiang YJ, Kane DA, Kelsh RN, Mullins MC, Nusslein-Volhard C: Mutations affecting the cardiovascular system and other internal organs in zebrafish. Development. 1996, 123: 293-302.

    CAS  Google Scholar 

  3. Warren KS, Wu JC, Pinet F, Fishman MC: The genetic basis of cardiac function: dissection by zebrafish (Danio rerio) screens. Philos Trans R Soc Lond B Biol Sci. 2000, 355: 939-944. 10.1098/rstb.2000.0629.

    Article  CAS  Google Scholar 

  4. Amsterdam A, Hopkins N: Mutagenesis strategies in zebrafish for identifying genes involved in development and disease. Trends Genet. 2006, 22: 473-478. 10.1016/j.tig.2006.06.011.

    Article  CAS  Google Scholar 

  5. Chen E, Ekker SC: Zebrafish as a genomics research model. Curr Pharm Biotechnol. 2004, 5: 409-413. 10.2174/1389201043376652.

    Article  CAS  Google Scholar 

  6. Liquori CL, Berg MJ, Siegel AM, Huang E, Zawistowski JS, Stoffer T, Verlaan D, Balogun F, Hughes L, Leedom TP, Plummer NW, Cannella M, Maglione V, Squitieri F, Johnson EW, Rouleau GA, Ptacek L, Marchuk DA: Mutations in a gene encoding a novel protein containing a phosphotyrosine-binding domain cause type 2 cerebral cavernous malformations. Am J Hum Genet. 2003, 73: 1459-1464. 10.1086/380314.

    Article  CAS  Google Scholar 

  7. Sabbah HN, Khaja F, Brymer JF, McFarland TM, Albert DE, Snyder JE, Goldstein S, Stein PD: Noninvasive evaluation of left ventricular performance based on peak aortic blood acceleration measured with a continuous-wave Doppler velocity meter. Circulation. 1986, 74: 323-329.

    Article  CAS  Google Scholar 

  8. DuBourg O, Jondeau G, Beauchet A, Hardy A, Bourdarias JP: Doppler-derived aortic maximal acceleration. A reliable index of left ventricular systolic function. Chest. 1993, 103: 1064-1067.

    Article  CAS  Google Scholar 

  9. Saeian K, Wann LS, Sagar KB: Doppler echocardiographic evaluation of left ventricular function. Echocardiography. 1990, 7: 21-25.

    Article  CAS  Google Scholar 

  10. Bennett ED, Barclay SA, Davis AL, Mannering D, Mehta N: Ascending aortic blood velocity and acceleration using Doppler ultrasound in the assessment of left ventricular function. Cardiovasc Res. 1984, 18: 632-638.

    Article  CAS  Google Scholar 

  11. Kopp R, Pelster B, Schwerte T: How does blood cell concentration modulate cardiovascular parameters in developing zebrafish (Danio rerio)?. Comp Biochem Physiol A Mol Integr Physiol. 2006

    Google Scholar 

  12. Bagatto B, Burggren W: A three-dimensional functional assessment of heart and vessel development in the larva of the zebrafish (Danio rerio). Physiol Biochem Zool. 2006, 79: 194-201. 10.1086/498185.

    Article  Google Scholar 

  13. Liebling M, Forouhar AS, Wolleschensky R, Zimmermann B, Ankerhold R, Fraser SE, Gharib M, Dickinson ME: Rapid three-dimensional imaging and analysis of the beating embryonic heart reveals functional changes during development. Dev Dyn. 2006, 235: 2940-2948. 10.1002/dvdy.20926.

    Article  Google Scholar 

  14. Whitehead KJ, Plummer NW, Adams JA, Marchuk DA, Li DY: Ccm1 is required for arterial morphogenesis: implications for the etiology of human cavernous malformations. Development. 2004, 131: 1437-1448. 10.1242/dev.01036.

    Article  CAS  Google Scholar 

  15. Schwerte T, Pelster B: Digital motion analysis as a tool for analysing the shape and performance of the circulatory system in transparent animals. J Exp Biol. 2000, 203: 1659-1669.

    CAS  Google Scholar 

  16. Kittleson MD, Rienle RD: Aortic Stenosis. Small Animal Cardiovascular Medicine. 1998, St. Louis, Mosby, OInc., 260-272.

    Google Scholar 

  17. Ho YL, Shau YW, Tsai HJ, Lin LC, Huang PJ, Hsieh FJ: Assessment of zebrafish cardiac performance using Doppler echocardiography and power angiography. Ultrasound Med Biol. 2002, 28: 1137-1143. 10.1016/S0301-5629(02)00564-1.

    Article  Google Scholar 

  18. Pan X, Yu H, Shi X, Korzh V, Wohland T: Characterization of flow direction in microchannels and zebrafish blood vessels by scanning fluorescence correlation spectroscopy. J Biomed Opt. 2007, 12: 14034-10.1117/1.2435173.

    Article  Google Scholar 

  19. Jones EA, Baron MH, Fraser SE, Dickinson ME: Measuring hemodynamic changes during mammalian development. Am J Physiol Heart Circ Physiol. 2004, 287: H1561-9. 10.1152/ajpheart.00081.2004.

    Article  CAS  Google Scholar 

  20. Osterloh K, Gaehtgens P, Pries AR: Determination of microvascular flow pattern formation in vivo. Am J Physiol Heart Circ Physiol. 2000, 278: H1142-52.

    CAS  Google Scholar 

  21. Zimmerhackl B, Tinsman J, Jamison RL, Robertson CR: Use of digital cross-correlation for on-line determination of single-vessel blood flow in the mammalian kidney. Microvascular Research. 1985, 30: 63-74. 10.1016/0026-2862(85)90038-X.

    Article  CAS  Google Scholar 

  22. Arfors KE, Bergqvist D, Intaglietta M, Westergren B: Measurements of blood flow velocity in the microcirculation. Ups J Med Sci. 1975, 80: 27-33.

    Article  CAS  Google Scholar 

  23. Jamison RL, Zimmerhackl B, Robertson CR: Videomicroscopic method for direct determination of blood flow to the papilla of the kidney. Ann Biomed Eng. 1985, 13: 273-280. 10.1007/BF02584244.

    Article  CAS  Google Scholar 

  24. Fritsche R, Schwerte T, Pelster B: Nitric oxide and vascular reactivity in developing zebrafish, Danio rerio. Am J Physiol Regul Integr Comp Physiol. 2000, 279: R2200-7.

    CAS  Google Scholar 

  25. Hou PC, Burggren WW: Cardiac output and peripheral resistance during larval development in the anuran amphibian Xenopus laevis. Am J Physiol. 1995, 269: R1126-32.

    CAS  Google Scholar 

  26. Burggren WW, Warburton SJ, Slivkoff MD: Interruption of cardiac output does not affect short-term growth and metabolic rate in day 3 and 4 chick embryos. J Exp Biol. 2000, 203: 3831-3838.

    CAS  Google Scholar 

  27. Lawson ND, Weinstein BM: In vivo imaging of embryonic vascular development using transgenic zebrafish. Dev Biol. 2002, 248: 307-318. 10.1006/dbio.2002.0711.

    Article  CAS  Google Scholar 

  28. Mably JD, Chuang LP, Serluca FC, Mohideen MA, Chen JN, Fishman MC: santa and valentine pattern concentric growth of cardiac myocardium in the zebrafish. Development. 2006, 133: 3139-3146. 10.1242/dev.02469.

    Article  CAS  Google Scholar 

  29. Zawistowski JS, Stalheim L, Uhlik MT, Abell AN, Ancrile BB, Johnson GL, Marchuk DA: CCM1 and CCM2 protein interactions in cell signaling: implications for cerebral cavernous malformations pathogenesis. Hum Mol Genet. 2005, 14: 2521-2531. 10.1093/hmg/ddi256.

    Article  CAS  Google Scholar 

  30. Zon LI, Peterson RT: In vivo drug discovery in the zebrafish. Nat Rev Drug Discov. 2005, 4: 35-44. 10.1038/nrd1606.

    Article  CAS  Google Scholar 

  31. Stainier DY, Fouquet B, Chen JN, Warren KS, Weinstein BM, Meiler SE, Mohideen MA, Neuhauss SC, Solnica-Krezel L, Schier AF, Zwartkruis F, Stemple DL, Malicki J, Driever W, Fishman MC: Mutations affecting the formation and function of the cardiovascular system in the zebrafish embryo. Development. 1996, 123: 285-292.

    CAS  Google Scholar 

  32. Meijering E, Jacob M, Sarria JC, Steiner P, Hirling H, Unser M: Design and validation of a tool for neurite tracing and analysis in fluorescence microscopy images. Cytometry A. 2004, 58: 167-176. 10.1002/cyto.a.20022.

    Article  CAS  Google Scholar 

  33. Uhlik MT, Abell AN, Johnson NL, Sun W, Cuevas BD, Lobel-Rice KE, Horne EA, Dell'Acqua ML, Johnson GL: Rac-MEKK3-MKK3 scaffolding for p38 MAPK activation during hyperosmotic shock. Nat Cell Biol. 2003, 5: 1104-1110. 10.1038/ncb1071.

    Article  CAS  Google Scholar 

Download references


The authors are grateful to Sue Donerly for technical assistance with zebrafish husbandry and embryo manipulation. Funding was generously provided by the National Institutes of Health to MHM (HL084971) and GLJ (GM30324 and GM68820).

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Michael H Malone or Gary L Johnson.

Additional information

Authors' contributions

MHM and NS conceived of the study. MHM performed zebrafish husbandry, morpholino knockdown and drug treatment experiments, acquired images, wrote the laser-scanning velocimetry data analyzer software, analyzed linescan images and drafted the manuscript. NS developed image acquisition protocols, processed images for linescans, performed in silico microangiography, and assisted in drafting the manuscript. LS performed zebrafish husbandry and morpholino knockdown experiments. KH assisted in developing image acquisition protocols. EL assisted with zebrafish husbandry and morpholino knockdown experiments. GLJ participated in the design of the study, coordinated essential collaborations, and assisted in drafting the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional File 1: position effects. Dependence of cardiac performance measurements on scan line position (PDF 98 KB)


Additional File 2: cardiovascular analysis of OSM morphant. Decreased circulation in OSM-deficient embryos only affects the aortic arch (PDF 129 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Malone, M.H., Sciaky, N., Stalheim, L. et al. Laser-scanning velocimetry: A confocal microscopy method for quantitative measurement of cardiovascular performance in zebrafish embryos and larvae. BMC Biotechnol 7, 40 (2007).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: