Alternative platelet differentiation pathways initiated by nonhierarchically related hematopoietic stem cells.

Rare multipotent stem cells replenish millions of blood cells per second through a time-consuming process, passing through multiple stages of increasingly lineage-restricted progenitors. Although insults to the blood-forming system highlight the need for more rapid blood replenishment from stem cells, established models of hematopoiesis implicate only one mandatory differentiation pathway for each blood cell lineage. Here, we establish a nonhierarchical relationship between distinct stem cells that replenish all blood cell lineages and stem cells that replenish almost exclusively platelets, a lineage essential for hemostasis and with important roles in both the innate and adaptive immune systems. These distinct stem cells use cellularly, molecularly and functionally separate pathways for the replenishment of molecularly distinct megakaryocyte-restricted progenitors: a slower steady-state multipotent pathway and a fast-track emergency-activated platelet-restricted pathway. These findings provide a framework for enhancing platelet replenishment in settings in which slow recovery of platelets remains a major clinical challenge.

Single-cell multi-omics identifies chronic inflammation as a driver of TP53-mutant leukemic evolution.

Understanding the genetic and nongenetic determinants of tumor protein 53 (TP53)-mutation-driven clonal evolution and subsequent transformation is a crucial step toward the design of rational therapeutic strategies. Here we carry out allelic resolution single-cell multi-omic analysis of hematopoietic stem/progenitor cells (HSPCs) from patients with a myeloproliferative neoplasm who transform to TP53-mutant secondary acute myeloid leukemia (sAML). All patients showed dominant TP53 'multihit' HSPC clones at transformation, with a leukemia stem cell transcriptional signature strongly predictive of adverse outcomes in independent cohorts, across both TP53-mutant and wild-type (WT) AML. Through analysis of serial samples, antecedent TP53-heterozygous clones and in vivo perturbations, we demonstrate a hitherto unrecognized effect of chronic inflammation, which suppressed TP53 WT HSPCs while enhancing the fitness advantage of TP53-mutant cells and promoted genetic evolution. Our findings will facilitate the development of risk-stratification, early detection and treatment strategies for TP53-mutant leukemia, and are of broad relevance to other cancer types.

Epigenetic programming defines haematopoietic stem cell fate restriction.

Haematopoietic stem cells (HSCs) are multipotent, but individual HSCs can show restricted lineage output in vivo. Currently, the molecular mechanisms and physiological role of HSC fate restriction remain unknown. Here we show that lymphoid fate is epigenetically but not transcriptionally primed in HSCs. In multi-lineage HSCs that produce lymphocytes, lymphoid-specific upstream regulatory elements (LymUREs) but not promoters are preferentially accessible compared with platelet-biased HSCs that do not produce lymphoid cell types, providing transcriptionally silent lymphoid lineage priming. Runx3 is preferentially expressed in multi-lineage HSCs, and reinstating Runx3 expression increases LymURE accessibility and lymphoid-primed multipotent progenitor 4 (MPP4) output in old, platelet-biased HSCs. In contrast, platelet-biased HSCs show elevated levels of epigenetic platelet-lineage priming and give rise to MPP2 progenitors with molecular platelet bias. These MPP2 progenitors generate platelets with faster kinetics and through a more direct cellular pathway compared with MPP2s derived from multi-lineage HSCs. Epigenetic programming therefore predicts both fate restriction and differentiation kinetics in HSCs.

CD34+CD19-CD22+ B-cell progenitors may underlie phenotypic escape in patients treated with CD19-directed therapies.

