Single-cell analysis uncovers that metabolic reprogramming by ErbB2 signaling is essential for cardiomyocyte proliferation in the regenerating heart

  1. Hessel Honkoop
  2. Dennis EM de Bakker
  3. Alla Aharonov
  4. Fabian Kruse
  5. Avraham Shakked
  6. Phong D Nguyen
  7. Cecilia de Heus
  8. Laurence Garric
  9. Mauro J Muraro
  10. Adam Shoffner
  11. Federico Tessadori
  12. Joshua Craiger Peterson
  13. Wendy Noort
  14. Alberto Bertozzi
  15. Gilbert Weidinger
  16. George Posthuma
  17. Dominic Grün
  18. Willem J van der Laarse
  19. Judith Klumperman
  20. Richard T Jaspers
  21. Kenneth D Poss
  22. Alexander van Oudenaarden
  23. Eldad Tzahor
  24. Jeroen Bakkers  Is a corresponding author
  1. Hubrecht Institute-KNAW and University Medical Center Utrecht, Netherlands
  2. Weizmann Institute of Science, Israel
  3. University Medical Center Utrecht, Utrecht University, Netherlands
  4. Duke University Medical Center, United States
  5. Vrije Universiteit Amsterdam, Netherlands
  6. Ulm University, Germany
  7. Max Planck Institute of Immunobiology and Epigenetics, Germany
  8. VU University Medical Center, Netherlands
  9. University Medical Center Utrecht, Netherlands
10 figures, 1 table and 1 additional file

Figures

Figure 1 with 4 supplements
Single-cell mRNA sequencing identifies different cardiomyocyte-populations in the injured zebrafish heart.

(a) Schematic of cryoinjury procedure on adult TgBAC(nppa:mCitrine) fish and immunohistochemistry on section of injured TgBAC(nppa:mCitrine) heart 7 dpi. Overview image on left and zoom-in of boxed …

Figure 1—source data 1

Single-cell mRNA sequencing data from cryoinjured TgBAC(nppa:mCitrine) zebrafish hearts at 7 dpi.

Raw reads are represented per gene (rows) and single cell (columns).

https://cdn.elifesciences.org/articles/50163/elife-50163-fig1-data1-v2.xlsx
Figure 1—source data 2

List of the clusters in the adult injured heart (raw data see SD1), as identified by StemID and the associated cells per cluster.

https://cdn.elifesciences.org/articles/50163/elife-50163-fig1-data2-v2.xlsx
Figure 1—source data 3

List of differentially expressed genes between cardiomyocytes clusters 2 and 7 of the adult heart dataset.

Only genes with a p-value<0.05 are listed. Per gene, the log2 fold change, adjusted p-value (padj) and associated gene name are given. GO-terms for genes upregulated between clusters 2 and 7 (p<0.01) are listed in separate excel sheets.

https://cdn.elifesciences.org/articles/50163/elife-50163-fig1-data3-v2.xlsx
Figure 1—figure supplement 1
TgBAC(nppa:mCitrine) expression recapitulates endogenous nppa gene expression.

(a) Design of the bacterial artificial chromosome (BAC) used for generation of the transgenic line TgBAC(nppa:mCitrine). (b) Transgenic mCitrine expression in the heart in relation to RFP expression …

Figure 1—figure supplement 2
Single-cell mRNA sequencing identifies different cell-populations in the injured zebrafish heart.

tSNE maps visualizing log2-transformed read-counts of genes with high expression in cardiomyocytes (mt-co1, tnnc1a, tnnt2a), in fibroblasts (fn1b), in endothelial cells (cdh5) and immune cells (ctsd)…

Figure 1—figure supplement 3
Cluster 7 cells display highest nppa:mCitrine expression.

(a) tSNE map of the cardiomyocyte clusters derived from the single-cell mRNA-sequencing. (b) mCitrine fluorescence levels of the nppa:mCitrine FACS sorted cells from injured hearts that were used …

Figure 1—figure supplement 4
Enhanced expression of genes elevated in cluster 7 versus cluster 2 is injury induced.

(a) In situ hybridizations for border zone genes nppb, mustn1, nppa and vmhc in zebrafish hearts at seven dpi compared to uninjured hearts. (b) In situ hybridizations for glycolysis genes hk1 and ldh…

Single-cell transcriptome of border zone cardiomyocytes resembles that of embryonic cardiomyocytes.

(a) Cartoon to illustrate the experimental procedure for single cell analysis of embryonic and adult cardiac cells. (b) tSNE map of combined adult (red) and embryonic datasets (light blue). (c) tSNE …

Figure 2—source data 1

Single-cell mRNA sequencing data from Tg(cmlc2:GFP) zebrafish hearts at 2 days post fertilization.

Raw reads are represented per gene (rows) and single cell (columns).

https://cdn.elifesciences.org/articles/50163/elife-50163-fig2-data1-v2.csv
Figure 2—source data 2

List of the clusters from the embryonic heart data (SD5) identified by StemID and the associated cells for each cluster.

