Single-cell analysis of immune recognition in chronic myeloid leukemia patients following tyrosine kinase inhibitor discontinuation.

Immunological control of residual leukemia cells is thought to occur in patients with chronic myeloid leukemia (CML) that maintain treatment-free remission (TFR) following tyrosine kinase inhibitor (TKI) discontinuation. To study this, we analyzed 55 single-cell RNA and T cell receptor (TCR) sequenced samples (scRNA+TCRαß-seq) from patients with CML (n = 13, N = 25), other cancers (n = 28), and healthy (n = 7). The high number and active phenotype of natural killer (NK) cells in CML separated them from healthy and other cancers. Most NK cells in CML belonged to the active CD56dim cluster with high expression of GZMA/B, PRF1, CCL3/4, and IFNG, with interactions with leukemic cells via inhibitory LGALS9-TIM3 and PVR-TIGIT interactions. Accordingly, upregulation of LGALS9 was observed in CML target cells and TIM3 in NK cells when co-cultured together. Additionally, we created a classifier to identify TCRs targeting leukemia-associated antigen PR1 and quantified anti-PR1 T cells in 90 CML and 786 healthy TCRß-sequenced samples. Anti-PR1 T cells were more prevalent in CML, enriched in bone marrow samples, and enriched in the mature, cytotoxic CD8 + TEMRA cluster, especially in a patient maintaining TFR. Our results highlight the role of NK cells and anti-PR1 T cells in anti-leukemic immune responses in CML.

Concurrent stem- and lineage-affiliated chromatin programs precede hematopoietic lineage restriction.

The emerging notion of hematopoietic stem and progenitor cells (HSPCs) as a low-primed cloud without sharply demarcated gene expression programs raises the question on how cellular-fate options emerge and at which stem-like stage lineage priming is initiated. Here, we investigate single-cell chromatin accessibility of Lineage-, cKit+, and Sca1+ (LSK) HSPCs spanning the early differentiation landscape. Application of a signal-processing algorithm to detect transition points corresponding to massive alterations in accessibility of 571 transcription factor motifs reveals a population of LSK FMS-like tyrosine kinase 3 (Flt3)intCD9high cells that concurrently display stem-like and lineage-affiliated chromatin signatures, pointing to a simultaneous gain of both lympho-myeloid and megakaryocyte-erythroid programs. Molecularly and functionally, these cells position between stem cells and committed progenitors and display multi-lineage capacity in vitro and in vivo but lack self-renewal activity. This integrative molecular analysis resolves the heterogeneity of cells along hematopoietic differentiation and permits investigation of chromatin-mediated transition between multipotency and lineage restriction.

A Combinatorial Single-cell Approach to Characterize the Molecular and Immunophenotypic Heterogeneity of Human Stem and Progenitor Populations.

Immunophenotypic characterization and molecular analysis have long been used to delineate heterogeneity and define distinct cell populations. FACS is inherently a single-cell assay, however prior to molecular analysis, the target cells are often prospectively isolated in bulk, thereby losing single-cell resolution. Single-cell gene expression analysis provides a means to understand molecular differences between individual cells in heterogeneous cell populations. In bulk cell analysis an overrepresentation of a distinct cell type results in biases and occlusions of signals from rare cells with biological importance. By utilizing FACS index sorting coupled to single-cell gene expression analysis, populations can be investigated without the loss of single-cell resolution while cells with intermediate cell surface marker expression are also captured, enabling evaluation of the relevance of continuous surface marker expression. Here, we describe an approach that combines single-cell reverse transcription quantitative PCR (RT-qPCR) and FACS index sorting to simultaneously characterize the molecular and immunophenotypic heterogeneity within cell populations. In contrast to single-cell RNA sequencing methods, the use of qPCR with specific target amplification allows for robust measurements of low-abundance transcripts with fewer dropouts, while it is not confounded by issues related to cell-to-cell variations in read depth. Moreover, by directly index-sorting single-cells into lysis buffer this method, allows for cDNA synthesis and specific target pre-amplification to be performed in one step as well as for correlation of subsequently derived molecular signatures with cell surface marker expression. The described approach has been developed to investigate hematopoietic single-cells, but have also been used successfully on other cell types. In conclusion, the approach described herein allows for sensitive measurement of mRNA expression for a panel of pre-selected genes with the possibility to develop protocols for subsequent prospective isolation of molecularly distinct subpopulations.

Single-cell molecular analysis defines therapy response and immunophenotype of stem cell subpopulations in CML.

Understanding leukemia heterogeneity is critical for the development of curative treatments as the failure to eliminate therapy-persistent leukemic stem cells (LSCs) may result in disease relapse. Here we have combined high-throughput immunophenotypic screens with large-scale single-cell gene expression analysis to define the heterogeneity within the LSC population in chronic phase chronic myeloid leukemia (CML) patients at diagnosis and following conventional tyrosine kinase inhibitor (TKI) treatment. Our results reveal substantial heterogeneity within the putative LSC population in CML at diagnosis and demonstrate differences in response to subsequent TKI treatment between distinct subpopulations. Importantly, LSC subpopulations with myeloid and proliferative molecular signatures are proportionally reduced at a higher extent in response to TKI therapy compared with subfractions displaying primitive and quiescent signatures. Additionally, cell surface expression of the CML stem cell markers CD25, CD26, and IL1RAP is high in all subpopulations at diagnosis but downregulated and unevenly distributed across subpopulations in response to TKI treatment. The most TKI-insensitive cells of the LSC compartment can be captured within the CD45RA- fraction and further defined as positive for CD26 in combination with an aberrant lack of cKIT expression. Together, our results expose a considerable heterogeneity of the CML stem cell population and propose a Lin-CD34+CD38-/lowCD45RA-cKIT-CD26+ population as a potential therapeutic target for improved therapy response.