CD19-directed immunotherapies have revolutionized the treatment of advanced B-cell acute lymphoblastic leukemia (B-ALL). Despite initial impressive rates of complete remission (CR) many patients ultimately relapse. Patients with B-ALL successfully treated with CD19-directed T cells eventually relapse, which, coupled with the early onset of CD22 expression during B-cell development, suggests that preexisting CD34+CD22+CD19- (pre)-leukemic cells represent an "early progenitor origin-related" mechanism underlying phenotypic escape to CD19-directed immunotherapies. We demonstrate that CD22 expression precedes CD19 expression during B-cell development. CD34+CD19-CD22+ cells are found in diagnostic and relapsed bone marrow samples of ∼70% of patients with B-ALL, and their frequency increases twofold in patients with B-ALL in CR after CD19 CAR T-cell therapy. The median of CD34+CD19-CD22+ cells before treatment was threefold higher in patients in whom B-ALL relapsed after CD19-directed immunotherapy (median follow-up, 24 months). Fluorescence in situ hybridization analysis in flow-sorted cell populations and xenograft modeling revealed that CD34+CD19-CD22+ cells harbor the genetic abnormalities present at diagnosis and initiate leukemogenesis in vivo. Our data suggest that preleukemic CD34+CD19-CD22+ progenitors underlie phenotypic escape after CD19-directed immunotherapies and reinforce ongoing clinical studies aimed at CD19/CD22 dual targeting as a strategy for reducing CD19- relapses. The implementation of CD34/CD19/CD22 immunophenotyping in clinical laboratories for initial diagnosis and subsequent monitoring of patients with B-ALL during CD19-targeted therapy is encouraged.

Processing single-cell RNA-seq datasets using SingCellaR.

Single-cell RNA sequencing has led to unprecedented levels of data complexity. Although several computational platforms are available, performing data analyses for multiple datasets remains a significant challenge. Here, we provide a comprehensive analytical protocol to interrogate multiple datasets on SingCellaR, an analysis package in R. This tool can be applied to general single-cell transcriptome analyses. We demonstrate steps for data analyses and visualization using bespoke pipelines, in conjunction with existing analysis tools to study human hematopoietic stem and progenitor cells. For complete details on the use and execution of this protocol, please refer to Roy et al. (2021).

Ezh2 is essential for the generation of functional yolk sac derived erythro-myeloid progenitors.

Yolk sac (YS) hematopoiesis is critical for the survival of the embryo and a major source of tissue-resident macrophages that persist into adulthood. Yet, the transcriptional and epigenetic regulation of YS hematopoiesis remains poorly characterized. Here we report that the epigenetic regulator Ezh2 is essential for YS hematopoiesis but dispensable for subsequent aorta-gonad-mesonephros (AGM) blood development. Loss of EZH2 activity in hemogenic endothelium (HE) leads to the generation of phenotypically intact but functionally deficient erythro-myeloid progenitors (EMPs), while the generation of primitive erythroid cells is not affected. EZH2 activity is critical for the generation of functional EMPs at the onset of the endothelial-to-hematopoietic transition but subsequently dispensable. We identify a lack of Wnt signaling downregulation as the primary reason for the production of non-functional EMPs. Together, our findings demonstrate a critical and stage-specific role of Ezh2 in modulating Wnt signaling during the generation of EMPs from YS HE.

Transitions in lineage specification and gene regulatory networks in hematopoietic stem/progenitor cells over human development.

Human hematopoiesis is a dynamic process that starts in utero 18-21 days post-conception. Understanding the site- and stage-specific variation in hematopoiesis is important if we are to understand the origin of hematological disorders, many of which occur at specific points in the human lifespan. To unravel how the hematopoietic stem/progenitor cell (HSPC) compartment changes during human ontogeny and the underlying gene regulatory mechanisms, we compare 57,489 HSPCs from 5 different tissues spanning 4 developmental stages through the human lifetime. Single-cell transcriptomic analysis identifies significant site- and developmental stage-specific transitions in cellular architecture and gene regulatory networks. Hematopoietic stem cells show progression from cycling to quiescence and increased inflammatory signaling during ontogeny. We demonstrate the utility of this dataset for understanding aberrant hematopoiesis through comparison to two cancers that present at distinct time points in postnatal life-juvenile myelomonocytic leukemia, a childhood cancer, and myelofibrosis, which classically presents in older adults.

Single-cell profiling of human bone marrow progenitors reveals mechanisms of failing erythropoiesis in Diamond-Blackfan anemia.

