Single-cell transcriptomics defines heterogeneity of epicardial cells and fibroblasts within the infarcted murine heart

  1. Julia Hesse
  2. Christoph Owenier
  3. Tobias Lautwein
  4. Ria Zalfen
  5. Jonas F Weber
  6. Zhaoping Ding
  7. Christina Alter
  8. Alexander Lang
  9. Maria Grandoch
  10. Norbert Gerdes
  11. Jens W Fischer
  12. Gunnar W Klau
  13. Christoph Dieterich
  14. Karl Köhrer
  15. Jürgen Schrader  Is a corresponding author
  1. Department of Molecular Cardiology, Medical Faculty and University Hospital Düsseldorf, Heinrich-Heine-University Düsseldorf, Germany
  2. Biologisch-Medizinisches-Forschungszentrum (BMFZ), Genomics & Transcriptomics Laboratory, Heinrich-Heine-University Düsseldorf, Germany
  3. Algorithmic Bioinformatics, Heinrich-Heine-University Düsseldorf, Germany
  4. Division of Cardiology, Pulmonology and Vascular Medicine, Medical Faculty and University Hospital Düsseldorf, Heinrich-Heine-University Düsseldorf, Germany
  5. Institute of Pharmacology and Clinical Pharmacology, Medical Faculty and University Hospital Düsseldorf, Heinrich-Heine-University Düsseldorf, Germany
  6. Section of Bioinformatics and Systems Cardiology, Klaus Tschira Institute for Integrative Computational Cardiology and Department of Internal Medicine III, University Hospital Heidelberg, Germany
4 figures, 1 table and 9 additional files

Figures

Figure 1 with 4 supplements
Cell populations in EpiSC from the infarcted heart.

(A) Schematic workflow. Epicardial stromal cells (EpiSC) and activated cardiac stromal cells (aCSC) were simultaneously collected from the surface and the myocardium of the isolated perfused heart …

Figure 1—source data 1

Source data for EpiSC population cell counts summarized in Figure 1C.

https://cdn.elifesciences.org/articles/65921/elife-65921-fig1-data1-v1.xlsx
Figure 1—figure supplement 1
Excluded hybrid and non-stromal cell populations in the EpiSC fraction.

(A) DoubletFinder tool was used to detect cell doublets with hybrid transcriptomes. (B) Heat map showing the expression of markers for cardiomyocytes (CM) and erythrocytes (Ery) used to identify …

Figure 1—figure supplement 2
Gene ontology terms typifying EpiSC populations.

Over-represented gene ontology (GO) terms in epicardial stromal cell (EpiSC) population marker genes as analyzed by Scillus.

Figure 1—figure supplement 3
Top 10 marker genes of EpiSC populations.

Heat map showing the expression of the top 10 most differentially expressed genes in identified epicardial stromal cell (EpiSC) populations (color-coded according to their identity).

Figure 1—figure supplement 4
Receptor–ligand pairs in EpiSC populations.

Feature plots and violin plots visualizing the expression levels of selected receptor–ligand pairs across different epicardial stromal cell (EpiSC) populations as identified by CellPhoneDB. (A) …

Figure 2 with 3 supplements
Molecular characterization and location of EpiSC populations.

(A) Expression of epicardial progenitor cell markers proposed in the literature as visualized in feature and violin plots. (B) Upper panel: RNA in situ hybridization of epicardial stromal cell …

Figure 2—figure supplement 1
Subclusters of Wt1-expressing EpiSC populations.

(A) UMAP plot of subclusters (Sub) of Wt1-expressing epicardial stromal cell (EpiSC)-7, -1, and -8. (B) UMAP plot of the subclustering with labeling according to previous identity. (C) Transfer of …

Figure 2—figure supplement 2
Expression of epicardial markers previously identified in the developing zebrafish heart in mouse post-MI EpiSC.

Expression of epicardial population markers identified in the developing zebrafish heart by Weinberger et al., 2020 visualized in the whole epicardial stromal cell (EpiSC) fraction (A and B) and in W…

Figure 2—figure supplement 3
Wt1-expressing subclusters in Wt1CreERT2RosatdTomato EpiSC.

