Spatially resolved RNA-sequencing of the embryonic heart identifies a role for Wnt/β-catenin signaling in autonomic control of heart rate
- Hubrecht Institute-KNAW, Netherlands
- University Medical Center Utrecht, Netherlands
Figures
Figure 1 with 1 supplement
Tomo-seq on isolated embryonic hearts reveals distinct clusters of gene expression.
(A) Hearts were isolated from 2-day-old embryos and 10 μm sections were made along the anterior-posterior axis, from outflow to inflow pole. Each section was collected in an individual tube followed by RNA isolation and cDNA transcription using section specific barcodes. After that, samples were pooled for linear amplification and sequence library preparation. (B) Pairwise correlation between individual sections across all genes detected at more than 20 reads in at least two sections. + 1 equals total positive correlation, 0 no correlation, −1 total negative correlation. Blocks of correlating sections can be observed. (C) Hierarchical cluster analysis of gene expression per section. Distinct gene expression clusters correspond to different regions of the heart shown in boxed areas: 1. arterial pole/OFT; 2. Ventricle; 3. AV canal; 4. Atrium. RNA sequencing reads per gene were normalized against the total read count per section as well as the total spike-in control RNA read counts. Genes expressed in min three sections z-value >1.25. See also Figure 1—figure supplement 1 and Figure 1—source data 1.
-
Figure 1—source data 1
Whole dataset table.
Spreadsheet containing the whole tomo-seq dataset of dissected heart of a 2-day-old zebrafish embryo. Sections #2-#41. Data normalized against total reads and total spike-in read count.
- https://doi.org/10.7554/eLife.31515.004
Figure 1—figure supplement 1
Tomo-seq statistics for the analyzed two dpf dataset.
(A) Total raw RNA-seq reads of the 2 dpf heart sample. Plotted are the total reads per 10 μm section. Sections #2 - #41 contained the tissue sample. (B) Total spike-in reads plotted per section. Spike-in expression per section was used for data normalization. (C) Basic statistics of RNA-seq results.
Figure 2
Transcriptome map of the embryonic heart with high spatial resolution.
(A–E) Tomo-seq expression traces and corresponding in situ hybridization for (A) myl7 (whole myocardium), (C) vmhc (ventricular myocardium), (D) myh6 (atrial myocardium), (E) has2 (endocardial cushions), (F) bmp4 (AVC myocardium, orange arrow; IFT myocardium, purple arrow) and (G) tbx2b (OFT myocardium, yellow arrow, AVC myocardium, orange arrow; IFT myocardium, purple arrow). Smoothening (LOESS) was applied to graphs A-E, span α = 0.3. Anterior up. Gray dashed line outlines the heart. White dashed line outlines the eyes. A, atrium; V, ventricle; AVC: atrioventricular canal; IFT, inflow tract; OFT, outflow tract; SA, sinoatrial region. Scale bars represent 50 μm.
Figure 3
Differential gene expression reveals distinct sub-compartments within the heart.
Superimposed expression profiles of has2, nppa, bmp4 and tbx2b. Expression peaks of the individual genes overlap with the position of the clusters (dashed lines). has2 and tbx2b in the OFT (yellow arrow). nppa in the ventricle. has2, bmp4 and tbx2b in the AVC (orange arrow). nppa in the atrium. bmp4 and tbx2b in the SA region (purple arrow). Smoothening (LOESS) was applied to has2 and nppa graphs, span α = 0.3. OFT, outflow tract; AVC: atrioventricular canal; SA, sinoatrial region.
Figure 4 with 1 supplement
Differentially expressed genes in cardiac sub-compartments.
(A–C) Volcano plots highlighting genes differentially expressed in the (A) atrium (n = 99) and ventricle (n = 196), (B) AV canal (n = 346) and (C) sinoatrial region (n = 502). Grey and red dots represent individual genes. Green lines indicate threshold of log2 fold change >2. Blue lines indicate threshold of p-value<0,05. (D–F) Expression traces and in situ hybridization analysis for representative example genes (gene names underlined in volcano plots) significantly upregulated in (D) atrium (smtnl1) and ventricle (mb), (E) AV canal (anxa5b and pdgfaa) and in the sinoatrial region (shox2 and fzd9b). Smoothening (LOESS) was applied to graphs D and E, span α = 0.2. Gray dashed line outlines the heart. White dashed line outlines the eyes. Anterior up. Scale bars represent 50 μm. See also Figure 4—figure supplement 1 and Figure 4—source data 1.
-
Figure 4—source data 1
Genes significantly upregulated in cardiac sub-compartments.
Ventricle (#12-#22), Atrium (#27- #37), AV canal (#19-#26) and sinoatrial region (#33-#39). Thresholds: logFC >2; p-value>0.05.
- https://doi.org/10.7554/eLife.31515.009
Figure 4—figure supplement 1
Identification of region-specific gene expression patterns.
(A–B) Additional gene expression plot and ISH analysis for representative example genes significantly upregulated in the ventricle (hspb11, rcan1a, csrp3) and atrium (apcdd1l, lmod2a). Gray dashed line: location of the heart. White dashed line outlines the eyes. Anterior to the top. (C–E) GO term enrichment for genes significantly enriched in the ventricle (C), atrium (D) and pacemaker domain (E). Smoothening (LOESS) was applied to graphs A and B, span α = 0.2. p-value<0.05. Scale bars represent 50 μm.
Figure 5 with 1 supplement
Wnt/β-catenin signaling in cardiac pacemaker cells.
3D reconstructions of confocal scans from whole mount embryos containing the tg(7xTCFXla.Siam:nlsmCherry) (reporting Wnt/β-catenin activity in red) and the tg(Isl1:GFF;UAS:GFP) (reporting isl1 expression and marking pacemaker cells in green) and stained for tropomyosin (myocardium in white). Anterior is up and posterior is down. Wild-type sibling (A-D’) and isl1K88X mutant (E-H’) embryos at 3 dpf. The blue dashed line boxes in A-D and E-H indicate the sinoatrial region shown enlarged in A’-D’ and E’-H’, respectively. In wild-type hearts, Isl+ pacemaker cells (blue arrowheads) co-expressed TCF-mCherry. In isl1-/- mutants, no TCF-mCherry expression was observed in the Isl1+ pacemaker cells (blue arrowheads). TCF-mCherry expression was also detected in the AV canal of wild type (D) and isl1-/- mutants (H) (arrows). Scale bars represent 50 μm. See also Figure 5—figure supplement 1.
Figure 5—figure supplement 1
Expression profiles of Wnt-signaling components.
(A) RNA expression plots of Wnt-signaling components expressed in the pacemaker domain (dark red). Superimposed expression profiles of ventricular marker gene vmhc (light blue) and atrial marker gene myh6 (light green) outline the position of the two chambers. Smoothening (LOESS) was applied to vmhc and myh6 graphs, span α = 0.3. SA, sinoatrial region.
Figure 6 with 5 supplements
Wnt/β-catenin signaling regulates parasympathetic control of heart rate.
(A) Experimental timeline. Embryos with the ubiquitous tg(TETRE:Mmu.Axin1-YFP)tud1 and cardiomyocyte-specific tg(myl7:TETAM2-2A-mCherry)ulm8 transgenes were incubated with doxycycline (25 μg/mL) plus dexamethasone (100 μM) to induce cardiomyocyte specific Axin1-YFP expression. Axin1-YFP was induced at different time points and the heart rate was measured at 3 dpf by high-speed video imaging and image analysis. (B) Atrial and corresponding ventricular kymographs from a representative mCherry-/YFP- sibling control and an embryo with Axin1-YFP overexpression (Axin1+). Note the shorter period per full heartbeat (indicated by white dotted lines) in the Axin1 expressing embryo. Movies are available as Figure 6—videos 1 and 2, respectively. (C) Result of heart rate analysis (ventricle and atrium) on sibling and Axin1+ embryos, showing a significant increase in heart rate in Axin1+ embryos. (D) Relative changes in heart rate measured at 3 dpf after induction of Axin1 expression at various time points. ∆heart rate was calculated by using the heart rates measured in the Axin1+ group and the control group using the following formula: (Axin1-YFP+ - mCherry-/YFP-) / (mCherry-/YFP-) x 100%. Induction of Axin1-YFP expression between 32 hpf and 56 hpf resulted in increased heart rates at 3 dpf. See also Figure 6—videos 1 and 2. n(control/Axin1+) 9/6 (12som); 20/15 (24hpf); 10/7 (28hpf); 18/16 (32hpf); 10/8 (36hpf); 30/21 (48hpf); 37/28 (52hpf); 27/23 (56hpf); 30/29 (72hpf). Column bar graph plotting mean with SEM. Atrial (green) and ventricular (blue) analysis (E–F) Heart rates at 3 dpf of control (mCherry-/YFP-) or Axin1-YFP+ embryos at baseline or after treatment with the parasympathetic agonist carbachol or the sympathetic agonist isoproterenol. Atrial (E) and ventricular (F) measurements. (left n = 27 (control), 26 (Axin+); right n = 20 (control), 14 (Axin+)) (G) Proposed model for Isl1/Wnt function in cardiac pacemaker cells. Isl1 is a central factor in pacemaker cell development. Isl1 expression is restricted to pacemaker cells in which it is required to establish the rhythmic cycle of membrane depolarizations and heart contractions and the activation of Wnt/β-catenin signaling. Wnt/β-catenin signaling is required to establish M2 cholinergic receptor signaling, which allows parasympathetic signals (e.g. acetylcholine) to adjust the rate of rhythmic membrane depolarizations and heart contractions. βAR, β-Adrenergic Receptor. Statistical significance (unpaired (D) and paired (E–F) t-test) *p<0.05; **p<0.01; ***p<0.001; ****p<0.0001. See also Figure 6—source data 1, Figure 6—figure supplement 1, Figure 6—figure supplement 2, Figure 6—figure supplement 3, Figure 6—video 1 and Figure 6—video 2.
-
Figure 6—source data 1
Heart rate measurements and detailed statistics shown in Figure 6C–F.
- https://doi.org/10.7554/eLife.31515.017
Figure 6—figure supplement 1
Heart rate variability in Axin1+ embryonic hearts.
Heart rate variability (HRV) in control and Axin1 overexpression embryos (n = 5 for each condition). Standard deviation of NN intervals (SDNN) between heart beats was measured on atrial and ventricular kymographs over a period of 10 s. Wild-type siblings in gray, Axin1+ embryos in orange. HRV in Axin1+ embryos is decreased. Statistical significance (unpaired t-test) *p<0.05; **p<0.01. See also Figure 6—figure supplement 1—source data 1.
-
Figure 6—figure supplement 1—source data 1
Heart rate measurements and detailed statistics shown in Figure 6—figure supplement 1.
- https://doi.org/10.7554/eLife.31515.014
Figure 6—figure supplement 2
No neuronal innervation was detected in the 3 dpf heart.
Confocal scans from whole mount embryos containing the tg(myl7-DsRed) (marking all cardiomyocytes) and the tg(Isl1:GFF;UAS:GFP) (reporting isl1 expression and marking pacemaker cells in green) and stained for acetylated tubulin (marking acetylated microtubules found in neurons and axons). Anterior is up and posterior is down. 3D reconstruction (A–D) and cross-sections (E–M). The skin immediately ventral to the heart is covered in axons. The Isl1-GFP-positive pacemaker cells (arrowhead) at the sinoatrial junction are not in direct contact with the acTubulin+ axons.
Figure 6—figure supplement 3
M2 muscarinic receptor is expressed at the venous pole.
(A–B) In situ hybridization for the muscarinic 2a cholinergic receptor (chrm2a) in wild-type embryos at 55 hpf (A) and 72hpf (B). (C) ISH for chrm2a at 72 hpf in Axin1-YFP+ embryo (Wnt inhibition started at 52 hpf). Arrowheads indicate expression at the venous pole. Gray dashed line: location of the heart; white dashed line: location of the eyes. Anterior up. Scale bars represent 50 μm.
Figure 6—video 1
Bright-field video recording of a mCherry-/YFP- sibling heart at 3dpf.
Lateral view (ventral to the left, anterior on the top). Note the regular, sequential beating of the heart at approximately 160 bpm. Also see Figure 6B. The playback speed is set to 150 fps, matching acquisition speed.
Figure 6—video 2
Bright field video recording of a Axin-YFP+ TetON heart at 3dpf.
Lateral view (ventral to the left, anterior on the top). The embryo was treated with doxycycline/dexamethasone from 52hpf onwards. Note the significantly elevated heart rate of approximately 250bpm. Also see Figure 6B. The playback speed is set to 150 fps, matching acquisition speed.
Tables
Key resources table
| Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |
|---|---|---|---|---|
| Gene (Danio rerio) | mb | NA | ZDB-GENE-040426–1430 | |
| Gene (D. rerio) | smtnl1 | NA | ZDB-GENE-050306–23 | |
| Gene (D. rerio) | anxa5b | NA | ZDB-GENE-030131–9076 | |
| Gene (D. rerio) | pdgfaa | NA | ZDB-GENE-030918–2 | |
| Gene (D. rerio) | fzd9b | NA | ZDB-GENE-000906–4 | |
| Gene (D. rerio) | shox2 | NA | ZDB-GENE-040426–1457 | |
| Strain, strain background (D. rerio) | Tupfel long fin (TL) | ZIRC,Eugene, OR | ZDB-GENO-990623–2 | |
| Genetic reagent (D. rerio) | Tg(UAS:GFP) | Kawakami lab Asakawa 2008 | ZDB-FISH-150901–15231 | |
| Genetic reagent (D. rerio) | Tg(Isl1:GFF;UAS:GFP)hu10018tg | this paper | see Material and methods section | |
| Genetic reagent (D. rerio) | Tg(7xTCFXla.Siam:nlsmCherry) | this paper | The tg(7xTCF-nlsmCherry) wnt-signaling reporter line was established in a TL wildtype background using the pDEST-7xTCF-nlsmCHERRY-polyA construct (Moro et al., 2012) received from the Smith lab (IMB, Brisbane, Australia). | |
| Genetic reagent (D. rerio) | isl1 mutant line. isl1sa29 | ZIRC, Eugene, OR | ZDB-GENO-990623–2 | |
| Genetic reagent (D. rerio) | Tg(myl7:EGFP): twu26tg | ZIRC, Eugene, OR | ZDB-FISH-150901–23018 | |
| Genetic reagent (D. rerio) | Tg(TETRE:Mmu.Axin1-YFP)tud1: tud1Tg | ZIRC, Eugene, OR, Gilbert Weidinger, Knopf 2012 | ZDB-FISH-150901–26117 | |
| Genetic reagent (D. rerio) | Tg(myl7:TETAM2-2A-mCherry)ulm8: ulm8Tg | ZIRC, Eugene, OR, Gilbert Weidinger, Knopf 2013 | ZDB-ALT-151214–5 | |
| Antibody | Anti-GFP (Chicken Antibodies, IgY Fraction) | Aves Labs, Tigard, OR | Cat# GFP-1010; RRID: AB_2307313 | 1:500 |
| Antibody | Living Colors DsRed Polyclonal Antibody | Clontech,Mountain View, CA | Cat# 632496; RRID: AB_10013483 | 1:200 |
| Antibody | Monoclonal Anti-Acetylated Tubulin antibody | Sigma-Aldrich | Cat# T7451; RRID: AB_609894 | |
| Antibody | Monoclonal Anti-Tropomyosin (Sarcomeric) antibody | Sigma-Aldrich | Cat# T9283; RRID: AB_261817 | 1:200 |
| Recombinant DNA reagent | BAC clone CH211-219F7 | BACPAC Resources Center (BPRC), Oakland, CA | ZDB-BAC-050218–781 | |
| Recombinant DNA reagent | pDEST-7xTCF-nlsmCHERRY-polyA | Moro et al. (2012) provided by Kelly Smith | ZDB-TGCONSTRCT-110113–2 | |
| Commercial assay or kit | TruSeq small RNA sample prep Kit | Illumina, San Diego, CA | RS-200–0012 or −0024,–0036, −0048 | |
| Commercial assay or kit | MessageAmpII Kit | Ambion, Foster City, CA | AM1751 | |
| Chemical compound, drug | Doxycycline hyclate | Sigma-Aldrich | D9891; CAS:24390-14-5 | |
| Chemical compound, drug | Dexamethasone | Sigma-Aldrich | D1756 CAS:50-02-2 | |
| Chemical compound, drug | Isoprenaline hydrochloride | Sigma-Aldrich | I5627; CAS:51-30-9 | |
| Chemical compound, drug | Carbamoylcholine chloride | Sigma-Aldrich | C4382; CAS:51-83-2 | |
| Software, algorithm | RStudio | https://www.rstudio.com/ RRID:SCR_000432 | ||
| Software, algorithm | Imaris data visualization software | Bitplane, Zürich, Switzerland | http://www.bitplane.com/imaris/imaris RRID:SCR_007370 | |
| Software, algorithm | HoKaWo Image acquisition module | Hamamatsu, Shizuoka, Japan | https://www.hamamatsu.com/eu/en/product/alpha/I/5006/U9304/index.html | |
| Other | Resource website with the complete 2dpf heart dataset | This paper | http://zebrafish.genomes.nl/tomoseq/ |
Additional files
-
Source Code 1
R source code
R script code used to process, analyze and visualize the tomo-seq datasets.
- https://doi.org/10.7554/eLife.31515.020
-
Transparent reporting form
- https://doi.org/10.7554/eLife.31515.021
Download links
A two-part list of links to download the article, or parts of the article, in various formats.
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Spatially resolved RNA-sequencing of the embryonic heart identifies a role for Wnt/β-catenin signaling in autonomic control of heart rate
eLife 7:e31515.
https://doi.org/10.7554/eLife.31515
