Dehao Huang1,2, Qianhao Zhao3, Mengyun Zhang4, Qitong Weng1,2, Qi Zhang1,2, Kaitao Wang5, Fang Dong6,7, Hui Cheng6,7, Fangxiao Hu1, Jinyong Wang1,2,8. 1. CAS Key Laboratory of Regenerative Biology, Guangzhou Institutes of Biomedicine and Health, Chinese Academy of Sciences, Guangzhou, China. 2. University of Chinese Academy of Sciences, Beijing, China. 3. Faculty of Forensic Medicine, Zhongshan School of Medicine, Sun Yat-Sen University, China. 4. GMU-GIBH Joint School of Life Sciences, Guangzhou Medical University, Guangzhou, China. 5. School of Biomedical Engineering, Sun Yat-Sen University, Shenzhen, China. 6. State Key Laboratory of Experimental Hematology & National Clinical Research Center for Blood Diseases, Institute of Hematology & Blood Diseases Hospital, Chinese Academy of Medical Sciences & Peking Union Medical College, Tianjin, China. 7. Center for Stem Cell Medicine & Department of Stem Cell and Regenerative Medicine, Chinese Academy of Medical Sciences & Peking Union Medical College, Tianjin, China. 8. State Key Laboratory of Stem Cell and Reproductive Biology, Institute for Stem Cell and Regeneration, Institute of Zoology, Chinese Academy of Sciences, Beijing, China.
Abstract
OBJECTIVES: The expression of transcription factor Hoxb5 specifically marks the functional haematopoietic stem cells (HSC) in mice. However, our recent work demonstrated that ectopic expression of Hoxb5 exerted little effect on HSC but could convert B-cell progenitors into functional T cells in vivo. Thus, cell type- and development stage-specific roles of Hoxb5 in haematopoietic hierarchy await more extensive exploration. In this study, we aim to investigate the effect of Hoxb5 expression in multipotent blood progenitor cells. MATERIALS AND METHODS: A Mx1cre/RosaLSL-Hoxb5-EGFP/+ mouse model was used to evaluate the effect of Hoxb5 expression in blood multipotent progenitor cells (MPP). Golden standard serial transplantation experiments were used to test the long-term haematopoiesis potential of Hoxb5-expressing MPP. Single-cell RNA-seq analysis was used to characterize the gained molecular features of Hoxb5-expressing MPP and to compare with the global transcriptome features of natural adult HSC and fetal liver HSC (FL HSC). RESULTS: Here, with a mouse strain engineered with conditional expression of Hoxb5, we unveiled that induced expression of Hoxb5 in MPP led to the generation of a de novo Sca1+ c-kit+ CD11b+ CD48+ (CD11b+ CD48+ SK) cell type, which can repopulate long-term multilineage haematopoiesis in serial transplantations. RNA-seq analysis showed that CD11b+ CD48+ SK cells exhibited acquired features of DNA replication and cell division. CONCLUSIONS: Our current study uncovers that Hoxb5 can empower MPP with self-renewal ability and indicates an alternative approach for generating HSC-like cells in vivo from blood lineage cells.
OBJECTIVES: The expression of transcription factor Hoxb5 specifically marks the functional haematopoietic stem cells (HSC) in mice. However, our recent work demonstrated that ectopic expression of Hoxb5 exerted little effect on HSC but could convert B-cell progenitors into functional T cells in vivo. Thus, cell type- and development stage-specific roles of Hoxb5 in haematopoietic hierarchy await more extensive exploration. In this study, we aim to investigate the effect of Hoxb5 expression in multipotent blood progenitor cells. MATERIALS AND METHODS: A Mx1cre/RosaLSL-Hoxb5-EGFP/+ mouse model was used to evaluate the effect of Hoxb5 expression in blood multipotent progenitor cells (MPP). Golden standard serial transplantation experiments were used to test the long-term haematopoiesis potential of Hoxb5-expressing MPP. Single-cell RNA-seq analysis was used to characterize the gained molecular features of Hoxb5-expressing MPP and to compare with the global transcriptome features of natural adult HSC and fetal liver HSC (FL HSC). RESULTS: Here, with a mouse strain engineered with conditional expression of Hoxb5, we unveiled that induced expression of Hoxb5 in MPP led to the generation of a de novo Sca1+ c-kit+ CD11b+ CD48+ (CD11b+ CD48+ SK) cell type, which can repopulate long-term multilineage haematopoiesis in serial transplantations. RNA-seq analysis showed that CD11b+ CD48+ SK cells exhibited acquired features of DNA replication and cell division. CONCLUSIONS: Our current study uncovers that Hoxb5 can empower MPP with self-renewal ability and indicates an alternative approach for generating HSC-like cells in vivo from blood lineage cells.
Haematopoietic stem cell (HSC) is the blood cell type that possesses dual features of self‐renewal and multilineage potential, which are critical for replenishing the entire haematopoietic system throughout an individual lifespan.
,
However, the absolute numbers of HSC in adults are extremely rare
,
and are not efficiently expanded in vitro.
,
Researchers have been attempting alternative approaches to generate engraftable blood progenitors by enforcing expressing those molecules highly expressed in HSC but absent in downstream progenies. Ectopic expression of Sox17 can confer self‐renewal potential on adult haematopoietic progenitors. However, this approach eventually led to leukemogenesis.
Likewise, miR‐125a is a non‐coding RNA gene preferentially expressed in HSC rather than blood progenies.
Ectopic expression of miR‐125a in mouse haematopoietic progenitors induced long‐term haematopoiesis, but the recipient mice suffered an MPN‐like disease after secondary transplantation.
,
,
Therefore, more extensive and innovative efforts are needed to develop safer approaches to convert blood progenitor cells into engraftable blood stem cells for ultimately therapeutic uses.Hoxb5, a member of HOX gene family, is preferentially expressed in HSC and uniquely marks the long‐term HSC.
,
Our recent study showed that the gain of function of Hoxb5 in pro‐pre‐B cells reprogrammed these cells into T lymphocytes in vivo.
Moreover, the latest research shows that exogeneous Hoxb5 expression confers protection against loss of self‐renewal to Hoxb5‐negative HSCs and can partially alter the cell fate of ST‐HSCs to that of LT‐HSCs.
Here, we further studied the potential role of Hoxb5 in the MPP cell context, an intermediate progeny of HSC without self‐renewal ability. Interestingly, conditional overexpression of Hoxb5 in MPP upon transplantation led to long‐term haematopoiesis in serially transplanted mice. More importantly, Hoxb5 resulted in a de novo cell type defined as CD11b+CD48+SK, which contributed to the sustainable long‐term haematopoiesis in serially transplanted recipients. CD11b+CD48+SK cells exhibited features related to DNA replication and cell division. This study reveals de novo evidence that Hoxb5 can efficiently reprogramme blood progenitors into engraftable blood stem cells.
MATERIALS AND METHODS
Mice
Animals were housed in the animal facility of the Guangzhou Institutes of Biomedicine and Health (GIBH). RosaLSL‐Hoxb5‐EGFP/+ mice were described as previously reported.
Mice of the CD45.1+ and Mx1‐cre strains were purchased from the Jackson laboratory. All the mouse lines were maintained on a pure C57BL/6 genetic background. All experiments were conducted in accordance with experimental protocols approved by the Animal Ethics Committee of GIBH.
Flow cytometry
Antibodies to CD2 (RM2‐5), CD3 (145‐2C11), CD4 (RM4‐5), CD8a (53–6.7), Gr1 (RB6‐8C5), CD11b (M1/70), Ter119 (TER‐119), B220 (6B2), c‐kit (2B8), Sca‐1 (E13‐161.7), CD135 (A2F10), CD150 (TC15‐12F12.2), CD19 (eBio1D3), CD48 (HM48‐1) ki‐67 (16A8), Fcγ RII/III (2.4G2), CD127 (SB/199), CD45.2 (104) and CD45.1(A20) were purchased from eBioscience or BioLegend. DAPI, 7‐AAD or PI was used to stain dead cells. Flow cytometry was performed on an LSR Fortessa (BD Biosciences) and data were processed by FlowJo software (Version: 10.4.0, Tree Star).
Cell sorting
Bone marrow cells used for transplantation or RNA‐seq were first incubated with the biotin‐conjugated antibody to Sca1 (anti‐Sca1 biotin) and then enriched using Anti‐Biotin MicroBeads by autoMACS Pro (Miltenyi Biotec). The enriched cells, stained with antibodies, were sorted by BD FACSAria III.
Transplantation
All recipients (CD45.1+, C57BL/6 background) were lethally irradiated (9 Gy, RS2000, Rad Source) at least 4 h, but less than 24 h before transplantation. MPP (400 cells/mouse) from RosaLSL‐Hoxb5‐EGFP/+ mouse or Mx1cre/RosaLSL‐Hoxb5‐EGFP/+ mouse for primary transplantation and donor‐derived CD48+CD11b+SK cells (2000 cells/mouse) from the primary recipients for secondary transplantation were retro‐orbitally transplanted into the recipients with the enriched Sca1− helper cells (0.25 million/mouse, CD45.1+). For third transplantation, total BM cells (10 million/mouse) from the secondary recipients were retro‐orbitally transplanted into the third recipients. To induce Hoxb5 expression, the primary recipients were intraperitoneally injected with polyinosinic‐polycytidylic acid (pIpC) (250 μg/mouse) every other day for six times starting from the 5th day before transplantation, until the 5th day after transplantation. Recipients were fed with the water added with trimethoprim‐sulfamethoxazole for 1 month after irradiation.
RNA‐seq and data analysis
cDNA of the single cell from adult wild‐type HSC (BM HSC, RosaLSL‐Hoxb5‐EGFP/+ mice, 8‐week old), fetal liver HSC (FL HSC, RosaLSL‐Hoxb5‐EGFP/+ mice, Day14.5, defined as CD45.2+Lin−Sca1+c‐kit+CD11b+CD150+), wild‐type MPP (WT‐MPP, RosaLSL‐Hoxb5‐EGFP/+ mice, 8‐week old) and CD48+CD11b+SK cells (primary recipients, Week 8 post‐transplantation) were extracted and amplified using Discover‐sc WTA Kit V2 (Vazyme). The expression level of B2m and Gapdh was used to assess the quality of the amplified cDNA by qPCR analysis. These qualified samples were used to prepare the sequencing library by the TruePrep DNA Library Prep Kit V2 (Vazyme), and the qualified libraries were sequenced by Illumina sequencer NextSeq 500. Raw data (FASTQ files) were generated using bcl2fastq software (version 2.16.0.10) and were uploaded to the Gene Expression Omnibus public database (GSE 183800). HISAT2 (version 2.1.0) was used to align the raw data, and the StringTie (version 1.3.4) was used to estimate the expression level in the transcripts per million (TPM) as previously reported.
,
Gene set enrichment analysis (GSEA) and gene‐ontology (GO)‐enrichment analysis (clusterProfiler package) were performed as described.
,
Spearman correlation coefficient between population was used for correlation analysis. Correlation analysis was performed by cor() function of R.
RESULTS
Enforced expression of Hoxb5 in MPP leads to long‐term haematopoiesis in transplantation setting
To evaluate the potential role of Hoxb5 in MPP, we chose Mx1cre/RosaLSL‐Hoxb5‐EGFP/+ mouse model in which ectopic Hoxb5 was inserted into Rosa26 in the form of Loxp‐Stop‐Loxp‐Hoxb5‐EGFP cassette as previously reported (Figure S1A).
Thus, the ectopic Hoxb5 can be conditionally induced in the presence of Cre recombinase, which is under the control of the Mx1 promoter, and can be induced by polyinosinic‐polycytidylic acid (pIpC). Upon transplantation into recipient animals, ectopic Hoxb5 expression would be turned on by injection with pIpC, and the GFP signal reports the expression of the ectopic Hoxb5 at single‐cell resolution. We sorted the conventional MPP (CD45.2+GFP−Lin [CD2, CD3, CD4, CD8, CD11b, Gr1, B220, Ter119]− CD48−Sca1+c‐kit+CD150−CD135+)
,
from the Sca1+ enriched bone marrow cells of Mx1cre/RosaLSL‐Hoxb5‐EGFP/+ mouse or RosaLSL‐Hoxb5‐EGFP/+ mouse. Four hundred sorted MPP along with 0.25 million Sca1− helper cells (CD45.1+) were retro‐orbitally transplanted into irradiated individual recipients (CD45.1+, C57BL/6 background) (Figure 1A,B). The recipients were intraperitoneally injected with pIpC (250 μg/mouse) every other day for six times starting from day 5 before transplantation. We assessed the reconstitution ability of the donor‐derived cells by analysing the peripheral blood (PB) chimeras every 4 weeks until Week 20 post‐transplantation. (Figure 1A). Amazingly, in the primary recipients transplanted with the MPP of the Mx1cre/RosaLSL‐Hoxb5‐EGFP/+ mouse, the ratio of the donor‐derived cells (CD45.2+GFP+) continuously increased and the minimum ratio was up to 62% at Week 20 post‐transplantation, whereas the control recipients transplanted with the RosaLSL‐Hoxb5‐EGFP/+ MPP showed a significantly low reconstitution ability, and the maximum donor‐derived cell (CD45.2+) ratio was 11.4% at the Week 20 post‐transplantation (Figure 1C). Furthermore, the contributions of donor‐derived cells in the spleen (SP) and bone marrow (BM) tissues of the Hoxb5‐expressing MPP recipients were significantly more than the control group (p < 0.001) (Figure S1B). In addition, multiple blood lineages including T cells (CD3+), B cells (CD19+) and myeloid (CD11b+or Gr1+) in the PB, SP and BM were also detected at Week 20 post‐transplantation in the primary recipients (Figure 1D). These results demonstrate that enforced expression of Hoxb5 in MPP leads to long‐term haematopoiesis.
FIGURE 1
Overexpression of Hoxb5 empowers long‐term reconstitution capacity on MPP. (A) The scheme for MPP transplantation. (B) Gating strategy for sorting the MPP. MPP were defined as CD45.2+GFP−Lin (CD2−CD3−CD4−CD8−CD11b−Gr1−B220−Ter119−) CD48−Sca1+c‐kit+CD150−CD135+ from either the Mx1cre/RosaLSL‐Hoxb5‐EGFP/+ mice or RosaLSL‐Hoxb5‐EGFP/+ control mice (8‐week old). The sorted MPP were transplanted into lethally irradiated (9.0 Gy) recipients (CD45.1+, C57BL/6 background, 400 cells/mouse) with Sca1− helper cells (CD45.1+, 0.25 million/mouse). Recipients were injected with pIpC (i.p. 250 μg/mouse) every other day for six times starting from day 5 before transplantation. (C) Contribution curves of the donor‐derived cells in peripheral blood cells (PB) of the primary recipients. The donor cells were defined as CD45.2+GFP+ (Mx1cre/RosaLSL‐Hoxb5‐EGFP/+ MPP recipients) or CD45.2+ (RosaLSL‐Hoxb5‐EGFP/+ recipients). The PB of the recipients transplanted with Mx1cre/RosaLSL‐Hoxb5‐EGFP/+ MPP (n = 6, as indicated by the red dot) or RosaLSL‐Hoxb5‐EGFP/+ MPP (n = 6, as indicated by the blue dot) was analysed every 4 weeks until the Week 20 after transplantation. Mean ± SD, ***p < 0.001 independent samples t‐test. (D) Lineage distribution of the recipients (n = 3) at Week 20 after Mx1cre/RosaLSL‐Hoxb5‐EGFP/+ MPP transplantation. Columns shown are percentages of donor‐derived T cells (CD3+), B cells (CD19+) and myeloid cells (CD11b+ or Gr1+) in PB, spleen (SP) and bone marrow (BM)
Overexpression of Hoxb5 empowers long‐term reconstitution capacity on MPP. (A) The scheme for MPP transplantation. (B) Gating strategy for sorting the MPP. MPP were defined as CD45.2+GFP−Lin (CD2−CD3−CD4−CD8−CD11b−Gr1−B220−Ter119−) CD48−Sca1+c‐kit+CD150−CD135+ from either the Mx1cre/RosaLSL‐Hoxb5‐EGFP/+ mice or RosaLSL‐Hoxb5‐EGFP/+ control mice (8‐week old). The sorted MPP were transplanted into lethally irradiated (9.0 Gy) recipients (CD45.1+, C57BL/6 background, 400 cells/mouse) with Sca1− helper cells (CD45.1+, 0.25 million/mouse). Recipients were injected with pIpC (i.p. 250 μg/mouse) every other day for six times starting from day 5 before transplantation. (C) Contribution curves of the donor‐derived cells in peripheral blood cells (PB) of the primary recipients. The donor cells were defined as CD45.2+GFP+ (Mx1cre/RosaLSL‐Hoxb5‐EGFP/+ MPP recipients) or CD45.2+ (RosaLSL‐Hoxb5‐EGFP/+ recipients). The PB of the recipients transplanted with Mx1cre/RosaLSL‐Hoxb5‐EGFP/+ MPP (n = 6, as indicated by the red dot) or RosaLSL‐Hoxb5‐EGFP/+ MPP (n = 6, as indicated by the blue dot) was analysed every 4 weeks until the Week 20 after transplantation. Mean ± SD, ***p < 0.001 independent samples t‐test. (D) Lineage distribution of the recipients (n = 3) at Week 20 after Mx1cre/RosaLSL‐Hoxb5‐EGFP/+ MPP transplantation. Columns shown are percentages of donor‐derived T cells (CD3+), B cells (CD19+) and myeloid cells (CD11b+ or Gr1+) in PB, spleen (SP) and bone marrow (BM)
Hoxb5 results in the occurrence of a de novo
CD11b
+
CD48
+
SK cell type associated with the long‐term engraftable feature
To investigate the cellular mechanism, we analysed the blood progenitor cells in the primary recipients at Week 20 post‐transplantation. We discovered a de novo donor‐derived Sca1+c‐kit+ population cells, which simultaneously expressed CD11b and CD48 surface markers. Certainly, this cell type is not identified in natural blood cells in the absence of Hoxb5 expression (Figure 2A). Consistent with previous reports, natural MPP transplantation cannot sustainably give rise to LSK cells in the bone marrow of recipient mice (Figure S1C). To further test whether the CD11b+CD48+SK cells are responsible for the long‐term repopulating feature in Hoxb5 expressing MPP, we sorted the GFP+CD11b+CD48+SK cells and transplanted them into secondary recipient mice (CD45.1+, C57BL/6 background, 2000 cells/mouse) with Sca1− helper cells (CD45.1+, 0.25 million/mouse). As expected, these CD11b+CD48+SK cells successfully reconstituted multilineage haematopoiesis in secondary recipients, as demonstrated by stable increases of donor‐derived cells (GFP+) in the PB after transplantation (Figure 2B). Of note, the donor chimeras achieved as high as 94.6% at Week 20 after transplantation and lineages of T, B and myeloid cells can be detected at Week 4, Week 12 and Week 20 post‐transplantation (Figure 2B,C). Moreover, the donor‐derived T, B and myeloid lineages in the PB, SP and BM also exhibited patterns resembling natural haematopoiesis at the Week 24 post‐transplantation (Figure 2D). More importantly, the donor‐derived CD11b+CD48+SK cells can still be detected in the BM of the secondary recipients (Figure 2E). These results indicate that the de novo CD11b+CD48+SK cell type is engraftable in the secondary recipients.
FIGURE 2
A de novo CD11b+CD48+SK cell population reconstitutes haematopoiesis in secondary recipients. (A) FACS analysis of the donor‐derived BM progenitors in the primary recipients at Week 20 post‐transplantation. Antibodies of lineages (CD2−CD3−CD4−CD8−Gr1−B220−Ter119−) (Lin), Sca1, c‐kit, CD11b and CD48 were stained the BM of Mx1cre/RosaLSL‐Hoxb5‐EGFP/+ MPP recipients. CD11b+CD48+SK cells were defined as CD45.2+GFP+Lin−Sca1+c‐kit+CD48+CD11b+ were sorted from the primary recipients at Week 20 post‐transplantation and retro‐orbitally transplanted into the lethally irradiated secondary recipients (CD45.1+, C57BL/6 background, 2000 cells/mouse). (B) Chimeras curves of the donor cells to the peripheral blood (PB) cells of the secondary recipients (n = 7, as indicated by the red dot). For the secondary transplantation, CD11b+CD48+SK cells were retro‐orbitally injected into the lethally irradiated recipients (9.0 Gy, 2000 cells/mouse). The donor‐derived cells (CD45.2+GFP+) in the PB were analysed every 4 weeks post‐transplantation. (C) Lineage distribution in PB of the secondary recipients (n = 4) at weeks 4, 12 and 20 post‐transplantation. Proportions of the CD3+ (T), CD19+ (B) and CD11b+ (myeloid) in donor‐derived cells were analysed. Each symbol represents an individual host mouse. (D) Lineage distribution of the recipients (n = 3) at Week 24 after CD11b+CD48+SK cell transplantation. Columns shown are percentages of donor‐derived T cells (CD3+), B cells (CD19+) and myeloid cells (CD11b+ or Gr1+) in PB, spleen (SP) and bone marrow (BM). (E) Immunophenotypes of the donor‐derived CD11b+CD48+SK cells in the bone marrow (BM) of the secondary recipients (two representative mice)
A de novo CD11b+CD48+SK cell population reconstitutes haematopoiesis in secondary recipients. (A) FACS analysis of the donor‐derived BM progenitors in the primary recipients at Week 20 post‐transplantation. Antibodies of lineages (CD2−CD3−CD4−CD8−Gr1−B220−Ter119−) (Lin), Sca1, c‐kit, CD11b and CD48 were stained the BM of Mx1cre/RosaLSL‐Hoxb5‐EGFP/+ MPP recipients. CD11b+CD48+SK cells were defined as CD45.2+GFP+Lin−Sca1+c‐kit+CD48+CD11b+ were sorted from the primary recipients at Week 20 post‐transplantation and retro‐orbitally transplanted into the lethally irradiated secondary recipients (CD45.1+, C57BL/6 background, 2000 cells/mouse). (B) Chimeras curves of the donor cells to the peripheral blood (PB) cells of the secondary recipients (n = 7, as indicated by the red dot). For the secondary transplantation, CD11b+CD48+SK cells were retro‐orbitally injected into the lethally irradiated recipients (9.0 Gy, 2000 cells/mouse). The donor‐derived cells (CD45.2+GFP+) in the PB were analysed every 4 weeks post‐transplantation. (C) Lineage distribution in PB of the secondary recipients (n = 4) at weeks 4, 12 and 20 post‐transplantation. Proportions of the CD3+ (T), CD19+ (B) and CD11b+ (myeloid) in donor‐derived cells were analysed. Each symbol represents an individual host mouse. (D) Lineage distribution of the recipients (n = 3) at Week 24 after CD11b+CD48+SK cell transplantation. Columns shown are percentages of donor‐derived T cells (CD3+), B cells (CD19+) and myeloid cells (CD11b+ or Gr1+) in PB, spleen (SP) and bone marrow (BM). (E) Immunophenotypes of the donor‐derived CD11b+CD48+SK cells in the bone marrow (BM) of the secondary recipients (two representative mice)To assess the long‐term haematopoiesis of the CD11b+CD48+SK cells, we performed a third round of transplantation using the total BM cells of the CD11b+CD48+SK cells secondary recipients (10 million/mouse, n = 6). The contribution of CD45.2+GFP+ donor cells to peripheral blood was measured at Weeks 8, 16, 20, 26 and 32 after transplantation. All the recipients were reconstituted with the CD45.2+GFP+ cells with a ratio range of 48.7%–74.2% at Week 8 post‐transplantation, and the average ratio was still as high as 49% at Week 32 after transplantation (Figure 3A). Moreover, the donor‐derived cell ratio has no significant difference (p = .075) at Week 32 compared with Week 8 post‐transplantation (Figure 3B). Meanwhile, the donor‐derived cells also showed multilineage distributions in PB at Weeks 8 and 20 post‐transplantation (Figure 3C). Furthermore, the average ratios of T cells (CD3+) were 7.1% (Week 8, n = 6), 17.4% (Week 20, n = 6) and 16.1% (Week 32, n = 6). The average ratios of B cells (CD19+) were 83.2% (Week 8, n = 6), 76.5% (Week 20, n = 6) and 82.8% (Week 32, n = 6). As for myeloid cells (CD11b+ or Gr1+), the average ratios were 11.1% (Week 8, n = 6), 10.6% (Week 20, n = 6) and 7.7% (Week 32, n = 6) post‐transplantation in PB, respectively (Figure 3D). Collectively, these results indicate that the CD11b+CD48+SK cells can sustain a long‐term haematopoiesis in serial transplantation settings.
FIGURE 3
CD11b+CD48+SK cells still maintained the repopulation capacity in the third transplantation recipients. (A) Chimeric curves of the donor cells in PB of the secondary recipients (n = 6, as indicated by the red dot). Donor‐derived cells (CD45.2+GFP+) in PB were analysed at weeks 8, 16, 20, 26 and 32 post‐transplantation. For the third transplantation, recipients (CD45.1+, C57BL/6 background) were lethally irradiated and then were retro‐orbitally injected with the nucleated BM cells (10 million/mouse) isolated from the secondary recipients. (B) Comparison of the donor‐derived cells (CD45.2+GFP+) at Weeks 8 and 32 post‐transplantation. Mean ± SD, no significant difference, p = 0.075, independent samples t‐test. (C) Representative FACS analysis (n = 3, as indicated by the red dot) of the PB from the third transplantation recipients (3rd) after transplanting with the total BM cells of the secondary recipients at Weeks 8 and 20 post‐transplantation. (D) Lineage distribution in PB of the third recipients (n = 6) at Weeks 8, 20 and 32 post‐transplantation. Proportions of CD3+ (T), CD19+ (B), CD11b+or Gr1+ (myeloid) in donor‐derived cells were analysed. Each symbol represents an individual host mouse
CD11b+CD48+SK cells still maintained the repopulation capacity in the third transplantation recipients. (A) Chimeric curves of the donor cells in PB of the secondary recipients (n = 6, as indicated by the red dot). Donor‐derived cells (CD45.2+GFP+) in PB were analysed at weeks 8, 16, 20, 26 and 32 post‐transplantation. For the third transplantation, recipients (CD45.1+, C57BL/6 background) were lethally irradiated and then were retro‐orbitally injected with the nucleated BM cells (10 million/mouse) isolated from the secondary recipients. (B) Comparison of the donor‐derived cells (CD45.2+GFP+) at Weeks 8 and 32 post‐transplantation. Mean ± SD, no significant difference, p = 0.075, independent samples t‐test. (C) Representative FACS analysis (n = 3, as indicated by the red dot) of the PB from the third transplantation recipients (3rd) after transplanting with the total BM cells of the secondary recipients at Weeks 8 and 20 post‐transplantation. (D) Lineage distribution in PB of the third recipients (n = 6) at Weeks 8, 20 and 32 post‐transplantation. Proportions of CD3+ (T), CD19+ (B), CD11b+or Gr1+ (myeloid) in donor‐derived cells were analysed. Each symbol represents an individual host mouse
Characterization of CD11b
+
CD48
+
SK cells at transcriptome level
To investigate the underlying molecular mechanisms, we characterized the CD11b+CD48+SK cells (n = 47) at transcriptome levels by single‐cell RNA‐seq analysis. Meanwhile, we also performed single‐cell RNA‐seq of the BM HSC (n = 36, RosaLSL‐Hoxb5‐EGFP/+, 8‐week old), WT‐MPP (n = 42, RosaLSL‐Hoxb5‐EGFP/+ mice, 8‐week old) and FL HSC (n = 56, RosaLSL‐Hoxb5‐EGFP/+ mice, Day 14.5 embryo) (Figure S1D) for comparisons. To unbiasedly cluster the four HSPC populations, UMAP analysis was performed and the results illustrated that CD11b+CD48+SK cells were artificially unique and distinct from natural FL HSC, BM HSC and WT‐MPP (Figure 4A). The co‐efficiency value between CD11b+CD48+SK cells and WT‐MPP is 0.988. And the one between CD11b+CD48+SK cells and FL HSC is 0.987. CD11b+CD48+SK cells are even closer to WT‐MPP rather than FL HSC (Figure 4B). To further investigate the underlying mechanisms of Hoxb5‐induced CD11b+CD48+SK cells in long‐term transplantation, we performed differential gene expression analysis of CD11b+CD48+SK cells with other three populations (FL HSC, BM HSC and WT‐MPP) (Figure 4C), and gene‐ontology (GO) analysis showed that the upregulated genes (adjusted p value <0.05) in CD11b+CD48+SK cells were related to translation‐regulation pathways, T‐cell activation and DNA replication pathways (Figure 4D). To further dissect the transcriptome signatures of CD11b+CD48+SK cells, we enriched the differentially expressed genes (adjusted p value <0.05) between FL HSC and WT‐MPP. The up‐ and downregulated differential expressed genes enriched above were respectively used as gene sets for Gene set enrichment analysis (GSEA) of CD11b+CD48+SK cells and WT‐MPP (Figure 4E,F). The results showed that the upregulated genes were enriched in the CD11b+CD48+SK cells (Figure 4E). Meanwhile, we also enriched the up‐ and downregulated differential expressed genes (adjusted p value <0.05) between FL HSC and BM HSC for GSEA between CD11b+CD48+SK cells and BM‐HSC. The GSEA suggested that the upregulated genes in FL HSC were also enriched in the CD11b+CD48+SK cells (Figure 4G,H). Moreover, we combined the leading edge genes upregulated in FL HSC compared with WT‐MPP or BM HSC and performed heatmap analysis for the four populations. The results showed that the expression levels of these genes in CD11b+CD48+SK cells were equivalent to those in FL HSC (Figure S1E). We further performed GO analysis using the leading edge genes and observed that they were involved in cell proliferation processes of DNA replication and chromosome segregation (Figure S1F). Besides the higher expression of Hoxb5 in CD11b+CD48+SK cells and FL HSC when compared with BM HSC and WT‐MPP, several genes regulating the cell cycle and haematopoiesis were also upregulated, including Birc5,
Gmnn,
Cdt1,
Cdc45
and Gins1
(Figure S1G). Furthermore, the essential genes in regulating HSC homeostasis or development, including Cdk6,
Satb1,
,
Runx3,
and Mybl2
,
, were exclusively upregulated in CD11b+CD48+SK cells (Figure S1H). Thus, despite CD11b+CD48+SK cells demonstrating coefficiency with natural MPP cells, they also share some common features with FL HSC, which might account for their acquired self‐renewal ability.
FIGURE 4
Characterization of CD11b+CD48+SK cells at single‐cell resolution. (A) UMAP analysis of the WT‐MPP (n = 42), BM HSC (n = 36), FL HSC (n = 56) and CD11b+CD48+SK cells (n = 47). Each colour represents one of the four cell populations. (B) Spearman correlation analysis of the WT‐MPP, BM HSC, FL HSC and CD11b+CD48+SK cells. (C) Heatmap analysis of CD11b+CD48+SK cells and other three populations (BM HSC, WT‐MPP and FL HSC). Upregulated differential expressed genes of CD11b+CD48+SK cells were used for plotting. (D) Gene ontology (GO) enrichment analysis of upregulated DEGs (adjusted p value <0.05) in (C) of CD11b+CD48+SK cells. (E) Gene set enrichment analysis of WT‐MPP (n = 42) and CD11b+CD48+SK cells (n = 47). The gene set used for analysis was from the upregulated genes in FL HSC (n = 56) versus WT‐MPP (n = 42) (adjusted p value <0.05). (F) Gene set enrichment analysis of WT‐MPP (n = 42) and CD11b+CD48+SK cells (n = 47). The gene set used for analysis was from the downregulated genes in FL HSC (n = 56) versus WT‐MPP (n = 42) (adjusted p value <0.05). (G) Gene set enrichment analysis of WT‐MPP (n = 42) and CD11b+CD48+SK cells (n = 47). The gene set used for analysis was from the upregulated genes in FL HSC (n = 56) versus BM HSC (n = 36) (adjusted p value <0.05). (H) Gene set enrichment analysis of WT‐MPP (n = 42) and CD11b+CD48+SK cells (n = 47). The gene set used for analysis was from the downregulated genes in FL HSC (n = 56) versus BM HSC (n = 36) (adjusted p value <0.05)
Characterization of CD11b+CD48+SK cells at single‐cell resolution. (A) UMAP analysis of the WT‐MPP (n = 42), BM HSC (n = 36), FL HSC (n = 56) and CD11b+CD48+SK cells (n = 47). Each colour represents one of the four cell populations. (B) Spearman correlation analysis of the WT‐MPP, BM HSC, FL HSC and CD11b+CD48+SK cells. (C) Heatmap analysis of CD11b+CD48+SK cells and other three populations (BM HSC, WT‐MPP and FL HSC). Upregulated differential expressed genes of CD11b+CD48+SK cells were used for plotting. (D) Gene ontology (GO) enrichment analysis of upregulated DEGs (adjusted p value <0.05) in (C) of CD11b+CD48+SK cells. (E) Gene set enrichment analysis of WT‐MPP (n = 42) and CD11b+CD48+SK cells (n = 47). The gene set used for analysis was from the upregulated genes in FL HSC (n = 56) versus WT‐MPP (n = 42) (adjusted p value <0.05). (F) Gene set enrichment analysis of WT‐MPP (n = 42) and CD11b+CD48+SK cells (n = 47). The gene set used for analysis was from the downregulated genes in FL HSC (n = 56) versus WT‐MPP (n = 42) (adjusted p value <0.05). (G) Gene set enrichment analysis of WT‐MPP (n = 42) and CD11b+CD48+SK cells (n = 47). The gene set used for analysis was from the upregulated genes in FL HSC (n = 56) versus BM HSC (n = 36) (adjusted p value <0.05). (H) Gene set enrichment analysis of WT‐MPP (n = 42) and CD11b+CD48+SK cells (n = 47). The gene set used for analysis was from the downregulated genes in FL HSC (n = 56) versus BM HSC (n = 36) (adjusted p value <0.05)
DISCUSSION
In this study, we explored the role of Hoxb5 in the MPP cell context. Hoxb5 expression leads to long‐term haematopoiesis of MPP in serial transplantation settings by generating a population of phenotypic CD11b+CD48+ SK cells. The Hoxb5‐induced CD11b+CD48+ SK cells shared mixed features of natural MPP and FL HSC, especially the molecular signatures of cell division and proliferation. The stemness feature is the only functional difference between HSC and their progeny MPP. However, the stemness‐losing mechanism along the differentiation path from HSC to MPP is unknown. The genes which are shut down or downregulated in MPP, such as Hoxb5, might be accountable for the loss of stemness from HSC to MPP. Moreover, this maybe indicated by the latest research, which reported that exogeneous Hoxb5 expression confers protection against loss of self‐renewal to Hoxb5‐negative HSC and can partially alter the cell fate of ST‐HSC to that of LT‐HSC.
Despite overexpressing Hoxb5 in MPP generated no phenotypic HSC, the de novo CD11b+CD48+SK cells can sustain long‐term engraftment with full blood lineage differentiation potential. Interestingly, FL HSC shared two features with Hoxb5‐expressing CD11b+CD48+SK cells, demonstrated by rapid proliferation and expressing CD11b surface marker.
,
But BM HSC lose the expression of CD11b,
which is consistent with their predominant dormancy under homeostasis. Therefore, the expression of CD11b is phenotypically associated with the fast‐expanding features of FL HSC and Hoxb5‐expressing CD11b+CD48+SK cells. Seemingly, the enforced expression of Hoxb5 in MPP activates a cell division machinery
,
without compromising their multilineage differentiation potential. However, the CD11b+CD48+SK cells have a differentiation bias towards lymphopoiesis, especially the B lymphopoiesis, which was consistent with the results of Hoxb5‐expressing total BM transplantation from the Vavcre/RosaLSL‐Hoxb5‐EGFP/+ as we previously reported,
which also showed no T lineage‐biased feature. Furthermore, there was also no T lineage‐biased differentiation in CD19cre/RosaLSL‐Hoxb5‐EGFP/+mouse, though T lymphocytes were abundantly produced in the recipients which were transplanted with pro‐pre‐B cells sorted from CD19cre/RosaLSL‐Hoxb5‐EGFP/+ mouse. Hence, these results suggested that Hoxb5 empowers different cell types with different cell fates.Reportedly, ectopic expression of either Sox17 or miR‐125a in MPP can confer a self‐renewal ability but eventually resulted in haematological malignancies. Sox17 eventually led to leukemogenesis within 374 days after transplantation and MiR‐125a‐induced MPN displayed a complex manner of oncogene dependency.
,
Interestingly, no haematological malignancies were found in the recipients transplanted with the Hoxb5‐expressing MPP. Thus, the self‐renewal feature activated by Hoxb5 might be insulated from oncogenesis.We also tested the engraftment potential of HOXB5‐expressing human MPP in immunodeficient animals. Unfortunately, the HOXB5‐expressing human MPP failed to recapitulate the long‐term engraftment phenotype of Hoxb5‐expressing murine MPP (data not shown). One possible reason is that the function of HOXB5 is not conservative between human and mouse species. However, we cannot exclude another possibility that HOXB5‐overexpressing human MPP need a humanized bone marrow micro‐environment for HOXB5‐reprgramming, which is not available in current immunodeficient animal models.In conclusion, our study reveals a rare role of Hoxb5 in empowering self‐renewal capacity on MPP, which provides insights into converting blood progenitors into alternative engraftable cell sources.
AUTHOR CONTRIBUTIONS
Dehao Huang and Qianhao Zhao are co‐first authors. Jinyong Wang, Fangxiao Hu and Qianhao Zhao designed the project. Fangxiao Hu, Qi Zhang and Dehao Huang conducted all the experiments and data analysis. Mengyun Zhang performed the RNA‐seq experiments and Qitong Weng analysed the RNA‐seq data. Qi Zhang performed a part of the mouse genotype experiments. Kaitao Wang and Qitong Weng performed the irradiation experiments. Hui Cheng, Fang Dong, Fangxiao Hu and Jinyong Wang discussed the data. Fangxiao Hu and Jinyong Wang wrote the manuscript and approved it.
CONFLICT OF INTEREST
The authors declare that there are no competing financial interests in relation to the work described.Figure S1Click here for additional data file.
Authors: Mark J Kiel; Omer H Yilmaz; Toshihide Iwashita; Osman H Yilmaz; Cox Terhorst; Sean J Morrison Journal: Cell Date: 2005-07-01 Impact factor: 41.582
Authors: Taro Sakamaki; Kevin S Kao; Katsuyuki Nishi; James Y Chen; Kay Sadaoka; Momo Fujii; Akifumi Takaori-Kondo; Irving L Weissman; Masanori Miyanishi Journal: Biochem Biophys Res Commun Date: 2021-01-05 Impact factor: 3.575
Authors: Shangqin Guo; Haitao Bai; Cynthia M Megyola; Stephanie Halene; Diane S Krause; David T Scadden; Jun Lu Journal: Proc Natl Acad Sci U S A Date: 2012-09-24 Impact factor: 11.205
Authors: Aravind Subramanian; Pablo Tamayo; Vamsi K Mootha; Sayan Mukherjee; Benjamin L Ebert; Michael A Gillette; Amanda Paulovich; Scott L Pomeroy; Todd R Golub; Eric S Lander; Jill P Mesirov Journal: Proc Natl Acad Sci U S A Date: 2005-09-30 Impact factor: 11.205
Authors: Edyta E Wojtowicz; Eric R Lechman; Karin G Hermans; Erwin M Schoof; Erno Wienholds; Ruth Isserlin; Peter A van Veelen; Mathilde J C Broekhuis; George M C Janssen; Aaron Trotman-Grant; Stephanie M Dobson; Gabriela Krivdova; Jantje Elzinga; James Kennedy; Olga I Gan; Ankit Sinha; Vladimir Ignatchenko; Thomas Kislinger; Bertien Dethmers-Ausema; Ellen Weersing; Mir Farshid Alemdehy; Hans W J de Looper; Gary D Bader; Martha Ritsema; Stefan J Erkeland; Leonid V Bystrykh; John E Dick; Gerald de Haan Journal: Cell Stem Cell Date: 2016-07-14 Impact factor: 24.633