Key Points
Substantial attrition of HSC function occurs within 24 hours of ex vivo culture independently of cell cycle progression.
Inhibition of JAK/STAT signaling during culture adaptation via ruxolitinib improves HSC function ex vivo.
Visual Abstract
Loss of long-term hematopoietic stem cell (LT-HSC) function ex vivo hampers the success of clinical protocols that rely on culture. However, the kinetics and mechanisms through which this occurs remain incompletely characterized. In this study, through time-resolved single-cell RNA sequencing, matched in vivo functional analysis, and the use of a reversible in vitro system of early G1 arrest, we defined the sequence of transcriptional and functional events that occur during the first ex vivo division of human LT-HSCs. We demonstrated that the sharpest loss in LT-HSC repopulation capacity happens early on, between 6 and 24 hours of culture, before LT-HSCs commit to cell cycle progression. During this time window, LT-HSCs adapt to the culture environment, limit the global variability in gene expression, and transiently upregulate gene networks involved in signaling and stress responses. From 24 hours, LT-HSC progression past early G1 contributes to the establishment of differentiation programs in culture. However, contrary to the current assumptions, we demonstrated that the loss of HSC function ex vivo is independent of cell cycle progression. Finally, we showed that targeting LT-HSC adaptation to culture by inhibiting the early activation of JAK/STAT signaling improves HSC long-term repopulating function ex vivo. Collectively, our study demonstrated that controlling early LT-HSC adaptation to ex vivo culture, for example, via JAK inhibition, is critically important to improve HSC gene therapy and expansion protocols.
Introduction
A trillion blood cells are produced daily in humans. This impressive output is achieved by rare hematopoietic stem cells (HSCs) that have the unique capacity of driving production of all blood cell types (differentiation), while maintaining a functional HSC pool (self-renewal). Owing to this extensive regenerative capacity, HSCs are the key functional units of HSC transplantation and HSC gene therapy (GT). HSC GT promises to be a curative, one-time treatment, and is under clinical investigation to treat more than 10 monogenic diseases.1 HSC GT currently requires an ex vivo culture step, which strongly reduces HSC functionality.2,3 Minimizing HSC attrition during ex vivo protocols would provide patients with larger numbers of functional HSCs, improve treatment efficacy and safety in current HSC GT indications, and pave the way for new therapeutic applications.
Long-term blood formation post-transplantation is driven by predominantly quiescent long-term HSCs (LT-HSCs).4,5 Quiescence (G0) is defined as the reversible absence of cell cycling. G0 is characterized by decreased cell size, reduced protein biosynthesis,6,7 high autophagic recycling,8,9 an elevated basal expression of cell stress response pathways,10-13 a relatively inactive glycolytic metabolism,14,15 and is maintained in vivo by the hypoxic bone marrow niche.16,17 Upon stimulation, HSCs exit quiescence (G0 to early G1) and progressively phosphorylate retinoblastoma (Rb) until the restriction point where commitment to division occurs and cell cycle progression (late G1 to S to G2 to M) ensues. CDK6 is a master regulator of HSC quiescence exit kinetics whose differential regulation within the HSC pool guarantees LT-HSC maintenance.18 An emerging body of work has demonstrated further heterogeneity within the quiescent state of LT-HSCs4 that is associated with distinct lineage preferences and kinetics of reconstitution after transplantation.19-22 Clinical GT protocols target the CD34+ fraction (a heterogeneous mixture of HSCs and progenitor cells) with a limited understanding of how ex vivo protocols impact the underlying biology unique to the LT-HSC population.
In current culture systems, HSCs inevitably exit quiescence and divide. Ex vivo culture increases protein synthesis rates,6,23 remodels the mitochondria to an activated state24 with increased oxidative metabolism and reactive oxygen species production,25 disrupts optimal proteostasis programs,23,26 and reduces the dependency on lysosomal recycling.27 These changes are all associated with a net decline in long-term repopulation capacity. Early studies showed that it is exclusively HSCs and progenitor cells in G0 before28 and after culture29,30 that engraft and not their cycling counterparts. From this and from an extensive body of literature on in vivo models,31-33 it has been assumed that cell cycle progression itself drives the loss of HSC function in culture. However, to date, no study has formally addressed this question in human LT-HSCs that undergo clinically relevant culture. Despite recent advances in strategies to expand HSCs ex vivo34 and the increased clinical adoption of HSC GT, we are still lacking a kinetic and mechanistic understanding of how HSCs lose self-renewal capabilities ex vivo. In this study, by pairing single-cell RNA sequencing (scRNA-seq) with in vivo functional analysis in a time-resolved manner, we dissect the transcriptional and functional changes that occur in highly purified human LT-HSCs during their first division in culture.
Methods
Human samples
Human biological samples were sourced ethically, and the research was conducted in accordance with the terms of the informed consents under an institutional review board/research ethics committee–approved protocol as specified below. Primary cord blood (CB) samples were obtained with informed consent from healthy donors from the Cambridge Blood and Stem Cell Biobank. CBs of both sexes were processed as a single sample, except for scRNA-seq for which only single sex CB pools were used. Mobilized peripheral blood (mPB) was obtained from healthy male donors, aged 25 to 28 years, by administering daily filgrastim (Neupogen; 10 μg/kg per day) for 5 days. Apheresis was performed on day 5 and 6 using the Optia Spectra (Terumo Blood and Cell Technologies). Sample preparation, flow cell sorting, and phenotyping are reported in the supplemental Methods, available on the Blood website.
Ex vivo culture
Two complementary culture systems were used, namely (1) experimental (EXPER) conditions comprising StemPro base media (Thermo Fisher Scientific) supplemented with nutrients (0.028%; Thermo Fisher Scientific), penicillin and streptomycin (1%), l-glutamine (1%), human low-density lipoprotein (50 ng/mL) (Stem Cell Technologies), and the following cytokines: stem cell factor (100 ng/mL), fms-like tyrosine kinase-3 ligand (FLT3L; 20 ng/mL), thrombopoietin (TPO; 100 ng/mL), erythropoietin (EPO; 3 units per mL), interleukin-6 (IL-6; 50 ng/mL), interleukin-3 (IL-3; 10 ng/mL), and granulocyte-macrophage colony-stimulating factor (GM-CSF; 20 ng/mL); and (2) gene therapy (GT) conditions comprising GMP Stem Cell Growth Medium (CellGenix) supplemented with l-glutamine (1%) (Thermo Fisher Scientific), penicillin and streptomycin (1%) (Thermo Fisher Scientific), and the following cytokines: SCF (300 ng/mL), FLT3L (300 ng/mL), IL-3 (60 ng/mL), and TPO (100 ng/mL). All cytokines were from Miltenyi/Peprotech, except EPO (Janssen). The CDK4/CDK6 inhibitor palbociclib (PD; PD0332991; Sigma, 200 nM), the JAK1/JAK2 inhibitor ruxolitinib (RUX; Selleckchem; 5-500 nM), UM171 (Stem Cell Technologies, 35 nM), or Z-VAD(OH)-FMK (Cayman Chemical, 100 nM) were added to GT, EXPER, or media supporting myeloid-erythroid-megakaryocytic (MEM) colony generation when indicated. The methodology for lentiviral vector production and all in vitro assays to assess HSC function, cell cycle, differentiation, and serial replating ability are described in the supplemental Methods.
Xenograft transplantation
NOD.Cg-PrkdcscidIl2rgtm1Wjl/SzJ (NSG) mice or NOD.Cg-Prkdcscid Il2rgtm1Wjl Tg(CMV-IL3,CSF2,KITLG)1Eav/ MloySzJ (NSG-SGM3) mice were bred in-house or obtained from Charles River. All animals were housed in a specific pathogen-free animal facility, and experiments were conducted according to the UK Home Office regulations or in accordance with the institutional guidelines approved by the University Health Network Animal Care Committee. The primary and secondary transplantation methods and analysis criteria are reported in the supplemental Methods.
scRNA-seq
scRNA-seq libraries were prepared using the Smart-Seq2 protocol,35 adapted as described in the supplemental Methods. Libraries were quantified using the KAPA library quantification kit (Roche) and were sequenced by paired-end sequencing (150 bp) using the Illumina HiSeq4000 (Illumina) and NovaSeq 6000 (Illumina) at the Cancer Research UK Cambridge Institute genomics core facility (Cambridge, United Kingdom). The scRNA-seq experimental design and bioinformatic methods are described in the supplemental Methods.
Statistical analysis
Appropriate statistical tests were performed using GraphPad Prism (v9.3), R Studio (v1.2), and Python (v3.8.6), further described in the supplemental Methods.
CB collection and experimental work on CB and mPB CD34+ cells were performed in accordance with regulated procedures approved under the 07/MRE05/44 and 18/EE/0199 research ethics committee research studies. Animal studies were ethically reviewed by the University of Cambridge Animal Welfare and Ethical Review Body and carried out in accordance with the Animals (Scientific Procedures) Act 1986, the GlaxoSmithKline policy on the care, welfare and treatment of animals, and institutional guidelines approved by the University Health Network Animal Care Committee.
Results
Loss of HSC function occurs early during ex vivo culture
HSC functional attrition in culture severely reduces the efficacy of HSC clinical approaches like ex vivo GT. We therefore sought to study the kinetics of this phenomenon using (1) purified human LT-HSCs (CD34+CD19−CD38−CD45RA−CD90+CD49f+),4 (2) early culture time points (≤72 hours) corresponding to key HSC cell cycle transitions,18 (3) 2 distinct ex vivo systems and HSC sources, namely “GT_mPB conditions” using mPB LT-HSCs and 2 hits of transduction with a lentiviral vector expressing GFP2,36,37; and “EXPER_CB conditions” using CB LT-HSCs (unless otherwise indicated) and promoting differentiation.18
We first characterized LT-HSC cell cycle kinetics in our ex vivo conditions. LT-HSCs in either the GT or EXPER conditions progressed past late G1 at ∼24 hours after culture initiation (Figure 1A-B) as assessed by both phosphorylation of Rb at the serine residue 807 to 811 and acquisition of the Ki-67 marker (supplemental Figure 1A-B). Concurrently, the proportion of cycling cells (S to G2 to M phases) increased within the 24- to 72-hour window (Figure 1C-D). LT-HSC time to first division was estimated to be 62.0 ± 6.5 hours for the GT_mPB conditions and 53.4 ± 6.5 hours for the EXPER_CB conditions (Figure 1E-G). Therefore, independent of HSC source and culture conditions, LT-HSCs progress to late G1 by ∼24 hours and their first division by 72 hours.
Next, we quantified the long-term repopulating cell frequency (%LTRC) at 0, 6, 24, and 62 or 72 hours of ex vivo culture and performed limiting dilution analysis (LDA) xenotransplantation in NSG mice. We transplanted fixed numbers of (1) GT_mPB CD34+ CD38− cells cultured for 62 hours (time of infusion into patients in HSC GT protocols) and (2) EXPER_CB LT-HSCs cultured for up to 72 hours. Human engraftment was analyzed 18 weeks after transplantation. Of note, in GT conditions, almost all engrafted animals displayed GFP+ cells (10/11; supplemental Table 1), indicating that transduced cells contributed to long-term engraftment. LDA showed that the %LTRC within the transplanted LT-HSC population was unchanged for the first 6 hours of culture. However, the %LTRC dropped conspicuously by 24 hours (GT, approximately threefold; P = .003; EXPER, approximately fivefold; P = .00015 when compared with 0 hours). An additional decrease in the %LTRC was observed past 24 hours in the EXPER conditions (approximately twofold; P = .0129 when compared to 24 hours) when LT-HSCs progressed through the cell cycle (Figure 1H-I; supplemental Tables 1 and 2). We next measured apoptosis using Annexin V/7-AAD staining. A decrease in viable cells was observed between 6 to 24 hours with no further decrease at the later time point in both the GT_mPB and EXPER_CB culture systems (Figure 1J-K). Interestingly, the loss in LTRC observed between 6 to 24 hours ex vivo largely outnumbered the loss in viability measured in the same time window. We therefore concluded that HSC loss ex vivo occurs predominantly before cells enter the late G1 phase of the cell cycle.
Transcriptome dynamics of LT-HSCs over the first cell cycle ex vivo
To better understand the kinetics of LTRC loss ex vivo, we investigated the transcriptional changes associated with the first cell division of LT-HSCs in culture. EXPER_CB LT-HSCs were cultured for 0, 6, 24, and 72 hours before scRNA-seq via an adapted Smart-Seq2 protocol.35 A total of 429 cells, collected over 2 independent experiments, passed quality control. Cells from both experiments were integrated through 2 distinct pipelines (Scanpy and Seurat 4) (supplemental Methods), which yielded highly concordant results. LT-HSCs grouped by culture duration in uniform manifold approximation and projection (UMAP) visualizations independently of the batch (Figure 2A; supplemental Figure 2A), of whether cell cycle regression was applied (supplemental Figure 2B), and of the integration pipeline used (supplemental Figure 2C).
At 0 hours, we identified 2 transcriptionally defined subsets of LT-HSCs in accordance with recent reports of transcriptional heterogeneity within the quiescent human LT-HSC fraction (supplemental Note 1; supplemental Figure 2D-H; supplemental Table 3). After 0 hours, pseudotime ordering largely recapitulated chronological time (Figure 2B; supplemental Figure 2I), and transcriptional assignment of cell cycle status agreed with the functional analysis (Figure 2C). LT-HSCs therefore progressed through distinct transcriptional states over time in culture. DESeq2 differential gene expression analysis of all pairwise comparisons in the data set (supplemental Table 4) showed that ex vivo culture has an extensive and rapid effect on the transcriptome. A total of 10 010 genes changed over the time course (hereafter termed “ex vivo modulated genes”; supplemental Table 5). Interestingly, of the gene expression changes that occurred between 0 and 72 hours (5980 genes), 75% (4460 genes) were already observed at 6 hours (Figure 2D) before the vast majority (>90%) of the cells entered late G1 in these culture conditions (Figure 1A).
We then analyzed ex vivo modulated genes using degPatterns38 to define gene (Figure 2E; supplemental Figure 3A; supplemental Table 6) and pathway (gene set variation analysis (GSVA)39; supplemental Figure 3B; supplemental Table 6) expression patterns over the time course. Similar dynamic patterns of upregulation and downregulation were observed at the gene and pathway level (supplemental Figure 3A-B). We found 5 expression patterns of interest, comprising (1) cumulative increases in the expression of pathways linked to transcription, RNA splicing, translation, and oxidative phosphorylation (Figure 2F); (2) progressive decreases in pathways linked to cell adhesion (Figure 2G); (3) early loss of the activator protein-1 transcription factors JUN and FOS, key to the HSC quiescence network (supplemental Figure 3C); and (4) transient patterns of gene (supplemental Figure 3D) or pathway (Figure 2H-I) expression indicating that HSCs initiate a stress response to adapt to culture signals. Particularly abundant were genes/pathways upregulated from 0 to 6 hours and subsequently decreased at 24 hours (27.3% of the total genes; 76.3% of “transient UP” genes). These included sphingolipid de novo biosynthesis (eg, DEGS1), apoptosis (eg, CFLAR and BIRC2), and stress pathways (eg, ATF4). The fifth pattern was upregulation of expression after the 6 hour time point of genes or gene sets related to proliferation and differentiation (Figure 2J-K). These included MYC (supplemental Figure 3E), MYC targets, cell cycle progression, signatures of active human LT-HSCs,22 and committed progenitors.40 Finally, we calculated scEntropy scores, which decrease during differentiation fate decisions.41 scEntropy scores significantly decreased only after 6 hours (Figure 2L), supporting the concept that cell-fate decisions occur at the 6 to 24 hour transition.
Dynamic changes in gene expression variability have been linked to cell-fate decisions in hematopoietic progenitors42 and T-cell activation.43 We therefore investigated gene expression variability over the time course, employing BASiCS 44,45 and determining maximally variable genes (MVGs) at each time point (supplemental Methods; supplemental Figure 3F). The number of MVGs was lowest at 6 hours (Figure 2M) showing that gene expression variability is restrained early in culture initiation. This is likely attributable to the adaptation of all LT-HSC subsets to the culture environment. MVGs at 6 hours were enriched for unfolded protein response regulatory genes (including DNAJC3 and XBP1; Figure 2N, supplemental Table 7), indicating that despite a global restriction in gene expression variability, LT-HSCs vary in their degree of activation of cell stress response genes at this critical time point.
In summary, our results indicate that, within the first 24 hours, LT-HSCs initially adapt to culture by restricting global gene variability and activating a transient stress response. Subsequently, MYC and cell cycle genes are upregulated and differentiation programs are initiated.
Cell cycle–dependent initiation of differentiation in cultured LT-HSCs
Ex vivo, changes to HSC metabolism and organelle biology occur concomitantly with cell cycle progression and loss of function.23-25,27,46 Our observation that the sharpest drop in LT-HSC repopulation capacity occurs before most LT-HSCs enter late G1 prompted us to formally test to what extent progression past early G1 contributes to HSC identity and function in culture. We took advantage of the CDK4/CDK6 inhibitor PD (PD033299), previously established to prevent division of CB LT-HSCs.18 First, we repeated this finding with LT-HSCs cultured in both the GT_mPB and EXPER_CB systems (Figure 3A-B), showing early G1 arrest of CB (Figure 3C-D) and mPB LT-HSCs (supplemental Figure 4A). PD treatment was reversible and did not compromise LT-HSC differentiation (supplemental Figure 4B-E), therefore allowing assessment of the effect of culture in the absence of cell cycle progression.
To disentangle which transcriptional differences in culture depend on cell cycle progression, scRNA-seq was performed on LT-HSCs treated with PD for 24 and 72 hours in the EXPER_CB system and for 62 hours in the GT_mPB system. These data were integrated with 0- and 6-hours cells in a unique embedding using Seurat4 (954 cells; n = 6 independent experiments). As expected, PD-treated LT-HSCs were transcriptionally allocated to G1 in both culture systems (supplemental Figure 4F-G). Cultured LT-HSCs pharmacologically arrested in G1, localized with the 24-hour untreated (UNTR) LT-HSCs in the UMAP (Figure 3E) and by pseudotime (Figure 3F) independently of the cell source, culture system, duration of PD treatment, and whether cell cycle regression was applied or not (supplemental Figure 4H-I). The impact of PD treatment on the transcriptome was relatively minimal with gene expression correlations being higher when comparing UNTR with PD conditions (Pearson correlation coefficient; EXPER_CB 24 hours, 0.96; EXPER_CB 72 hours, 0.871; GT_mPB 62 hours, 0.881) than when comparing any UNTR time point comparisons (Figure 3G-H). Furthermore, most pathways that were significantly enriched by PD treatment either at 24 hours or at 72 hours were related to cell cycle progression (Figure 3I; supplemental Table 8), demonstrating the low off-target effects of PD.
It remains unclear if progression past early G1 contributes to the initiation of differentiation in HSCs as it does in other stem cells.47-49 Transcriptional priming of differentiation occurs from 24 hours of culture onward as HSCs begin cycling (Figure 2K). Interestingly, when compared with UNTR conditions, several myelo-erythroid progenitor signatures were significantly downregulated in PD-treated LT-HSCs (Figure 4A) with genes linked to neutropoiesis (PRTN3, ELANE, and CTSG50) (supplemental Figure 5H) and erythropoiesis (TPI151) (supplemental Table 8) being significantly reduced. Given that temporary cell cycle blockade with PD did not impact in vitro terminal differentiation outcomes (supplemental Figure 4B-E), we interpreted these data to be a delay in lineage specification rather than a change in lineage instruction. We concluded that S to G2 to M progression facilitates the transcriptional onset of differentiation in ex vivo HSCs.
Loss of long-term repopulation capacity ex vivo is independent of cell cycle progression
Next, we tested whether the inhibition of cell cycle progression impacts HSC functional hallmarks. An increase in cell size (Figure 4B) and mitochondrial activity (Figure 4C), both of which are HSC activation hallmarks,15,25,46 were observed during EXPER_CB LT-HSC culture and were unchanged upon PD-induced G1 arrest, demonstrating that cell growth and mitochondrial metabolism are independent of cell cycle progression.
To formally determine if progression past early G1 contributed to the loss of long-term repopulation capacity, we calculated the %LTRC in ex vivo culture with PD, using a similar LDA setup as in Figure 1 (Figure 4D; supplemental Tables 1 and 9). LT-HSCs treated with PD exhibited loss of function that was indistinguishable from UNTR controls. At 18 weeks after transplantation, we observed no difference in the graft size, nor in the lineage output, in mice transplanted with EXPER_CB cells cultured with PD for 24 hours (supplemental Figure 5B-C), 72 hours (Figure 4E; supplemental Figure 5D), or GT_mPB cells cultured with PD for 62 hours (Figure 4F; supplemental Figure 5E) with respect to UNTR control cells cultured in the same conditions. LDA found no difference in the %LTRC observed between PD-treated and UNTR cells cultured for any duration in either of the EXPER_CB conditions (24 hours P = .405; 72 hours P = .426; Figure 4G; supplemental Tables 1 and 9). Secondary transplantation experiments revealed that the frequency of serially transplantable HSCs in the injected population was comparable between PD-treated and UNTR cells for both the EXPER_CB (Figure 4H; P = .190) and GT_mPB (Figure 4I; P = .860) conditions (supplemental Table 10) with no significant difference in the graft size observed (supplemental Figure 5F-G). Because PD does not compromise LT-HSC function (supplemental Figure 4B-E), our transplantation experiments demonstrated that cell cycle progression past early G1 was not the main driver of loss of human LT-HSC function in culture.
Inhibition of JAK/STAT signaling during adaptation improves LT-HSC self-renewal
We next hypothesized that inhibition of signaling pathways induced at adaptation may reduce ex vivo HSC functional attrition. JAK/STAT signaling stood out as a potential candidate, because JAK/STAT target genes were transiently upregulated 6 hours after culture initiation (Figure 5A). Consistently, at 30 minutes of GT culture, mPB HSC/multipotent progenitors (HSC/MPP; CD34+CD38−CD45RA−) significantly increased the phosphorylation of STAT5, STAT3, and, to a lesser extent, STAT1 (supplemental Figure 6A-B; supplemental Table 11). STAT phosphorylation could be decreased by the addition of RUX, a US Food and Drug Administration (FDA) approved JAK inhibitor, to the culture (supplemental Figure 6B). To determine the effect of RUX treatment in HSC culture systems, we assessed the differentiative and proliferative output of RUX-treated mPB LT-HSCs cultured in MEM medium for 14 days. The absolute cell numbers per well decreased significantly with increasing doses of RUX (supplemental Figure 6C), but the proportions of erythroid progenitors or myeloid cells did not change significantly, particularly at concentrations ≤50 nM (supplemental Figure 6D). RUX treatment therefore decreased differentiation globally with no major impact on lineage commitment, suggesting it may maintain LT-HSCs ex vivo.
To test whether RUX addition ex vivo could improve the regenerative function of cultured HSCs, we performed serial replating of mPB HSC/MPPs in both GT and EXPER conditions as an in vitro surrogate for self-renewal. By tertiary replating, in which colonies are largely generated by LT-HSCs, substantial increases in the number of colonies were observed at low doses of RUX in both media and across all donors (GT, 3-fold to 38-fold increase; EXPER, 3-fold to 37-fold; Figure 5B-D; supplemental Table 12). At 500 nM RUX, the positive effect was dampened or abrogated, likely because of decreased proliferation (supplemental Figure 6C) and potential off-target inhibition of tropomyosin receptor kinase B (TRK-β) activity.52,53 The number of colonies produced by 10 nM RUX-treated HSC/MPPs was still significantly lower than those generated by uncultured (0 hour) HSC/MPPs (Figure 5E; supplemental Table 12), indicating that RUX reduces the loss of functional HSCs induced during ex vivo adaptation instead of expanding the absolute number of initial HSCs.
Given the clinical need to develop robust human HSC expansion methods ex vivo, we performed serial replating to benchmark RUX action against previously reported strategies. In this assay, RUX (10 nM) treatment generated significantly higher numbers of tertiary colonies than treatment with UM171 (35 nM), a pyrimidoindole derivate in clinical trials for HSC expansion applications54-56 (RUX, 5-fold to 16-fold when compared with DMSO; UM171, 0.7-fold to 8-fold when compared with DMSO) (Figure 5F; supplemental Table 12). Interestingly, RUX and UM171 did not synergize in this assay (Figure 5F). RUX (10 nM) treatment also outperformed lowering TPO concentration, a strategy beneficial in human HSC expansion cultures57 (0.3-fold to 8-fold when compared with DMSO; Figure 5F). Finally, a pan-caspase inhibitor (Z-VAD(OH)-FMK, 100 nM), previously used for CD34+ expansion,58 only led to a modest increase in the tertiary replating ability of mPB HSC/MPPs (1.5-fold to 3-fold when compared with DMSO; Figure 5G; supplemental Table 12).
To formally assess whether ex vivo RUX treatment can increase the long-term in vivo repopulation ability of cultured LT-HSCs, we treated mPB CD34+CD38− with 10 nM RUX for 62 hours and performed serial transplantations with an LDA design as in Figure 4 (Figure 5H). No difference in graft size was observed at 8 weeks in PB (supplemental Figure 6E). However, at 18 weeks after the transplant, the grafts of mice transplanted with high cell doses of RUX-treated CD34+CD38− were significantly larger (P = .0059; Figure 5I; supplemental Figure 6F) than those of control mice. Although the %LTRC determined by LDA statistics were similar after primary transplantation (supplemental Table 13), we observed an approximately threefold increase in the %LTRC in RUX-treated CD34+CD38− when compared with DMSO treatment after secondary transplantation (P = .0491; Figure 5J; supplemental Figure 6G; supplemental Table 14). These data collectively demonstrate that ex vivo JAK/STAT signaling inhibition via RUX improves HSC regenerative function.
Discussion
Our study establishes the sequence of transcriptional and functional events that occur in human LT-HSCs during their first ex vivo division. Our data demonstrate that the early phases of LT-HSC adaptation to culture are critically detrimental to their function, irrespective of cell cycle progression. Our data also indicate the potential of targeting key pathways induced at adaptation, such as JAK/STAT signaling, to improve the safety and efficacy of any clinical protocol involving HSC culture.
It has long been assumed that division causes loss of LT-HSC function ex vivo. In vivo, there is a causative relationship between excessive proliferation of LT-HSCs and loss of self-renewal (reviewed in Johnson et al59). In addition, self-renewal capacity is enriched in vivo in G0 HSCs when compared with those in G1 or other cell cycle phases.28,59 Ex vivo, the G0/G1 fraction of cultured CD34+ is enriched in phenotypic HSC and contains the highest frequency of repopulating cells.29,30 Using a unique experimental system that efficiently and reversibly arrests LT-HSC progression past early G1 ex vivo, we demonstrate that cell cycle progression and division do not drive the sizable loss of LT-HSC repopulation capacity observed ex vivo over 62 or 72 hours. These findings allay concerns that gene editing strategies that require progression through the S phase may be associated with poor engraftment.60
Differentiation is tightly connected to cell cycle in most stem cell types, including HSCs,47-49,59,61-63 in which CDKs and cyclins control HSC fate regulators through direct phosphorylation or chromatin binding.64-67 We observe that temporary cell cycle blockade dampens or delays the upregulation of myelo-erythroid lineage specification gene signatures, indicating that cell cycle progression facilitates the establishment of LT-HSC differentiation programs. HSC self-renewal and lineage commitment are controlled independently in vitro and in vivo in mice.68-70 Our study of cultured human HSCs further supports this concept by temporally uncoupling loss of repopulation capacity (occurring pre-Rb phosphorylation) from establishment of differentiation (occurring post-Rb phosphorylation).
Given the extensive transcriptional rewiring that occurs within 24 hours of culture, HSC loss of function in this time window is certainly multifactorial. Adaptation, through its characteristic transient expression dynamics and high intercellular variability of stress-related genes, likely purges functional HSCs through differentiation or death, akin to the preferential culling of unfit HSCs that occurs in vivo in response to oxidative, genotoxic, and proteostatic damage.10-13,26,71,72 Early metabolic adaptations also likely contribute to later LT-HSC functional loss. In this study, we observed transient upregulation of the ceramide biosynthesis enzyme, DEGS1, the inhibition of which preserves LT-HSC repopulation capacity after ex vivo culture.26 Finally, we demonstrate that early induction of JAK/STAT signaling is detrimental to HSC function during culture.
The availability of large numbers of functionally fit HSCs remains a major limiting factor for most HSC transplantations and GT. Our data identify the FDA–approved JAK1/2 inhibitor, as an effective compound to mitigate loss of HSC function during ex vivo culture. Independent work from our groups demonstrate that tyrosine-unphosphorylated STAT5 (uSTAT5) is a critical regulator of self-renewal vs differentiation decisions in mouse HSCs.73,74 RUX shifts the uSTAT5/pSTAT5 ratio toward increased uSTAT5, promoting mouse HSC self-renewal.74 We propose this same mechanism underlies the beneficial effect of RUX on cultured human HSCs. This could be complemented by direct or indirect action of RUX on the MAPK,75 PI3K,57,76 NF-κB77 and TGF-β78 pathways. More broadly, by pinpointing irreversible HSC fate decisions to the first 24 hours of culture, our work identifies an early functional bottleneck and an untapped window of opportunity to preserve HSC function in any application that requires an ex vivo step. Although shortening HSC GT protocols may be useful,79 our data indicate that it should be complemented with strategies to minimize early HSC functional attrition. Similarly, targeting early HSC adaptation could further improve HSC expansion, complementing current protocols optimized to expand HSC over prolonged culture periods.55-57,80
Acknowledgments
The authors thank the cord blood (CB) and mobilized peripheral blood donors for their kind donation of human tissue, the Cambridge Blood and Stem Cell Biobank, specifically Joanna Baxter and the nurses who obtained consent and collected CB samples; the National Institute for Health Research Cambridge Biomedical Research Centre Cell Phenotyping Hub for their flow cytometry services and advice, Reiner Schulte and Gabriela Grondys-Kotarba at the Cambridge Institute for Medical Research Flow Cytometry Core Facility for their assistance with cell sorting, and the Cancer Research UK Cambridge Institute genomics center for sequencing. The authors also thank Diego Lopez, Julia Sanchez Viladevall, and Benjy Tan Jek Yang for proofreading the manuscript.
This work was supported by a Wellcome Trust and Royal Society Sir Henry Dale Fellowship (number 107630/Z/15/Z [E.L.]), the European Hematology Association Non-Clinical Research Fellowship Award (RG20 [E.L.]), the UK Research and Innovation Frontier Research Grant (EP/Y026586/1 [E.L.]), core support grants from Wellcome and the Medical Research Council (MRC) to the Wellcome-MRC Cambridge Stem Cell Institute (203151/Z/16/Z [E.L., B.G., and A.R.G.]), GlaxoSmithKline-Varsity Alliance (N.F., E.L., and B.G.), the MRC iCASE PhD studentship (1942750 [C.S.J.]), the Cancer Research UK Cambridge Cancer Centre (CTRQQR-2021∖100012 [S. Belluschi and G.M.]), the Deutsche Forschungsgemeinschaft Research Fellowship (ME 5209/1-1 [N.M.]), the Cancer Research UK (C1163/A21762 [B.G. and A.R.G.]), Blood Cancer UK (18002 [B.G.]), Gates Cambridge Scholarship, Gates Cambridge Trust (M.D.), Alborada Trust (A.R.G.), William B. Harrison Foundation (A.R.G.), and the University of Toronto Medicine by Design initiative with funding from the Canada First Research Excellence Fund (J.E.D.). This research was funded in whole, or in part, by the Wellcome Trust.
Authorship
Contribution: E.L., N.F., C.S.J., M.W., S. Belluschi, and B.G. conceptualized the study; C.S.J., S. Belluschi, W.W.Y.L., and N.K.W. developed the methodology; C.S.J., M.W., S. Belluschi, W.M., W.W.Y.L., K.B.K., G.K., N.M., E.F.C., J.M., G.M., C.G.-W., S.S., S. Basheer, J.L., M.D., and N.F. performed the investigation; K.S., X.W., C.S.J., S. Belluschi, and E.L. were responsible for visualization; C.S.J., M.W., G.M., K.S., X.W., E.D., and E.L. conducted formal analysis and statistics; K.S. and X.W. performed data curation; E.L., K.S., and C.S.J. were responsible for validation; N.F., S.J.H., and C.B. obtained resources; C.S.J., K.S., and E.L. wrote the original draft; C.S.J., M.W., W.M., K.S., E.L., and A.R.G. participated in reviewing and editing the manuscript; E.L., N.F., B.G., J.E.D., C.B., N.K.W., and A.R.G. supervised the project; E.L. and N.F. managed the project; and E.L., N.F., and A.R.G. were responsible for funding acquisition.
Conflicts-of-interest disclosure: E.L. reports receiving research funds from GlaxoSmithKline. A.R.G. and J.L. report serving as a consulting for Incyte. The remaining authors declare that they have no competing interests.
Correspondence: Elisa Laurenti, Department of Haematology, University of Cambridge, Cambridge CB2 0AW, United Kingdom; email: [email protected].
References
Author notes
C.S.J., M.W., K.S., and S. Belluschi contributed equally to this study.
Sequencing files and metadata associated to Smart-Seq2 data sets are deposited at Gene Expression Omnibus (accession numbers GSE213365 [time-course of EXPER_CB LT-HSCs] and GSE213370 [time-course of GT_mPB LT-HSCs] under the super series GSE213372). All codes are publicly available on GitHub at https://github.com/elisa-laurenti/LT-HSC_exvivo.
The data sets presented here are also browsable at the following URLs: https://bioinf.stemcells.cam.ac.uk/shiny/laurenti/LT_HSC_exvivo_CB/; https://bioinf.stemcells.cam.ac.uk/shiny/laurenti/LT_HSC_exvivo_CB_PD/; and https://bioinf.stemcells.cam.ac.uk/shiny/laurenti/LT_HSC_exvivo_mPB_PD/.
All data needed to evaluate the presented conclusions are in the main text/supplemental Materials. All other data/materials will be distributed by the authors upon request ([email protected]).
The online version of this article contains a data supplement.
There is a Blood Commentary on this article in this issue.
The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.
This feature is available to Subscribers Only
Sign In or Create an Account Close Modal