Differential regulation by CD47 and thrombospondin-1 of extramedullary erythropoiesis in mouse spleen

  1. Rajdeep Banerjee
  2. Thomas J Meyer
  3. Margaret C Cam
  4. Sukhbir Kaur
  5. David D Roberts  Is a corresponding author
  1. Laboratory of Pathology, Center for Cancer Research, National Cancer Institute, National Institutes of Health, United States
  2. CCR Collaborative Bioinformatics Resource, Office of Science and Technology Resources, National Cancer Institute, National Institutes of Health, United States
7 figures, 7 tables and 1 additional file

Figures

Figure 1 with 2 supplements
Effects of Cd47 or Thbs1 gene disruption on spleen cell numbers and content of cells expressing erythropoietic precursor markers or the proliferation marker Ki67.

(A) Total spleen cell numbers in WT, Cd47−/− and Thbs1−/− C57BL/6 mice determined after lysis of RBC (mean ± SEM, n=3). Flow cytometry was performed to analyze gated singlet spleen cells stained with Ter119 antibody (B), CD34 antibody (C), Sca1 antibody (D), Ki67 antibody with the indicated gating for high and low expression (E), cKit antibody (F), Epor antibody (G), Ermap antibody (H), Gypa antibody (I), or Aqp1 antibody (J). Further analysis of Epor (K, O), Ermap (L, P), Gypa (M, Q), and Aqp1 expression (N, R) was performed after gating for Ter119 expression. The percentages of cells positive for the indicated surface markers are presented (n=3 or 4 mice of each genotype). p-values were determined using a two-tailed t test for two-samples assuming equal variances in GraphPad Prism. *=p < 0.05, **=p < 0.01, ***=p < 0.001.

Figure 1—figure supplement 1
Enlargement of spleen in the absence of CD47 and bulk RNAseq analysis of lin− Cd47−/− vs WT spleen cells.

(A) Representative images of spleens from WT, Cd47−/−, and Thbs1−/− mice. (B) Volcano plot for differentially expressed genes between lineage-depleted Naive Cd47−/− vs WT CD8+-enriched spleen cells. Red dots represent significant genes with >2 fold change and p<0.001, and the top 30 genes were labelled. (C) GSEA plot showing RBC Heme Metabolism gene set enrichment. (D) Heat map visualization of the top differentially expressed genes in the RBC Heme Metabolism gene set.

Figure 1—figure supplement 2
Flow cytometry analysis strategy.

The sequential flow cytometry gating strategy is illustrated.

Figure 2 with 1 supplement
Effects of Cd47 or Thbs1 gene disruption on the percentage of Ter119+ and Ter119 spleen cells expressing markers of multipotent and committed erythroid precursors and cell proliferation.

Spleen cells isolated from WT, Cd47−/− and Thbs1−/− mice (2 male and 2 female of each genotype) were costained with Ter119 antibody along with cKit, Ki67, Sca1 Ermap, Gypa, Epor, or Aqp1 antibodies and acquired on an LSRFortessa SORP. After gating for singlet cells, the percentages of Ter119+ and Ter119 cells positive for stem cell markers cKit (A) and Sca1 (B), high or low levels of the proliferation marker Ki67 (C, D), were compared among WT, Cd47−/− and Thbs1−/− mouse spleens (n=3–4). The proliferation of CD34+ and CD34 populations of Ter119+ spleen cells from WT, Cd47−/−, Thbs1−/− mice was evaluated by staining with CD34, Ter119 and Ki67 antibodies. Ter119+CD34+ cells (E), Ter119+CD34 cells (F), and Ter119CD34+ cells (G) were also quantified (left panels) and analyzed for the proliferation marker Ki67 (center and right panels). p-values were determined using a two-tailed t test for two-samples assuming equal variances in GraphPad Prism. *=p < 0.05, **=p < 0.01, ***=p < 0.001.

Figure 2—figure supplement 1
Flow cytometry analysis strategy.

The sequential flow cytometry gating strategy is illustrated.

Figure 3 with 1 supplement
Effects of Cd47 and Thbs1 gene deletion on stem cell and erythroid precursor populations in mouse spleen identified using single cell RNA sequence analysis.

(A) tSNE clustering analysis of lineage-depleted spleen cells from WT, Cd47−/− and Thbs1−/−mice (n = 3). The encircled area contains erythroid cells (clusters 12) and stem cells (cluster 14). (B) Cell type analysis using Immgen and Mouse RNAseq and SingleR (v.1.0) databases. (C) Distribution of WT, Cd47−/− and Thbs1−/− spleen cells in each cluster of the tSNE plot. (D) Enlarged plots of clusters 12 and 14 and cells in these clusters colored by genotype.

