Dynamic changes of ribosome biosynthesis during HSC development
Our prior studies indicated that ribosome biogenesis (RiBi) process was enriched both in human and mouse HECs25, 37. To substantiate these findings, we utilized our single-cell RNA transcriptome data38, 39 and conducted gene set enrichment analysis (GSEA), the results reaffirmed the substantial enrichment of RiBi-related terms in both mouse and human HECs when compared to AECs (Fig S1A). Subsequently, employing a total of 529 genes associated with RiBi, we investigated their expression profiles throughout entire HSC ontogeny in mice. We found that the RiBi score was dramatically elevated during the differentiation of AEC into T1 pre-HSC, then gradually decreased (Fig. 1A). In contrast to the random selection of 529 genes, it is noteworthy that these RiBi genes were able to effectively distinguish consecutive developmental stages. Additionally, several of these genes, including ribosomal subunit genes and assembly factors displayed developmental stage-specific expression patterns, indicating the acquirement of distinct RiBi features during HSC development (Fig. 1B, Fig S1B-D). Further functional annotation of these RiBi genes unveiled that massive rRNA biosynthesis-related terms were prominently enriched in HECs and pre-HSCs, encompassing rRNA transcription, rRNA modification and Ribosome assembly (Fig. 1C, Fig S1E-H). Meanwhile, there was a notable upregulation in the expression levels of RNA polymerase I constituents responsible for rRNA synthesis in both HECs and pre-HSCs, which was in harmony with the amplified transcription of the 47S rRNA precursor (47S pre-rRNA), as evidenced by quantitative PCR assays (Fig. 1D, E).
Moreover, the expression of these RiBi-related genes positively correlated with hematopoietic feature genes, whereas it negatively correlated with genes featuring arterial endothelial identity during AEC to pre-HSCs transition (Fig. 1F, Fig S1I). Collectively, these results indicated that RiBi is notably activated and may play a crucial role during HSC development.
Disruption of RiBi impairs the generation of HSCs
To determine the role of RiBi in HSC generation, we sorted CD31+CD45− cells from E10.0 AGM, which encompass the majority of hematopoietic precursors. Subsequently, these cells were treated with CX-5461, a specific inhibitor of RNA polymerase I that has been previously reported to disturb rRNA transcription and ultimately RiBi33, 36. CX-5461 treatment could effectively restrain the expression level of 47S pre-rRNA in CD31+CD45− cells (Fig S2C). Of note, following the co-culture assay, the proportion of hematopoietic cells (CD45+Kit+) exhibited a significant decrease concomitant with increasing concentrations of CX-5461 (Fig S2A, B). This observation was further substantiated by explant culture assay, which demonstrated a notable decline in both the number and proportion of CD45+ hematopoietic cells as well as immunophenotypic hematopoietic stem and progenitor cells (HSPCs, CD45+Kit+Sca1+CD201+) after treatment of CX5461 (Fig. 2A-F). Importantly, the transplantation experiments utilizing the derivatives of explant cultures showed significantly lower chimerism in the peripheral blood at 16 weeks post-transplantation in CX5461 treatment groups when compared to controls (Fig. 2G), suggesting blocking RiBi dramatically hindered the production of adult repopulating HSCs. Taken together, these results validated that RiBi was pivotal for the generation of HSCs during embryonic development.
Interfering RiBi prevents the generation of HSCs from HECs rather than pre-HSCs
The generation of HSCs is a complex and tightly regulated process that primarily involves initial HEC specification, followed by a fate transition to pre-HSCs, which ultimately acquire definitive hematopoietic function. We then separately sorted HSC-primed HECs (CD41−CD43−CD45−CD31+CD201+Kit+CD44+, PK44) at E10.0 and the functional T1 pre-HSCs at E11.0 as previously identified25, 39, and combined flow cytometric analysis and transplantation assays to determine the stages at which the defect occurred due to the inhibition of RiBi (Fig. 3A). Following the co-culture of HSC-primed HECs, the number of CD45+ hematopoietic cells was significantly reduced in groups treated with CX-5461 compared to the controls (Fig. 3B). Particularly, both the quantity and proportion of the immunophenotypic HSPCs (CD45+Kit+Sca1+CD201+) also substantially diminished upon treatment with CX-5461 (Fig. 3C-E). This finding was further validated by transplantation experiments, which demonstrated a failure in reconstitution in CX-5461 treatment groups at 16 weeks post-transplantation, in contrast to the controls (Fig. 3F). PK44-represented HECs exhibit a continuum of cellular states of EHT, encompassing the endothelial, dual endothelial and hematopoietic, and hematopoietic potential25, 40. We further found CX-5461 treatment only led to diminishing the hemogenic capacity of PK44 cells, while leaving endothelial potential unaffected (Fig S3A-C). Nevertheless, upon conducting co-culture experiments with T1 pre-HSCs (CD31+CD41lowCD45−), we observed that although the progenies of T1 pre-HSCs presented a slight reduction in CX-5461 treatment groups, inhibition of RiBi in T1 pre-HSCs did not compromise their ability to generate functional HSCs (Fig. 3G-K). Overall, these results suggest that RiBi acts at HEC stage but is not required for subsequent pre-HSC maturation to form functional HSCs.
Perturbation of RiBi in HECs specifically induces cell cycle arrest
Previous study indicated the activation of cell cycle is essential for EHT and the ensuing subsequent HSPCs generation24. When compared to AECs and T1 pre-HSCs, HECs showed a more active cycling status, characterized by a greater proportion of cells in the G1/S phase (Fig. 4A, B). Meanwhile, cells with higher RiBi level were more likely to enter the G1/S phase, and the expression of genes related to RiBi exhibited a positive correlation with those associated with the S/G2/M phases (Fig. 4C-E). As several studies have clarified that RiBi might be involved in regulating cell cycle progress 32, 36, we therefore hypothesized that inhibition of RiBi in HECs could induce cell cycle arrest, subsequently disrupting their hemogenic capacity. Of note, when treated HECs with CX-5461 for a 48-hour short duration, there was a marked increase in the proportion of cells in the G0 phase and a corresponding decrease in the S/G2/M phase, indicating the occurrence of cell cycle arrest in HECs upon RiBi inhibition (Fig. 4F, G). To exclude the effect of CX-5461 on cell apoptosis, additional experiments were conducted to assess the cell death status, which revealing that there was no significant difference in HECs for two groups (Fig S3D, E). Intriguingly, we also found that the distribution of each cell cycle phase in the progenies of T1 pre-HSCs remained unchanged after treatment with CX-5461 (Fig. 4H, I). In summary, these findings implied that RiBi might preserve the ability of HECs to generate HSCs though modulating progression of the cell cycle.
Runx1 pre-configurates RiBi signature in AECs prior to HECs
A previous study has elucidated that Runx1 appears to directly regulate RiBi in mouse adult bone marrow HSPCs33. To explore whether RiBi is also regulated by Runx1 during EHT, we employed and re-analyzed Runx1 ChIP-seq data from AECs to HECs in our previous report41. First, the number of Runx1 binding peaks in the genome are similar between AECs and HECs (Fig. 5A, C; Fig S4A). Consistent with the observation that the amount of open genomic regions tending to be occupied by Runx1 was also similar between AECs and HECs as evidenced by ATAC-seq data (Fig. 5B, C). Second, we observed that over half of the Runx1-binding peaks were annotated to promoter regions and exhibited three major enrichment patterns, comprising AEC-specific (C1), HEC-specific (C3) and common shared in AECs and HECs (C2) (Fig. 5C, D). Furthermore, we conducted Gene Ontology (GO) analysis on the genes that exhibit Runx1 binding within their promoter regions accompanying with upregulated expression, which indicated enrichment of several RiBi-related terms as early as in C1 including ribonucleoprotein complex export from nucleus and ribonucleoprotein complex localization, while more specific terms, such as ribonucleoprotein complex biogenesis, regulation of transcription by RNA polymerase I, and ribosome biogenesis were significantly enriched in C2 (Fig. 5E, Fig S4B-D). Accurately, we found that Runx1 exhibited a preference for binding to promoters rather than introns of RiBi-related genes as well as the genes encoding ribosomal proteins in both AECs and HECs, and the expression level of these genes was up-regulated in HECs compared with AECs (Fig. 5F-H). Collectively, these results suggested that Runx1 already binds to RiBi-related genes in AECs to facilitate the subsequent up-regulation of these genes in HECs.
Next, we utilized and deep-analyzed the published scRNA-seq data26 of EHT-related cells from Runx1+/+ and Runx1+/− littermates, in which the identified pre-HEC transcriptomically corresponds to the AEC we annotated25, 38. We indeed observed a continuous increase in RiBi level along the development path from conflux AECs to IACs in Runx1+/+ (wild type, WT) embryos (Fig. 5I). Notably, Runx1 haploinsufficiency resulted in the increase of the proportion of pre-HEC/AEC compared to control26. Strikingly, the expression levels of several ribosomal protein genes in pre-HEC/AEC were indeed found to be reduced in Runx1+/− compared with in Runx1+/+ (Fig. 5J), further confirming the notion that Runx1 was involved in modulating RiBi gene expression as early as in AECs. In summary, these results together demonstrated that the binding of Runx1 onto the genome participating in regulating RiBi gene expression might contribute to HSC development.