Ribosome dysfunction underlies the pathogenesis of many cancers and heritable ribosomopathies. Here, we investigate how mutations in either ribosomal protein large (RPL) or ribosomal protein small (RPS) subunit genes selectively affect erythroid progenitor development and clinical phenotypes in Diamond-Blackfan anemia (DBA), a rare ribosomopathy with limited therapeutic options. Using single-cell assays of patient-derived bone marrow, we delineated two distinct cellular trajectories segregating with ribosomal protein genotypes. Almost complete loss of erythroid specification was observed in RPS-DBA. In contrast, we observed relative preservation of qualitatively abnormal erythroid progenitors and precursors in RPL-DBA. Although both DBA genotypes exhibited a proinflammatory bone marrow milieu, RPS-DBA was characterized by erythroid differentiation arrest, whereas RPL-DBA was characterized by preserved GATA1 expression and activity. Compensatory stress erythropoiesis in RPL-DBA exhibited disordered differentiation underpinned by an altered glucocorticoid molecular signature, including reduced ZFP36L2 expression, leading to milder anemia and improved corticosteroid response. This integrative analysis approach identified distinct pathways of erythroid failure and defined genotype-phenotype correlations in DBA. These findings may help facilitate therapeutic target discovery.

The bromodomain inhibitor JQ1+ reduces calcium-sensing receptor activity in pituitary cell lines.

Corticotrophinomas represent 10% of all surgically removed pituitary adenomas, however, current treatment options are often not effective, and there is a need for improved pharmacological treatments. Recently, JQ1+, a bromodomain inhibitor that promotes gene transcription by binding acetylated histone residues and recruiting transcriptional machinery, has been shown to reduce proliferation in a murine corticotroph cell line, AtT20. RNA-Seq analysis of AtT20 cells following treatment with JQ1+ identified the calcium-sensing receptor (CaSR) gene as significantly downregulated, which was subsequently confirmed using real-time PCR and Western blot analysis. CaSR is a G protein-coupled receptor that plays a central role in calcium homeostasis but can elicit non-calcitropic effects in multiple tissues, including the anterior pituitary where it helps regulate hormone secretion. However, in AtT20 cells, CaSR activates a tumour-specific cAMP pathway that promotes ACTH and PTHrP hypersecretion. We hypothesised that the Casr promoter may harbour binding sites for BET proteins, and using chromatin immunoprecipitation (ChIP)-sequencing demonstrated that the BET protein Brd3 binds to the promoter of the Casr gene. Assessment of CaSR signalling showed that JQ1+ significantly reduced Ca2+e-mediated increases in intracellular calcium (Ca2+i) mobilisation and cAMP signalling. However, the CaSR-negative allosteric modulator, NPS-2143, was unable to reduce AtT20 cell proliferation, indicating that reducing CaSR expression rather than activity is likely required to reduce pituitary cell proliferation. Thus, these studies demonstrate that reducing CaSR expression may be a viable option in the treatment of pituitary tumours. Moreover, current strategies to reduce CaSR activity, rather than protein expression for cancer treatments, may be ineffective.

Heterogeneous disease-propagating stem cells in juvenile myelomonocytic leukemia.

Juvenile myelomonocytic leukemia (JMML) is a poor-prognosis childhood leukemia usually caused by RAS-pathway mutations. The cellular hierarchy in JMML is poorly characterized, including the identity of leukemia stem cells (LSCs). FACS and single-cell RNA sequencing reveal marked heterogeneity of JMML hematopoietic stem/progenitor cells (HSPCs), including an aberrant Lin-CD34+CD38-CD90+CD45RA+ population. Single-cell HSPC index-sorting and clonogenic assays show that (1) all somatic mutations can be backtracked to the phenotypic HSC compartment, with RAS-pathway mutations as a "first hit," (2) mutations are acquired with both linear and branching patterns of clonal evolution, and (3) mutant HSPCs are present after allogeneic HSC transplant before molecular/clinical evidence of relapse. Stem cell assays reveal interpatient heterogeneity of JMML LSCs, which are present in, but not confined to, the phenotypic HSC compartment. RNA sequencing of JMML LSC reveals up-regulation of stem cell and fetal genes (HLF, MEIS1, CNN3, VNN2, and HMGA2) and candidate therapeutic targets/biomarkers (MTOR, SLC2A1, and CD96), paving the way for LSC-directed disease monitoring and therapy in this disease.

