Research Papers

Quantifying Function in the Early Embryonic Heart

[+] Author and Article Information
Brennan M. Johnson

School of Biomedical Engineering,
Colorado State University,
Fort Collins, CO 80523

Deborah M. Garrity

Department of Biology,
Colorado State University,
Fort Collins, CO 80523

Lakshmi Prasad Dasi

School of Biomedical Engineering,
Colorado State University,
Fort Collins, CO 80523
e-mail: lakshmi.dasi@colostate.edu

1Corresponding author. Present address: Department of Mechanical Engineering, Colorado State University, Room A103D Engineering, 1374 Campus Delivery, Fort Collins, CO 80523-1374.

Contributed by the Bioengineering Division of ASME for publication in the Journal of Biomechanical Engineering. Manuscript received September 12, 2012; final manuscript received January 31, 2013; accepted manuscript posted February 19, 2013; published online April 2, 2013. Assoc. Editor: Naomi Chesler.

J Biomech Eng 135(4), 041006 (Apr 02, 2013) (11 pages) Paper No: BIO-12-1402; doi: 10.1115/1.4023701 History: Received September 12, 2012; Revised January 31, 2013; Accepted February 19, 2013

Congenital heart defects arise during the early stages of development, and studies have linked abnormal blood flow and irregular cardiac function to improper cardiac morphogenesis. The embryonic zebrafish offers superb optical access for live imaging of heart development. Here, we build upon previously used techniques to develop a methodology for quantifying cardiac function in the embryonic zebrafish model. Imaging was performed using bright field microscopy at 1500 frames/s at 0.76 μm/pixel. Heart function was manipulated in a wild-type zebrafish at ∼55 h post fertilization (hpf). Blood velocity and luminal diameter were measured at the atrial inlet and atrioventricular junction (AVJ) by analyzing spatiotemporal plots. Control volume analysis was used to estimate the flow rate waveform, retrograde fractions, stroke volume, and cardiac output. The diameter and flow waveforms at the inlet and AVJ are highly repeatable between heart beats. We have developed a methodology for quantifying overall heart function, which can be applied to early stages of zebrafish development.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Bartman, T., Walsh, E. C., Wen, K.-K., McKane, M., Ren, J., Alexander, J., Rubenstein, P. A., and Stainier, D. Y. R., 2004, “Early Myocardial Function Affects Endocardial Cushion Development in Zebrafish,” PLoS Biol., 2(5), pp. 673–681. [CrossRef]
Broekhuizen, M. L. A., Hogers, B., DeRuiter, M. C., Poelmann, R. E., Gittenberger-de Groot, A. C., and Wladimiroff, J. W., 1999, “Altered Hemodynamics in Chick Embryos After Extraembryonic Venous Obstruction,” Ultrasound Obstet. Gynecol., 13(6), pp. 437–445. [CrossRef] [PubMed]
Hogers, B., DeRuiter, M. C., Gittenberger-de Groot, A. C., and Poelmann, R. E., 1999, “Extraembryonic Venous Obstructions Lead to Cardiovascular Malformations and can be Embryolethal,” Cardiovasc. Res., 41(1), pp. 87–99. [CrossRef] [PubMed]
Hogers, B., DeRuiter, M. C., Gittenberger-de Groot, A. C., and Poelmann, R. E., 1997, “Unilateral Vitelline Vein Ligation Alters Intracardiac Blood Flow Patterns and Morphogenesis in the Chick Embryo,” Circ. Res., 80(4), pp. 473–481. [CrossRef] [PubMed]
Hove, J. R., Köster, R. W., Forouhar, A. S., Acevedo-Bolton, G., Fraser, S. E., and Gharib, M., 2003, “Intracardiac Fluid Forces are an Essential Epigenetic Factor for Embryonic Cardiogenesis,” Nature, 421(6919), pp. 172–177. [CrossRef] [PubMed]
Kopp, R., Pelster, B., and Schwerte, T., 2007, “How Does Blood Cell Concentration Modulate Cardiovascular Parameters in Developing Zebrafish (Danio Rerio)?,” Comp. Biochem. Physiol., Part A: Mol. Integr. Physiol., 146(3), pp. 400–407. [CrossRef] [PubMed]
Filas, B. A., Efimov, I. R., and Taber, L. A., 2007, “Optical Coherence Tomography as a Tool for Measuring Morphogenetic Deformation of the Looping Heart,” Anat. Rec., 290(9), pp. 1057–1068. [CrossRef]
Auman, H. J., Coleman, H., Riley, H. E., Olale, F., Tsai, H. J., and Yelon, D., 2007, “Functional Modulation of Cardiac Form Through Regionally Confined Cell Shape Changes,” PLoS Biol., 5(3), pp. 604–615. [CrossRef]
Stainier, D. Y., Lee, R. K., and Fishman, M. C., 1993, “Cardiovascular Development in the Zebrafish. I. Myocardial Fate Map and Heart Tube Formation,” Development, 119(1), pp. 31–40. [PubMed]
Thisse, C., and Zon, L. I., 2002, “Organogenesis—Heart and Wood Formation From the Zebrafish Point of View,” Science, 295(5554), pp. 457–462. [CrossRef] [PubMed]
Bakkers, J., 2011, “Zebrafish as a Model to Study Cardiac Development and Human Cardiac Disease,” Cardiovasc. Res., 91(2), pp. 279–288. [CrossRef] [PubMed]
Chico, T. J. A., Ingham, P. W., and Crossman, D. C., 2008, “Modeling Cardiovascular Disease in the Zebrafish,” Trends Cardiovasc. Med., 18(4), pp. 150–155. [CrossRef] [PubMed]
Dahme, T., Katus, H. A., and Rottbauer, W., 2009, “Fishing for the Genetic Basis of Cardiovascular Disease,” Dis. Model. Mech., 2(1–2), pp. 18–22. [CrossRef] [PubMed]
Stainier, D. Y. R., and Fishman, M. C., 1994, “The Zebrafish as a Model System to Study Cardiovascular Development,” Trends Cardiovasc. Med., 4(5), pp. 207–212. [CrossRef] [PubMed]
Beis, D., Bartman, T., Jin, S. W., Scott, I. C., D'Amico, L. A., Ober, E. A., Verkade, H., Frantsve, J., Field, H. A., Wehman, A., Baier, H., Tallafuss, A., Bally-Cuif, L., Chen, J. N., Stainier, D. Y., and Jungblut, B., 2005, “Genetic and Cellular Analyses of Zebrafish Atrioventricular Cushion and Valve Development,” Development, 132(18), pp. 4193–4204. [CrossRef] [PubMed]
Beis, D., and Stainier, D. Y. R., 2006, “In Vivo Cell Biology: Following the Zebrafish Trend,” Trends Cell Biol., 16(2), pp. 105–112. [CrossRef] [PubMed]
Harvey, R. P., 2002, “Patterning the Vertebrate Heart,” Nat. Rev. Genet., 3(7), pp. 544–556. [CrossRef] [PubMed]
Lieschke, G. J., and Currie, P. D., 2007, “Animal Models of Human Disease: Zebrafish Swim Into View,” Nat. Rev. Genet., 8(5), pp. 353–367. [CrossRef] [PubMed]
Roman, B. L., Pham, V. N., Lawson, N. D., Kulik, M., Childs, S., Lekven, A. C., Garrity, D. M., Moon, R. T., Fishman, M. C., Lechleider, R. J., and Weinstein, B. M., 2002, “Disruption of Acvrl1 Increases Endothelial Cell Number in Zebrafish Cranial Vessels,” Development, 129(12), pp. 3009–3019. [PubMed]
Butcher, J. T., McQuinn, T. C., Sedmera, D., Turner, D., and Markwald, R. R., 2007, “Transitions in Early Embryonic Atrioventricular Valvular Function Correspond With Changes in Cushion Biomechanics That are Predictable by Tissue Composition,” Circ. Res., 100(10), pp. 1503–1511. [CrossRef] [PubMed]
Liu, A. P., Nickerson, A., Troyer, A., Yin, X., Cary, R., Thornburg, K., Wang, R., and Rugonyi, S., 2011, “Quantifying Blood Flow and Wall Shear Stresses in the Outflow Tract of Chick Embryonic Hearts,” Comput. Struct., 89(11–12), pp. 855–867. [CrossRef] [PubMed]
Liu, A.P., Yin, X., Shi, L., Li, P., Thornburg, K. L., Wang, R., and Rugonyi, S., 2012, “Biomechanics of the Chick Embryonic Heart Outflow Tract at HH18 Using 4D Optical Coherence Tomography Imaging and Computational Modeling,” PLoS One, 7(7), e040869. [CrossRef]
Johnson, B. M., Garrity, D. M., and Dasi, L. P., 2013, “The Transitional Cardiac Pumping Mechanics in the Embryonic Heart,” J. Cardiovasc. Eng. Technol., (in press). [CrossRef]
Forouhar, A. S., Liebling, M., Hickerson, A., Nasiraei-Moghaddam, A., Tsai, H.-J., Hove, J. R., Fraser, S. E., Dickinson, M. E., and Gharib, M., 2006, “The Embryonic Vertebrate Heart Tube is a Dynamic Suction Pump,” Science, 312(5774), pp. 751–753. [CrossRef] [PubMed]
Bagatto, B., and Burggren, W., 2006, “A Three-Dimensional Functional Assessment of Heart and Vessel Development in the Larva of the Zebrafish (Danio Rerio),” Physiol. Biochem. Zool., 79(1), pp. 194–201. [CrossRef] [PubMed]
Westerfield, M., 1995, The Zebrafish Book, University of Oregon Press, Eugene, OR.
Schwerte, T., and Pelster, B., 2000, “Digital Motion Analysis as a Tool for Analysing the Shape and Performance of the Circulatory System in Transparent Animals,” J. Exp. Biol., 203(11), pp. 1659–1669. [PubMed]
Tam, J., and Roorda, A., 2010, “Enhanced Detection of Cell Paths in Spatiotemporal Plots for Noninvasive Microscopy of the Human Retina,” Proceedings of the 7th IEEE International Symposium on Biomedical Imaging: From Nano to Macro, pp. 584–587.
Drew, P. J., Blinder, P., Cauwenberghs, G., Shih, A. Y., and Kleinfeld, D., 2010, “Rapid Determination of Particle Velocity From Space-Time Images Using the Radon Transform,” J. Comput. Neurosci., 29(1–2), pp. 5–11. [CrossRef] [PubMed]
Duncan, D. D., Lemaillet, P., Ibrahim, M., Nguyen, Q. D., Hiller, M., and Ramella-Roman, J., 2010, “Absolute Blood Velocity Measured With a Modified Fundus Camera,” J. Biomed. Opt., 15(5), p. 056014. [CrossRef] [PubMed]
Japee, S. A., Pittman, R. N., and Ellis, C. G., 2005, “A New Video Image Analysis System to Study Red Blood Cell Dynamics and Oxygenation in Capillary Networks,” Microcirculation (N.Y.), 12(6), pp. 489–506. [CrossRef]
Japee, S. A., Pittman, R. N., and Ellis, C. G., 2005, “Automated Method for Tracking Individual Red Blood Cells Within Capillaries to Compute Velocity and Oxygen Saturation,” Microcirculation (N.Y.), 12(6), pp. 507–515. [CrossRef]
Malone, M. H., Sciaky, N., Stalheim, L., Hahn, K. M., Linney, E., and Johnson, G. L., 2007, “Laser-Scanning Velocimetry: A Confocal Microscopy Method for Quantitative Measurement of Cardiovascular Performance in Zebrafish Embryos and Larvae,” BMC Biotechnol., 7(40). [CrossRef]
Pan, X. T., Yu, H., Shi, X., Korzh, V., and Wohland, T., 2007, “Characterization of Flow Direction in Microchannels and Zebrafish Blood Vessels by Scanning Fluorescence Correlation Spectroscopy,” J. Biomed. Opt., 12(1), 014034. [CrossRef] [PubMed]
Tam, J., Martin, J. A., and Roorda, A., 2010, “Noninvasive Visualization and Analysis of Parafoveal Capillaries in Humans,” Invest. Ophthalmol. Visual Sci., 51(3), pp. 1691–1698. [CrossRef]
Tam, J., Tiruveedhula, P., and Roorda, A., 2011, “Characterization of Single-File Flow Through Human Retinal Parafoveal Capillaries Using an Adaptive Optics Scanning Laser Ophthalmoscope,” Biomed. Opt. Express, 2(4), pp. 781–793. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 3

