The goal of this protocol is to detail a novel method for the assessment of segmental cardiac function in embryonic zebrafish under both physiological and pathological conditions.
Zebrafish are increasingly utilized as a model organism for cardiomyopathies and regeneration. Current methods evaluating cardiac function fail to reliably detect segmental mechanics and are not readily feasible in zebrafish. Here we present a semiautomated, open-source method for the quantitative assessment of four-dimensional (4D) segmental cardiac function: displacement analysis of myocardial mechanical deformation (DIAMOND). Transgenic embryonic zebrafish were imaged in vivo using a light-sheet fluorescence microscopy system with 4D cardiac motion synchronization. Acquired 3D digital hearts were reconstructed at end-systole and end-diastole, and the ventricle was manually segmented into binary datasets. Then, the heart was reoriented and isotropically resampled along the true short axis, and the ventricle was evenly divided into eight portions (I–VIII) along the short axis. Due to the different resampling planes and matrices at end-systole and end-diastole, a transformation matrix was applied for image registration to restore the original spatial relationship between the resampled systolic and diastolic image matrices. After image registration, the displacement vector of each segment from end-systole to end-diastole was calculated based on the displacement of mass centroids in three dimensions (3D). DIAMOND shows that basal myocardial segments adjacent to the atrioventricular canal undergo the highest mechanical deformation and are the most susceptible to doxorubicin-induced cardiac injury. Overall, DIAMOND provides novel insights into segmental cardiac mechanics in zebrafish embryos beyond traditional ejection fraction (EF) under both physiological and pathological conditions.
Chemotherapy-induced cardiac toxicity and ensuing heart failure are one of the main reasons for chemotherapy discontinuation1. Therefore, cardiac functional assessment plays a crucial role in the identification of cardiac toxicity and, more importantly, in the prediction of early cardiac injury following chemotherapy2. However, current approaches for cardiac functional assessment encounter limitations. Methods such as left ventricular ejection fraction (LVEF) provide only global and often delayed cardiac mechanics after injury3,4. Tissue Doppler imaging provides segmental myocardial deformation information but suffers from significant intraobserver and interobserver variability, in part due to ultrasound beam angle dependency5. Two-dimensional (2D) speckle tracking utilizes the B-mode of echocardiography, which theoretically eliminates the angle dependency, but its accuracy is limited by out-of-plane motion6. Therefore, a rigorous approach for quantifying segmental cardiac function is lacking in both research and clinical settings.
In this context, we developed a 4D quantification method for the analysis of segmental cardiac function that we named displacement analysis of myocardial mechanical deformation (DIAMOND), to determine the displacement vectors of myocardial mass centroids in 3D space. We applied DIAMOND for the in vivo assessment of cardiac function and doxorubicin-induced cardiac toxicity with zebrafish (Danio rerio) as the animal model, chosen due to their regenerating myocardium and highly conserved developmental genes7. We further compared segmental DIAMOND displacement with global ejection fraction (EF) determination and 2D strain following doxorubicin treatment. By integrating DIAMOND displacement with 4D light-sheet fluorescent microscopy (LSFM) acquired rendering of embryonic zebrafish hearts, DIAMOND shows that the basal myocardial segments adjacent to the atrioventricular canal undergo the highest mechanical deformation and are the most susceptible to acute doxorubicin cardiac injury8.
All methods described here have been approved by the UCLA Institutional Animal Care and Use Committee (IACUC), and experiments were performed in compliance with protocols approved by the UCLA Office of Animal Research.
1. Breeding Tg(cmlc2:mCherry) zebrafish and collection of embryos
2. Doxorubicin treatment to induce cardiac injury
3. Notch pathway modulation
4. LSFM imaging and post-imaging synchronization
5. Reconstruction of the 3D systolic and diastolic heart
6. Segmentation of the ventricle
7. Creation of rectangular parallelepipeds for image registration
8. Resample systolic and diastolic 3D hearts along the short axis plane
9. Division of the resampled heart
10. Registration of systolic and diastolic image matrices
11. Output of the displacement vectors
The process by which DIAMOND was developed to assess 3D segmental cardiac function is presented in Figure 1. Following LSFM image acquisition and reconstruction in 3D of the embryonic zebrafish heart (Figure 1A), the true short axis plane was determined as the plane perpendicular to the vertical and horizontal long axes, both of which are determined in a multiplane viewer (Figure 1B). The heart was then resampled along the short axis plane (Figure 1C), and divided into eight equal segments constituted by even angles according to a virtual division line (red dotted line) connecting the center of the endocardial ventricular cavity to the center of the atrioventricular canal (Figure 1E). The 3D depictions of the identified segments are illustrated in a cross-sectional view (Figure 1F) and in comparison with the raw data (Figure 2). Segments VII and VIII were removed from the analysis because they encompass the atrioventricular canal and thus contain less myocardium compared with other segments. The different resampling planes for end-systole (HS) and end-diastole (HD) lead to distinct coordinate systems for end-systolic and end-diastolic matrixes, which need to be registered to restore their original spatial relationship (Figure 1G). The coordinate system of the end-systolic matrix was chosen as the reference for consistency. To determine the transformation matrix (Tm) from the end-diastolic matrix to the end-systolic matrix, a matrix of three parallelepipeds, which is asymmetrical in 3D and has the same dimension as the original image matrix, was virtually created. The parallelepipeds were resampled twice, first in the short axis plane of the end-systole matrix, and then in the short axis plane of the end-diastole matrix, leading to different transformed parallelepipeds for end-systole (green) and end-diastole (red) (Figure 1H).
The green and red parallelepipeds were then registered together by a rigid body registration algorithm and Tm was calculated and applied to the end-diastole matrix to restore the coordinates (Figure 1I). This process permits subsequent tracking in the 3D space of the displacement vectors of mass centroids from any segment of the ventricle during the cardiac cycle (Figure 1J). DIAMOND displacement of ventricular segments I–VI can be tracked during multiple time points in the cardiac cycle (Figure 1K), which can be simplified for quantitative analysis to two time points ranging from end-systole to end-diastole (Figure 1L). The segments generated by DIAMOND can be visualized in Figure 2, where each color represent one cardiac segment.
With DIAMOND, we uncovered segmental heterogeneity of cardiac function and susceptibility to doxorubicin-induced myocardial injury in zebrafish. Following a 24 h treatment with 10 µM doxorubicin from 3–4 dpf (Figure 3A), we compared DIAMOND displacement of ventricular segments between control and chemotherapy-treated groups (Figure 3B) and 48 h after treatment (Figure 3C). All DIAMOND figures follow the same graphical pattern as the resampled ventricles along the short axis (Figure 1E). The data are presented as percentages by normalizing the L2-norm of the displacement vector to the inner perimeter of the heart, with the X (green), Y (blue), and Z components (orange) illustrated as their weighted contributions. At 4 dpf, the average L2-norm of the segmental displacement vectors in control fish ranged from 6.6–11.3 µm, or 3.8–6.6% after normalization. Our results indicate that under control conditions, the basal segments I and VI undergo the largest displacements and are the ones most susceptible to doxorubicin induced cardiac injury (Figure 3B, 29% decrease from 6.6–4.7%, n = 10 control and n = 8 doxorubicin, p < 0.01). At 6 dpf, the average L2-norm of the segmental displacement vectors in the control fish ranged from 6.8–14 µm, or 3.9–8% after normalization. At 6 dpf, the basal segments I and VI recovered DIAMOND displacement to control levels, suggesting segmental regeneration (Figure 3C, n = 10 control and n = 8 doxorubicin). In parallel, a worsening in 2D basal strain from -53 to -38% was observed at 4 dpf following doxorubicin treatment, followed by a return to control levels at 6 dpf, corroborating the DIAMOND displacement results (Figure 3D, 3E). A parallel decrease in global ejection fraction in response to doxorubicin at 4 dpf with recovery at 6 dpf was also observed (Figure 3F, 3G).
We next applied DIAMOND during doxorubicin treatment and Notch pathway modulation using the Notch inhibitor DAPT and rescue using Notch downstream effectors NICD and NRG1 mRNA (Figure 4A). NICD and NRG1 mRNA microinjection rescued the decrease in DIAMOND displacement and EF after acute chemotherapy-induced injury at 4 dpf (Figure 4B, 4D). Exposure to the Notch inhibitor DAPT together with doxorubicin led to a more diffuse decrease in DIAMOND displacement, in addition to the basal segments I and VI (Figure 4B). Moreover, inhibition of the Notch pathway after chemo-induced injury further hindered the recovery of DIAMOND displacement of the basal segments and EF at 6 dpf. The inhibition was rescued by the Notch downstream effectors NICD and NRG1 (Figure 4C, 4E).
Figure 1: 4D DIAMOND displacement development. (A) Raw images were captured by light-sheet fluorescent microscopy. (B and C) Reconstructed 3D heart was resampled along the true short axis plane view. (D) Schematic illustration of embryonic zebrafish heart. (E and F) 2D and 3D illustrations of the division of the ventricle into eight segments excluding segments VII and VIII. (G) The different coordinate systems of end-systole and end-diastole after resampling. (H) A group of rectangular parallelepipeds was created for the generation of a transformation matrix (Tm). (I) Registered end-systolic and end-diastolic coordinate systems by applying Tm. (J) Displacement vector of the segmental mass centroid from end-systole to end-diastole. (K) DIAMOND displacement of ventricular segments I–VI tracked during multiple time points in the cardiac cycle. (L) DIAMOND displacement of ventricular segments I–VI from end-systole to end-diastole. This figure from Chen et al.8 is reproduced with permission from the American Society for Clinical Investigation (ASCI). Please click here to view a larger version of this figure.
Figure 2: DIAMOND segmentation of the embryonic zebrafish heart compared with raw data in 3D. The embryonic zebrafish heart was divided into six segments (volumes) depicted here in different colors for the calculation of DIAMOND displacements (left). The displacement vector of each segment computed by DIAMOND represents its segmental cardiac function. The atrium and outflow tract were removed during segmentation. Scale bar = 50 µm. Please click here to view a larger version of this figure.
Figure 3: DIAMOND unravels the segmental heterogeneity in cardiac function and susceptibility to chemotherapy-induced injury. (A) Experimental schedule of doxorubicin treatment. (B and C) Segmental comparison of DIAMOND displacement vectors normalized to the inner myocardial perimeter between control and doxorubicin-treated groups at 4 and 6 dpf (t tests, **p < 0.01, n = 8–10 per group). (D and E) Assessment of strain in the ventricular base depicting a similar injury and regeneration pattern as the DIAMOND displacement vectors (*p < 0.05, n = 6–8 per group). (F and G) Decrease in ejection fraction in response to doxorubicin at 4 dpf with recovery at 6 dpf, following a pattern similar to segmental DIAMOND displacements at the global ventricular level (t tests, **p < 0.01, error bars SEM, n = 6–10 per group). This figure from Chen et al.8 is reproduced with permission from the ASCI. Please click here to view a larger version of this figure.
Figure 4: DIAMOND mechanics for assessment of Notch-mediated myocardial recovery following doxorubicin-induced injury. (A) Experimental schedule. (B and C) NICD and NRG1 Notch downstream effectors rescued the reduction of DIAMOND displacement in segments I and VI at 4 dpf. At 6 dpf, inhibition of Notch signaling by DAPT impaired the restoration of segmental cardiac function (ANOVA, **p < 0.01 Dox vs. control; †p < 0.05, ††p < 0.01, Dox + DAPT vs. control, n = 6–10 per group). (D and E) Ejection fraction corroborates DIAMOND mechanics at the global level (ANOVA, *p < 0.05, **p < 0.01, error bars SEM, n = 5–11 per group). This figure from Chen et al.8 is reproduced with permission from the ASCI. Please click here to view a larger version of this figure.
A rigorous strategy for quantification of segmental myocardial function is critical to assess cardiac mechanics beyond traditional EF, known to be an insensitive and delayed indicator of myocardial injury1,4,12. Hence, there has been a growing interest in markers of early myocardial changes, and a growing body of literature supports myocardial deformation parameters as an early indicator to forecast ventricular dysfunction4,13. Echocardiographic measurement of left ventricular (LV) strain provides an established method of myocardial deformation measurement13. However, tissue Doppler-based strain imaging suffers from a number of shortcomings due to angle dependency and intraobserver and interobserver variability14. Speckle tracking echocardiography (STE) can resolve angle-independent 2D and 3D tissue deformation, but the accuracy of 2D speckle tracking is affected by through-plane motion6, while 3D speckle tracking requires superior spatial resolution to resolve the positive ultrasound interference patterns (speckles) in 3D and high temporal resolution to track the speckles between frames15. In the present protocol, we describe DIAMOND displacement as a novel myocardial deformation parameter for in vivo quantification of 4D segmental cardiac function in zebrafish. Compared with EF and 2D strain as reference standards, DIAMOND provides additional segmental deformation information without being affected by through-plane motion. By integrating DIAMOND with 4D LSFM, our technique can assess the displacement vector of a heart segment 20–30 µm in width, which is currently impossible for even the most advanced 3D STE system, which has millimeter-range resolution16.
To apply DIAMOND, it is critical to have a comprehensive understanding of the anatomical structure of the embryonic zebrafish heart. During image segmentation, it is essential that the atrioventricular canal and the outflow tract are correctly identified and segmented out from the rest of the myocardium when the user is performing step 6 in the protocol. Furthermore, the horizontal and the vertical long axes of the ventricle must be accurately determined in order to derive the true short axis plane for image resampling in step 8.
The major rate limiting factor of applying DIAMOND is the manual segmentation of the ventricle, which becomes time-consuming when multiple phases during the cardiac cycle may need to be evaluated. With the advancement of machine learning and neural networks, an automated cardiac segmentation method17,18,19,20 could be integrated with DIAMOND to provide monitoring of segmental cardiac function throughout the entire cardiac cycle. Further applications of DIAMOND also include the integration with echocardiography, micro-CT, or micro-MRI, suitable in larger animal models for the multiscale assessment of cardiac injury and regeneration21. However, the method will first require adaptation to the presence of myocardial fibers leading to more complex cardiac deformation including torsion in mammals22,23.
Overall, DIAMOND provides a novel method to evaluate segmental cardiac function in embryonic zebrafish under both physiological and pathological conditions and may be used as a platform for high-throughput in vivo screening of pathways associated with chemotherapy-induced cardiac toxicity.
The authors have nothing to disclose.
The present work was funded by American Heart Association grants 16SDG30910007 and 18CDA34110338, and by National Institutes of Health grants HL083015, HL111437, HL118650, and HL129727.
Amira6 | FEI | Image analyzing software | |
DAPT | Millipore Sigma | D5942-5MG | |
Doxorubicin hydrochloride | Millipore Sigma | D1515-10MG | |
Ethyl 3-aminobenzoate methanesulfonate | Millipore Sigma | E10521-10G | Tricaine |
MATLAB | MathWorks | Programming environment | |
MATLAB Image Processing Toolbox | MathWorks | Image processing toolbox |