Skip to main content

Single-cell map of diverse immune phenotypes in the acute myeloid leukemia microenvironment

Abstract

Background

Knowledge of immune cell phenotypes, function, and developmental trajectory in acute myeloid leukemia (AML) microenvironment is essential for understanding mechanisms of evading immune surveillance and immunotherapy response of targeting special microenvironment components.

Methods

Using a single-cell RNA sequencing (scRNA-seq) dataset, we analyzed the immune cell phenotypes, function, and developmental trajectory of bone marrow (BM) samples from 16 AML patients and 4 healthy donors, but not AML blasts.

Results

We observed a significant difference between normal and AML BM immune cells. Here, we defined the diversity of dendritic cells (DC) and macrophages in different AML patients. We also identified several unique immune cell types including T helper cell 17 (TH17)-like intermediate population, cytotoxic CD4+ T subset, T cell: erythrocyte complexes, activated regulatory T cells (Treg), and CD8+ memory-like subset. Emerging AML cells remodels the BM immune microenvironment powerfully, leads to immunosuppression by accumulating exhausted/dysfunctional immune effectors, expending immune-activated types, and promoting the formation of suppressive subsets.

Conclusion

Our results provide a comprehensive AML BM immune cell census, which can help to select pinpoint targeted drug and predict efficacy of immunotherapy.

Background

Acute myeloid leukemia (AML), as a heterogeneous disease caused by several specific mutant genes (such as FLT3, NPM1, DNMT3A, IDH1, IDH2, TET2, fusion genes, and so on), is characterized by increased proliferation of abnormal myeloid progenitors with blocking terminal differentiation in BM and other tissues [1]. Rapid clonal expansion of malignant blasts harbors within the BM microenvironment, replaces heterogeneous hematopoietic cells and stromal cells, and impairs normal hematopoiesis and immune cell development [2]. The expanding malignant cells not only impairs stromal cells and hematopoiesis, but also remodels BM immune microenvironment [3]. The malignant myeloid cells can impair osteogenesis [4], HSC-progenitor transition [5], myelo-erythropoiesis [6], erythroid differentiation [7], macrophage phagocytosis [8], dendritic cell differentiation [9], T cell anti-tumor function [10], and Natural killer (NK) cell immune surveillance [11]. However, compared with other cancer types, especially solid tumor, the immune cell types, immune status, and molecular mechanisms of AML patient BM microenvironment are poorly understood.

AML microenvironment includes complex interactions between immunosuppressive cell types, cytokines, and surface stimulatory molecules. AML cells can decrease MHC-I/II expression, produce reactive oxygen species (ROS) and indolamine-2,3-dioxygenase (IDO), and increase inhibitory ligands PDL1, B7-H3 (CD276), and Galectin 9 (Gal), which lead to escape immune surveillance and T cell exhaustion [12]. Antibody-based immune checkpoint blockade for AML patients does not seem ideal [13, 14], and underlying mechanism is not very clear. NK cells, as central players of innate immune system, owns robust anti-tumor effects. Under AML stress, several mechanisms involved impairing NK cell function [11]: 1) decreased expression of IFN-γ, TNF-α, NKp30, NKp44, and NKp46, and increased inhibitory NKG2A and KIR2DL2 in NK cells [15,16,17]; 2) increased AML cell resistance to NK cell-mediated cytotoxicity [18,19,20]; 3) suppressed by immunosuppressive cell types, such as DC and Treg [21,22,23]. Treg is also a critical player of immune response, which can limit activation and proliferation of cytotoxic lymphocytes by secreting anti-inflammatory cytokines, competing cytokines/costimulators, and contact-dependent suppression [24]. What cannot be ignored is that roles of myeloid lineage in AML microenvironment, such as myeloid-derived suppressor cells (MDSC), macrophages, and DC. MDSC is defined as innate myeloid cells with immunosuppressive function during cancer, which can be divided into several subsets with explicit developmental trajectory [25]. MDSC-like cells also expanded in AML patients [26], but its heterogeneity and features of developmental trajectory do not benefit for screening effective specific targets, which should be identified by scRNA-seq analysis and meticulous functional assay. And malignant cells polarized macrophages towards a tumor supporting status [8, 27], which have been reported in many other tumor types [28]. In addition, DC dysfunction impairs the immune response of AML patients [29, 30]. In summary, AML blasts can systematically change the immune status of BM microenvironment to support malignant cell growth and resist immune surveillance.

In this study, we characterized the immune components of different AML patient-derived BM cells by sc-RNAseq data analysis. Different from concentrating on analysis of clonal heterogeneity and hierarchies of AML cells [31, 32], we focused on the mature immune cell types and immune status under the stress of AML and chemotherapy. To assess the composition, function, and status of immune cells, we performed scRNA-seq analysis using publicly-available datasets (GSE116256) [32]. Our analysis revealed different AML patients owned unique immune profiles, diversity of immunosuppressive DC and macrophage subsets, exhausted and dysfunctional T/NK subpopulations, and suppressive T cells with unusual developmental trajectories. Finally, we investigated several new immune cell types or status, such as TH17/Treg intermediate population, cytotoxic CD4+ effectors, T cell: erythrocyte complexes, activated Tregs, and CD8+ memory-like cells. These results have deepened our understanding of immune cell components and status, which also noted us the choice of immunotherapy strategy should be customized according to the AML patient’s specific BM immune microenvironment, rather than using trial-and-error strategies. In the foreseeable future, the most optimized and efficient immunotherapy choice will be solved by high-throughput scRNA-seq technology, based on further comprehensive investigations of tumor cell: immune cell and/or immune cell: immune cell interaction.

Methods

scRNA-seq datasets

The scRNA-seq datasets of AML BM cells and healthy donor BM cells were acquired from the Gene Expression Omnibus (GEO) database (GSE116256). The scRNA-seq data was acquired from BM cells of 16 AML patients and 4 healthy donors. And the information about cell preparation and single-cell transcriptome profiling can be got from the paper of Bradley E. Bernstein and his colleagues [32]. The GSM numbers of all these samples with other detailed information (days from diagnosis, gender, age, mutations, and so on) are listed in Supplementary Table 1.

Quality control and data processing

Single-cell datasets of AML patients and healthy donors were integrated using “merge” function in version 3.2.2 of Seurat R package [33]. We filtered cells that have unique feature counts over 3000, less than 200, and ≥ 10% mitochondrial counts. The merged dataset was normalized using Seurat “NormalizedData” function with a global-scaling normalization method “LogNormalize”, and multiplied this by a scale factor (10,000 by default). And then scaled by performing Seurat “ScaleData” function with regression of the variation of “nCount_RNA” and “percent.mt”. Performing Seurat “JackStrawPlot” function and “ElbowPlot” function helped to select suitable dimensionality. Dimension reduction analysis was performed by Seurat “RunPCA” function, and non-linear dimensional reduction was performed by Seurat “RunUMAP” function.

Reconstructing cell development trajectories

To explore the developmental progression of naïve CD4+ T cells to TH17-like cells and/or Treg cells, we used Monocle package (version 2.14.0) for reconstructing their development trajectories [34]. We extracted the dataset of naïve CD4+ cluster, TH17-like cluster, and Treg cluster, and then selected the cluster feature genes for the trajectory reconstruction.

Survival analysis