The quiescent fraction of chronic myeloid leukemic stem cells depends on BMPR1B, Stat3 and BMP4-niche signals to persist in patients in remission.

Chronic myelogenous leukemia arises from the transformation of hematopoietic stem cells by the BCR-ABL oncogene. Though transformed cells are predominantly BCR-ABL-dependent and sensitive to tyrosine kinase inhibitor treatment, some BMPR1B+ leukemic stem cells are treatment-insensitive and rely, among others, on the bone morphogenetic protein (BMP) pathway for their survival via a BMP4 autocrine loop. Here, we further studied the involvement of BMP signaling in favoring residual leukemic stem cell persistence in the bone marrow of patients having achieved remission under treatment. We demonstrate by single-cell RNA-Seq analysis that a sub-fraction of surviving BMPR1B+ leukemic stem cells are co-enriched in BMP signaling, quiescence and stem cell signatures, without modulation of the canonical BMP target genes, but enrichment in actors of the Jak2/Stat3 signaling pathway. Indeed, based on a new model of persisting CD34+CD38- leukemic stem cells, we show that BMPR1B+ cells display co-activated Smad1/5/8 and Stat3 pathways. Interestingly, we reveal that only the BMPR1B+ cells adhering to stromal cells display a quiescent status. Surprisingly, this quiescence is induced by treatment, while non-adherent BMPR1B+ cells treated with tyrosine kinase inhibitors continued to proliferate. The subsequent targeting of BMPR1B and Jak2 pathways decreased quiescent leukemic stem cells by promoting their cell cycle re-entry and differentiation. Moreover, while Jak2-inhibitors alone increased BMP4 production by mesenchymal cells, the addition of the newly described BMPR1B inhibitor (E6201) impaired BMP4-mediated production by stromal cells. Altogether, our data demonstrate that targeting both BMPR1B and Jak2/Stat3 efficiently impacts persisting and dormant leukemic stem cells hidden in their bone marrow microenvironment.

Single-Cell Analyses Reveal Megakaryocyte-Biased Hematopoiesis in Myelofibrosis and Identify Mutant Clone-Specific Targets.

Myelofibrosis is a severe myeloproliferative neoplasm characterized by increased numbers of abnormal bone marrow megakaryocytes that induce fibrosis, destroying the hematopoietic microenvironment. To determine the cellular and molecular basis for aberrant megakaryopoiesis in myelofibrosis, we performed single-cell transcriptome profiling of 135,929 CD34+ lineage- hematopoietic stem and progenitor cells (HSPCs), single-cell proteomics, genomics, and functional assays. We identified a bias toward megakaryocyte differentiation apparent from early multipotent stem cells in myelofibrosis and associated aberrant molecular signatures. A sub-fraction of myelofibrosis megakaryocyte progenitors (MkPs) are transcriptionally similar to healthy-donor MkPs, but the majority are disease specific, with distinct populations expressing fibrosis- and proliferation-associated genes. Mutant-clone HSPCs have increased expression of megakaryocyte-associated genes compared to wild-type HSPCs, and we provide early validation of G6B as a potential immunotherapy target. Our study paves the way for selective targeting of the myelofibrosis clone and illustrates the power of single-cell multi-omics to discover tumor-specific therapeutic targets and mediators of tissue fibrosis.

C/EBPα and GATA-2 Mutations Induce Bilineage Acute Erythroid Leukemia through Transformation of a Neomorphic Neutrophil-Erythroid Progenitor.