Figure 3—figure supplement 1
Single cell RNA sequence post-filter QC plots.

The scatter plots, histograms and violin plots of three samples lineage-depleted of WT, Cd47−/− and Thbs1−/− spleen cells are shown. The number of UMI per cell filter was set to exclude cells with <2000 UMI in nCount RNA. The Feature RNA plot shows the number of genes with non-zero expression detected per cell. Filtering was set to exclude cells having <15% mitochondrial gene expression as presented in the percent mt plot. The log10Genes per UMI plot represents the scores for complexity of the RNA library found in each cell. No filter was set in this row.

Figure 4 with 2 supplements
Effects of Cd47 and Thbs1 gene deletion on gene expression in stem cell and erythroid precursor clusters.

High resolution tSNE plots showing the distribution of mRNAs encoding the multipotent stem cell markers CD34 and Ly6a (Sca1) and Gata2, the erythropoietic markers Kit and Epor, the proliferation marker Mki67, erythroid differentiation transcription factors Klf1 and Gata1, and erythroid differentiation and extramedullary erythropoiesis markers Ermap, Aqp1, Tmem56, Trim10, Gypa1, Spta1, Sptb, Ebp42, Xpo1, and RanBP2 in clusters 12 and 14. Expression levels were normalized to maximum expression of each mRNA in these clusters.

Figure 4—source data 1

Differential expression of erythropoietic, stem cell, and proliferation associated markers in cell clusters 12 and 14.

https://cdn.elifesciences.org/articles/92679/elife-92679-fig4-data1-v1.docx
Figure 4—source data 2

Differential mRNA expression of the nuclear export protein Xpo1 and nuclear pore protein synthesis instructor Ranbp2 in cluster12 and CD34+ and CD34 subsets of cluster 12 cells.

https://cdn.elifesciences.org/articles/92679/elife-92679-fig4-data2-v1.docx
Figure 4—figure supplement 1
The distributions of mRNA expression of the indicated genes related to stem cells and erythropoietic cells are shown throughout the 18 clusters as a tSNE projection.
Figure 4—figure supplement 2
Violin plots of erythropoietic, stem cell, and proliferation associated marker mRNAs expressed in 18 cell clusters.

Data are from two female and one male mouse of each genotype.

Differential effects of Cd47 and Thbs1 gene deletion on mRNA expression levels in erythroid precursor and stem cell clusters 12 and 14.

(A) Violin plots comparing mRNA expression levels of the indicated genes in Cd47−/−(red), Thbs1−/−(blue) and WT spleen cells (green) in cluster 12. (B) Violin plots comparing mRNA expression levels of the indicated genes in Cd47−/−, Thbs1−/−, and WT spleen cells in cluster 14. *=p < 0.05 relative to WT cells in the respective cluster. n = 2 female and 1 male mice of each genotype; *=p < 0.05.

Figure 6 with 1 supplement
Re-clustering of lineage-depleted spleen cells selected for expression of erythroid signature genes.

(A) tSNE plot showing the distribution of cells selected for expressing threshold levels of Gypa, Ermap, Klf1, Gata1, and/or Aqp1. (B) Re-clustered cells expressing the erythroid progenitor signature are displayed in an UMAP projection. (C) Immgen and mouse RNAseq main cell type annotation of reclustered cells expressing the erythroid gene signature. (D) Distribution of WT (purple), Cd47−/− (green), and Thbs1−/− cells (orange) in the erythroid precursor and T cell clusters. Data are from two female and one male mouse of each genotype.

Figure 6—figure supplement 1
Strategy for reclustering spleen cells that express a gene signature for committed erythroid precursors.

(A) The left panels show the distribution of cells expressing threshold levels of the 5 gene signature in the original TSNE projection. (B) The RBC progenitor module scoring was performed, as module Scores (a.k.a. Signature Scores) calculated for each cell using expression of five genes (Gypa, Ermap, Klf1, Gata1, and Aqp1). The threshold was set manually (red dashed line on subplots) to separate high-scoring cells from low-scoring cells and expressed by the color on the TSNE plot (top left, high score cells are dark red and low scores cells are represented as light red). RBC Progenitor Module Scoring Plots show the distributions of module scores in lineage-depleted spleen cells from each mouse, with the threshold score indicated by the dotted red lines. (C) The volcano plot presents differentially expressed genes contrasting the reclustered erythropoietic progenitor cell clusters and the T cell cluster.

Gene expression in reclustered lineage-depleted spleen cells selected for expression of erythroid signature genes.