https://cdn.elifesciences.org/articles/50163/elife-50163-fig2-data2-v2.xls
Figure 2—source data 3

Combined single-cell mRNA sequencing data from embryonic and adult hearts.

Raw reads represented per gene (rows) and single cell (columns).

https://cdn.elifesciences.org/articles/50163/elife-50163-fig2-data3-v2.csv
Figure 2—source data 4

List of pairwise differentially expressed genes between all cardiomyocyte clusters from the combined embryonic and adult datasets.

Only genes with a p-value<0.01 are listed. Per gene, the log2 fold change, adjusted p-value (padj) and associated gene name are given. GO-terms for genes upregulated between clusters 5 and 6 (p<0.01) are listed in separate excel sheets.

https://cdn.elifesciences.org/articles/50163/elife-50163-fig2-data4-v2.xls
Figure 3 with 1 supplement
Pseudo time analysis reveals dedifferentiation and metabolic changes in border zone cardiomyocytes.

(a) Bar plot of StemID scores for the cardiomyocyte clusters (clusters #2, 1, 4 and 7) calculated by the formula: number of significant links for each cluster multiplied by the median transcriptome …

Figure 3—source data 1

List of genes that are differentially upregulated in the respective modules identified in the pseudo timeline analysis (related to Figure 3c).

https://cdn.elifesciences.org/articles/50163/elife-50163-fig3-data1-v2.xlsx
Figure 3—figure supplement 1
Induction of hexokinase1 expression precedes expression of myomesin1b.

In situ hybridization for the glycolysis gene hexokinase1 (a,b) and the embryonic cardiac gene myomesin1b (c,d) on sections of injured zebrafish hearts at 3 dpi (a,c) and 7 dpi (b,d). While hk1

Figure 4 with 1 supplement
Border zone cardiomyocytes undergo a metabolic switch from mitochondrial OXPHOS to glycolysis.

(a) Succinate dehydrogenase (SDH) enzyme staining on a seven dpi heart section with injury area separated by dashed line. Quantification of SDH activity in remote zone (RZ), border zone (BZ) and …

Figure 4—figure supplement 1
Energy metabolism genes are differentially expressed between cluster 7 and cluster 2 cells.

(a) Plot showing differentially expressed genes between cluster 7 versus cluster 2 cells. Differentially expressed genes (p-value<0.05) are highlighted in red (upregulated in cluster 7 cells) and …

2-Deoxy glucose impairs cardiomyocyte proliferation.

(a) Experimental design for the 2-DG injections to inhibit glycolysis in injured zebrafish hearts. (b) Confocal image of injured zebrafish hearts at seven dpi either injected with PBS or 2-DG …

Figure 6 with 1 supplement
Nrg1/ErbB2 signaling induces glycolysis genes in zebrafish and mouse.

(a) Cartoon showing experimental procedure to induce cardiomyocyte specific Nrg1 expression in zebrafish. Panels show in situ hybridization for hexokinase 1 (hk1) and lactate dehydrogenase a (ldha) …

Figure 6—source data 1

Oligonucleotide sequences for real-time PCR analysis.

https://cdn.elifesciences.org/articles/50163/elife-50163-fig6-data1-v2.xlsx
Figure 6—figure supplement 1
Hypoxia sensor phd3:GFP expression does not correlate with expression of the glycolytic gene ldha.

(a) Overview picture of phd3:GFP heart 7dpi co-stained for expression of ldha by fluorescent ISH. (b) Zoom-ins of trabecular border zone cardiomyocytes. Note high ldha expression and low expression …

Figure 7 with 1 supplement
Glycolysis inhibitors impair mitogenic effect of ErbB2 activation in cardiomyocytes.

(a) Cartoon showing the experimental procedure to analyse the effects of glycolysis inhibitors (2-DG and lonidamine) on cardiomyocyte proliferation. (b) Immunofluorescence analysis on P7 cardiac …

Figure 7—figure supplement 1
Activation of ErbB2 signaling in murine MI model induces glycolytic gene expression while repressing mitochondrial genes.

(a) Cartoon showing experimental procedure. Myocardial infarction (MI) was induced by left anterior descending coronary artery ligation at the same day when Dox was removed to induce cardiomyocyte …

Author response image 1
Different read cut-offs do not alter clustering pattern, while altering number of included cells.

(A) Distribution of reads per cell. Green line represents 3,500 reads, red line indicates the mean number of reads per cell (10,443). (B) tSNE maps generated with different read cutoffs. Read …

Author response image 2
GO-term analysis on differentially expressed genes in cardiomyocyte clusters.

Representative GO-terms for cardiomyocyte clusters. GO-terms in upper table are obtained from differentially expressed genes between adult clusters 2 and 7 (p<0.01). Lower table GO-terms are …

Author response image 3
Glycolytic gene expression in the embryonic zebrafish heart.

Expression of eno1a (a, b, c, d), hk1 (e, f, g, h) and ldha (i, j, k, l) during zebrafish development at 28hpf (a, e, i), 48 hpf (b, f, j) and 5dpf (c, d, g, h, k, l) respectively. Outline of the …

Tables

Key resources table
Reagent type
(species) or
resource
DesignationSource or
reference
IdentifiersAdditional
information
Strain, strain background (Danio rerio)Tupfel Long Fin (TL)ZIRCZDB-GENO-990623–2
Genetic reagent (Danio rerio)Tg(myl7:dsRED)s879TgChi et al., 2008ZDB-FISH-150901–3078
Genetic reagent (Danio rerio)Tg(myl7:GFP)twu34TgHuang et al., 2003ZDB-FISH-150901–212
Genetic reagent (Danio rerio)Tg(phd3:GFP)sh144Santhakumar et al., 2012ZDB-FISH-150901–26851
Genetic reagent (Danio rerio)Tg(gata4:EGFP)ae1Heicklen-Klein and Evans, 2004ZDB-FISH-150901–14762
Genetic reagent (Danio rerio)TgBAC(nppa:mCitrine)This paperMore info on generation of this line can be obtained from the Materials and methods section ‘Transgenic zebrafish lines and cryoinjury’.
Cell line (Mus musculus)TetRE-caErbb2 X MYH6-tTA X MYH6-cre ROSA26-tdTomatoD'Uva et al., 2015More info on how to obtain primary cell line can be obtained from Materials and methods sections ‘Transgenic mouse lines and animal procedures’ and ‘Pharmacological inhibition of glycolysis’.
Biological sample (Mus musculus)TetRE-caErbb2 X MYH6-tTA heartsD'Uva et al., 2015 from:
Xie et al., 1999 and
Yu et al., 1996.
The Jackson Laboratory, stock no. 010577;
Biological sample (Danio rerio)Tg(cmlc2:CreER)pd10 x Tg(β-act2:BSNrg1)pd107heartsGemberling et al., 2015ZDB-FISH-150901–25249
x
ZDB-FISH-150901–25354
Biological sample (Danio rerio)Tg(gata4:EGFP)ae1heartsLin et al. 2009ZDB-FISH-150901–14762
Biological sample (Danio rerio)Tg(phd3:GFP)sh144heartsSanthakumar et al., 2012ZDB-FISH-150901–26851
AntibodyAnti-AuroraB kinase (mouse monoclonal)BD Transduction Laboratories#611082, RRID:AB_22277081:200
AntibodyAnti-Ki67 (Rabbit monoclonal)Cell Marque#275R1:200
AntibodyAnti-MCT4 (rabbit polyclonal)Santa Cruz#SC50329, RRID:AB_21893331:200
AntibodyAnti-PCNA (mouse monoclonal)DAKO#M0879, RRID:AB_21606511:800
AntibodyAnti-GFP (chicken polyclonal)Aves#GFP-1010, RRID:AB_23073131:1000
AntibodyAnti-Mef2C (rabbit polyclonal)Santa Cruz/Biorbyt#SC313, RRID:AB_631920 / #orb256682Both 1:1000
Chemical compound, drug2-DeoxyglucoseSigma-Aldrich#D61341 mg/g or 3 mM
Chemical compound, drugLonidamineSigma-AldrichL490080 uM
Chemical compound, drug2NBDGCaymanchem#11046400 uM
Chemical compound, drugAG1478Selleck ChemicalS27285M, from 10 mM stock in DMSO
Chemical compound, drugDoxycyclineHarlan LaboratoriesTD02503
Chemical compound, drugTrizolLife technologies#15596026
Chemical compound, drugFast SYBR Green Master MixApplied Biosystems#4385612
Chemical compound, drugCollagenase type IIGibco171010150.1%
Chemical compound, drugTrypLE Express Enzyme (1x), phenol redGibco12605036
Software, algorithmFIJISchindelin et al., 2012Version 2.0.0, RRID:SCR_002285
Software, algorithmRaceID2/StemIDGrün, D. et al. De Novo Prediction of Stem Cell Identity using Single-Cell Transcriptome Data. Cell Stem Cell 19, 266–277 (2016).RRID:SCR_017242
Software, algorithmRstudioRstudioRRID:SCR_000432Version 1.2.1335
Software, algorithmImarisBitplaneRRID:SCR_0007370V9.3.1
Software, algorithmGene Set Enrichment AnalysisGenepattern, Broad InstituteRRID:SCR_003199# of permutations = 1000
Commercial assay, kitSuperscript III First strand synthesis systemThermo Fisher Scientific#18080051Input 200 ng RNA
Commercial assay, kitMiRNeasyQiagen217004
Commercial assay, kitHigh capacity cDNA Reverse Transcription kitApplied Biosystems4374966Input 1 ug RNA
Commercial assay, kitNeonatal Dissociation kitMiltenyi Biotec130-098-373
Commercial assay, kitTruSeq small RNA primersIllumina20005613

Additional files

Download links