The TCGA AML data (file “TCGA-LAML.htseq_fpkm.tsv”, file “TCGA-LAML.survival.tsv”, and file “gencode.v22.annotation.gene.probeMap”) were download from UCSC Xena (http://xena.ucsc.edu/) [35] and used to assess the prognostic effect of single functional genes, preference gene sets, and gene sets from cluster biomarkers. Cluster biomarkers were got through performing Seurat “FindAllMarkers” function and reporting only the positive ones. We used package “survival” and “survminer” packages to get the survival curve.

Results

A scRNA-seq census of AML BM immune cells pre- and post-treatment

We hypothesized the immune phenotypes and status were remodeled by uncontrollable AML blasts, it might be identifiable in data generated from recent efforts to distinguish AML hierarchies [32]. Bernstein and his colleagues showed an atlas of AML cell states by scRNA-seq, and found monocyte-like AML cells suppressed T cell activity by expressing immunomodulatory genes [32]. To characterize the dynamic changes of mature hematopoietic cell lineages’ states at more refined levels, we first downloaded and explored the relevant datasets from GSE116256 (Supplementary Table 1) [32]. The scRNA-seq data from BM cells of 16 AML patients and 4 healthy donors was performed uniform manifold approximation and projection (UMAP) analysis [33]. These 36,477 BM-derived cells segregated into 22 populations (Fig. 1a). These populations were identified based on the expression of canonical marker genes for mature terminal lineages (Fig. 1b) and remarkable genes for hematopoietic stem/progenitor cells (HSPCs) or leukemia stem/progenitor cells (LSPC) (dubbed “SP-like cells”) (Supplementary Figure 1A-D).

Fig. 1
figure1

Dissection and clustering of AML patient BM cells and healthy donor BM cells. a The UMAP projection of BM cells from 16 AML patients and 4 healthy donors, showing the formation of 22 main clusters, including 3 for T/NK cells (Cluster 0, Cluster 3, and Cluster 17), 3 for B lineage (Cluster 12, Cluster 14, and Cluster 16), 3 for mature myeloid lineage (Cluster 1, Cluster 11, and Cluster 20), 3 for erythroid lineage (Cluster 8, Cluster 18, and Cluster 19), 9 for SP-like cells (Cluster 2, Cluster 4, Cluster 5, Cluster 6, Cluster 7, Cluster 9, Cluster 10, Cluster 13, and Cluster 15), and 1 for non-hematopoietic stromal cells (Cluster 21). Each dot represents to one single cell, and is colored according to cell cluster. b CD3E, CD4, CD8A, CD1C, CD14, NCAM1, and CD79B-expressing (expression value > 0) of across 36,477 BM-derived cells illustrated in UMAP plots. c Histogram of cell-type fractions for each AML patient and healthy donors’ BM cells, colored based on cell type

Then we analyzed the compositions of mature hematopoietic lineages in AML samples during treatment (Fig. 1c). Compared with healthy BM samples, some AML samples (AML420B.D0, AML420B.D15, AML420B.D31; AML328.D113, AML328.D171; AML707B.D18, AML707B.D41, AML707B.D97, AML707B.D113; AML870.D14) have higher proportions of CD4+ T cells (Cluster 0), which indicated that chemotherapy increased the frequency of CD4+ T cells with wiping out malignant cells. Most untreated AML samples had a low frequency of CD8+-T/NKT/NK subset (Cluster 3). And the chemotherapy can increase the frequency of this cluster, especially in the second to third weeks after treatment (Supplementary Figure 1E). This trend also occurred within MKI67+ T cell population (Cluster 17) (Supplementary Figure 1B), LYZhiEREGhi monocyte precursor population (Cluster 7), and erythroid lineage 1 (Cluster 8) (Supplementary Figure 1E) in some patients. Interestingly, Mono/Mac subpopulation (Cluster 1) and erythroid lineage 2 (Cluster 18) showed two different treatment response models (Supplementary Figure 1E), which might be associated with the functional heterogeneity of mature cells, needing more detailed analysis. And the treatment could decrease the frequency of B cells (Cluster 16) initially, and increase subsequently. Of note, some AML patients had disproportional number of mast cells (Cluster 20), and treatment can decrease this population (Supplementary Figure 1E). Interestingly, analyses of AML samples from TCGA indicated that patients with high expression of signature genes of MSI2hiMAP3K8hi SP-like population and mast cells, low LYZhiEREGhi monocyte precursor, DC, and MKI67+ T showed significantly better overall survival (Supplementary Figure 1F). And specific-gene signature of other clusters, such as CD4+ cluster and CD8+-T/NKT/NK cluster did not yield significant survival association, which is not consistent with perception of T/NK cell anti-tumor function. And these results foreboded the immune cell function might degrade or inverse in a broader sense.

AML is a high-risk hematological malignancy and show high heterogeneity with complex mutant and/or fusion gene combinations. AML cells occupies the niche of hematopoiesis, leads to ineffective hematopoiesis, induces immune dysfunction. Periphery blood (PB) samples of AML patients showed us the immune dysfunction of lymphocytes based on the analysis of flow cytometry and immune-related risk of factors, indicated that the immune signatures corrected clinical outcomes [36]. The immune status of BM microenvironment are more sensitive to reflect immune response for clearing AML blast cells compared to PB [37]. So dissecting the immune landscape of AML is important for predicting the immune status and screening suitable immunoregulatory drugs. Our analysis based on scRNA-seq of AML BM cells can detect much more details to dissect the immune landscape, and help to find new immune cell subsets and treatment strategies.

The unique composition of DC subtypes in AML patients

To uncover the spectrum of DC heterogeneity and states, we utilized UMAP to re-cluster 1293 DC-like cells (Fig. 2a), and got 5 clusters (Table 1). Consistent with previous studies, the DC-like population derived from AML patients and healthy donors, expresses CD11c (ITGAX), CD18 (ITGB2), and MHC II molecules (HLA-DRB5, HLA-DRB1, and HLA-DRA) at high levels, but CD11b (ITGAM) at a low level (Fig. 2b), which identified the DC population [38]. Cluster 0 express CD36, FLT3, LY86, PTPRCAP, GLIPR2, IL3RA (CD123),and SEMA4D at high levels, but LYZ at a low level, and was defined as plasmacytoid dendritic cells (pDC) [39]. Compared with other 4 clusters, CD206+ DC subpopulation (cluster 1) express higher levels of MRC1 (CD206), FCER1A, IFI30, IL18, and Fc gamma receptor II (FCGR2A, FCGR2B, and FCGR2C), with an obvious CD4-expressing feature (Supplementary Figure 2A). Cluster 2 (dubbed “CLEC7A+ DCs”) is identified by high expression of CLEC7A and ITGAE (CD103) (Supplementary Figure 2A).

Fig. 2
figure2

Dissection and clustering of DC-like cells in AML patients and healthy donors. a UMAP plot of DC-like cells from Fig. 1a-represented DC cluster. These DC-like cells can be divided into 5 subsets. b Expression of DC-related genes across the transcriptionally defined DC clusters. c Dot plot of select canonical DC-related genes (cDC1: CLEC9A, ANPEP, and FBXO27; cDC2: CLEC12A, CLEC10A, SIRPA, and DENND3; pDC: IL3RA and JCHAIN) differentially expressed between different DC subsets. d Proportion of 5 DC subsets in total DC cells in each AML patient or healthy donor. Only samples of containing ≥10 DC-like cells were represented. e Dot plot of differentially expressed co-inhibitory and co-stimulation molecule genes. f Dot plot of differentially expressed cytokine genes. g The Kaplan-Meier overall survival curves of TCGA AML patients grouped by specific DC subset (CD206+ DC and CX3CR1+ DC) gene sets and DC-related genes (TGFB1, MRC1, CLEC7A, ITGAX, ITGB2, CX3CR1, CCL22, and TNFSF8). + represents censored observations, and P value was calculated by multivariate Cox regression

Table 1 Immune cell types and functions in AML BM microenvironment

CD1C+ DC (cluster 3) express AREG, S100A8, S100A9, CD86, and CD1C at high levels. And Cluster 4 was identified as CX3CR1+ DC with CD8A expression. CD1C+ DC shows the cDC1 gene expression pattern, while cDC2 pattern in Cluster 1–4, and pDC pattern in cluster 0 (Fig. 2c). cDC1s are critical for eliciting anti-tumor CD8+ T cell responses and T helper 1 (TH1) responses, cDC2 for CD4+ T cell response, and pDC for producing large amounts of type I IFN [40]. These results indicated that T cell response biased to CD4+ components, but not CD8+ components. The 5 clusters also show obvious differential expression pattern of transcription factor (Supplementary Figure 2B), pattern recognition receptors (Supplementary Figure 2C), cell adhesion/migration molecules (Supplementary Figure 2D), cytokines (Supplementary Figure 2E), and chemokine receptors (Supplementary Figure 2F).

As Fig. 2d shown, healthy donor BM-derived DC cells mainly located on the CD1C+ DC population. Compared with healthy states, AML states downregulated the expression of SELL (CD62L), CD44, and CD2, which indicates that the cell adhesion/migration function of AML patient BM-derived DC is changed (Supplementary Figure 2D). As known, DC plays critical roles in T cell response, so we illuminated the expression level of T cell function-related costimulatory and coinhibitory molecules (Fig. 2e). CD1C+ subset express many functional molecules, such as CD80, CD86, VEGFA, CXCL8, CXCL2, IL1B, IL1RN, IL10, IL15, IL23A, IL7, OSM, VEGFB, and GDF9 (Fig. 2f). And found that CLEC7A+ DC cluster showed low levels of these molecules except TNFSF9, CX3CR1+ DC cluster with high levels of costimulatory and coinhibitory molecules, and other clusters with high levels of costimulatory molecules. CLEC7A is a functional receptor in DC/Macrophage to enhance NK cells-mediated tumoricidal activity [41]. And high expression of TNFSF9 (4-1BBL) also indicated this population might enhance T cell-mediated tumoricidal activity, as previously reported [42], but significantly corrected with poor prognosis puzzlingly (Fig. 2g). CLEC7A+ DC appeared in a few AML samples (AML475.D0, AML371.D0, and AML556.D31). pDC subset express much more TNFSF4, TNFSF13, LTB, HGF, and CCL28 (Fig. 2f), which appeared in many AML samples (AML1012.D0, AML210A.D0, AML329.D0, AML419A.D0, AML556.D31, AML707B.D41, and AML921A.D0). TNFSF4 (OX40L) provides co-stimulatory signals to enhance T cell function through TRAF2 and TRAF5 [43]. DC-derived TNFSF13 (APRIL) showed antitumor potential by upregulating proliferation and survival of T cells [44]. HGF is a negative cytokine of cancer immunotherapies involved in reactive recruitment of neutrophils, which can impair T cell expansion and effector function [45], but significantly corrected with good prognosis (Fig. 2g). And CCL18 is associated with recruitment of tumor-associated macrophages (TAMs), Treg, tumor-associated dendritic cells (TADCs), and cancer-associated fibroblast (CAFs). Even so, the direct costimulation of TNFSF4 and TNFSF13 is the mainstream power to drive the antitumor potential. FCGR2A (CD32a) and FCGR2C (CD32c) are critical for phagocytosis and cross-presentation of antibody-coated antigens, but FCGR2B (CD32b) as an inhibitory FcγR [46]. CD206+ DC subset with high expression of both activate and inhibitory FcγR, might involve a specific dysfunction of immune checkpoint. We also found TNFSF12, TNFSF8 (CD30L), and IL18, expressed at high levels in CD206+ DC cells (Fig. 2g). TNFSF8 plays both positive and negative roles in T cell-mediated immune function [43], but it has more important roles in regulating Treg function under AML status based on the high CD30 expression level in Treg and dysfunctional proliferation T subset (Fig. 4d). IL18 combining different cytokines can induce different immune responses, such as type 1 response, type 2 response, and innate-type allergic inflammation [47]. IL16, IL32, CCL5, CCL23, PDGFA, PDGFC, TNFSF14, TNFSF10, XCL1 and XCL2 were highly enriched in the CX3CR1+ DC subset. TNFSF10 (TRAIL), as a proapoptotic ligand, involved the CD4+ T cell apoptosis mediating by TRAIL/DR5 signaling [48]. Proinflammatory TNFSF14 (LIGHT) expresses lower than TNFSF10. And PDGF was expressed by TADC [49], involved formation of tumor microenvironment [50] and Treg induction [51]. The proportions of this CX3CR1+ DC subset are increased in most AML samples (AML1012.D0, AML210A.D0, AML329.D0, AML329.D20, AML419A.D0, AML707B.D41, and AML707B.D112). The CX3CR1+ DC subset expressed CD274 and PDCD1LG2, as the ligands of PDCD1 mediating T-cell suppression. Taken together, immunosuppression-related DC cells increased in most AML samples, especially Treg-related CD206+ DC and T cell suppression-related CX3CR1+ DC, which infirmed by the survival analysis of TCGA AML data (Fig. 2g and Supplementary Figure 2F). And expression of several DC-related genes (TGFB1, MRC1, CLEC7A, ITGAX, ITGB2, CX3CR1, CCL22 and TNFSF8) significantly corrected with poor prognosis of TCGA AML patients (Fig. 2g).

Dysfunctional DC impairs innate and adaptive immune response. Educated DC, as mentor of T/NK cells, has immense potential in tumor immunotherapy [52]. Our results showed that AML patient-derived DC compositions changed whatever treated or not. And two immunosuppressive DC subset, CD206+ DC and CX3CR1+ DC became the mainstream, which might influence immunotherapy. These results note us to identify the major DC subset, before adopting specific DC-based tumor immunotherapy. Several clinical studies based on DC vaccines showed favorable outcomes in AML patients [53, 54]. But it just worked in some patients, which might remind us that the subsets and states of AML patient autologous DC influenced the DC-based tumor immunotherapy. We can harvest more powerful antitumor-DC by metabolic status reprogramming [55] and DC regeneration from induced pluripotent stem cell (iPSC)/other somatic cells [56, 57].

The diverse immune phenotypes of macrophages in different AML patients

In order to better understand the heterogeneity of monocytes and macrophages within and across AML patients and healthy donors, we extracted and clustered this population (4487 cells) by UMAP (Supplementary Figure 3A), and got 10 subsets (Cluster 0:MS4A6Ahigh subset, Cluster 1:CD163high subset, Cluster 2:MKI67+ subset, Cluster 3:MARCOhigh subset, Cluster 4:FCGR3A+ subset, Cluster 5:MRC1+ subset, Cluster 6:CD14highITGAMhigh subset, Cluster 7:CXCL8high subset, Cluster 8:IL10highTGFB1high subset, and Cluster 9:IFNGhigh subset) (Table 1). Cluster 9:IFNGhigh subset (119 cells) express both CD14 and CD3D (Supplementary Figure 3B), consisted of T cells bound to monocytes or monocyte debris [58], and were excluded in subsequent analysis (Fig. 3a).

Fig. 3
figure3

Dissection and clustering of mature myeloid lineages in AML patients and healthy donors. a UMAP plot of Monocyte/Macrophages from Fig. 1a-represented Mono/Mac cluster. These mature myeloid cells can be divided into 9 subsets without CD14+CD3D+ subset. b Dot plot of differentially key cell-type marker genes. c Heatmap of mean expression of selected cytokine genes in each cell subset. d Proportion of 9 mature myeloid subsets in myeloid population subsets in each AML patient or healthy donor. Only samples of containing ≥20 myeloid cells were represented. e The Kaplan-Meier overall survival curves of TCGA AML patients grouped by specific MACRO+ subset gene sets and myeloid-related genes (CD163, ITGAM, MMP9, and CCL5). + represents censored observations, and P value was calculated by multivariate Cox regression

Interestingly, MKI67+ subset is the major composition with VEGFA-expressing at high level in healthy donors BM-derived monocytes/macrophages (Fig. 3b). MS4A6Ahigh subset is a common population in most AML samples and part of normal healthy donor samples. This subset showed an M2-like pattern (CSF1R, and TGFB1) and expressed a Treg-attracting chemokine CCL22 (Fig. 3c). CD163high subset with M2 phenotype (CD163, CXCR4, IL4R, TGFB1, and CCL18) is a novel population, and as a dominant population in AML419 but not other patients. CX3CR1 and CSF1R are highly expressed in FCGR3A+ (CD16+) subset, which can be identified as an inflammatory phenotype in anti-tumor immune response [59]. Although in preclinical models of glioblastoma, targeting TAMs and microglia using CSF1R inhibitor combined with radiotherapy could enhance survival [60], might be not suitable for AML patients and impairs anti-tumor immune response [61]. This reminds us CSF1R inhibitors (such as BLZ945 and PLX3397) might not bring benefit to some AML patients (AML210A, AML329, and AML707B) without identifying major groups of CSF1R-expressing immune cells. In solid tumor models (melanoma, breast cancer, and colon carcinoma), pattern recognition scavenger receptor MARCO defined immunosuppressive TAMs. And antibody targeting this subset can enhance the effect of anti-CTLA4 checkpoint therapy to block the tumor growth and metastasis [62]. Cluster 3:MARCOhigh subset showed higher expression of several typic genes (M1-like gene: IL1B; M2-like genes: CD163, IL4R, CSF1R, CX3CR1,and TGFB1) than MKI67+ subset, which was positively correlated with M2-like phenotype. Anti-MARCO immunotherapy might especially benefit AML patient (AML210A, AML314, AML328, AML329, and AML707B) with high frequencies of MARCOhigh subset. MRC1+ subset also showed a M2-like phenotype with high expression of MRC1 (CD206), CSF1R, and TGFB1, and enriched in some AML patients, such as AML210A, AML707B, and AML921. And these patients might benefit from innate defense regulators RP-182, which can trigger a conformational switch of CD206 and enforce the TAMs to M1-like cells [63]. CXCL8high subset with high expression of HIF1A, also expressed other important cytokines or regulators, such as MMP9, TNF, IL1B, VEGFA, CCL5, TNFSF8, and only enriched in AML371 patient (Fig. 5d). For AML475, chemotherapy can decrease the frequency of a high immunosuppressive IL10highTGFB1high populations (Fig. 3d). What’s more, MARCOhigh subset-gene signature and expression of several macrophages-related genes (CD163, ITGAM and CCL5) significantly corrected with poor prognosis of TCGA AML patients, but MMP9 with good prognosis (Fig. 3e and Supplementary Figure 3C).

Many macrophage-based therapies had been developed, such as anti-SIRPα antibody [64], anti-MARCO antibody [62], chimeric antigen receptor macrophages (CAR-M) [65,66,67], and adoptively transferred macrophages with IFN-γ backpacks [68]. Mono/Mac compositions of AML patients showed substantial variation. The large proportional differences of Mono/Mac subsets indicated the immune status are diverse in different AML patients. This might note us simplex immunotherapy, by driving the directed differentiation of Mono/Mac, is limited and variable for AML treatment.

Burnt conventional T/NK lineage and expanding suppressive subsets in AML BM microenvironment

To comprehensively analyze the functional change of T/NK lineage in AML states, we utilized UMAP for identifying T subsets with more details (Fig. 4a). Thus, the T/NK lineage was divided into 10 clusters (Table 1). CD4+ naïve T cells (Cluster 0) was identified by the high expression of naïve/TCM-state-related genes (CCR7, LEF1, TCF7, and SELL) (Fig. 4b). CD69high CD4+ T cells (Cluster 1) was identified by high expressing CD69, and CD69low CD4+ T cells (Cluster 5) accordingly. Of note, Cluster 4 (dubbed “TH17-like cells”) express RORC, RORA, STAT3,BATF, AHR, IRF4, and MAF. Cluster 6 (dubbed “Treg cells”) uniquely express IL2RA and FOXP3 at high levels; Proliferating T cell cluster (Cluster 7) express canonical proliferation markers MKI67 and PCNA. IFN-CD4+ subpopulation (Cluster 8) express interferon-stimulated genes (ISGs) (IFIT3, ISG15, ISG20, MX1, and IFNAR1). CD8+-T/NK/NKT cells comprised 3 clusters distinct from CD4+ T cells and included: NK/NKT-like cells (Cluster 2) expressing cytotoxic genes (GZMA, PRF1, GNLY, NKG7, GZMB, and GZMH) and NK-related genes (TYROBP, KLRF1, FCGR3A, NCAM1 (CD56), and CD160); Cluster 3 (dubbed “GNLY+GZMH+CD8+ cells”) expressing CCL5, GNLY, and GZMH at high levels; Cluster 9 (dubbed “GZMK+GZMA+CD8+ cells”) expressing CD160, GZMK, and GZMA.

Fig. 4
figure4

Diversity of T/NK subsets revealed by scRNA-seq analysis. a UMAP plot of sc-RNAseq data (n = 10,096 cells) showed 10 distinct clusters. b Dot plot of differently key cell-type marker genes. c Histogram showed the fractions of different cell-type in T/NK populations for each AML patient and healthy donors’ BM cells, colored based on cell type. d Dot plot showed the transcript expression pattern of stimulation molecules and their receptors. e The Kaplan-Meier overall survival curves of TCGA AML patients grouped by specific Treg gene sets, dysfunctional/exhausted-gene set (LAG3, TIGIT, CTLA4, HAVCR2, TOX, PDCD1, CD274, PDCD1LG2), and several genes (CD274, PDCD1LG2, and BATF). + represents censored observations, and P value was calculated by multivariate Cox regression

To uncover the T/NK functional states in AML state before and after treatment, we assessed the proportion of different subsets to total T/NK lymphocyte in both AML patients and healthy individuals. The proportions of CD4+ naïve T cells to total T/NK lymphocytes in most AML patients are much more than healthy individuals, whatever before or after treatment (Fig. 4c). Moreover, chemotherapy could promote the terminal differentiation of naïve T helper cells except AML328 and AML722B. The results profiled that AML-states inhibited the terminal differentiation of CD4+ naïve T cells, and eliminating AML cells can reverse this state in most cases, which was reported by other researchers [32]. CD69highCD4+ T cells accounted for 60% of T/NK lymphocytes in the healthy individuals’ BM. Only parts of AML patients (such as AML707B and AML420B) showed the increasing proportions of CD69highCD4+ T cells after treatment obviously. The chemotherapy cannot drive the naïve T helpers towards to CD69highCD4+ activating T subset in most cases. And CD69lowCD4+ T cells appeared in some samples (such as AML1012.D0, AML314, AML722B, and AML371) at high frequencies, which indicated these CD4+ T cells cannot be activated for performing their immune-activated function. The proportions of TH17-like cells to total T/NK lymphocytes in most AML patients are much more than healthy individuals, whatever before or after treatment (Fig. 4c, Supplementary Figure 4A), which is positively correlated with increasing IL6 secreted from malignant cells [7]. As shown, TH17-like T cells is rare in BM cells of healthy individuals. On the previous reports, TH17 cells can be regarded as a more favorable outcome in AML patients, controversially [69, 70]. We then analyzed the proportion of Treg cells by total T/NK lymphocytes in AML patients. Parts of AML samples, but not all, showed significantly higher levels of Treg subpopulation. Interestingly, the proportions of Treg cells increased within 4 weeks after treatment and then decreased. The increased Treg cells might create an immunosuppressive niche and impair the immune activation for eliminating malignant cells. And this indicated that immunotherapy of inhibiting Treg expansion and function within 4 weeks is the most effective and economical treatment strategy to enhance the treatment effect of conventional chemotherapy. Interestingly, some AML samples (AML556, AML420B, AML475, and AML707B) showed higher proportions of proliferation T cell subsets. As previously reported, dysfunctional T cells are the major intratumoral proliferating T cell compartment with dysfunctional signature [71]. Although the proliferation T cluster showed the proliferation-related markers (MKI67 and PCNA), this cluster acquired a dysfunctional signature of expressing LAG3, TIGIT, CTLA4, HAVCR2, and TOX at high levels, similar to the immunosuppressive signature of Treg. We observed a strong ISGs pattern in cluster 8 in AML patients (AML328, AML329, and AML707B). The IFN response module can mark the activated CD4+ T cells, but acquired a dysfunctional signature of expressing LAG3, PRDM1, and TIGIT.

AML samples showed higher proportions of GNLY+GZMH+CD8+ T cells and lower proportions of GZMK+GZMA+CD8+ T cells compared to healthy donors. Interestingly, NKT/NK-like cluster and GZMK+GZMA+CD8+ T cluster showed significant exhaustion-related genes (LAG3, TIGIT, HAVCR2, and TOX in NKT/NK-like cluster; PDCD1, LAG3, TIGIT, CTLA4, and TOX in GZMK+GZMA+CD8+ T cluster) (Fig. 4d) identified in many kinds of tumors [72, 73].

In conclusion, under the AML states, T/NK cells showed a diversity of subpopulations and function, and transited into immunosuppressive states. Our analysis of TCGA AML data showed the Treg and dysfunctional/exhausted T/NK subsets might represent poor prognosis (Fig. 4e and Supplementary Figure 4B). Immunotherapies have changed the methods of cancer treatment subversively, such as the application of checkpoint inhibitors targeting PD-1 and CTLA-4 [74, 75]. But these strategies are not always efficacious, which is partly due to the difference of heterogeneity and/or ability of tumor-infiltrating T lymphocytes. Our results identified distinct states of T/NK cell lineages, which confirmed that different response of checkpoint inhibitors in different AML patients with different outcomes. Importantly, increasing frequencies of TH17-like cluster and Treg cluster, as common cell signature, are obvious in AML patients. And the dynamic changes of T/NK subsets also might be associated with the outcomes, and help to adjust the immunotherapies.

Novel T clusters and functional states in AML patients

To characterize the T/NK cellular diversity in AML BM microenvironment, we further analyzed specific clusters for detecting more function and state changes. We acquired the expression data of CD69highCD4+ T population and CD69lowCD4+ T population for UMAP (Fig. 5a), and found a unique cytotoxic CD4+ T subset existed in both AML patients and healthy donors, and CD69low LTBhighCD4+ T subset in AML patients but rarely in healthy donors (Fig. 5b and Fig. 5c) (Table 1). Cytotoxic CD4+ effectors, with anti-tumor activity, was identified by other researchers in other tumors [76,77,78], and might lead to therapeutic benefit. And we first identified a cytotoxic CD4+ T subset under AML stress by scRNA-seq analysis. Although CD69lowLTBhighCD4+ T subset is approximated T cell: erythrocyte complexes (Supplementary Figure 5A), which might involve cell:cell communication events, but not technical artefact or random association [58].

Fig. 5
figure5

Unique CD4+ subsets revealed by scRNA-seq analysis. a UMAP plot of CD69highCD4+ population and CD69lowCD4+ population. These 2515 cells can be divided into 4 subsets. b Heatmap showing average expression level of cell-type genes expressed by the 4 clusters. c Histogram showed the fractions of different cell-type in CD69highCD4+ population and CD69lowCD4+ population for each AML patient and healthy donors’ BM cells, colored based on cell type. d The branched trajectory of state transition of naïve CD4+ T cells, TH17-like cells, and Treg cells in a two-dimensional state-space inferred by Monocle (version 2.14.0). Each dot corresponded to one single cell, colored according to its cluster label. e Expression maps showing log-normalized expression of typical markers (IL7R, IL2RA, MKI67, KLRB1, FOXP3, and CTLA4) in the differentiation of Naïve CD4+ T to TH17-like cells and/or Treg cells. Data are shown as log-normalized expression. Yellow indicates high expression, dark blue indicates low expression. f Violin plot showing the expression levels of functional genes (IGF1R, RORC, KLRB1, AHNAK, and IL10RA)

As known, TH17 lineage and Treg lineage strongly linked, and differentiation imbalance leads to abnormal immune states [79]. So we then analyzed the developmental trace and function of naïve CD4+ T cells, TH17-like cells, and Treg cells (Fig. 5d). The expression of signature genes and known functional markers suggested the direction of naïve CD4+ T population to TH17-like population, nearly TH17/Treg intermediate population, and then terminal Treg population (Supplementary Figure 5B). These 3 populations have different functional gene expression pattern, such as high expression level of IGF1R in Naïve CD4+ T cell, RORC and KLRB1 in TH17/Treg intermediate population, and IL10RA in Treg population (Fig. 5e, Supplementary Figure 5C). Signaling of IGF-IGF1R drive the favorite of naïve T to TH17 but not Treg, and prime the TH17 cell fate [80]. KLRB1 is a remarkable proinflammatory marker of TH17 cells and proinflammatory FoxP3+ cells [81, 82]. RORC, as a faithfully TH17-specific transcript factor, balances the generation of TH17 and Treg subsets, and expresses at TH0-like intermediate population toward to TH17 subset and Treg subset [83]. Triggered signaling of IL10RA is important for Treg-mediated suppression of TH17 proinflammatory response [84]. To illustrate the state transition of these 3 populations, we performed branched trajectory analysis by Monocle (Fig. 5d). As shown in Fig. 5d, TH17/Treg intermediate population is halfway between naïve CD4+ T cells and Treg cells. As expected, Treg-specific genes (IL2RA, FOXP3, and CTLA4) express at high levels in left terminal branches, which enriched Treg cells. And KLRB1 express at high level on the left side of the above branch. IL7R can be used as a marker to exclude non-Treg cells, which located on the non-Treg positions. Interestingly, Treg population was divided into MKI67+ activated Treg subset and MKI67 resting Treg subset. In the healthy donor, resting Treg cells and naïve CD4+ T cells are the major population in state 1 and 6. But in the AML state, Naïve CD4+, TH17/Treg intermediate cells, and MKI67+ Treg cells are increasing before and after treatment (Supplementary Figure 5D). This reminds us that arresting the directed differentiation of TH17/Treg intermediate cells to Treg cells, or transdifferentiation of TH17 to Treg, might be beneficial for therapeutic effect combing conventional chemotherapy.

To feature the changes of CD8+-T/NKT/NK cells, we analyzed three related populations (NK/NKT-like cluster, GNLY+GZMH+CD8+ cluster, and GZMK+GZMA+CD8+ cluster) by UMAP method with more details (Fig. 6a). These clusters can be divided into 6 small subsets (Cluster 0:GZMAlowGNLY+CD8+ subset, Cluster 1:TIGIT+CD8+ T subset, Cluster 2:Naïve CD8+ T subset, Cluster 3:NK subset, Cluster 4:GZMAlowGNLYlowCD8+ T subset, and Cluster 5:GZMA+GNLYlowCD8+ T subset) (Fig. 6b) (Table 1). Interestingly, the proportions of GZMA+GNLYlowCD8+ T subset to total CD8+-T/NKT/NK in most AML patients, were decreased. And the GZMA+GNLYlowCD8+ T subset showed exhausting gene expression pattern (TIGIT, PDCD1, and CTLA4). Meanwhile, TIGIT+CD8+ T population was increased in most AML patients obviously, and GZMAlowGNLY+CD8+ T population in part of AML patients (AML556) (Fig. 6c). TIGIT is an inhibitory receptor expressed on dysfunctional T cells, as previously reported [85]. And blockade of TIGIT targets CD8+ CTL or NK cells, prevents exhaustion and promotes target [85, 86]. TIGIT-blocking strategies might enhance the AML treatment mediated by CD8+ CTL and NK cells in most cases (except AML314.D0, AML556.D15, and AML870.D14). Likewise, blockade of PDCD1 and CTLA4 also enhance the AML treatment. GZMAlowGNLYlowCD8+ T subset expressed low levels of granzyme genes (GZMA, GZMB, GZMH, and GZMK), GNLY, KLRG1, ITGAE, B3GAT1, and PRF1, but high levels of TCF7, RUNX3, CD69, and IL7R (Fig. 6b and d), which is consisted with phenotype of CD8+ memory T cells [87, 88]. This subset plays positive roles in improving cancer immunotherapy [89]. This CD8+ memory-like cells have a high proportion in part of samples (AML314.D0, AML556.D0, AML556.D15, AML556.D31, AML707B.D41, and AML870.D14). In most AML patients, the proportion of transitional GZMA+GNLYlowGZMK+CD8+ effector T subset (except AML556.D0, AML556.D31,) and GZMAlowGNLY+ cytotoxic T subset (AML314.D0, AML556.D0, AML556.D15, AML556.D31, and AML870.D14), are decreased. And naïve-like CD8+ T subset with inhibitory molecules pattern (LAG3, TIGIT, and CTLA4), seems to be repulsive to CD8+ memory-like subset.

Fig. 6
figure6

Dysfunctional/exhausted CD8+-T/NK subsets revealed by scRNA-seq analysis. a UMAP plot of NK/NKT-like population, GNLY+GZMH+CD8+ T population, and GZMK+GZMA+CD8+ T population. These 2901 cells can be divided into 6 subsets. b Dot plot showed the transcript expression pattern of cell-type genes. c Histogram showed the fractions of different cell-type in NK/NKT-like population, GNLY+GZMH+CD8+ T population, and GZMK+GZMA+CD8+ T population for each AML patient and healthy donors’ BM cells, colored based on cell type. d Dot plot showed the transcript expression pattern of memory-like CD8+ T-related genes

Detailed classification of tumor-infiltrating lymphocytes (TILs) and identification of cellular crosstalk in AML BM niche will provide important means of assistance based on existed mature chemotherapies [90]. By examining immune cell types can help to formulate personalized high-efficiency immunotherapy trials. Our results showed AML-derived TILs exhibited heterogenous, and combined immunotherapy strategies based on heterogenous cell-types (such as Treg, TH17-like cells, exhausted NK cell, dysfunctional CTLs, and cytotoxic CD4+ T cells) provide assistance for inhibiting relapse after chemotherapy.

Discussion

Macrophages is a plastic heterogeneous population, which involved the survival and drug resistance of AML. Neonatal macrophages can mature into classically activated (M1) and alternatively activated (M2) macrophages depending on the microenvironment. Meanwhile, M1 macrophages and M2 macrophages can translate into each other with changes of immune status. Repolarization of M2 macrophages toward M1 macrophages is helpful to establish new therapeutic targets in AML [91]. Our results showed the diversity of macrophages in AML BM niche. As reported in Renca tumor, single anti-CSF1R treatment decreased part of TAMs, which suggested differential sensitivity of distinct macrophage subsets for specifical antibody-drug [92]. But anti-CSF1R treatment for AML should consider the major CSF1R-expressing population in BM microenvironment, and avoid to eliminate anti-tumor population and impair T cell response [61]. Yu and his colleagues also showed us anti-CD40 agonist therapy could amplify immune-activated cDC1 subset, increase effector memory CD8+ CTLs, and induce the activation and expansion of TH1-like CD4+ T cells in MC38 tumor model [92]. Targeting mature myeloid lineages is an attractive therapeutic approach for AML, but firstly need to identify the mainstream of myeloid subset in BM microenvironment. The scRNA-seq analysis can help to identify the tumor microenvironment information in detail, and implement more effective immunotherapy approaches combined with conventional chemotherapy.

Many animal models showed us the evidence of antileukemic T-cell immunity with exhausted program or antigen-specific T cell tolerance [93, 94]. Our results identified the abnormal T-cell subset alterations in AML before and after treatment by scRNA-seq analysis. Augmented Treg is an obvious feature of newly diagnosed AML patients, and the proportion of Treg cells increased within 4 weeks after treatment, which are at odds with other research based on other research approaches [95,96,97,98]. Tumor microenvironment components secrete immunosuppressive chemokines, such as CCL17 and CCL22, to attract Treg cells. And dysfunctional DC populations expand Treg population mediated by cytokines and costimulatory [99]. As previously reported, activated Treg, with antigen stimulation and oligoclonal skewing, presented during early lymphocyte recovery, which is consistent with our results of scRNA-seq analysis. Furthermore, our data also indicated that there exist two Treg subsets with different states, function, or development path, in AML BM microenvironment. Meanwhile, Treg cells not only impairs conventional T-cell function, but also induces DC apoptosis and dysfunction [99]. We also found several dysfunctional DC subsets involved in the AML process, especially CX3CR1+ DC subset with high expression level of CD274 and PDCD1LG2. These dysfunctional DC subsets, conspire with Treg cells, contributed to the T cell exhaustion and anergy. And increased TH17-like cells, as an abnormal population in BM microenvironment, closed to a TH17/Treg intermediate state, might illustrate the contradictory of TH17 phenotype and function in AML patients [100, 101]. A tumor-bearing mouse model showed tumor-infiltrating Treg cells can be converted from IL17A+FoxP3neg cells fostered by TGFβ and PGE2 [102], which is consistent with our results about two Treg states and development pathways of naïve CD4+ T cell to TH17-like subsets and then Treg subsets in AML patients.

Conclusions

In summary, we analyzed scRNA-seq dataset of AML patient-derived BM cells and characterized their immune cell landscape with more details (Table 1). What’s more, these results can help us to predict the prognosis of AML (Table 2). We also identified several infrequently reported immune cell types in AML patients, such as TH17/Treg intermediate population, CD8+ memory T cells, different types of macrophages, and dysfunctional DC subsets. Exhausted conventional T cells and immunosuppressive T cells (Treg and other T subsets) can be used as targets of anti-CTLA4, anti-PD1, and anti-CD25 therapies. But these don’t always work because of the diversity of T/NK cells and other immune cell types. So the explosive emergence of immune-regulated drugs emphasizes the strong need for identification of predictive biomarkers, which help to illustrate which cell populations are the most critical targets in AML. Additionally, inhibition of TH17/Treg intermediate cells toward to Treg direction should be noted, which can be developed as a new immunotherapy strategy. What’s more, we found that mature myeloid lineages exist extremely high diversity of monocyte/macrophages and DC. Targeting single or small macrophage subset and/or DC subset, don’t seem to work for most AML patients. But as our scRNA-seq analysis of AML BM cells, this approach might be as a mean of diagnosis to help identification of effective immunotherapy strategies by targeting macrophages, DC, and TILs [103]. With the maturation and popularization of scRNA-seq technology, this technology will provide more and more details about occurrence of disease, selection and development of treatment and prognosis approaches, and predicting disease risk, as it is used in SARS-CoV-2 prevention and treatment [104,105,106].

Table 2 Summarization of major immune cell subset in AML BM microenvironment and prognosis

Availability of data and materials

All data generated or analyzed in this study are included in this article. Other data that are relevant to this article are available from the corresponding author upon reasonable request.

Abbreviations

AML:

Acute myeloid leukemia

scRNA-seq:

single-cell RNA sequencing

BM:

Bone marrow

DC:

Dendritic Cells

Treg:

Regulatory T

TH17:

T helper cell 17

ROS:

Reactive oxygen species

IDO:

Indolamine-2,3-dioxygenase

NK:

Natural killer

MDSCs:

Myeloid-derived suppressor cells

UMAP:

Uniform manifold approximation and projection

HSPCs:

Hematopoietic stem/progenitor cells

LSPC:

Leukemia stem/progenitor cells

PB:

Periphery blood

pDC:

plasmacytoid dendritic cells

TAMs:

Tumor-associated macrophages

TADCs:

Tumor-associated dendritic cells

CAFs:

Cancer-associated fibroblast

iPSC:

induced pluripotent stem cell

CAR-M:

Chimeric antigen receptor macrophages

References

  1. 1.

    Gebru MT, Wang H-G. Therapeutic targeting of FLT3 and associated drug resistance in acute myeloid leukemia. J Hematol Oncol. 2020;13(1):155.

    PubMed  PubMed Central  Article  Google Scholar 

  2. 2.

    Shafat MS, Gnaneswaran B, Bowles KM, Rushworth SA. The bone marrow microenvironment - home of the leukemic blasts. Blood Rev. 2017;31(5):277–86.

    PubMed  Article  Google Scholar 

  3. 3.

    Kokkaliaris KD, Scadden DT. Cell interactions in the bone marrow microenvironment affecting myeloid malignancies. Blood Adv. 2020;4(15):3795–803.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  4. 4.

    Chen Y, Hoffmeister LM, Zaun Y, Arnold L, Schmid KW, Giebel B, Klein-Hitpass L, Hanenberg H, Squire A, Reinhardt HC, Dührsen U, Bertram S, Hanoun M. Acute myeloid leukemia–induced remodeling of the human bone marrow niche predicts clinical outcome. Blood Adv. 2020;4(20):5257–68.

    PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Miraki-Moud F, Anjos-Afonso F, Hodby KA, Griessinger E, Rosignoli G, Lillington D, Jia L, Davies JK, Cavenagh J, Smith M, Oakervee H, Agrawal S, Gribben JG, Bonnet D, Taussig DC. Acute myeloid leukemia does not deplete normal hematopoietic stem cells but induces cytopenias by impeding their differentiation. Proc Natl Acad Sci U S A. 2013;110(33):13576–81.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    Boyd AL, Reid JC, Salci KR, Aslostovar L, Benoit YD, Shapovalova Z, Nakanishi M, Porras DP, Almakadi M, Campbell CJV, Jackson MF, Ross CA, Foley R, Leber B, Allan DS, Sabloff M, Xenocostas A, Collins TJ, Bhatia M. Acute myeloid leukaemia disrupts endogenous myelo-erythropoiesis by compromising the adipocyte bone marrow niche. Nat Cell Biol. 2017;19(11):1336–47.

    CAS  PubMed  Article  Google Scholar 

  7. 7.

    Zhang TY, Dutta R, Benard B, Zhao F, Yin R, Majeti R. IL-6 blockade reverses bone marrow failure induced by human acute myeloid leukemia. Sci Transl Med. 2020;12:538.

    Google Scholar 

  8. 8.

    Yahya SA-M, Lacramioara B, Bertram O, Judith MH, Robert FL, Aniththa T, Judith S, Renata K, Klaus L, Thomas S, Rainer H, Ulrich D, Cyrus K. Acute myeloid leukemia cells polarize macrophages towards a leukemia supporting state in a growth factor independence 1 dependent manner. Haematologica. 2016;101(10):1216–27.

    Article  CAS  Google Scholar 

  9. 9.

    Valencia J, Fernandez-Sevilla LM, Fraile-Ramos A, Sacedón R, Jiménez E, Vicente A, Varas A. Acute Lymphoblastic Leukaemia Cells Impair Dendritic Cell and Macrophage Differentiation: Role of BMP4. Cells. 2019;8(7):722.

    CAS  PubMed Central  Article  PubMed  Google Scholar 

  10. 10.

    Lamble AJ, Kosaka Y, Laderas T, Maffit A, Kaempf A, Brady LK, Wang W, Long N, Saultz JN, Mori M, Soong D, LeFave CV, Huang F, Adams H, Loriaux MM, Tognon CE, Lo P, Tyner JW, Fan G, McWeeney SK, Druker BJ, Lind EF. Reversible suppression of T cell function in the bone marrow microenvironment of acute myeloid leukemia. Proc Natl Acad Sci. 2020;117(25):14331–41.

    CAS  PubMed  Article  Google Scholar 

  11. 11.

    Lion E, Willemen Y, Berneman ZN, Van Tendeloo VF, Smits EL. Natural killer cell immune escape in acute myeloid leukemia. Leukemia. 2012;26(9):2019–26.

    CAS  PubMed  Article  Google Scholar 

  12. 12.

    Epperly R, Gottschalk S, Velasquez MP. A bump in the road: how the hostile aml microenvironment affects car t cell therapy. Front Oncol. 2020;10:262.

    PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Daver N, Basu S, Garcia-Manero G, Cortes JE, Ravandi F, Jabbour EJ, Hendrickson S, Pierce S, Ning J, Konopleva M, Andreeff M, Kornblau SM, Pemmaraju N, Bueso-Ramos CE, Blando J, Hidalgo Lopez JE, Allison J, Kantarjian HM, Sharma P. Phase IB/II study of Nivolumab in combination with Azacytidine (AZA) in patients (pts) with relapsed acute myeloid leukemia (AML). Blood. 2016;128(22):763.

    Article  Google Scholar 

  14. 14.

    Davids MS, Kim HT, Bachireddy P, Costello C, Liguori R, Savell A, Lukez AP, Avigan D, Chen YB, McSweeney P, LeBoeuf NR, Rooney MS, Bowden M, Zhou CW, Granter SR, Hornick JL, Rodig SJ, Hirakawa M, Severgnini M, Hodi FS, Wu CJ, Ho VT, Cutler C, Koreth J, Alyea EP, Antin JH, Armand P, Streicher H, Ball ED, Ritz J, et al. Ipilimumab for patients with relapse after allogeneic transplantation. N Engl J Med. 2016;375(2):143–53.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Stringaris K, Sekine T, Khoder A, Alsuliman A, Razzaghi B, Sargeant R, Pavlu J, Brisley G, de Lavallade H, Sarvaria A, Marin D, Mielke S, Apperley JF, Shpall EJ, Barrett AJ, Rezvani K. Leukemia-induced phenotypic and functional defects in natural killer cells predict failure to achieve remission in acute myeloid leukemia. Haematologica. 2014;99(5):836–47.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  16. 16.

    Costello RGT, Sivori S, Marcenaro E, Lafage-Pochitaloff M, Mozziconacci M-J, Reviron D, Gastaut JA, Pende D, Olive D, Moretta A. Defective expression and function of natural killer cell–triggering receptors in patients with acute myeloid leukemia. Blood. 2002;99(10):3661–7.

    CAS  PubMed  Article  Google Scholar 

  17. 17.

    Fauriat C, Just-Landi S, Mallet F, Arnoulet C, Sainty D, Olive D, Costello RT. Deficient expression of NCR in NK cells from acute myeloid leukemia: evolution during leukemia treatment and impact of leukemia cells in NCRdull phenotype induction. Blood. 2007;109(1):323–30.

    CAS  PubMed  Article  Google Scholar 

  18. 18.

    Nanbakhsh A, Pochon C, Mallavialle A, Amsellem S, Bourhis JH, Chouaib S. C-Myc regulates expression of NKG2D ligands ULBP1/2/3 in AML and modulates their susceptibility to NK-mediated lysis. Blood. 2014;123(23):3585–95.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. 19.

    Kearney CJ, Ramsbottom KM, Voskoboinik I, Darcy PK, Oliaro J. Loss of DNAM-1 ligand expression by acute myeloid leukemia cells renders them resistant to NK cell killing. Oncoimmunology. 2016;5(8):e1196308.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  20. 20.

    Paczulla AM, Rothfelder K, Raffel S, Konantz M, Steinbacher J, Wang H, Tandler C, Mbarga M, Schaefer T, Falcone M, Nievergall E, Dörfel D, Hanns P, Passweg JR, Lutz C, Schwaller J, Zeiser R, Blazar BR, Caligiuri MA, Dirnhofer S, Lundberg P, Kanz L, Quintanilla-Martinez L, Steinle A, Trumpp A, Salih HR, Lengerke C. Absence of NKG2D ligands defines leukaemia stem cells and mediates their immune evasion. Nature. 2019;572(7768):254–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Fauriat C, Moretta A, Olive D, Costello RG. Defective killing of dendritic cells by autologous natural killer cells from acute myeloid leukemia patients. Blood. 2005;106(6):2186–8.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    Ebata K, Shimizu Y, Nakayama Y, Minemura M, Murakami J, Kato T, Yasumura S, Takahara T, Sugiyama T, Saito S. Immature NK cells suppress dendritic cell functions during the development of leukemia in a mouse model. J Immunol. 2006;176(7):4113–24.

    CAS  PubMed  Article  Google Scholar 

  23. 23.

    Ustun C, Miller JS, Munn DH, Weisdorf DJ, Blazar BR. Regulatory T cells in acute myelogenous leukemia: is it time for immunomodulation? Blood. 2011;118(19):5084–95.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  24. 24.

    Sojka DK, Huang Y-H, Fowell DJ. Mechanisms of regulatory T-cell suppression - a diverse arsenal for a moving target. Immunology. 2008;124(1):13–22.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  25. 25.

    Alshetaiwi H, Pervolarakis N, McIntyre LL, Ma D, Nguyen Q, Rath JA, Nee K, Hernandez G, Evans K, Torosian L, Silva A, Walsh C, Kessenbrock K. Defining the emergence of myeloid-derived suppressor cells in breast cancer using single-cell transcriptomics. Sci Immunol. 2020;5:44.

    Article  CAS  Google Scholar 

  26. 26.

    Pyzer AR, Stroopinsky D, Rajabi H, Washington A, Tagde A, Coll M, Fung J, Bryant MP, Cole L, Palmer K, Somaiya P, Karp Leaf R, Nahas M, Apel A, Jain S, McMasters M, Mendez L, Levine J, Joyce R, Arnason J, Pandolfi PP, Kufe D, Rosenblatt J, Avigan D. MUC1-mediated induction of myeloid-derived suppressor cells in patients with acute myeloid leukemia. Blood. 2017;129(13):1791–801.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  27. 27.

    Wang L, Zheng G. Macrophages in leukemia microenvironment. Blood Science. 2019;1(1):29–33.

    Article  Google Scholar 

  28. 28.

    Pathria P, Louis TL, Varner JA. Targeting tumor-associated macrophages in Cancer. Trends Immunol. 2019;40(4):310–27.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Rickmann M, Macke L, Sundarasetty BS, Stamer K, Figueiredo C, Blasczyk R, Heuser M, Krauter J, Ganser A, Stripecke R. Monitoring dendritic cell and cytokine biomarkers during remission prior to relapse in patients with FLT3-ITD acute myeloid leukemia. Ann Hematol. 2013;92(8):1079–90.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. 30.

    Lau CM, Nish SA, Yogev N, Waisman A, Reiner SL, Reizis B. Leukemia-associated activating mutation of Flt3 expands dendritic cells and alters T cell responses. J Exp Med. 2016;213(3):415–31.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Wu J, Xiao Y, Sun J, Sun H, Chen H, Zhu Y, Fu H, Yu C, Weigao E, Lai S, Ma L, Li J, Fei L, Jiang M, Wang J, Ye F, Wang R, Zhou Z, Zhang G, Zhang T, Ding Q, Wang Z, Hao S, Liu L, Zheng W, He J, Huang W, Wang Y, Xie J, Li T, et al. A single-cell survey of cellular hierarchy in acute myeloid leukemia. J Hematol Oncol. 2020;13(1):128.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  32. 32.

    van Galen P, Hovestadt V, Wadsworth Ii MH, Hughes TK, Griffin GK, Battaglia S, Verga JA, Stephansky J, Pastika TJ, Lombardi Story J, Pinkus GS, Pozdnyakova O, Galinsky I, Stone RM, Graubert TA, Shalek AK, Aster JC, Lane AA, Bernstein BE. Single-cell RNA-Seq reveals AML hierarchies relevant to disease progression and immunity. Cell. 2019;176(6):1265–1281 e1224.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  33. 33.

    Stuart T, Butler A, Hoffman P, Hafemeister C, Papalexi E, Mauck WM III, Hao Y, Stoeckius M, Smibert P, Satija R. Comprehensive Integration of Single-Cell Data. Cell. 2019;177(7):1888–1902.e1821.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Qiu X, Hill A, Packer J, Lin D, Ma Y-A, Trapnell C. Single-cell mRNA quantification and differential analysis with census. Nat Methods. 2017;14(3):309–15.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Goldman MJ, Craft B, Hastie M, Repečka K, McDade F, Kamath A, Banerjee A, Luo Y, Rogers D, Brooks AN, Zhu J, Haussler D. Visualizing and interpreting cancer genomics data via the Xena platform. Nat Biotechnol. 2020;38(6):675–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Tang L, Wu J, Li C-G, Jiang H-W, Xu M, Du M, Yin Z, Mei H, Hu Y. Characterization of immune dysfunction and identification of prognostic immune-related risk factors in acute myeloid leukemia. Clin Cancer Res. 2020;26(7):1763–72.

    CAS  PubMed  Article  Google Scholar 

  37. 37.

    Tan J, Yu Z, Huang J, Chen Y, Huang S, Yao D, Xu L, Lu Y, Chen S, Li Y. Increased PD-1+Tim-3+ exhausted T cells in bone marrow may influence the clinical outcome of patients with AML. Biomark Res. 2020;8(1):6.

    PubMed  PubMed Central  Article  Google Scholar 

  38. 38.

    Collin M, McGovern N, Haniffa M. Human dendritic cell subsets. Immunology. 2013;140(1):22–30.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  39. 39.

    Villani AC, Satija R, Reynolds G, Sarkizova S, Shekhar K, Fletcher J, Griesbeck M, Butler A, Zheng S, Lazo S, Jardine L, Dixon D, Stephenson E, Nilsson E, Grundberg I, McDonald D, Filby A, Li W, De Jager PL, Rozenblatt-Rosen O, Lane AA, Haniffa M, Regev A, Hacohen N. Single-cell RNA-seq reveals new types of human blood dendritic cells, monocytes, and progenitors. Science. 2017;356:6335.

    Article  CAS  Google Scholar 

  40. 40.

    Collin M, Bigley V. Human dendritic cell subsets: an update. Immunology. 2018;154(1):3–20.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  41. 41.

    Chiba S, Ikushima H, Ueki H, Yanai H, Kimura Y, Hangai S, Nishio J, Negishi H, Tamura T, Saijo S, Iwakura Y, Taniguchi T. Recognition of tumor cells by Dectin-1 orchestrates innate immune cells for anti-tumor responses. Elife. 2014;3:e04177.

    PubMed  PubMed Central  Article  Google Scholar 

  42. 42.

    Huang R, Li X, He Y, Zhu W, Gao L, Liu Y, Gao L, Wen Q, Zhong JF, Zhang C, Zhang X. Recent advances in CAR-T cell engineering. J Hematol Oncol. 2020;13(1):86.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  43. 43.

    Ward-Kavanagh LK, Lin WW, Sedy JR, Ware CF. The TNF receptor superfamily in co-stimulating and co-inhibitory responses. Immunity. 2016;44(5):1005–19.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    Shurin MR, Ma Y, Keskinov AA, Zhao R, Lokshin A, Agassandian M, Shurin GV. BAFF and APRIL from Activin A–treated dendritic cells Upregulate the antitumor efficacy of dendritic cells <em>in vivo</em>. Cancer Res. 2016;76(17):4959–69.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  45. 45.

    Glodde N, Bald T, van den Boorn-Konijnenberg D, Nakamura K, O’Donnell JS, Szczepanski S, Brandes M, Eickhoff S, Das I, Shridhar N, Hinze D, Rogava M, van der Sluis TC, Ruotsalainen JJ, Gaffal E, Landsberg J, Ludwig KU, Wilhelm C, Riek-Burchardt M, Müller AJ, Gebhardt C, Scolyer RA, Long GV, Janzen V, MWL T, Kastenmüller W, Mazzone M, Smyth MJ, Tüting T, Hölzel M. Reactive Neutrophil Responses Dependent on the Receptor Tyrosine Kinase c-MET Limit Cancer Immunotherapy. Immunity. 2017;47(4):789–802.e789.

    CAS  PubMed  Article  Google Scholar 

  46. 46.

    Boruchov AM, Heller G, Veri MC, Bonvini E, Ravetch JV, Young JW. Activating and inhibitory IgG fc receptors on human DCs mediate opposing functions. J Clin Invest. 2005;115(10):2914–23.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  47. 47.

    Nakanishi K. Unique action of Interleukin-18 on T Cells and other immune cells. Front Immunol. 2018;9:763.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  48. 48.

    Cheng L, Yu H, Wrobel JA, Li G, Liu P, Hu Z, Xu XN, Su L. Identification of pathogenic TRAIL-expressing innate immune cells during HIV-1 infection in humanized mice by scRNA-Seq. JCI Insight. 2020;5(11):e135344.

    PubMed Central  PubMed  Google Scholar 

  49. 49.

    Triplett TA, Cardenas KT, Lancaster JN, Hu Z, Selden HJ, Jasso GJ, Balasubramanyam S, Chan K, Li L, Chen X, Marcogliese AN, Davé UP, Love PE, Ehrlich LI. Endogenous dendritic cells from the tumor microenvironment support T-ALL growth via IGF1R activation. Proc Natl Acad Sci U S A. 2016;113(8):E1016–25.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  50. 50.

    López de Andrés J, Griñán-Lisón C, Jiménez G, Marchal JA. Cancer stem cell secretome in the tumor microenvironment: a key point for an effective personalized cancer treatment. J Hematol Oncol. 2020;13(1):136.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  51. 51.

    Agrawal S, Ganguly S, Hajian P, Cao J-N, Agrawal A. PDGF upregulates CLEC-2 to induce T regulatory cells. Oncotarget. 2015;6(30):28621–32.

    PubMed  PubMed Central  Article  Google Scholar 

  52. 52.

    Wang Y, Xiang Y, Xin VW, Wang X-W, Peng X-C, Liu X-Q, Wang D, Li N, Cheng J-T, Lyv Y-N, Cui S-Z, Ma Z, Zhang Q, Xin H-W. Dendritic cell biology and its role in tumor immunotherapy. J Hematol Oncol. 2020;13(1):107.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  53. 53.

    Khoury HJ, Collins RH Jr, Blum W, Stiff PS, Elias L, Lebkowski JS, Reddy A, Nishimoto KP, Sen D, Wirth ED 3rd, Case CC, DiPersio JF. Immune responses and long-term disease recurrence status after telomerase-based dendritic cell immunotherapy in patients with acute myeloid leukemia. Cancer. 2017;123(16):3061–72.

    CAS  PubMed  Article  Google Scholar 

  54. 54.

    Anguille S, Van de Velde AL, Smits EL, Van Tendeloo VF, Juliusson G, Cools N, Nijs G, Stein B, Lion E, Van Driessche A, Vandenbosch I, Verlinden A, Gadisseur AP, Schroyens WA, Muylle L, Vermeulen K, Maes MB, Deiteren K, Malfait R, Gostick E, Lammens M, Couttenye MM, Jorens P, Goossens H, Price DA, Ladell K, Oka Y, Fujiki F, Oji Y, Sugiyama H, et al. Dendritic cell vaccination as postremission treatment to prevent or delay relapse in acute myeloid leukemia. Blood. 2017;130(15):1713–21.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  55. 55.

    Du X, Wen J, Wang Y, Karmaus PWF, Khatamian A, Tan H, Li Y, Guy C, Nguyen T-LM, Dhungana Y, Neale G, Peng J, Yu J, Chi H. Hippo/Mst signalling couples metabolic state and immune function of CD8α+ dendritic cells. Nature. 2018;558(7708):141–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  56. 56.

    Rosa FF, Pires CF, Kurochkin I, Ferreira AG, Gomes AM, Palma LG, Shaiv K, Solanas L, Azenha C, Papatsenko D, Schulz O, Reis Sousa C, Pereira C-F. Direct reprogramming of fibroblasts into antigen-presenting dendritic cells. Science Immunology. 2018;3(30):eaau4292.

    PubMed  Article  Google Scholar 

  57. 57.

    Silk KM, Silk JD, Ichiryu N, Davies TJ, Nolan KF, Leishman AJ, Carpenter L, Watt SM, Cerundolo V, Fairchild PJ. Cross-presentation of tumour antigens by human induced pluripotent stem cell-derived CD141(+)XCR1+ dendritic cells. Gene Ther. 2012;19(10):1035–40.

    CAS  PubMed  Article  Google Scholar 

  58. 58.

    Burel JG, Pomaznoy M, Lindestam Arlehamn CS, Weiskopf D, da Silva AR, Jung Y, Babor M, Schulten V, Seumois G, Greenbaum JA, Premawansa S, Premawansa G, Wijewickrama A, Vidanagama D, Gunasena B, Tippalagama R, de Silva AD, Gilman RH, Saito M, Taplitz R, Ley K, Vijayanand P, Sette A, Peters B. Circulating T cell-monocyte complexes are markers of immune perturbations. Elife. 2019;8:e46045.

    PubMed  PubMed Central  Article  Google Scholar 

  59. 59.

    Olingy CE, Dinh HQ, Hedrick CC. Monocyte heterogeneity and functions in cancer. J Leukoc Biol. 2019;106(2):309–22.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  60. 60.

    Akkari L, Bowman RL, Tessier J, Klemm F, Handgraaf SM, de Groot M, Quail DF, Tillard L, Gadiot J, Huse JT, Brandsma D, Westerga J, Watts C, Joyce JA. Dynamic changes in glioma macrophage populations after radiotherapy reveal CSF-1R inhibition as a strategy to overcome resistance. Sci Transl Med. 2020;12:552.

    Article  CAS  Google Scholar 

  61. 61.

    Lei F, Cui N, Zhou C, Chodosh J, Vavvas DG, Paschalis EI. CSF1R inhibition by a small-molecule inhibitor is not microglia specific; affecting hematopoiesis and the function of macrophages. Proc Natl Acad Sci. 2020;117(38):23336–8.

    CAS  PubMed  Article  Google Scholar 

  62. 62.

    Georgoudaki AM, Prokopec KE, Boura VF, Hellqvist E, Sohn S, Östling J, Dahan R, Harris RA, Rantalainen M, Klevebring D, Sund M, Brage SE, Fuxe J, Rolny C, Li F, Ravetch JV, Karlsson MC. Reprogramming tumor-associated macrophages by antibody targeting inhibits Cancer progression and metastasis. Cell Rep. 2016;15(9):2000–11.

    CAS  PubMed  Article  Google Scholar 

  63. 63.

    Jaynes JM, Sable R, Ronzetti M, Bautista W, Knotts Z, Abisoye-Ogunniyan A, Li D, Calvo R, Dashnyam M, Singh A, Guerin T, White J, Ravichandran S, Kumar P, Talsania K, Chen V, Ghebremedhin A, Karanam B, Bin Salam A, Amin R, Odzorig T, Aiken T, Nguyen V, Bian Y, Zarif JC, de Groot AE, Mehta M, Fan L, Hu X, Simeonov A, et al. Mannose receptor (CD206) activation in tumor-associated macrophages enhances adaptive and innate antitumor immune responses. Sci Transl Med. 2020;12:530.

    Article  CAS  Google Scholar 

  64. 64.

    Yanagita T, Murata Y, Tanaka D, Motegi SI, Arai E, Daniwijaya EW, Hazama D, Washio K, Saito Y, Kotani T, Ohnishi H, Oldenborg PA, Garcia NV, Miyasaka M, Ishikawa O, Kanai Y, Komori T, Matozaki T. Anti-SIRPα antibodies as a potential new tool for cancer immunotherapy. JCI Insight. 2017;2(1):e89140.

    PubMed  PubMed Central  Article  Google Scholar 

  65. 65.

    Klichinsky M, Ruella M, Shestova O, Lu XM, Best A, Zeeman M, Schmierer M, Gabrusiewicz K, Anderson NR, Petty NE, Cummins KD, Shen F, Shan X, Veliz K, Blouch K, Yashiro-Ohtani Y, Kenderian SS, Kim MY, O’Connor RS, Wallace SR, Kozlowski MS, Marchione DM, Shestov M, Garcia BA, June CH, Gill S. Human chimeric antigen receptor macrophages for cancer immunotherapy. Nat Biotechnol. 2020;38:947–53.

    CAS  PubMed  Article  Google Scholar 

  66. 66.

    Zhang L, Tian L, Dai X, Yu H, Wang J, Lei A, Zhu M, Xu J, Zhao W, Zhu Y, Sun Z, Zhang H, Hu Y, Wang Y, Xu Y, Church GM, Huang H, Weng Q, Zhang J. Pluripotent stem cell-derived CAR-macrophage cells with antigen-dependent anti-cancer cell functions. J Hematol Oncol. 2020;13(1):153.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  67. 67.

    Morrissey MA, Williamson AP, Steinbach AM, Roberts EW, Kern N, Headley MB, Vale RD. Chimeric antigen receptors that trigger phagocytosis. Elife. 2018;7:e36688.

    PubMed  PubMed Central  Article  Google Scholar 

  68. 68.

    Shields CW, Evans MA, Wang LL-W, Baugh N, Iyer S, Wu D, Zhao Z, Pusuluri A, Ukidve A, Pan DC, Mitragotri S. Cellular backpacks for macrophage immunotherapy. Sci Adv. 2020;6(18):eaaz6579.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  69. 69.

    Abousamra NK, Salah El-Din M, Helal R. Prognostic value of Th17 cells in acute leukemia. Med Oncol. 2013;30(4):732.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  70. 70.

    Han Y, Ye A, Bi L, Wu J, Yu K, Zhang S. Th17 cells and interleukin-17 increase with poor prognosis in patients with acute myeloid leukemia. Cancer Sci. 2014;105(8):933–42.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  71. 71.

    Li H, van der Leun AM, Yofe I, Lubling Y, Gelbard-Solodkin D, van Akkooi ACJ, van den Braber M, Rozeman EA, Haanen J, Blank CU, Horlings HM, David E, Baran Y, Bercovich A, Lifshitz A, Schumacher TN, Tanay A, Amit I. Dysfunctional CD8 T cells form a proliferative, dynamically regulated compartment within human melanoma. Cell. 2019;176(4):775–789 e718.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  72. 72.

    Sade-Feldman M, Yizhak K, Bjorgaard SL, Ray JP, de Boer CG, Jenkins RW, Lieb DJ, Chen JH, Frederick DT, Barzily-Rokni M, Freeman SS, Reuben A, Hoover PJ, Villani AC, Ivanova E, Portell A, Lizotte PH, Aref AR, Eliane JP, Hammond MR, Vitzthum H, Blackmon SM, Li B, Gopalakrishnan V, Reddy SM, Cooper ZA, Paweletz CP, Barbie DA, Stemmer-Rachamimov A, Flaherty KT, et al. Defining T cell states associated with response to checkpoint immunotherapy in melanoma. Cell. 2018;175(4):998–1013 e1020.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  73. 73.

    Kurachi M. CD8+ T cell exhaustion. Semin Immunopathol. 2019;41(3):327–37.

    PubMed  Article  Google Scholar 

  74. 74.

    Salik B, Smyth MJ, Nakamura K. Targeting immune checkpoints in hematological malignancies. J Hematol Oncol. 2020;13(1):111.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  75. 75.

    Wu K, Yi M, Qin S, Chu Q, Zheng X, Wu K. The efficacy and safety of combination of PD-1 and CTLA-4 inhibitors: a meta-analysis. Exp Hematol Oncol. 2019;8(1):26.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  76. 76.

    Oh DY, Kwek SS, Raju SS, Li T, McCarthy E, Chow E, Aran D, Ilano A, Pai C-CS, Rancan C, Allaire K, Burra A, Sun Y, Spitzer MH, Mangul S, Porten S, Meng MV, Friedlander TW, Ye CJ, Fong L. Intratumoral CD4+ T Cells Mediate Anti-tumor Cytotoxicity in Human Bladder Cancer. Cell. 2020;181(7):1612–1625.e1613.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  77. 77.

    Quezada SA, Simpson TR, Peggs KS, Merghoub T, Vider J, Fan X, Blasberg R, Yagita H, Muranski P, Antony PA, Restifo NP, Allison JP. Tumor-reactive CD4(+) T cells develop cytotoxic activity and eradicate large established melanoma after transfer into lymphopenic hosts. J Exp Med. 2010;207(3):637–50.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  78. 78.

    García-Guerrero E, Sánchez-Abarca LI, Domingo E, Ramos TL, Bejarano-García JA, Gonzalez-Campos JA, Caballero-Velázquez T, Pérez-Simón JA. Selection of tumor-specific cytotoxic T lymphocytes in acute myeloid leukemia patients through the identification of T-cells capable to establish stable interactions with the leukemic cells: “doublet technology”. Front Immunol. 2018;9:1971.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  79. 79.

    Weaver CT, Harrington LE, Mangan PR, Gavrieli M, Murphy KM. Th17: an effector CD4 T cell lineage with regulatory T cell ties. Immunity. 2006;24(6):677–88.

    CAS  PubMed  Article  Google Scholar 

  80. 80.

    DiToro D, Harbour SN, Bando JK, Benavides G, Witte S, Laufer VA, Moseley C, Singer JR, Frey B, Turner H, Bruning J, Darley-Usmar V, Gao M, Conover C, Hatton RD, Frank S, Colonna M, Weaver CT. Insulin-Like Growth Factors Are Key Regulators of T Helper 17 Regulatory T Cell Balance in Autoimmunity. Immunity. 2020;52(4):650–667.e610.

    CAS  PubMed  Article  Google Scholar 

  81. 81.

    Pesenacker AM, Bending D, Ursu S, Wu Q, Nistala K, Wedderburn LR. CD161 defines the subset of FoxP3+ T cells capable of producing proinflammatory cytokines. Blood. 2013;121(14):2647–58.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  82. 82.

    Mousset CM, Hobo W, Woestenenk R, Preijers F, Dolstra H, van der Waart AB. Comprehensive Phenotyping of T cells using flow Cytometry. Cytometry Part A. 2019;95(6):647–54.

    Article  Google Scholar 

  83. 83.

    Noack M, Miossec P. Th17 and regulatory T cell balance in autoimmune and inflammatory diseases. Autoimmun Rev. 2014;13(6):668–77.

    CAS  PubMed  Article  Google Scholar 

  84. 84.

    Chaudhry A, Samstein RM, Treuting P, Liang Y, Pils MC, Heinrich JM, Jack RS, Wunderlich FT, Brüning JC, Müller W, Rudensky AY. Interleukin-10 signaling in regulatory T cells is required for suppression of Th17 cell-mediated inflammation. Immunity. 2011;34(4):566–78.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  85. 85.

    Guillerey C, Harjunpää H, Carrié N, Kassem S, Teo T, Miles K, Krumeich S, Weulersse M, Cuisinier M, Stannard K, Yu Y, Minnie SA, Hill GR, Dougall WC, Avet-Loiseau H, Teng MWL, Nakamura K, Martinet L, Smyth MJ. TIGIT immune checkpoint blockade restores CD8(+) T-cell immunity against multiple myeloma. Blood. 2018;132(16):1689–94.

    CAS  PubMed  Article  Google Scholar 

  86. 86.

    Zhang Q, Bi J, Zheng X, Chen Y, Wang H, Wu W, Wang Z, Wu Q, Peng H, Wei H, Sun R, Tian Z. Blockade of the checkpoint receptor TIGIT prevents NK cell exhaustion and elicits potent anti-tumor immunity. Nat Immunol. 2018;19(7):723–32.

    CAS  PubMed  Article  Google Scholar 

  87. 87.

    McCully ML, Ladell K, Andrews R, Jones RE, Miners KL, Roger L, Baird DM, Cameron MJ, Jessop ZM, Whitaker IS, Davies EL, Price DA, Moser B. CCR8 expression defines tissue-resident memory T cells in human skin. J Immunol. 2018;200(5):1639–50.

    CAS  PubMed  PubMed Central  Google Scholar 

  88. 88.

    Bai R, Lv Z, Xu D, Cui J. Predictive biomarkers for cancer immunotherapy with immune checkpoint inhibitors. Biomark Res. 2020;8(1):34.

    PubMed  PubMed Central  Article  Google Scholar 

  89. 89.

    Gabriel SS, Kallies A. Tissue-resident memory T cells keep cancer dormant. Cell Res. 2019;29(5):341–2.

    PubMed  PubMed Central  Article  Google Scholar 

  90. 90.

    Chen Y, Wang Y, Luo H, Meng X, Zhu W, Wang D, Zeng H, Zhang H. The frequency and inter-relationship of PD-L1 expression and tumour mutational burden across multiple types of advanced solid tumours in China. Exp Hematol Oncol. 2020;9(1):17.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  91. 91.

    Li Y, You MJ, Yang Y, Hu D, Tian C. The role of tumor-associated macrophages in leukemia. Acta Haematol. 2020;143(2):112–7.

    CAS  PubMed  Article  Google Scholar 

  92. 92.

    Zhang L, Li Z, Skrzypczynska KM, Fang Q, Zhang W, O'Brien SA, He Y, Wang L, Zhang Q, Kim A, Gao R, Orf J, Wang T, Sawant D, Kang J, Bhatt D, Lu D, Li CM, Rapaport AS, Perez K, Ye Y, Wang S, Hu X, Ren X, Ouyang W, Shen Z, Egen JG, Zhang Z, Yu X. Single-cell analyses inform mechanisms of myeloid-targeted therapies in Colon Cancer. Cell. 2020;181(2):442–459 e429.

    CAS  Article  Google Scholar 

  93. 93.

    Hasegawa K, Tanaka S, Fujiki F, Morimoto S, Nakajima H, Tatsumi N, Nakata J, Takashima S, Nishida S, Tsuboi A, Oka Y, Oji Y, Kumanogoh A, Sugiyama H, Hosen N. An Immunocompetent mouse model for MLL/AF9 leukemia reveals the potential of spontaneous cytotoxic T-cell response to an antigen expressed in leukemia cells. PLoS One. 2015;10(12):e0144594.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  94. 94.

    Zhang L, Chen X, Liu X, Kline DE, Teague RM, Gajewski TF, Kline J. CD40 ligation reverses T cell tolerance in acute myeloid leukemia. J Clin Invest. 2013;123(5):1999–2010.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  95. 95.

    Li Z, Philip M, Ferrell PB. Alterations of T-cell-mediated immunity in acute myeloid leukemia. Oncogene. 2020;39(18):3611–9.

    CAS  PubMed  Article  Google Scholar 

  96. 96.

    Wang X, Zheng J, Liu J, Yao J, He Y, Li X, Yu J, Yang J, Liu Z, Huang S. Increased population of CD4(+)CD25(high), regulatory T cells with their higher apoptotic and proliferating status in peripheral blood of acute myeloid leukemia patients. Eur J Haematol. 2005;75(6):468–76.

    PubMed  Article  Google Scholar 

  97. 97.

    Williams P, Basu S, Garcia-Manero G, Hourigan CS, Oetjen KA, Cortes JE, Ravandi F, Jabbour EJ, Al-Hamal Z, Konopleva M, Ning J, Xiao L, Hidalgo Lopez J, Kornblau SM, Andreeff M, Flores W, Bueso-Ramos C, Blando J, Galera P, Calvo KR, Al-Atrash G, Allison JP, Kantarjian HM, Sharma P, Daver NG. The distribution of T-cell subsets and the expression of immune checkpoint receptors and ligands in patients with newly diagnosed and relapsed acute myeloid leukemia. Cancer. 2019;125(9):1470–81.

    CAS  PubMed  Article  Google Scholar 

  98. 98.

    Kanakry CG, Hess AD, Gocke CD, Thoburn C, Kos F, Meyer C, Briel J, Luznik L, Smith BD, Levitsky H, Karp JE. Early lymphocyte recovery after intensive timed sequential chemotherapy for acute myelogenous leukemia: peripheral oligoclonal expansion of regulatory T cells. Blood. 2011;117(2):608–17.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  99. 99.

    Zou W. Regulatory T cells, tumour immunity and immunotherapy. Nat Rev Immunol. 2006;6(4):295–307.

    CAS  PubMed  Article  Google Scholar 

  100. 100.

    Yu S, Liu C, Zhang L, Shan B, Tian T, Hu Y, Shao L, Sun Y, Ji C, Ma D. Elevated Th22 cells correlated with Th17 cells in peripheral blood of patients with acute myeloid leukemia. Int J Mol Sci. 2014;15(2):1927–45.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  101. 101.

    Ersvaer E, Liseth K, Skavland J, Gjertsen BT, Bruserud Ø. Intensive chemotherapy for acute myeloid leukemia differentially affects circulating TC1, TH1, TH17 and TREG cells. BMC Immunol. 2010;11(1):38.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  102. 102.

    Downs-Canner S, Berkey S, Delgoffe GM, Edwards RP, Curiel T, Odunsi K, Bartlett DL, Obermajer N. Suppressive IL-17A+Foxp3+ and ex-Th17 IL-17AnegFoxp3+ Treg cells are a source of tumour-associated Treg cells. Nat Commun. 2017;8(1):14649.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  103. 103.

    Liu D. Cancer biomarkers for targeted therapy. Biomark Res. 2019;7(1):25.

    PubMed  PubMed Central  Article  Google Scholar 

  104. 104.

    Lü M, Qiu L, Jia G, Guo R, Leng Q. Single-cell expression profiles of ACE2 and TMPRSS2 reveals potential vertical transmission and fetus infection of SARS-CoV-2. Aging (Albany NY). 2020;12(20):19880–97.

    Article  Google Scholar 

  105. 105.

    Cao Y, Su B, Guo X, Sun W, Deng Y, Bao L, Zhu Q, Zhang X, Zheng Y, Geng C, Chai X, He R, Li X, Lv Q, Zhu H, Deng W, Xu Y, Wang Y, Qiao L, Tan Y, Song L, Wang G, Du X, Gao N, Liu J, Xiao J, Su X-D, Du Z, Feng Y, Qin C, et al. Potent Neutralizing Antibodies against SARS-CoV-2 Identified by High-Throughput Single-Cell Sequencing of Convalescent Patients&#x2019; B Cells. Cell. 2020;182(1):73–84.e16.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  106. 106.

    Ziegler CG, Allon SJ, Nyquist SK, Mbano IM, Miao VN, Tzouanas CN, Cao Y, Yousif AS, Bals J, Hauser BM, Feldman J, Muus C, Wadsworth MH, Kazer SW II, Hughes TK, Doran B, Gatter GJ, Vukovic M, Taliaferro F, Mead BE, Guo Z, Wang JP, Gras D, Plaisant M, Ansari M, Angelidis I, Adler H, JMS S, Taylor CJ, Lin B, et al. SARS-COV-2 Receptor Ace2 is an interferon-stimulated gene in human airway epithelial cells and is detected in specific cell subsets across tissues. Cell. 2020;181(5):1016–1035.e1019.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

Download references

Acknowledgments

We thank Yimeng Du, Linna Chen, and Lingling Shang for assistance with thoughtful discussion.

Funding

This work was supported by the Postdoctoral Research Start-up Funding of the First Affiliated Hospital of Zhengzhou University (RQ. Guo), Postdoctoral Research Start-up Funding of Henan Province (RQ. Guo), Key Research and Development and Promotion Project of Henan province (RQ. Guo), National Natural Science Foundation of China (No. 82000112, W. Li; No. U1804192, YM. Li), Key scientific research projects of colleges and universities in Henan Province (No. 18B310026, ZX. Jiang), National Natural Science Foundation of Henan Province (No. 182300410301, R. Guo), and The medical science and technology research project of Henan province (No. 201701004, ZX. Jiang).

Author information

Affiliations

Authors

Contributions

RQ. Guo, DL. Liu, ZX. Jiang, and YP. Song conceived and supervised the project. CY. Du, YM. Li, R. Guo, and XS. Xie supervised the project. RQ. Guo, W. Li, MD. LV, FJ. Cao, GH. Wu, YY. Zhang, FC Gao, HL. Pang, YD. Li, HZ. Xing, CY. Liang, and TX. Lyu collected and analyzed the data. RQ. Guo wrote the manuscript. The authors read and approved the final manuscript.

Corresponding authors

Correspondence to Wei Li or Delong Liu or Yongping Song or Zhongxing Jiang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

No individual data were used in this study.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Table.

.

Additional file 2: Supplementary Figure 1.

Analysis of differences in gene expression between BM-derived cells in AML patients and healthy donors, dynamic changes of cell-type proportion, and survival curves of TCGA AML patients. Expression of Hallmark signatures: top genes (A), lineage marker genes (B), HSPC pattern genes (C) and monocyte pattern genes (D). E, dynamic changes of cell-type (CD8+-T/NK/NK, MKI67+ T, LYZhiEREGhi monocyte precursor, Mono/Mac, Erythroid lineage 1, Erythroid lineage 2, B cell, and Mast cell) proportion before and after treatment, and healthy donor-derived BM cells are represented at the end of plots. F, The Kaplan-Meier overall survival curves of TCGA AML patients grouped by the cluster-specific gene sets. + represents censored observations, and P value was calculated by multivariate Cox regression.

Additional file 3: Supplementary Figure 2.

Dot plot of differentially surface markers (A), transcription factors (B), pattern recognition receptors (C), cell adhesion/migration molecules (D), and chemokine receptors (E). F, The Kaplan-Meier overall survival curves of TCGA AML patients grouped by specific DC subset (pDC, CLEC7A+ DC, and CD1C+ DC) gene sets. + represents censored observations, and P value was calculated by multivariate Cox regression.

Additional file 4: Supplementary Figure 3. A,

UMAP plot of Monocyte/Macrophages from Fig. 1a-represented Mono/Mac cluster. These mature myeloid cells can be divided into 10 subsets before filtering possible cell-cell complexes. B, Expression levels of CD14 and CD3D across Mono/Mac population illustrated in UMAP plots. C, The Kaplan-Meier overall survival curves of TCGA AML patients grouped by specific subset gene sets. + represents censored observations, and P value was calculated by multivariate Cox regression.

Additional file 5: Supplementary Figure 4. A,

dynamic changes of proportion of distinct cell-types in total T/NK cells before and after treatment, and healthy donor-derived BM cells, as control, are represented at the end of plots. B, The Kaplan-Meier overall survival curves of TCGA AML patients grouped by specific NK/NKT-like gene set and IFN-CD4+ gene set. + represents censored observations, and P value was calculated by multivariate Cox regression.

Additional file 6: Supplementary Figure 5. A,

Violin plot showing the expression levels of HBA2, HBB, and LTB in 4 clusters (CD69highGZMA-CD4+ T, CD69highGZMA+CD4+ T, CD4+ Cytotoxic T, CD69lowLTBhighCD4+ T) from Fig. 5a-represented cells. B, The state-space of Naïve CD4+ T cluster, TH17-like cluster, and Treg cluster. Each dot corresponded to one single cell, colored according to its state (total 6 states). C, Expression maps showing log-normalized expression of typical markers (ITGB1, GZMA, IGF1R, and IL10RA) in the differentiation of Naïve CD4+ T to TH17-like cells and/or Treg cells. Data are shown as log-normalized expression. Yellow indicates high expression, dark blue indicates low expression. D, Typic state-space of some AML samples (AML420B.D0 and AML420B.D14; AML328.D0, AML328.D29, AML328.D113, and AML328.D171) are represented, and BM4 as healthy control. Each dot corresponded to one single cell, colored according to its state (total 6 states).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Guo, R., Lü, M., Cao, F. et al. Single-cell map of diverse immune phenotypes in the acute myeloid leukemia microenvironment. Biomark Res 9, 15 (2021). https://doi.org/10.1186/s40364-021-00265-0

Download citation

Keywords

  • Acute myeloid leukemia
  • Microenvironment
  • Single-cell RNA sequencing
  • Immune phenotypes
  • Bone marrow
  • Immune cells
  • Myeloid cells
  • T lymphocytes
\