Distribution of the multipotent stem cell markers CD34 and Ly6a (Sca1), erythropoietic markers Gata2, Kit, and Epor, the erythroid differentiation transcription factors Klf1 and Gata1, and erythroid differentiation and extramedullary erythropoiesis marker the proliferation marker Mik67, and the erythroid markers Ermap, Tfrc, Aqp1, Tfrc, Tmem56, Trim10, Gypa1, Ebp42, Spta1, Sptb, Xpo1, and Ranbp2 in the erythroid lineage cluster (left) and T cell cluster (right). Expression levels were normalized to maximum expression of each mRNA in these clusters. Data are from two female and one male mouse of each genotype.

Figure 7—source data 1

Differential expression of erythropoietic, stem cell, and proliferation associated markers in reclustered erythroid and T cell clusters.

https://cdn.elifesciences.org/articles/92679/elife-92679-fig7-data1-v1.docx

Tables

Table 1
Expression of erythropoiesis-associated genes in lineage-depleted Cd47−/− versus WT spleen cells.
GeneErythropoiesis expression/functionFold change Cd47-/-/WT*T statisticp-value
Ermapextramedullary erythropoiesis marker21.56.049.35x10–5
Tal1extramedullary erythropoiesis marker13.95.133.55x10–4
Gypaextramedullary erythropoiesis marker89.33.673.86x10–3
Gata1extramedullary erythropoiesis marker7.694.865.38x10–4
Kelextramedullary erythropoiesis marker29.04.618.04x10–4
Slc4a1extramedullary erythropoiesis marker151.44.905.05x10–4
Klf1extramedullary erythropoiesis marker20.65.063.95x10–4
Cldn13extramedullary erythropoiesis marker49.24.421.09x10–3
Trim10extramedullary erythropoiesis marker49.74.111.83x10–3
Eporextramedullary erythropoiesis marker12.34.915.01x10–4
Sptbextramedullary erythropoiesis marker36.55.213.16x10–4
Rhagextramedullary erythropoiesis marker33.04.647.68x10–4
Hba-a1erythroblasts63.86.465.23x10–5
Hbb-bserythroblasts59.26.604.34x10–5
Gata1BFU-E through erythroblasts7.694.865.38x10–4
Tfrc (CD71)CFU-E through erythroblasts3.126.197.58x10–5
KitProgenitors through CFU-E1.726.247.02x10–5
Sox6Adult definitive erythropoiesis35.53.560.0046
Aqp1Adult definitive erythropoiesis26.96.514.91x10–5
Nr3c1Adult definitive erythropoiesis1.122.800.018
Mki67Proliferation marker4.437.651.14x10–5
Cd34Multipotent progenitors through CFU-E1.601.590.141
Ly6a (Sca1)Multipotent progenitors1.171.020.331
Anpep (CD13)Multipotent progenitors1.161.140.28
Cd33Multipotent progenitors–1.06–0.370.71
Gata2Multipotent progenitors1.080.660.52
Xpo1Stability of nuclear Gata11.233.862.7x10–3
  1. *

    Gene enrichment in naïve Cd47−/− vs WT spleen cells depleted for CD4, CD11b, CD11c, CD19, CD45R (B220), CD49b (DX5), CD105, MHC Class II, Ter-119, and TCRγ/δ.

  2. Reported markers of stress-induced extramedullary erythropoiesis (Delic et al., 2020; Thompson et al., 2010).

