Erythropoietin directly remodels the clonal composition of murine hematopoietic multipotent progenitor cells
Abstract
The cytokine erythropoietin (EPO) is a potent inducer of erythrocyte development and one of the most prescribed biopharmaceuticals. The action of EPO on erythroid progenitor cells is well established, but its direct action on hematopoietic stem and progenitor cells (HSPCs) is still debated. Here, using cellular barcoding, we traced the differentiation of hundreds of single murine HSPCs, after ex vivo EPO exposure and transplantation, in five different hematopoietic cell lineages, and observed the transient occurrence of high-output myeloid-erythroid-megakaryocyte-biased and myeloid-B-cell-dendritic cell-biased clones. Single-cell RNA sequencing analysis of ex vivo EPO-exposed HSPCs revealed that EPO induced the upregulation of erythroid associated genes in a subset of HSPCs, overlapping with multipotent progenitor (MPP) 1 and MPP2. Transplantation of barcoded EPO-exposed MPP2 confirmed their enrichment in myeloid-erythroid-biased clones. Collectively, our data show that EPO does act directly on MPP independent of the niche and modulates fate by remodeling the clonal composition of the MPP pool.
Editor's evaluation
This paper will be of broad interest to readers in the field of cytokine signaling, experimental hematology, and clinical hematology. Erythropoietin is one of the most widely used cytokines clinically, but the cells it exerts its effects on has been debated. This study has combined clonal lineage tracing and single-cell sequencing to understand the cell population that responds to erythropoietin and indicates that erythropoietin acts directly on multipotent progenitors to transiently modulate their output.
https://doi.org/10.7554/eLife.66922.sa0Introduction
Erythrocytes are the most numerous hematopoietic cells in our body and are constantly renewed (Sender et al., 2016). The major inducer of erythroid cell development in steady state and anemic conditions is the cytokine erythropoietin (EPO) (Richmond et al., 2005). Recombinant EPO is widely used to treat anemia and is one of the most sold biopharmaceuticals (Walsh, 2014). Previously, EPO was thought to solely target erythroid-committed progenitors and induce their increased proliferation and survival via the EPO receptor (EPOR) (Koury, 2016). Recently, EPO has also been suggested to act on hematopoietic stem and progenitor cells (HSPCs) (Cheshier et al., 2007; Shiozawa et al., 2010; Grover et al., 2014; Giladi et al., 2018; Yang et al., 2017; Tusi et al., 2018; Singh et al., 2018; Dubart et al., 1994), but the nature of EPO’s effect on HSPC fate remains unresolved despite potential adverse side effects during long-term EPO usage in the clinics and associations of high EPO levels with leukemias (Ma et al., 2010; Weinreb et al., 2020; Li et al., 2015).
It is well established that EPO can induce HSPCs to cycle, as evidenced by a number of bulk and single-cell studies in vitro and in vivo (Cheshier et al., 2007; Giladi et al., 2018; Yang et al., 2017; Singh et al., 2018; Dubart et al., 1994). Its role in modulating HSPC fate is less clear, however, with a lack of studies that functionally assess HSPC fate at the single-cell level in vivo and analyze the direct effect of EPO on HSPCs, and not the effect of the surrounding niche. More specifically, the upregulation of erythroid associated genes in HSPC (LSK CD150+ Flt3- CD48-) in response to in vivo EPO has been observed with bulk transcriptomics (Singh et al., 2018), suggesting that HSPCs deviate their fate toward erythroid production. Recent single-cell RNA sequencing (scRNAseq) analysis observed different changes of lineage-associated gene expressions after in vivo EPO exposure of HSPCs (Yang et al., 2017; Giladi et al., 2018). As changes of gene expression do not necessarily result in cell fate modification (Weinreb et al., 2020), functional validations in vivo are necessary. In one study, such functional validation was performed using bulk transplantation of in vivo EPO-exposed HSPCs (LSK CD150+ Flt3-) (Grover et al., 2014). This study showed increased erythroid production and decreased myeloid cell production, concluding that EPO deviates the fate of HSPCs in favor of erythroid production. As EPO has also been shown to target hematopoietic niche cells (osteoblasts and osteocytes [Li et al., 2015; Singbrant et al., 2011], endothelial cells [Ito et al., 2017; Singbrant et al., 2011], adipocytes [Alvarez and Noguchi, 2017; Zhang et al., 2014], and mesenchymal stem cells [Shiozawa et al., 2010; Tari et al., 2017]), it remains however still unclear whether EPO acts directly on HSPCs, via their environment, or both.
It is now established that HSPCs encompass cells with different long-term reconstitution capacity after transplantation, as well as a heterogeneous output in terms of quantity (lineage bias) and type of cells (lineage restriction) (Dykstra et al., 2007; Morita et al., 2010; Oguro et al., 2013; Sanjuan-Pla et al., 2013; Carrelha et al., 2018; Yamamoto et al., 2013; Lu et al., 2011; Verovskaya et al., 2013). Recently, the HSPC compartment was subdivided into long-term hematopoietic stem cells (LT-HSC) and different multipotent progenitors (MPPs) (MPP1–4; Wilson et al., 2008; Cabezas-Wallscheid et al., 2014), with variation around the phenotypic definition (Pietras et al., 2015). Interestingly, HSPC composition responds to irradiation with HSC transiently self-renewing less and increasing their production of MPP2–3 (Pietras et al., 2015). In the case of EPO, conflicting results suggest that either HSC or MPP respond to EPO (Giladi et al., 2018; Yang et al., 2017; Singh et al., 2018), but the difference in HSPC definition and the lack of functional validation make it difficult to compare these studies. As HSPCs are functionally heterogeneous, and the current phenotypic definition partially captures this heterogeneity, a single-cell in vivo lineage-tracing approach is needed to assess whether EPO can influence HSPC fate decisions.
To analyze the functional effect of EPO on the differentiation of individual HSPCs (C-Kit+ Sca1+ CD150+ Flt3-) removing the effect of the niche, we here utilized cellular barcoding technology that allowed us to trace the progeny of hundreds of single HSPCs in vivo. By analyzing cellular barcodes in five mature hematopoietic lineages and HSPCs, we observed transient induction of high-output myeloid-erythroid-(MEK)-biased barcode clones compensated by myeloid-B-cell-dendritic cell (MBDC)-biased clones after ex vivo EPO exposure. scRNAseq of ex vivo EPO-exposed HSPCs revealed upregulation of erythroid associated genes in a subset of the compartment with overlap to gene signatures of MPP1 (C-Kit+ Sca1+ Flt3- CD150+ CD48- CD34+) and MPP2 (C-Kit+ Sca1+ Flt3- CD150+ CD48+) (Cabezas-Wallscheid et al., 2014; Wilson et al., 2008) and not of LT-HSCs (C-Kit+ Sca1+ Flt3- CD150+ CD48- CD34-). Transplantation of barcoded MPP2 confirmed their enrichment in ME-biased clones in response to EPO. Moreover, the increased contribution of biased HSPC clones to the mature cell lineages after EPO exposure did not match their frequency in HSPCs, indicating that they were differentiating more than self-renewing, a property associated to MPPs. The transient effect of EPO on HSPCs further corroborates an action of EPO on MPP1/2 rather than LT-HSCs. Altogether, our results are consistent with a model in which perturbations induce clonal remodeling of HSPC contributing to hematopoiesis, with biased MPPs transiently contributing more than LT-HSC. They also demonstrate a direct effect of EPO on MPPs after transplantation with implications for basic HSC research and therapeutic applications in the clinic.
Results
EPO exposure induces biases in single HSPCs
Given the debate surrounding which HSPC subset is responding to EPO, we decided to analyze the direct effect of EPO on the differentiation of HSPCs defined as C-Kit+ Sca1+ Flt3- CD150+ (encompassing LT-HSC [C-Kit+ Sca1+ Flt3- CD150+ CD48- CD34-], MPP1 [C-Kit+ Sca1+ Flt3- CD150+ CD48- CD34+], and MPP2 [C-Kit+ Sca1+ Flt3- CD150+ CD48+]; Cabezas-Wallscheid et al., 2014; Wilson et al., 2008) at the single-cell level by cellular barcoding. To this purpose, we generated a new high-diversity lentiviral barcode library (LG2.2, 18,026 barcodes in reference list), consisting of random 20 nucleotides sequences positioned adjacent to the green fluorescent protein (GFP) gene, enabling the tracking of many individual cells in parallel. Using this LG2.2 library, we labeled single HSPCs (Figure 1—figure supplement 1a) with unique genetic barcodes as previously described (Naik et al., 2013), exposed them to EPO (1000 ng/ml) or PBS for 16 hr ex vivo, and transplanted around 2600 cells (mean 2684 cells ± 175 cells) of which around 10% barcoded cells into irradiated mice (Figure 1a). Note that HSPCs kept their sorting phenotype after ex vivo culture albeit a slight downregulation of C-Kit (Matsuoka et al., 2011) and upregulation of Flt3 (Figure 1—figure supplement 1f). At day 30 after transplantation, the earliest timepoint at which HSPCs produce simultaneously erythroid, myeloid, and lymphoid cells (Boyer et al., 2019), barcoded (GFP+) erythroblasts (E; Ter119+ CD44+34), myeloid cells (M; Ter119- CD19- CD11c- CD11b+), and B-cells (B; Ter119- CD19+) (Figure 1—figure supplement 1b,c,e) were sorted from the spleen and their barcode identity assessed through PCR and deep sequencing. Note that bone and spleen had similar barcoding profiles (Figure 1—figure supplement 2). No difference in chimerism was observed between the EPO and control group in the spleen and blood, even when mTdTomato/mGFP donor mice were used to better assess the erythroid lineage (Figure 1b and c). On average, we detected around 80 barcodes per mouse, of which most were detected in several lineages (Figure 1d). Comparison of the numbers of barcodes producing each lineage showed that EPO exposure resulted in the same number of engrafting and differentiating cells as in control (Figure 1d). Notably, the number of erythroid restricted cells remained stable in the EPO group as compared to control (Figure 1e), indicating that the response to EPO is more complex than a direct instruction of erythroid-restricted HSPCs.
To quantify the effect of EPO on HSPC lineage biases, barcode-labeled HSPCs were classified based on the balance of their cellular output in the M, B, and E lineages. With this classification using a 10% threshold, cells classify, for example, as ME-biased if they have above 10% of their output in the M and E lineage and under 10% of their output in the B lineage (Figure 1e, other thresholds in Figure 1—figure supplement 3b). Interestingly, application of this classification revealed that although the proportion of lineage-biased HSPCs in the control and EPO groups was similar (Figure 1f), their contribution to the different lineages was increased by EPO exposure (Figure 1g). In the control group, balanced HSPCs (MBE) produced the majority of all lineages, as previously published (Perié et al., 2015). In the EPO group, ME- and MB-biased clones produced most cells of the analyzed lineages (Figure 1g). ME-biased HSPCs produced the majority of erythroid cells (57% ± 10%), MB-biased HSPCs produced the majority of B-cells (58% ± 36%), and ME- and MB-biased clones contributed the majority of myeloid cells (MB-biased 45% ± 38% and ME-biased 20% ± 13%, together 65% ± 25%). To test the significance of this effect, we used a permutation test that compares the effect size between the control and EPO groups to the one of all random groupings of mice (Tak et al., 2019). The contributions of the ME- and of the MB-biased HSPC classes to the different lineages were significantly different in the EPO and control groups (Table 1). These results were reproduced in an additional experiment (Figure 1—figure supplement 3c–f, Supplementary file 1). A lower EPO concentration (160 ng/ml) as well as an additional single injection of EPO (133 µg/kg) during transplantation gave similar results (Figure 2, Figure 2—figure supplement 1a, Table 1). Also at 6 weeks post-transplantation, similar results were obtained (Figure 1—figure supplement 4, Supplementary file 1). In summary, ex vivo EPO priming of HSPCs modified the output balance of HSPCs rather than the number of lineage-restricted and -biased cells. Balanced clones produced a smaller percentage of the mature cells; ME-biased HSPCs produced most of the erythroid cells and MB-biased HSPCs produced most of the B cells.
Contribution of ME- and MB-biased HSPCs to the DC and MkP lineage
To further characterize the cells produced by the ME-biased and MB-biased HSPCs, we repeated our experimental setup including the analysis of the megakaryocyte and dendritic cell (DC) lineages (Figure 3). Megakaryocyte progenitors (MkP) were chosen as proxy for the production of platelets that are not suitable for barcode analysis. Barcoded (GFP+) DCs (DC; Donor Ter119- CD19- CD11c+ CD11b-) and MkP (MkP; C-Kit+ Sca-1- CD150+ CD41+) (Figure 1—figure supplement 1c–e) were sorted together with M, E, and B cells, 4 weeks after transplantation of control or EPO-exposed HSPCs (1000 ng/ml). In both groups, the majority of clones produced also DCs (Figure 3a). In the control group, balanced HSPCs produced the majority of DCs (65% ± 9%) (Figure 3b). However, in the EPO group, balanced HSPCs decreased their contribution to the DC lineage (36% ± 25%) and MB-biased HSPCs significantly increased their contribution (86% ± 43% EPO vs. 22% ± 11% control group) (Figure 3b, Table 1), thus, they were MBDC-biased HSPCs. In contrast, the ME-biased HSPCs produced few DCs in both groups (Figure 3b), indicating that ME-biased HSPCs are restricted both in their B and DC production compared to the M and E production.
The majority of the MkP production came from the ME-biased HSPCs in both groups (58% ± 21% control and 55% ± 14% EPO group, Figure 3c and d), indicating that ME-biased HSPCs were also MkP-biased HSPCs (thus MEK-biased). We did not detect a high contribution of MkP-restricted HSPCs (Carrelha et al., 2018; Sanjuan-Pla et al., 2013; Rodriguez-Fraticelli et al., 2018) to the MkP lineage (Figure 3d). Finally, as high EPO exposure has been linked to changes in macrophage numbers (Theurl et al., 2016; Wang et al., 2018; Gilboa et al., 2017; Kuzmac et al., 2014; Ulyanova et al., 2016; Mausberg et al., 2011; Liao et al., 2018; Bretz et al., 2018), we analyzed the contribution of control and EPO-exposed HSPCs to the myeloid lineage in more detail, but could not detect changes in the percentage of the different myeloid subsets produced (Figure 3—figure supplement 1a–c).
Effect of EPO on short-term HSPC self-renewal
In light of previous studies that suggested changes in HSPC proliferation after in vivo EPO exposure (Cheshier et al., 2007; Dubart et al., 1994; Giladi et al., 2018; Yang et al., 2017; Singh et al., 2018), we next explored if the short-term self-renewal capacity of HSPCs was impacted. To this end, we analyzed barcodes in bone marrow HSPCs in addition to the spleen E, M, and B lineages at week 4 after transplantation of control or EPO-exposed HSPCs (160 and 1000 ng/ml) (Figure 4). We reasoned that barcodes of HSPCs differentiating and short-term self-renewing (dividing to give rise to other HSPCs) after transplantation are detected in both compartments, while detection in only HSPCs or mature lineages indicates a prevalence of short-term self-renewal or differentiation respectively. Most of the barcodes detected in HSPCs overlapped with barcodes in the mature cells (Figure 4b, left) in both the control and two EPO groups, showing that most of the transplanted cells had given rise to other HSPCs and differentiated irrespective of the treatment. Some barcodes were only detected in mature cells (Figure 4b, right), indicating that some HSPCs had only differentiated or were below the limit of detection. These HSPCs were equally abundant in the control and two EPO groups (Figure 4b, right).
To analyze if different lineage biases correlated to different short-term self-renewal capacity, we analyzed the proportion of biased HSPC classes, as previously defined, within the HSPC compartment (Figure 4c). In the control group, balanced and ME-biased HSPCs contributed most to the HSPC reads (34% ± 36% MBE and 37% ± 34% ME-biased HSPCs), while barcodes of MB-biased HSPCs contributed less (15% ± 32%) (Figure 4c), a trend that has been previously described (Kim et al., 2014; Aiuti et al., 2013; Carrelha et al., 2018). Surprisingly, the pattern of contributions of different biased HSPC subsets to HSPC reads was unchanged in the EPO groups (Figure 4c), implying that the extent of short-term self-renewal was unchanged after ex vivo EPO exposure.
To study if the increased production of cells by the ME- and MB-biased HSPCs to the mature cells observed after ex vivo EPO exposure (Figures 1—3) correlated with short-term self-renewal capacity of HSPCs, we analyzed the contribution of barcodes detected or not in HSPCs to the E, M, and B lineages (Figure 4d). In the control group, the majority of mature cells were derived from barcodes also present in HSPCs. However, in both EPO groups, the contribution of barcodes detected in HSPCs to mature cells was significantly lower (Figure 4d, Table 1), implying that the increased contribution of biased HSPC classes to the mature cell lineages after ex vivo EPO exposure was most likely caused by cells differentiating more than short-term self-renewing.
EPO exposure induces an erythroid program in a subgroup of HSPCs
To further characterize the effect of EPO exposure on HSPCs, we performed scRNAseq of barcoded C-Kit+ Sca1+ Flt3- CD150+ cells after ex vivo culture in medium supplemented with EPO or PBS using the 10X Genomics Chromium platform. 1706 cells from control and 1595 cells from the EPO group passed our quality control. To compare the HSPCs injected with noncultured hematopoietic cells, we generated a reference map of 44,802 C-kit+ cells from Dahlin et al., 2018 and used published signatures as detailed in Materials and methods (Pietras et al., 2015; Wolf et al., 2019) to annotate this map (Figure 5—figure supplement 1a, b and e, f). Projection of our single-cell data on this map showed that both the control and the EPO-exposed HSPCs similarly overlapped with non-MPP4 LSK cells, according to their sorting phenotype (Figure 5—figure supplement 1e and f). These results indicate that neither the ex vivo culture itself nor the EPO treatment dramatically affected the global identity of the sorted HSPCs.
When comparing the EPO and control groups, we found 1176 differentially expressed genes (Figure 5a) and this number was significantly higher than the number expected due to chance (p-value=0.01) as assessed by permutation testing. Among the most upregulated genes in the EPO-exposed HSPCs were genes with erythroid association as Hbb-bs, Erdr1, Wtap, Kmt2d, or Nfia (Starnes et al., 2009), and GATA1 targets (Abhd2, Cbx3, Kdelr2, Pfas), cell cycle-related genes (Tubb5, Hist1h2ap), as well as genes previously described to be induced in HSPCs after in vivo EPO exposure, such as Bmp2k (Shiozawa et al., 2010) and Ifitm1 (Giladi et al., 2018); (Figure 5a). Genes involved in stem cell maintenance, such as Serpina3g, Mecom, Txnip, Meis1, Pdzk1ip1(Giladi et al., 2018), Sqstm1 (Meenhuis et al., 2011), Smad7 (Blank et al., 2006), and Aes (Steffen et al., 2011), were among the most downregulated genes in the EPO-exposed HSPCs (Figure 5a).
As our cellular barcoding data suggests that single HSPCs differ in their response to EPO, we assessed the heterogeneity of EPO responses at the transcriptomic level. UMAP-based visualization of the data suggested that a subgroup of EPO-exposed cells was transcriptomically distinct (Figure 5b), independently of the number of principal component analysis (PCA) components and genes used in the analysis (Figure 5—figure supplement 1c). To test this observation, we defined an EPO response signature based on differentially expressed genes between the EPO and control group. Plotting the expression of the EPO response signature at the single-cell level showed that the majority of the transcriptomic differences between the control and EPO group were indeed driven by this small subgroup of cells (Figure 5c). Reasoning that this subgroup contains the cells directly responding to EPO, we defined as EPO responders, cells in the 90th percentile of EPO response signature expression (Figure 5c) for subsequent analysis. Importantly, unsupervised clustering analysis of the data (Figure 5—figure supplement 2a and b) showed similar results. The genes encoding EPOR, as well as the alternative EPORs EphB4, CD131, and CRFL3, were equally expressed between the EPO responders, nonresponders, and control groups (Figure 5—figure supplement 1d). Reasoning that the EPO responders correspond to MEK-biased HSPCs, we also looked for potential MBDC-biased HSPCs but could not detect a subgroup of cells with upregulation of lymphoid-associated genes, suggesting that the MBDC bias is not a direct effect of EPO exposure but more an indirect effect. In summary, the scRNAseq analysis corroborated our functional barcoding data, showing that a subset of HSPCs can respond directly to EPO stimulation.
EPO responder HSPCs overlap with MPP1 and MPP2 signatures
As our barcode analysis suggested that the effect of direct EPO exposure on HSPCs is caused by cells differentiating more than self-renewing, we next wanted to assess which of the HSPC subsets are the EPO responders in our scRNAseq dataset. We annotated the UMAP-based visualization of our data with published signatures of the HSC (dormant HSCs [Cabezas-Wallscheid et al., 2017] and LT-HSC [Wilson et al., 2015]), MPP1 (Cabezas-Wallscheid et al., 2017), and MPP2 (Pietras et al., 2015) subsets included in our HSPC gate and analyzed its overlap with the previously defined EPO responder and nonresponder cells (Figure 6a and b). Relative to the control group and nonresponders of the EPO group, the EPO responders had a reduced expression of HSC gene signatures and increased expression of MPP1 and MPP2 signatures (Figure 6a and b). An annotation of the reference map generated from data of Dahlin et al., 2018 likewise showed a low overlap of EPO responders with the most quiescent HSC subsets (Figure 6c and d). The independent analysis using unsupervised clustering further supported this result (Figure 5—figure supplement 2c and d). All in all, our scRNAseq analysis implied that, in line with our barcoding results, the HSPCs directly reacting to EPO are most likely MPP cells of the MPP1 and MPP2 subsets.
EPO exposure induces ME biases in single MPP2
To confirm that MPP2 are a subset within HSPCs reacting directly to EPO as predicted by the scRNAseq analysis, we transplanted barcoded control or EPO-exposed (1000 ng/ml) MPP2 together with unbarcoded CD48- HSPCS (C-Kit+ Sca1+ Flt3- CD150+ CD48-) (Figure 7—figure supplement 1a) and analyzed their barcoded progeny in the E, M, and B lineages of the spleen at week 4 after transplantation (Figure 7). We found an equivalent engraftment as for the entire HSPC compartment and no difference between the EPO-treated and the control group (Figure 7a and b). Applying the same classification as in Figure 1 to quantify the effect of EPO on MPP2 lineage biases, we observed that, as for the whole HSPC compartment, ME-biased cells contributed more to the M and E lineages (Figure 7d and e, other threshold in Figure 7—figure supplement 1c). Similarly to our data on whole HSPC compartment (Figure 1g), the proportion of the differently biased MPP2 was similar between the control and EPO groups (Figure 7c). This data confirms that the MPP2 population is enriched in HSPCs responding to EPO.
Transient effect of EPO exposure
Finally, we reasoned that if EPO directly acts on MPPs 1/2 with short reconstitution capacity after transplantation rather than long-term repopulating HSC, then the EPO effect should be transient. To test this hypothesis, we repeated the experiment and analyzed barcodes in the E, M, and B lineages at 4 months after transplantation of control or EPO-exposed HSPCs (160 and 1000 ng/ml) (Figure 8). In the control group, as reported before (Wu et al., 2018), the chimerism at 4 months was higher and the number of barcodes detected was lower than at 1 month post-transplantation (Figures 8c and d and 1b–d). While ME-biased HSPCs still have significant higher contribution to the M and E lineages (Figure 8a and b, p-values in Table 1), the majority of cells in all lineages were produced by balanced HSPCs (Figure 8a and b), implying that the effect of direct EPO exposure on HSPCs is fading away. This confirms that the effect of direct EPO exposure on HSPCs is likely caused by MPP cells with a short reconstitution capacity after transplantation.
Discussion
EPO is a key regulator of hematopoiesis and is classically considered to support the proliferation and survival of erythroid-committed progenitors. By analyzing the in vivo fate of hundreds of EPO-stimulated vs. untreated transplanted (c-Kit+ Sca1+ Flt3- CD150+) HSPCs at the single-cell level, we established that EPO can change HSPC differentiation in the absence of an EPO-stimulated bone marrow microenvironment. Collectively our results yield two important conclusions: (i) EPO has a direct effect on HSPCs, that is, not solely due to the effects of EPO on the surrounding niche; and (ii) EPO directly remodels the clonal composition of HSPC by inducing fate-biased MPP and reducing the output of HSC.
Specifically, we observe that EPO induced MEK-biased (ME) and MBDC-biased (MB) HSPCs that produced the majority (>60%) of mature hematopoietic cells at 4 and 6 weeks after transplantation. In contrast, balanced HSPCs (MBE) had a reduced output of mature cells in response to EPO. The increased erythroid-associated gene signature in a subset of HSPCs after ex vivo EPO exposure suggests that EPO directly induces high-output MEK-biased HSPCs, which is indirectly compensated for by the occurrence of high output of MBDC-biased HSPCs to maintain a balanced production of hematopoietic cells.
These biased clones had a higher propensity to differentiate than to self-renew, and their response to EPO was transient, suggesting that EPO-responsive cells are multipotent progenitors, and not LT-HSCs. This is supported by transcriptomic analysis showing that EPO responders express the MPP1/MPP2 gene signatures. Transplantation of barcoded EPO-exposed MPP2 confirmed their enrichment in ME-biased clones in response to EPO.
Similar to studies that assessed the effect of high systemic EPO exposure (Tusi et al., 2018; Yang et al., 2017; Singh et al., 2018; Giladi et al., 2018), we found that MPPs, not HSCs, are responding to EPO. The occurrence of myeloid, megakaryocytic, and erythroid gene expression in MPP1 after bleeding (Yang et al., 2017) is in line with our findings. Previously, long-term EPO exposure in Tg6 transgenic mice did however not change the in vitro differentiation outcome of MPP2 (Singh et al., 2018). Furthermore, we did not detect a fate deviation toward erythroid production at the expense of myeloid production as seen for in vivo EPO-exposed HSPCs after transplantation (Grover et al., 2014). These differences could be due both to the duration and route of EPO exposure, as well as the indirect effects of systemic EPO exposure through other cells, for example, from the bone marrow niche.
In addition, our data shows that direct cytokine stimulation leads to a clonal remodeling of the HSPC compartment, with a transient increase in the contribution of fate-biased MPPs. Without longitudinal barcoding data within the same animal, we cannot distinguish if EPO is transiently changing the fate and outcome of the same HSPCs over time or if EPO is pushing the differentiation of some HSPCs that will be replaced by more balanced and stable HSPCs.
Different studies have suggested that the behavior of transplanted HSPC differs from native HSPC (Sun et al., 2014; Busch et al., 2015). Transplantation seems to favor the long-term output from HSC, whereas steady-state hematopoiesis is maintained more by MPPs’ contribution than LT-HSC (Sun et al., 2014; Busch et al., 2015; Schoedel et al., 2016). Our work is in line with a model in which MPPs are a highly malleable cell population that can rapidly respond to changing demands for new cells, such as transplantation or infection (Pietras et al., 2015).
The direct effect of EPO on MPPs we described here could be one of the factors underlying the development of adverse side effects and comorbidities during long-term EPO use in the clinics and associations of high EPO levels with leukemias (Rainville et al., 2016; Ma et al., 2010; Halawi et al., 2017). To translate these results to the clinics and understand the side effect of EPO treatment, further work is required to determine if HSPCs and erythroid progenitors like CFU-E are responding to the same dose and duration of EPO exposure.
Materials and methods
Mice
Male C57BL/6J CD45.1+, C57BL/6J CD45.2+, and Rosa26CreERT2;mT/mG mice from Jackson Laboratory or bred at Institute Curie aged between 7 and 13 weeks were used in all experiments. All procedures were approved by the responsible national ethics committee (APAFIS# 10955–201708171446318v1).
Barcode library, barcode reference list, and lentivirus production
Request a detailed protocolThe LG2.2 barcode library is composed of a DNA stretch of 180 bp with a 20 bp ‘N’-stretch. DsDNA was generated by 10 PCR rounds and cloned into the XhoI-EcoRI site of the lentiviral pRRL-CMV-GFP plasmid (Dull et al., 1998). Subsequently, ElectroMaxStbl4 cells were transformed, and >10,000 colonies picked for amplification by Maxiprep. To create the barcode reference list (https://github.com/PerieTeam/Eisele-et-al.-; Eisele et al., 2022), barcode plasmids were PCR amplified twice in duplicate and sequenced as described below. Sequencing results were filtered for barcode reference list generation as previously described in Naik et al., 2013. Lentiviruses were produced by transfecting the barcode plasmids and p8.9-QV and pVSVG into HEK293T cells in DMEM-GlutaMAX (Gibco) supplemented with 10% FCS (Eurobio), 1% MEM NEAA (Sigma), and 1% sodium pyruvate (Gibco) using polyethyleneimine (Polysciences). Supernatant was 0.45 µm filtered, concentrated by 1.5 hr ultracentrifugation at 31,000 × g, and frozen at –80°C. HEK293T were tested for mycoplasma contamination every 6 months.
HSPC and MPP2 isolation, barcoding, EPO treatment, and transplantation
View detailed protocolIsolation and labeling of cells with the barcoding library were performed as described in Naik et al., 2013. Briefly, after isofluorane anesthesia and cervical dislocation, bone marrow cells were isolated from femur, tibia, and iliac bones of mice by flushing, and C-Kit+ cells were enriched with anti-CD117 magnetic beads on the MACS column system (Miltenyi). Cells were stained for C-Kit, Flt3, CD150, Sca-1 (Key resources table) propidium iodide (PI) (Sigma) (1:5000), and if appropriate CD48. Lineage staining was not performed after C-Kit+ MACS enrichment for transplantation cohorts. For HSPC cohorts, HSPCs (Figure 1—figure supplement 1a) were sorted and transduced with the barcode library in StemSpanMedium SFEM (STEMCELL Technologies) with 50 ng/ml mSCF (STEMCELL Technologies) through 1.5 hr of centrifugation at 300 × g and 4.5 hr incubation at 37°C to obtain 10% barcoded cells. After transduction, cells were incubated with human recombinant EPO (Eprex, erythropoietin alpha, Janssen) at a final concentration of 1000 or 160 ng/ml or PBS for 16 hr at 37°C. After the incubation, the cells were transplanted by tail vein injection in recipient mice 6 Gy sublethally irradiated 3 hr before on a CIXD irradiator. Mice were allocated to groups of 4–5 mice for each condition randomly without masking. When indicated, cells were injected together with additional EPO (133 µg/kg). On average, 2600 cells (mean 2684 cells ± 175 cells) were injected in the tail vein of each mouse. For the MPP2 cohort, MPP2 and CD48- HSPCs (Figure 7—figure supplement 1a) were sorted. Both populations were cultured alike, but only MPP2 were transduced with the barcode library and treated with 1000 ng/ml recombinant EPO as described above. After the culture, barcoded MPP2 and unbarcoded CD48- HSPCs were mixed at a ratio of 32/45 (to be as close as possible to the original ratio of both populations in the HSPCs) and transplanted as described above. A FACSAria (BD Biosciences) was used for sorting. FACSDiva software (BD Biosciences) was used for measurements and FlowJo v.10 (TreeStar) for analysis.
Cell progeny isolation for barcode analysis
Request a detailed protocolSpleens were mashed and both blood and spleen cells were separated based on Ter119 using a biotinylated anti-Ter119 antibody (Key resources table) and anti-biotinylated beads on the MACS column system (Miltenyi). Ter119+ cells were stained for Ter119 and CD44 (Chen et al., 2009). Ter119- cells were stained for CD45.1 CD11b, CD11c, CD19, and, if appropriate, CD115, Siglec-F, and Ly6G (Key resources table). Bone marrow cells were flushed from bones and enriched for C-Kit+ cells as above. When appropriate, the C-Kit- fraction was further separated based on Ter119 and stained as above. C-Kit+ cells were stained for C-Kit, Flt3, CD150, Sca-1, and, if appropriate, CD41 (Key resources table), and PI (1:5000) as described above. For analyzed and/or sorted populations, see Figure 1—figure supplement 1. Populations were only sorted for mice with an engraftment (donor cells percentage) of above 5% in spleen, bone, and blood.
Lysis, barcode amplification, and sequencing
Request a detailed protocolSorted cells were lysed in 40 μl Viagen Direct PCR Lysis Reagent (cell) (Euromedex) supplemented with 0.5 mg/ml Proteinase K Solution RNA grade (Invitrogen) at 55°C for 120 min, 85°C for 30 min, and 95°C for 5 min. Samples were then split into two replicates, and a three-step nested PCR was performed to amplify barcodes and prepare for sequencing. The first step amplifies barcodes (top-LIB [5′TGCTGCCGTCAACTAGA ACA-3′] and bot-LIB [5′GATCTCGAATCAGGCGCTTA-3′]). A second step adds unique 4 bp plate indices as well as Read 1 and 2 Illumina sequences (PCR2-Read1-plate-index-forward 5′ACACTCTTTCCCTACACGACGCTCTTCCGATCTNNNNCTAGAACACTCGAGATCAG3′ and PCR2-Read2-reverse 5′GTGACTGGAGTTCAGACGTGTGCTCTTCCGAT CGATCTCGAATCAGGCGCTTA3′). In a third step, P5 and P7 flow cell attachment sequences and one of 96 sample indices of 7 bp are added (PCR3-P5-forward 5′AATGATA CGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT3′ and PCR7-P7-sample-index-reverse 5′CAAGCAGAAGACGGCATACGAGANNNNNNNGTGACTGGAGTTCAGA CGTGCTCTTCCGATC3′) (PCR program: hot start 5 min 95°C, 15 s at 95°C; 30 s at 57.2°C; 30 s at 72°C, 5 min 72°C, 30 [PCR1-2] or 15 cycles [PCR 3]). Both index sequences (sample and plate) were designed based on Faircloth and Glenn, 2012 such that sequences differed by at least 2 bp (https://github.com/PerieTeam/Eisele-et-al.-). To avoid lack of diversity at the beginning of the reads, at least four different plate indices were used for each sequencing run. Primers were ordered desalted as high-performance liquid chromatography purified. During lysis and each PCR, a mock control was added. The DNA amplification by the three PCRs was monitored by the run on a large 2% agarose gel. PCR3 products for each sample and replicate were pooled, purified with the Agencourt AMPure XP system (Beckman Coulter), diluted to 5 nM, and sequenced on a HiSeq system (Illumina) (SR-65bp) at the Institute Curie facility with 10% of PhiX spike-in.
Barcode sequence analysis
Request a detailed protocolSequencing results were filtered, and barcodes were categorized in progenitor classes as in Naik et al., 2013 and further explained on GitHub (https://github.com/PerieTeam/Eisele-et-al.-). In brief, sequencing results were analyzed using R-3.4.0 (R Development Core Team 2019; http://wwwR-project.org.), Excel, and GraphPad Prism version 8.0 for Mac (GraphPad Software, La Jolla, CA, https://www.graphpad.com). Reads were first filtered for perfect match to the input index- and common-sequences using XCALIBR (https://github.com/NKI-GCF/xcalibr) and filtered against the barcode reference list. Samples were then filtered for containing at least 5000 reads and normalized to 105 per sample. Samples with a Pearson correlation between duplicates below 0.9 were discarded, and barcodes present in one of the two replicates were set to zero. Samples with under 10 barcodes were filtered out, unless indicated in the figure legend. The mean of the replicates was used for further processing. When the mean percentage of barcodes shared between different sequencing runs was higher than within the same sequencing run for mice of a same transduction batch, reads below the read quartile of the mean percentage of barcodes shared between mice of a same transduction batch but sequenced on different sequencing runs were set to zero in order to equalize the barcode sharing between mice transplanted from a same transduction batch in different sequencing runs to the barcode sharing between mice within each sequencing run. After filtering, read counts of each barcode in the different cell lineages were normalized enabling categorization into classes of biased output toward the analyzed lineages using a threshold of 10% of barcode reads (other thresholds in Figure 1—figure supplement 3c and Figure 7—figure supplement 1c). Statistics on barcoding results were performed using a permutation test as in Tak et al., 2021. Significance of flow cytometry results was assessed using Student’s t-test. Some mice were excluded from the analysis due to death before readout or due to a donor cell engraftment <5%, as well as the filtering out of mice for which one or more cell subset samples did not pass the barcode data filtering steps as detailed above.
scRNAseq and analysis
Request a detailed protocolscRNAseq was performed using the 10X Genomics platform on one pool of HSPCs isolated from eight mice, barcoded and culture with or without EPO for 16 hr in vitro as described above. Sequencing libraries were prepared using the Chromium Single Cell 3′ v2 kit and sequenced on a HiSeq system (Illumina) at the Institut Curie NGS facility. Data was analyzed using Cellranger v3 (10X Genomics), Seurat v3 (Satija et al., 2015), and customized scripts. Raw sequencing reads were processed using Cellranger. To obtain a reads/cell/gene count table, reads were mapped to the mouse GRCm38.84 reference genome. scRNAseq analysis was performed using Seurat (Satija et al., 2015). During filtering, Gm, Rik, and Rp genes were discarded as noninformative genes. Cells with less than 1000 genes per cell and with a high percentage of mitochondrial genes were removed from downstream analyses. Following our filtering procedures, the average UMI count per cell was 5157, with mitochondrial genes accounting for 5% of this. The average number of genes detected per cell was 2337. Cell cycle annotation using the cyclone method from the scran R package showed that 2938 cells were in G1 phase, 233 cells were in G2M phase, and 127 cells were in S phase. No batch effect was detected between the EPO and no-EPO group; therefore, no batch correction was applied. Data normalization was performed using the default Seurat approach, and differentially expressed genes were determined using a logistic regression in Seurat. Unsupervised clustering was performed on the significant variable genes using the 10 first PCA followed by the nonlinear dimensionality reduction technique UMAP (McInnes et al., 2018; Figure 5—figure supplement 2). Unsupervised Louvain clustering of the data was performed across a range of resolution parameters, and the resolution value that led to the most stable clustering profiles was chosen (Blondel et al., 2008; Figure 5—figure supplement 2). Annotation of the clusters was obtained by mapping published signatures using the AddModuleScore method of Seurat. The signatures are defined in the following publications: dHSC and MPP1 signatures were obtained from Cabezas-Wallscheid et al., 2017. The MolO LT-HSC signature was taken from Wilson et al., 2015, and the MPP2 and 4 signature was taken from Pietras et al., 2015. An Excel file listing the genes in these signatures is available on GitHub (https://github.com/PerieTeam/Eisele-et-al.-). To identify EPO responder cells in the EPO group, differential expression analysis was performed between the control and EPO groups (lists of DEGs are available at https://github.com/PerieTeam/Eisele-et-al.-). Subsequently, genes that were differentially expressed (adjusted p-value<0.05) between the EPO and control groups were transformed into an EPO response signature that when overlaid onto the UMAP-based visualization was enriched only in a subset of the EPO group cells. Briefly this signature was obtained by taking the background-corrected mean expression values of both the up- and downregulated genes per cell as implemented in the AddModuleScore method of Seurat. Within each cell, these two signature scores were used to create a composite EPO response score by subtracting the downregulated response from the upregulated response signature. Cells in the upper 90th percentile with regards to the expression of the EPO response signature were labeled EPO responders.
To perform supervised cell-type annotation, a reference map was generated from a published single-cell sequencing dataset of 44,802 C-Kit+ cells from Dahlin et al., 2018. Preprocessing was performed using a scanpy pipeline (Wolf et al., 2019). Data was then visualized using the nonlinear nondimensionality reduction technique UMAP (McInnes et al., 2018). Annotation of the reference map was obtained by overlaying published signatures as above using the AddModuleScore method of Seurat and also known markers as Flt3, slamf1, and Gata1 (Figure 5—figure supplement 1b). For the erythroid progenitors, these markers are Gata1, Klf1, Epor, Gypa, Hba-a2, Hba-a1 (Figure 5—figure supplement 1a). Cells were mapped onto the reference map using a k-nearest-neighbors mapping approach. Briefly, for each cell in the query dataset, the nearest neighbors in the PCA space of the reference dataset were determined using the nn2 function of the RANN package, and the mean UMAP 1 and 2 coordinates of the 10 nearest neighbors were taken as the reference point for the new cell of interest. To benchmark our mapping approach, cells from an independent dataset of erythroid progenitors Tusi et al., 2018 were used without additional preprocessing (Figure 5—figure supplement 1a).
Data availability statement
Request a detailed protocolRaw data are available at zenodo doi:10.5281/zenodo.5645045. All codes to filter and process raw data, as well as filtered data, are available at https://github.com/PerieTeam/Eisele-et-al.-. Contact author is leila.perie@curie.fr.
Data availability
all data and scripts are available on the github of the Perie lab https://github.com/PerieTeam/Eisele-et-al.- (copy archived at swh:1:rev:ff1da6c9b3ec2b8e14e5921aeb2ac70fa2bcced0) .
-
GitHubID GitHub. Erythropoietin directly remodels the clonal composition of murine hematopoietic multipotent progenitor cells.
References
-
Lentiviral hematopoietic stem cell gene therapy in patients with Wiskott-Aldrich syndromeScience (New York, N.Y.) 341:1233151.https://doi.org/10.1126/science.1233151
-
Fast unfolding of communities in large networksJournal of Statistical Mechanics 2008:10008.https://doi.org/10.1088/1742-5468/2008/10/P10008
-
The effect of bleeding on hematopoietic stem cell cycling and self-renewalStem Cells and Development 16:707–717.https://doi.org/10.1089/scd.2007.0017
-
A Third-Generation Lentivirus Vector with A Conditional Packaging SystemJournal of Virology 72:8463–8471.https://doi.org/10.1128/JVI.72.11.8463-8471.1998
-
SoftwareEisele et al, version swh:1:rev:ff1da6c9b3ec2b8e14e5921aeb2ac70fa2bcced0Software Heritage.
-
Erythropoietin guides multipotent hematopoietic progenitor cells toward an erythroid fateThe Journal of Experimental Medicine 211:181–188.https://doi.org/10.1084/jem.20131189
-
A higher prevalence of hematologic malignancies in patients with thalassemia: Background and culpritsAmerican Journal of Hematology 92:414–416.https://doi.org/10.1002/ajh.24682
-
Bone Marrow Endothelial Cells Induce Immature and Mature B Cell Egress in Response to ErythropoietinCell Structure and Function 42:149–157.https://doi.org/10.1247/csf.17018
-
Tracking erythroid progenitor cells in times of need and times of plentyExperimental Hematology 44:653–663.https://doi.org/10.1016/j.exphem.2015.10.007
-
Erythropoietin promotes bone formation through EphrinB2/EphB4 signalingJournal of Dental Research 94:455–463.https://doi.org/10.1177/0022034514566431
-
Low level of C-kit expression marks deeply quiescent murine hematopoietic stem cellsStem Cells (Dayton, Ohio) 29:1783–1791.https://doi.org/10.1002/stem.721
-
undefined Uniform Manifold Approximation and Projection.Journal of Open Source Software 3:861.https://doi.org/10.21105/joss.00861
-
Heterogeneity and hierarchy within the most primitive hematopoietic stem cell compartmentThe Journal of Experimental Medicine 207:1173–1182.https://doi.org/10.1084/jem.20091318
-
Targeting EPO and EPO receptor pathways in anemia and dysregulated erythropoiesisExpert.Opinion on Therapeutic Targets 20:287–301.https://doi.org/10.1517/14728222.2016.1090975
-
Turning cells red: signal transduction mediated by erythropoietinTrends in Cell Biology 15:146–155.https://doi.org/10.1016/j.tcb.2005.01.007
-
Spatial reconstruction of single-cell gene expression dataNature Biotechnology 33:495–502.https://doi.org/10.1038/nbt.3192
-
Biopharmaceutical benchmarks 2014Nature Biotechnology 32:992–1000.https://doi.org/10.1038/nbt.3040
-
Lineage tracing on transcriptional landscapes links state to fate during differentiationScience (New York, N.Y.) 367:eaaw3381.https://doi.org/10.1126/science.aaw3381
-
Geographic clonal tracking in macaques provides insights into HSPC migration and differentiationThe Journal of Experimental Medicine 215:217–232.https://doi.org/10.1084/jem.20171341
-
Single cell transcriptomics reveals unanticipated features of early hematopoietic precursorsNucleic Acids Research 45:1281–1296.https://doi.org/10.1093/nar/gkw1214
-
Erythropoietin action in stress response, tissue maintenance and metabolismInternational Journal of Molecular Sciences 15:10296–10333.https://doi.org/10.3390/ijms150610296
Article and author information
Author details
Funding
ATIP Avenir CNRS
- Leïla Perié
Labex Cell(n)scale (ANR-10-LABX-31)
- Leïla Perié
Idex Paris-Science-Lettres (ANR-10-IDEX-0001-02 PSL)
- Leïla Perié
H2020 European Research Council (758170-Microbar)
- Leïla Perié
H2020 Marie Skłodowska-Curie Actions (666003)
- Almut S Eisele
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Dr. T Schumacher for discussion and lentiviral library production and Dr. K Duffy for advices on permutation testing. We thank the Institute Curie flow cytometry, next-generation sequencing, and animal facility. We thank Fahima Di Federico from the UMR168 BMBC facility for amplifying the barcode plasmid pool.
Ethics
All procedures were approved by the responsible national ethics committee (APAFIS#10955-201708171446318 v1).
Copyright
© 2022, Eisele et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,486
- views
-
- 242
- downloads
-
- 14
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
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)
Further reading
-
- Immunology and Inflammation
Adipose tissue inflammation is now considered to be a key process underlying metabolic diseases in obese individuals. However, it remains unclear how adipose inflammation is initiated and maintained or the mechanism by which inflammation develops. We found that microRNA-802 (Mir802) expression in adipose tissue is progressively increased with the development of dietary obesity in obese mice and humans. The increasing trend of Mir802 preceded the accumulation of macrophages. Adipose tissue-specific knockout of Mir802 lowered macrophage infiltration and ameliorated systemic insulin resistance. Conversely, the specific overexpression of Mir802 in adipose tissue aggravated adipose inflammation in mice fed a high-fat diet. Mechanistically, Mir802 activates noncanonical and canonical NF-κB pathways by targeting its negative regulator, TRAF3. Next, NF-κB orchestrated the expression of chemokines and SREBP1, leading to strong recruitment and M1-like polarization of macrophages. Our findings indicate that Mir802 endows adipose tissue with the ability to recruit and polarize macrophages, which underscores Mir802 as an innovative and attractive candidate for miRNA-based immune therapy for adipose inflammation.
-
- Immunology and Inflammation
Chronic antigenic stimulation can trigger the formation of interleukin 10 (IL-10)-producing T-regulatory type 1 (TR1) cells in vivo. We have recently shown that murine T-follicular helper (TFH) cells are precursors of TR1 cells and that the TFH-to-TR1 cell transdifferentiation process is characterized by the progressive loss and acquisition of opposing transcription factor gene expression programs that evolve through at least one transitional cell stage. Here, we use a broad range of bulk and single-cell transcriptional and epigenetic tools to investigate the epigenetic underpinnings of this process. At the single-cell level, the TFH-to-TR1 cell transition is accompanied by both, downregulation of TFH cell-specific gene expression due to loss of chromatin accessibility, and upregulation of TR1 cell-specific genes linked to chromatin regions that remain accessible throughout the transdifferentiation process, with minimal generation of new open chromatin regions. By interrogating the epigenetic status of accessible TR1 genes on purified TFH and conventional T-cells, we find that most of these genes, including Il10, are already poised for expression at the TFH cell stage. Whereas these genes are closed and hypermethylated in Tconv cells, they are accessible, hypomethylated, and enriched for H3K27ac-marked and hypomethylated active enhancers in TFH cells. These enhancers are enriched for binding sites for the TFH and TR1-associated transcription factors TOX-2, IRF4, and c-MAF. Together, these data suggest that the TR1 gene expression program is genetically imprinted at the TFH cell stage.