Image Processing. A comparison is made between the unprocessed data (top row) and the processed data (bottom row). For each frame, a sliding average of the intensity of the previous frames is subtracted from the current frame. This removes static portions of the image and creates a black background. Moving portions of the image are easily seen as gray or white pixels. The upper ST plot segment (top middle) reveals notable background artifacts, which result in multiple inaccurate variance peaks (top right). Background artifacts are greatly reduced in the processed images (bottom left and middle), resulting in significantly more accurate variance peaks (bottom right) and, consequently, more accurate velocity estimation.

Grahic Jump Location
Fig. 2

Zebrafish embryos and imaging equipment. (a) Zebrafish embryos next to a United States penny for scale. (b) Ventral view of embryonic zebrafish heart. Blood cells can easily be seen through the transparent tissues surrounding the heart. As blood flows through the heart, it proceeds through the inlet (i), atrium (a), atrioventricular junction (AVJ), and finally through the ventricle (v). Scale bar indicates 100 μm. (c) Experimental setup, consisting of bright field stereomicroscope, high speed camera, and computer.

Grahic Jump Location
Fig. 1

Flow chart of methodology

Grahic Jump Location
Fig. 5

ST plot analysis with Radon transforms. Part (a) shows a cartoon of an embryonic heart. ST plot reference lines are placed as shown at the atrial inlet and AVJ. A sample section of a typical ST plot is shown in (b). The ST plot must be analyzed in smaller sections, called bins. One such bin is outlined in (b). Part (c) shows the Radon transform analysis at 0 deg, 90 deg, and 138 deg. The projection of this bin is summed along each of these projection angles (θ) to produce a corresponding plot. When the projection angle matches the streak angle (138 degree sign (°) in this case), the resulting plot contains prominent peaks. A complete Radon transform contains an analysis of the ST plot bin from every angle, where Radon transform values are represented by intensity, as seen in (d). Part (e) shows a plot of the Radon transform variance for each angle. The projection angle with the highest variance is perpendicular to the streak angle in the ST plot. The velocity (pixels/frame) is then tan(θ–90 deg).