Table 2
Differential mRNA expression of erythropoietic, stem cell, and proliferation associated markers in WT, Cd47−/−, and Thbs1−/− cells in clusters 12 and 14.
Cd47−/− vs WTThbs1−/− vs WT
ClusterGenep-valueAvg log2 FC p-valueAvg log2 FC
12Klf10.463–0.1150.0023–0.510
12Aqp10.331–0.1300.0011–0.472
12Tfrc0.0170.5310.930.009
14Tfrc0.640.0710.9220.016
12Epor0.506–0.0784.75x10–5–0.576
12Ermap5.1x10–30.3090.65–0.073
12Gata10.570.0070.0011–0.377
12Mki679.16x10–60.9970.220.456
14Mki670.00890.7980.0017–0.062
12Kit0.00150.3310.0570.264
14Kit0.0200.4090.580.222
12Xpo13.31x10–80.4720.00690.323
14Xpo10.0790.2110.190.148
12Ranbp10.240.1080.130.079
14Ranbp10.910.0690.076–0.486
12Ranbp21.98x10–140.7540.0920.269
14Ranbp20.00560.3650.880.042
12Nr3c18.5x10–40.3208.8x10–40.430
14Nr3c10.0120.1538.6x10–50.428
12Ddx463.04x10–80.5302.68x10–100.779
14Ddx460.00140.3850.00230.490
Table 3
Co-expression of the indicated erythropoiesis related genes in WT, Cd47−/− and Thbs1−/− spleen cells in cluster 12 was quantified and expressed as a percentage of the total cell number of each genotype in cluster 12.
Cluster 12 gene coexpressionWTCd47−/−Thbs1−/−
Cd34_Ermap0.0%1.8%1.6%
Cd34_Klf16.5%8.1%3.1%
Cd34_Gata16.5%8.4%0.0%
Mki67_Cd342.2%6.3%0.0%
Mki67_Ermap55.4%63.0%54.7%
Mki67_Kit48.9%67.2%57.8%
Mki67_Ly6a37.0%48.2%43.8%
Mki67_Klf162.0%74.6%57.8%
Mki67_Epor54.4%62.7%53.1%
Klf1_Ermap67.4%68.0%62.5%
Klf1_Gata178.3%82.4%59.4%
Ermap_Gata162.0%65.1%56.2%
Kit_Ermap54.4%56.7%59.4%
Kit_Klf171.7%77.1%78.1%
Ly6a_Ermap42.4%37.7%46.9%
Ly6a_Klf164.1%59.2%67.2%
Table 4
Differential expression of mitochondrial encoded genes in cluster 12 cells from WT, Cd47−/−, and Thbs1−/− spleens.
Genep-val Cd47−/−vs WTAvg log2FC Cd47−/−vs_WTp-val Thbs1-/- vs_WTAvg log2FC Thbs1-/- vs_WT
Mt-Atp60.887–0.0275.43x10–11–0.812
Mt-Atp87.23x10–200.980.01680.377
Mt-Co10.0280.1352.60x10–7–0.532
Mt-Co20.2260.0683.37x10–10–0.693
Mt-Co30.2070.0971.02x10–8–0.671
Mt-Cytb0.338–0.0069.35x10–10–0.790
Mt-Nd10.0740.156.78x10–6–0.636
Mt-Nd20.1350.1131.00x10–7–0.748
Mt-Nd37.67x10–90.5311.36x10–7–0.813
Mt-Nd40.1520.0935.84x10–6–0.558
Mt-Nd4l3.50x10–100.6586.49x10–80.787
Mt-Nd54.63x10–80.6120.01540.375
Mt-Nd61.65x10–70.4164.39x10–60.524
Table 5
Differential mRNA expression of erythropoietic, stem cell, and proliferation associated markers in reclustered WT, Cd47−/−, and Thbs1−/− erythroid and T cell clusters.
Cd47−/− vs WTThbs1−/− vs WT
ClusterGenep-valueAvg log2 FCp-valueAvg log2 FC
ErythroidKlf10.507–0.0950.0016–0.526
ErythroidAqp10.196–0.0888.1x10–4–0.390
T cellsAqp10.6620.0270.350–0.092
ErythroidTfrc0.0150.5490.9360.054
T cellsTfrc----
ErythroidEpor0.767–0.0460.026–0.219
ErythroidErmap0.00640.3260.496–0.010
T cellsErmap0.3350.0760.6970.063
ErythroidGata10.380.0090.0040–0.337
T cellsGata10.0044–.0.1620.3020.093
ErythroidMki674.6x10–61.0150.1180.545
T cellsMki67----
ErythroidKit0.00660.2930.0590.275
T cellsKit----
ErythroidXpo12.65x10–90.5140.00250.373
T cellsXpo10.2380.0470.1370.066
ErythroidRanbp10.0990.1330.410.078
T cellsRanbp10.278–0.1240.9700.019
ErythroidRanbp23.2x10–140.7550.0580.298
T cellsRanbp24.5x10–40.2750.2100.094
ErythroidNr3c16.5x10–40.3190.00110.447
T cellsNr3c10.00180.2740.00180.233
ErythroidDdx461.65x10–90.5302.98x10–90.775
T cellsDdx462.3x10–50.3487.5x10–60.446
ErythroidHba-a10.9060.0340.0088–2.168
Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Commercial assay or kitCD8a+T Cell Isolation KitMiltenyi BiotecCat#: 130–104–075
Commercial assay or kitCD8a (Ly-2) MicroBeads mouseMiltenyi BiotecCat#: 130-117-044
Chemical compound, drugEDTA solutionSigma-AldrichCat#: E80082 mM
Peptide, recombinant proteinbovine serum albumin (BSA)Sigma-AldrichCat#: A79060.5%
Chemical compound, drugACK Lysing Buffer, 100 mLQuality BiologicalsCat#: 118-156-721
Strain, strain background (Mus musculus, C57BL/6)WT miceJackson LaboratoriesWT C57BL/6
Strain, strain background (M. musculus, C57BL/6)Cd47-/- miceJackson LaboratoriesB6.129S7-Cd47tm1Fpl/J; Strain:003173Lindberg et al., 1996;
PMID:8864123
Strain, strain background (M. musculus, C57BL/6)Thbs1-/- miceJackson LaboratoriesB6.129S2-Thbs1tm1Hyn/J; Strain:006141Lawler et al., 1998;
PMID:9486968
AntibodyAnti-mouse Ter119-APC (Rat monoclonal)BiolegendCat#: 116212;
RRID:AB_313713
IgG2b
FACS (1 μg/100 μl/1 million cells)
Antibodyanti-mouse CD34-Percp/Cy5.5 (Rat monoclonal)BiolegendCat#: 119327;
RRID:AB_2728136
IgG2a
FACS (1 μg/100 μl/1 million cells)
Antibodyanti mouse Sca1 PE/Cy7 (Rat monoclonal)BiolegendCat#: 108114;
RRID:AB_493596
IgG2a
FACS (0.5 μg/100 μl/1 million cells)
Antibodyanti-mouse cKit PE (Rat monoclonal)BiolegendCat#: 105808;
RRID:AB_313217
IgG2b
FACS (0.1 μg/100 μl/1 million cells)
Antibodyanti-mouse Ki67- PE/Cy7 (Rat monoclonal)BiolegendCat#: 652425;
RRID:AB_2632693
IgG2a
FACS (0.5 μg/100 μl/1 million cells)
AntibodyIgG2b, κ APC Isotype Control Antibody (Rat monoclonal)BiolegendCat#: 400611;
RRID:AB_326555
IgG2b
FACS (1 μg/100 μl/1 million cells)
AntibodyIgG2a, κ PerCP/Cyanine5.5 Isotype Control Antibody (Rat monoclonal)BiolegendCat#: 400531;
RRID:AB_2864286
FACS (1 μg/100 μl/1 million cells)
AntibodyIgG2a, κ PE/Cyanine7 Isotype Control Antibody (Rat monoclonal)BiolegendCat#: 400521;
RRID:AB_326542
FACS (0.25 μg/100 μl/1 million cells)
AntibodyIgG2b, κ PE Isotype Control Antibody (Rat monoclonal)BiolegendCat#: 400608;
RRID:AB_326552
FACS (0.1 μg/100 μl/1 million cells)
AntibodyIgG2a, κ Alexa Fluor 488 Isotype Control Antibody (Rat monoclonal)BiolegendCat#: 400525;
RRID:AB_2864283
FACS (0.25 μg/100 μl/1 million cells)
Antibodyanti-Rabbit IgG (H+L) Cross-Adsorbed, Alexa Fluor 594 (Goat polyclonal)Thermo FisherCat#: A-11012FACS (0.2 μg/100 μl/1 million cells)
AntibodyAnti-ERMAP (Rabbit polyclonal)Thermo FisherCat#: BS-12333RFACS (1:100 dilution)
AntibodyAnti-GYPA (Rabbit polyclonal)Thermo FisherCat#: BS-2575RFACS (1:100 dilution)
AntibodyAnti-Aquaporin 1 (Rabbit polyclonal)Thermo FisherCat#: PA5-78806FACS (1 μg/100 μl/1 million cells)
AntibodyAnti-EPOR (Rabbit polyclonal)Bioss; Thermo FisherCat#: BS-1424RFACS (1 μg/100 μl/1 million cells)
Author response table 1
ExtramedullaryCD47KO/WT_FCCD47KO/WT_pvalDEG rank
erythropoiesis Gene
Ermap3.2029418680.0359388562087
Gypa23.092892990.000169485163
Gata12.8091056640.0448971972273
Slc4a14.3528122480.0057327231641
Klf130.095438120.00044213976
Trim1010.89360420.003179996672
Sptb4.4430918770.000325321607
Rhag25.908694368.45 E-05121

Additional files

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)

  1. Rajdeep Banerjee
  2. Thomas J Meyer
  3. Margaret C Cam
  4. Sukhbir Kaur
  5. David D Roberts
(2024)
Differential regulation by CD47 and thrombospondin-1 of extramedullary erythropoiesis in mouse spleen
eLife 12:RP92679.
https://doi.org/10.7554/eLife.92679.3