Acute erythroid leukemia (AEL) commonly involves both myeloid and erythroid lineage transformation. However, the mutations that cause AEL and the cell(s) that sustain the bilineage leukemia phenotype remain unknown. We here show that combined biallelic Cebpa and Gata2 zinc finger-1 (ZnF1) mutations cooperatively induce bilineage AEL, and that the major leukemia-initiating cell (LIC) population has a neutrophil-monocyte progenitor (NMP) phenotype. In pre-leukemic NMPs Cebpa and Gata2 mutations synergize by increasing erythroid transcription factor (TF) expression and erythroid TF chromatin access, respectively, thereby installing ectopic erythroid potential. This erythroid-permissive chromatin conformation is retained in bilineage LICs. These results demonstrate that synergistic transcriptional and epigenetic reprogramming by leukemia-initiating mutations can generate neomorphic pre-leukemic progenitors, defining the lineage identity of the resulting leukemia.

Cell-intrinsic depletion of Aml1-ETO-expressing pre-leukemic hematopoietic stem cells by K-Ras activating mutation.

Somatic mutations in acute myeloid leukemia are acquired sequentially and hierarchically. First, pre-leukemic mutations, such as t(8;21) that encodes AML1-ETO, are acquired within the hematopoietic stem cell (HSC) compartment, while signaling pathway mutations, including KRAS activating mutations, are late events acquired during transformation of leukemic progenitor cells and are rarely detectable in HSC. This raises the possibility that signaling pathway mutations are detrimental to clonal expansion of pre-leukemic HSC. To address this hypothesis, we used conditional genetics to introduce Aml1-ETO and K-RasG12D into murine HSC, either individually or in combination. In the absence of activated Ras, Aml1-ETO-expressing HSC conferred a competitive advantage. However, activated K-Ras had a marked detrimental effect on Aml1-ETO-expressing HSC, leading to loss of both phenotypic and functional HSC. Cell cycle analysis revealed a loss of quiescence in HSC co-expressing Aml1-ETO and K-RasG12D, accompanied by an enrichment in E2F and Myc target gene expression and depletion of HSC self-renewal-associated gene expression. These findings provide a mechanistic basis for the observed absence of KRAS signaling mutations in the pre-malignant HSC compartment.

Unravelling Intratumoral Heterogeneity through High-Sensitivity Single-Cell Mutational Analysis and Parallel RNA Sequencing.

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for resolving transcriptional heterogeneity. However, its application to studying cancerous tissues is currently hampered by the lack of coverage across key mutation hotspots in the vast majority of cells; this lack of coverage prevents the correlation of genetic and transcriptional readouts from the same single cell. To overcome this, we developed TARGET-seq, a method for the high-sensitivity detection of multiple mutations within single cells from both genomic and coding DNA, in parallel with unbiased whole-transcriptome analysis. Applying TARGET-seq to 4,559 single cells, we demonstrate how this technique uniquely resolves transcriptional and genetic tumor heterogeneity in myeloproliferative neoplasms (MPN) stem and progenitor cells, providing insights into deregulated pathways of mutant and non-mutant cells. TARGET-seq is a powerful tool for resolving the molecular signatures of genetically distinct subclones of cancer cells.

SCL/TAL1 cooperates with Polycomb RYBP-PRC1 to suppress alternative lineages in blood-fated cells.

During development, it is unclear if lineage-fated cells derive from multilineage-primed progenitors and whether active mechanisms operate to restrict cell fate. Here we investigate how mesoderm specifies into blood-fated cells. We document temporally restricted co-expression of blood (Scl/Tal1), cardiac (Mesp1) and paraxial (Tbx6) lineage-affiliated transcription factors in single cells, at the onset of blood specification, supporting the existence of common progenitors. At the same time-restricted stage, absence of SCL results in expansion of cardiac/paraxial cell populations and increased cardiac/paraxial gene expression, suggesting active suppression of alternative fates. Indeed, SCL normally activates expression of co-repressor ETO2 and Polycomb-PRC1 subunits (RYBP, PCGF5) and maintains levels of Polycomb-associated histone marks (H2AK119ub/H3K27me3). Genome-wide analyses reveal ETO2 and RYBP co-occupy most SCL target genes, including cardiac/paraxial loci. Reduction of Eto2 or Rybp expression mimics Scl-null cardiac phenotype. Therefore, SCL-mediated transcriptional repression prevents mis-specification of blood-fated cells, establishing active repression as central to fate determination processes.