Grahic Jump Location
Fig. 4

Creation of spatiotemporal plots. Part (a) shows a series of three frames with mock cells moving in the vertical direction. The vertical gray column marks the pixels that lie along an arbitrary reference line. The pixels along this line are plotted side by side to create the spatiotemporal plot in (b). As cells move along the reference line, their changing position is seen as an angled streak in the ST plot. The vertical axis of the ST plot represents position along the reference line, while the horizontal axis represents time. The velocity of the cells along the reference line can be determined by calculating the slope (dy/dt) of the streaks in the ST plot.

Grahic Jump Location
Fig. 6

Measurement of lumen diameter. ST plots were created from reference lines which spanned the width of the lumen at the atrial inlet (a) and AVJ. The lines in (b) mark the flow boundary used for measuring orifice diameter, ‘D’. Measurements were used to plot diameter versus time. Time was normalized such that one complete cardiac cycle spans from 0 to 1.

Grahic Jump Location
Fig. 7

Cardiac physiological parameters. Plots were created for orifice diameter (a), blood velocity (b), flow rate (c), and accumulated blood volume (d) versus time at each orifice. Results are shown for C2. Time was normalized such that one complete cardiac cycle spans from 0 to 1.

Grahic Jump Location
Fig. 8

Changes in flow rate and volume at the atrial inlet and AVJ. Plots of the estimated flow rate versus time (a) and (c) and accumulated orifice volume versus time (b) and (d) for each experimental case. Time was normalized such that one complete cardiac cycle spans from 0 to 1.

Grahic Jump Location
Fig. 9

Bin size optimization. Bin size was optimized based on an ST plot segment with changing velocity. The local ST plot streak angle was manually measured as a baseline. Bin size was varied and results were compared to the baseline. Results are shown (from left to right) of the ST plot bin, the resulting Radon transform, and the variance of the derivative of the intensity of the Radon transform.

Grahic Jump Location
Fig. 10

Span size optimization. Span size was optimized at low‐moderate velocities (top) and peak velocities (bottom). Span size was varied and results were shown of (from left to right) the processed image, the resulting ST plot, the resulting Radon transform, and the variance of the derivative of the intensity of the Radon transform.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In