UMAP plot of the epicardial stromal cell (EpiSC) fraction from lineage tracing experiments using Wt1CreERT2RosatdTomato mice with assignment of previously defined subclusters of Wt1-expressing EpiSC.

Functional characterization of EpiSC populations.

(A, B) Heat maps (A) and violin plots (B) of selected genes showing the expression of cardiogenic factors and Notch target genes. (C, D) Heat maps (C) and violin plots (D) of genes associated with …

Figure 4 with 7 supplements
Comparison of EpiSC, aCSC, and CSC.

(A) Canonical correlation analysis (CCA) space alignment of epicardial stromal cells (EpiSC), activated cardiac stromal cells (aCSC), and CSC single-cell RNA sequencing (scRNAseq) data in one UMAP …

Figure 4—figure supplement 1
Cell populations in aCSC from infarcted myocardium.

(A) UMAP plot of clustered single-cell RNA sequencing (scRNAseq) data of the pooled activated cardiac stromal cell (aCSC) fraction (n = 24,470 single cells). Each point represents a single cell, and …

Figure 4—figure supplement 1—source data 1

Source data for aCSC population cell counts summarized in Figure 4—figure supplement 1B.

https://cdn.elifesciences.org/articles/65921/elife-65921-fig4-figsupp1-data1-v1.xlsx
Figure 4—figure supplement 2
Excluded hybrid and non-stromal cell populations in the aCSC fraction.

(A) DoubletFinder tool was used to detect cell doublets with hybrid transcriptomes. (B) Heat map showing the expression of markers for mural cells, endothelial cells (EC), immune cells, …

Figure 4—figure supplement 3
Top 10 marker genes of aCSC populations.

Heat map showing the expression of the top 10 most differentially expressed genes in identified activated cardiac stromal cell (aCSC) populations (color-coded according to their identity).

Figure 4—figure supplement 4
Cell populations in CSC from uninjured hearts.

(A) UMAP plot of clustered single-cell RNA sequencing (scRNAseq) data of the pooled cardiac stromal cell (CSC) fraction (n = 24,781 single cells). Each point represents a single cell and identified …

Figure 4—figure supplement 4—source data 1

Source data for CSC population cell counts summarized in Figure 4—figure supplement 4B.

https://cdn.elifesciences.org/articles/65921/elife-65921-fig4-figsupp4-data1-v1.xlsx
Figure 4—figure supplement 5
Excluded hybrid and non-stromal cell populations in the CSC fraction.

(A) DoubletFinder tool was used to detect cell doublets with hybrid transcriptomes. (B) Heat map showing the expression of markers for mural cells, cardiomyocytes (CM), endothelial cells (EC), glial …

Figure 4—figure supplement 6
Top 10 marker genes of CSC populations.

Heat map showing the expression of the top 10 most differentially expressed genes in identified cardiac stromal cell (CSC) populations (color-coded according to their identity).

Figure 4—figure supplement 7
Comparison of EpiSC to CSC.

(A–C) Relative expression of epicardial and stromal cell markers (A), genes associated with HIF-1 signaling (B), and genes associated with cell proliferation (C) in cardiac stromal cell (CSC) and …

Figure 4—figure supplement 7—source data 1

Source data for differentially expressed genes among EpiSC and CSC in CCA cluster A displayed in Figure 4—figure supplement 7D.

https://cdn.elifesciences.org/articles/65921/elife-65921-fig4-figsupp7-data1-v1.xlsx

Tables