Niche-mediated depletion of the normal hematopoietic stem cell reservoir by Flt3-ITD-induced myeloproliferation.

Although previous studies suggested that the expression of FMS-like tyrosine kinase 3 (Flt3) initiates downstream of mouse hematopoietic stem cells (HSCs), FLT3 internal tandem duplications (FLT3 ITDs) have recently been suggested to intrinsically suppress HSCs. Herein, single-cell interrogation found Flt3 mRNA expression to be absent in the large majority of phenotypic HSCs, with a strong negative correlation between Flt3 and HSC-associated gene expression. Flt3-ITD knock-in mice showed reduced numbers of phenotypic HSCs, with an even more severe loss of long-term repopulating HSCs, likely reflecting the presence of non-HSCs within the phenotypic HSC compartment. Competitive transplantation experiments established that Flt3-ITD compromises HSCs through an extrinsically mediated mechanism of disrupting HSC-supporting bone marrow stromal cells, with reduced numbers of endothelial and mesenchymal stromal cells showing increased inflammation-associated gene expression. Tumor necrosis factor (TNF), a cell-extrinsic potent negative regulator of HSCs, was overexpressed in bone marrow niche cells from FLT3-ITD mice, and anti-TNF treatment partially rescued the HSC phenotype. These findings, which establish that Flt3-ITD-driven myeloproliferation results in cell-extrinsic suppression of the normal HSC reservoir, are of relevance for several aspects of acute myeloid leukemia biology.

Single-cell transcriptomics uncovers distinct molecular signatures of stem cells in chronic myeloid leukemia.

Recent advances in single-cell transcriptomics are ideally placed to unravel intratumoral heterogeneity and selective resistance of cancer stem cell (SC) subpopulations to molecularly targeted cancer therapies. However, current single-cell RNA-sequencing approaches lack the sensitivity required to reliably detect somatic mutations. We developed a method that combines high-sensitivity mutation detection with whole-transcriptome analysis of the same single cell. We applied this technique to analyze more than 2,000 SCs from patients with chronic myeloid leukemia (CML) throughout the disease course, revealing heterogeneity of CML-SCs, including the identification of a subgroup of CML-SCs with a distinct molecular signature that selectively persisted during prolonged therapy. Analysis of nonleukemic SCs from patients with CML also provided new insights into cell-extrinsic disruption of hematopoiesis in CML associated with clinical outcome. Furthermore, we used this single-cell approach to identify a blast-crisis-specific SC population, which was also present in a subclone of CML-SCs during the chronic phase in a patient who subsequently developed blast crisis. This approach, which might be broadly applied to any malignancy, illustrates how single-cell analysis can identify subpopulations of therapy-resistant SCs that are not apparent through cell-population analysis.

Distinct myeloid progenitor-differentiation pathways identified through single-cell RNA sequencing.

According to current models of hematopoiesis, lymphoid-primed multi-potent progenitors (LMPPs) (Lin(-)Sca-1(+)c-Kit(+)CD34(+)Flt3(hi)) and common myeloid progenitors (CMPs) (Lin(-)Sca-1(+)c-Kit(+)CD34(+)CD41(hi)) establish an early branch point for separate lineage-commitment pathways from hematopoietic stem cells, with the notable exception that both pathways are proposed to generate all myeloid innate immune cell types through the same myeloid-restricted pre-granulocyte-macrophage progenitor (pre-GM) (Lin(-)Sca-1(-)c-Kit(+)CD41(-)FcγRII/III(-)CD150(-)CD105(-)). By single-cell transcriptome profiling of pre-GMs, we identified distinct myeloid differentiation pathways: a pathway expressing the gene encoding the transcription factor GATA-1 generated mast cells, eosinophils, megakaryocytes and erythroid cells, and a pathway lacking expression of that gene generated monocytes, neutrophils and lymphocytes. These results identify an early hematopoietic-lineage bifurcation that separates the myeloid lineages before their segregation from other hematopoietic-lineage potential.