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

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

Involvement of MAF/SPP1 axis in the development of bone marrow fibrosis in PMF patients

Primary myelofibrosis (PMF) is a myeloproliferative neoplasm characterized by hyperplastic megakaryopoiesis and myelofibrosis. We recently described the upregulation of MAF (v-maf avian musculoaponeurotic fibrosarcoma oncogene homolog) in PMF CD34+ hematopoietic progenitor cells (HPCs) compared to healthy donor. Here we demonstrated that MAF is also upregulated in PMF compared with the essential thrombocytemia (ET) and polycytemia vera (PV) HPCs. MAF overexpression and knockdown experiments shed some light into the role of MAF in PMF pathogenesis, by demonstrating that MAF favors the megakaryocyte and monocyte/macrophage commitment of HPCs and leads to the increased expression of proinflammatory and profibrotic mediators. Among them, we focused our further studies on SPP1 and LGALS3. We assessed SPP1 and LGALS3 protein levels in 115 PMF, 47 ET and 24 PV patients plasma samples and we found that SPP1 plasma levels are significantly higher in PMF compared with ET and PV patients. Furthermore, in vitro assays demonstrated that SPP1 promotes fibroblasts and mesenchymal stromal cells proliferation and collagen production. Strikingly, clinical correlation analyses uncovered that higher SPP1 plasma levels in PMF patients correlate with a more severe fibrosis degree and a shorter overall survival. Collectively our data unveil that MAF overexpression contributes to PMF pathogenesis by driving the deranged production of the profibrotic mediator SPP1

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 upregulation 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

Characterisation and discovery of novel miRNAs and moRNAs in JAK2V617F mutated SET2 cells.

In order to gain insights into a possible role of microRNAs in myeloproliferative neoplasms, we performed short RNAs massive sequencing and extensive bioinformatic analysis in the JAK2V617F mutated SET2 cell line. Overall, 652 known mature miRNAs were detected, of which 21 were highly expressed, thus being responsible of most of miRNA-mediated gene repression. microRNA putative targets were enriched in specific signaling pathways, providing information about cell activities under massive post-transcriptional regulation. The majority of miRNAs were mixtures of sequence variants, called isomiRs, mainly due to alternative, non-canonical processing of hairpin precursors. We also identified 78 novel miRNAs (miRNA*) derived from known hairpin precursors. Both major and minor (*) form of miRNAs were expressed concurrently from half of expressed hairpins, highlighting the relevance of miRNA* and the complexity of strand selection bias regulation. Finally, we discovered that SET2 cells express a number of miRNA-offset RNAs (moRNAs), short RNAs derived from genomic regions flanking mature miRNAs. We provide novel data about the possible origin of moRNAs, while their functional role remains to be elucidated. Overall, this study shed light on the complexity of microRNA-mediated gene regulation in SET2 cells and represents the basis for future studies in JAK2V617F mutated cellular models