Key resources table
Reagent type
(species) or resource
DesignationSource or referenceIdentifiersAdditional information
Strain, strain background (Mus musculus)C57BL/6JJanvier LabsStrain name: C57BL/6JRj
Strain, strain background (M. musculus)STOCK Wt1tm2(cre/ERT2)Wtp/JThe Jackson Laboratory;
PMID:18568026
Stock no. 010912; RRID:IMSR_JAX:010912
Strain, strain background (M. musculus)B6;129S6-Gt(ROSA)26Sortm14(CAG-tdTomato)Hze/JThe Jackson LaboratoryStock no. 007908; RRID:IMSR_JAX:007908
Chemical compound, drugTamoxifenSigma-AldrichT5648
Chemical compound, drug7AADBioLegend420404(1:100)
Chemical compound, drugFixable Viability Dye eFluor 780eBioscience65–0865(1:100)
OtherCollagenase CLS IIBiochromC2-22Enzyme blend; now available from Sigma-Aldrich, C2-BIOC
AntibodyAPC-conjugated rat monoclonal anti-CD31BioLegend102410, clone 390; RRID:AB_312905(1:100)
AntibodyPE-Cy7 rat monoclonal anti-CD45BD Biosciences552848, clone 30-F11; RRID:AB_394489(1:400)
Commercial assay or kitRNAscope 2.5 HD Reagent Kit-REDAdvanced Cell Diagnostics322350
Sequence-based reagentRNAscope Probe Mm-MslnAdvanced Cell Diagnostics443241
Sequence-based reagentRNAscope Probe Mm-Wt1Advanced Cell Diagnostics432711
Sequence-based reagentRNAscope Probe Mm-Cd44Advanced Cell Diagnostics479191
Sequence-based reagentRNAscope Probe Mm-Ifit3Advanced Cell Diagnostics508251
Sequence-based reagentRNAscope Probe Mm-Sfrp2Advanced Cell Diagnostics400381
Sequence-based reagentRNAscope Probe Mm-Pcsk6Advanced Cell Diagnostics593361
Sequence-based reagentRNAscope Probe Mm-Top2aAdvanced Cell Diagnostics491221
Sequence-based reagentRNAscope Probe Mm-Dkk3Advanced Cell Diagnostics400931
Software, algorithmCellRanger10x Genomics
Software, algorithmSeuratPMID:29608179RRID:SCR_007322
Software, algorithmDoubletFinderPMID:30954475RRID:SCR_018771
Software, algorithmScillushttps://github.com/xmc811/
Scillus Xu, 2021
Software, algorithmGOrillaPMID:19192299RRID:SCR_006848
Software, algorithmVelocyto.pyPMID:30089906RRID:SCR_018167
Software, algorithmCellPhoneDBPMID:32103204RRID:SCR_017054
Software, algorithmSCENICPMID:28991892RRID:SCR_017247

Additional files

Supplementary file 1

Average gene expression levels in EpiSC, aCSC, and CSC populations.

https://cdn.elifesciences.org/articles/65921/elife-65921-supp1-v1.xlsx
Supplementary file 2

Genes with significantly enriched expression among EpiSC, aCSC, and CSC populations.

https://cdn.elifesciences.org/articles/65921/elife-65921-supp2-v1.xlsx
Supplementary file 3

Selected ligand–receptor interactions between EpiSC populations as predicted by CellPhoneDB.

https://cdn.elifesciences.org/articles/65921/elife-65921-supp3-v1.xlsx
Supplementary file 4

Average gene expression levels in subclusters of Wt1-expressing EpiSC populations.

https://cdn.elifesciences.org/articles/65921/elife-65921-supp4-v1.xlsx
Supplementary file 5

Genes with significantly enriched expression among subclusters of Wt1-expressing EpiSC populations.

https://cdn.elifesciences.org/articles/65921/elife-65921-supp5-v1.xlsx
Supplementary file 6

Genes with significantly enriched expression in group I and group II populations vs group III populations.

https://cdn.elifesciences.org/articles/65921/elife-65921-supp6-v1.xlsx
Supplementary file 7

Average gene expression levels in clusters from CCA space alignment of EpiSC, aCSC, and CSC with separation in EpiSC, aCSC, and CSC, including p-values of differential expression.

https://cdn.elifesciences.org/articles/65921/elife-65921-supp7-v1.xlsx
Supplementary file 8

Genes with significantly enriched expression among clusters from CCA space alignment of EpiSC, aCSC, and CSC.

https://cdn.elifesciences.org/articles/65921/elife-65921-supp8-v1.xlsx
Transparent reporting form
https://cdn.elifesciences.org/articles/65921/elife-65921-transrepform-v1.docx

Download links