- Title
-
Effects of Extended Pharmacological Disruption of Zebrafish Embryonic Heart Biomechanical Environment on Cardiac Function, Morphology and Gene Expression
- Authors
- Foo, Y.Y., Motakis, E., Tiang, Z., Shen, S., Lai, J.K.H., Chan, W.X., Wiputra, H., Chen, N., Chen, C.K., Winkler, C., Foo, R.S.Y., Yap, C.H.
- Source
- Full text @ Dev. Dyn.
|
|
|
Controls (n = 5) and treated hearts before (n = 5) and after (n = 5) washout were quantified for (A–C) ventricular outlet velocity, (D) heart rate, (E) ejected stroke volume, and (F) ejected cardiac output. For (D–F), mean and standard deviations were presented and a two-sample t-test was used, where *P < 0.05. Note that the heart rate plot (D) differs from Figure 3B as the chosen embryos were specifically assessed for other parameters and the sample size was limited to five |
(A) LSFMM image and spatial WSS contour of 5 dpf treated heart I (pre-washout) during systole and diastole, respectively. A, atrium; V, ventricle. Spatial–temporal WSS assessment of controls at 5 dpf (hearts I, II, and III), and treated hearts at 5 dpf before and after BDM washout. Treated hearts I and II were collected before BDM washout and treated hearts III and IV were analyzed after BDM removal. Area-averaged WSS experienced at the inlet, mid-ventricular, and outlet were measured over a course of cardiac cycle for controls and treated hearts (B–G), respectively. (H) Time-averaged WSS (TAWSS) of controls and treated hearts. (I) Oscillatory shear index (OSI) of treated hearts I, II, III, and IV at inlet and outlet regions. Note that OSI of controls was nearly zero based on the calculated oscillatory velocity index |
(A) Left: PCA plot of biological replicates sequenced by GIS and NovogeneALT and two clusters of different biological conditions (green markers were controls and red markers were treated samples) was distinguishable. Right: PCA plot of the expressions after removal of batch effects accounting for differences of the two service providers. (B) Number of detected genes in each sample from five independent experiments (N = 5) where biological replicates 1–3 were sequenced by GIS and biological replicates 4–5 were sequenced by NovogeneALT. (C) Number of significantly up-regulated (FDR ≤0.05 and log FC ≤−1) and down-regulated genes (FDR ≤0.05 and log FC ≥1) |
Real-time RT-qPCR validation of a subset of differentially expressed genes between control and BDM-treated hearts. Expression of genes of interest were first normalized to a housekeeping gene (actb) and then plotted as relative to expression levels in control. Data were presented as mean ± SE. Sample size was three for (A–I) but limited to two only for (J, K) as the expression level in some treated samples for col15a1b and erbb4a were too low to be detected. For (A–I), statistical test was one-way ANOVA test, where *P < 0.05 |
Correlation of log2 FC of nine genes (Figure 7A–I) between RNA-seq and quantitative real-time PCR (RT-qPCR). A strong, statistically significant Pearson correlation (R = 0.9269) was shown between the expression levels detected using RNA-seq and real-time PCR. Each data point represented the mean value of three biological replicates |
Confocal images of GFP+ myocardium. Extensive trabeculae formation (white arrowheads) was evident in control hearts (A). Inner linings of the BDM (B) and blebbistatin-treated hearts (C) remained smooth, which is a common sign of an under-developed myocardium. Scale bars: 50 μm |
Apoptosis detection in BDM- and blebbistatin-treated zebrafish embryonic hearts using TUNEL assay. Apoptotic cells (red) were detected in the endothelium and myocardium (G, K) of treated hearts. Scale bars: 50 μm |
|
Schematic diagram showing several possible mechano-responsive pathways of a malformed heart after a 3-day BDM treatment |