Abstract
There is a strong correlation between myeloid-derived suppressor cells (MDSC) and resistance to immune checkpoint blockade (ICB), but the detailed mechanisms underlying this correlation are largely unknown. Using single-cell RNA sequencing analysis in a bilateral tumor model, we found that immunosuppressive myeloid cells with characteristics of fatty acid oxidative metabolism dominate the immune-cell landscape in ICB-resistant subjects. In addition, we uncovered a previously underappreciated role of a serine/threonine kinase, PIM1, in regulating lipid oxidative metabolism via PPARγ-mediated activities. Enforced PPARγ expression sufficiently rescued metabolic and functional defects of Pim1−/− MDSCs. Consistent with this, pharmacologic inhibition of PIM kinase by AZD1208 treatment significantly disrupted the myeloid cell–mediated immunosuppressive microenvironment and unleashed CD8+ T-cell–mediated antitumor immunity, which enhanced PD-L1 blockade in preclinical cancer models. PIM kinase inhibition also sensitized nonresponders to PD-L1 blockade by selectively targeting suppressive myeloid cells. Overall, we have identified PIM1 as a metabolic modulator in MDSCs that is associated with ICB resistance and can be therapeutically targeted to overcome ICB resistance.
Introduction
Immune checkpoint blockade (ICB) has generated remarkable responses in treating cancer (1, 2). However, this success has been limited to a small proportion of patients due to multiple resistance mechanisms, including T-cell exhaustion and an immunosuppressive tumor microenvironment (TME; refs. 3–5). Among these mechanisms, myeloid-derived suppressor cells (MDSC), which constitute the largest population of suppressor cells in the TME, are considered a major barrier to successful ICB therapy (4, 6). MDSCs adapt to the glucose-limited and lipid-rich TME by changing their energy consumption toward fatty acid oxidative (FAO) metabolism (7, 8). Although it is known that MDSCs in the TME increase expression of key FAO enzymes, mitochondrial mass and function, and oxygen consumption rate (OCR), which are all associated with elevated suppressive activity (7), the detailed mechanisms by which oxidative metabolism is regulated in MDSCs remain unclear.
To identify the cellular and molecular mechanisms that cause ICB resistance, we performed single-cell RNA sequencing (scRNA-seq) on immune cells from responder and nonresponder tumors following PD-L1 blockade. This led us to discover that the expression of Pim1 strongly correlated with increased FAO and PPARγ-driven lipid metabolism in MDSCs. PIM1 is a member of a serine/threonine protein kinase family that functions downstream of the JAK–STAT pathway (9–11). PIM kinases are constitutively active when expressed in cells, and their activity is directly correlated with their expression levels (9, 10). Previous research has shown that PIM kinases, predominately PIM1, regulate various aspects of cellular metabolism in myocardiocytes, adipocytes, and T cells (10, 12–14). However, the function of PIM1 in myeloid cells remains largely unknown. Furthermore, it has been shown that PIM kinases function as an emergency backup of the AKT–mTOR pathway in effector T cells and thereby are dispensable for normal T-cell activation, expansion, and function (10, 13). Taken together, these data suggest that targeting PIM-dependent metabolic pathways in immunosuppressive myeloid cells would have minimal effects on T-cell antitumor activity.
In this study, we have demonstrated that PIM1 was highly expressed in MDSCs with enhanced FAO metabolism, and this was strongly associated with ICB resistance in both mice and humans. Genetic deletion of Pim1 in MDSCs impaired fatty acid uptake and their immunosuppressive function by reducing the expression of PPARγ, a key regulator of lipid uptake and oxidative metabolism. Furthermore, pharmacologic inhibition of PIM kinases selectively diminished numbers of MDSCs in the TME, which in turn promoted the antitumor function of cytotoxic T cells and tumor control. Finally, PIM kinase inhibition not only improved the efficacy of PD-L1 blockade but also overcame ICB resistance in nonresponders.
Materials and Methods
Mice
Six- to eight-week-old male and female mice were used for tumor studies. C57BL/6 and C57BL/6-CD45.1 mice were purchased from Charles River. Pim1−/− mice were originally obtained from Dr. Robert Woodland (University of Massachusetts) and were backcrossed to C57BL/6 background for at least 12 generations in our animal facility prior to use. All mice were bred and maintained under the guideline of Institutional Animal Care and Use Committee of the Medical College of Wisconsin.
Tumor inoculation and treatments
MC38, B16-GMCSF (Gvax), and B16-Bl6 tumor cell lines were originally purchased from the ATCC and kindly provided by Dr. Lily Wang (Cleveland Clinic) in 2016. These cell lines have not been further authenticated or tested for mycoplasma since. All tumor cell lines were cultured in DMEM (Cat. #12-64-F, Lonza) supplemented with 10% FCS (Cat. #SH3007103, Hyclone), 2 mmol/L glutamine (Cat. #25-005-CI, Corning), and 100 U/mL penicillin/streptomycin (Cat. #30-002-CI, Corning) and passaged fewer than 8 times prior to inoculation. To establish the bilateral tumor model, 1 × 106 MC38 cells were s.c. injected into the left and right flanks of mice. Once the size of the tumors reached 200 mm3, mice were randomized and received the following treatments: (i) vehicle control, (ii) oral administration of AZD1208 (Selleckchem) for 10 days, (iii) three doses of i.p. injection with 200 μg anti–PD-L1 (clone #10F.9G2, BioXcell) every other day, or (iv) combination of AZD1208 and α–PD-L1 mAb. For the scRNA-seq study, tumors on the left side were surgically removed one day after the completion of anti–PD-L1 treatments to isolate immune cells, whereas tumors from the right side were kept intact to determine responder status.
For the melanoma model, 5 × 104 B16-BL6 cells were injected intradermally in the left flank. The tumor-bearing mice were treated on days 3, 6, and 9, with the following treatments: (i) injection of 106 irradiated (16,000 rads) B16-GMCSF cells (Gvax) on the right flank, (ii) Gvax treatment and oral administration of AZD1208 (30 mg/kg), (iii) Gvax treatment and i.p. injection with 200 μg anti–PD-L1, or (iv) combination of Gvax, AZD1208, and anti–PD-L1. Tumor volumes were measured 3 times a week and calculated as [longest dimension × (perpendicular dimension) 2]/2. Mice were euthanized when the tumor was greater than 2,000 mm3. Mice were considered cured when tumor size was lower than 10 mm3.
Immune-cell isolation
Tumors were harvested from mice and minced into small pieces, followed by a 1-hour incubation at 37°C with RPMI supplemented with 1% FBS, 2 mg/mL Collagenase Type I (Cat. #LS004196, Worthington Biochemical Corporation) and Collagenase Type XI (Cat. #LS004188, Worthington Biochemical Corporation), and 30 mg/mL DNase (Cat. #11284932001, Sigma-Aldrich). After incubation, tumor samples were mashed against a 70 μm cell strainer to harvest immune cells, which were subsequently enriched by Lymphocyte Cell Separation Medium (Cat. #25-072-CV, Corning) and lysed in ACK lysis buffer (Cat. #BW10548E, Lonza). Single-cell suspensions were then used for flow cytometry staining and FACS.
Flow cytometry
Single-cell suspensions were blocked with CD16/32 Fc blocking antibody (Cat. #101320, BioLegend) for 15 minutes and incubated with antibodies against surface markers for 30 minutes at 4°C, followed by three washes in FACS buffer. These samples were then run on an LSR II Green flow cytometer (BD Biosciences) and analyzed by FlowJo software (BD Biosciences).
Intracellular transcription factor staining was performed after surface staining. The cells were fixed with buffer from the True-Nuclear Transcription Factor Buffer Set (Cat. #424401, BioLegend) for 1 hour according to the manufacturer's protocol. Cells were then washed with permeabilization buffer and stained with antibodies against transcription factors in permeabilization buffer.
For intracellular cytokine staining, single-cell suspensions were stimulated with anti-CD3/CD28 (Cat. #100340 and 102116, respectively, BioLegend) in the presence of brefeldin A (Cat. #420601, BioLegend) and 10 ng/mL IL2 (Cat. #212-12, PeproTech) for 6 hours at 37°C. Subsequently, cells were collected for staining of surface markers and fixed by fixation buffer. Cells were then washed with permeabilization buffer and stained for cytokine-specific antibodies in permeabilization buffer.
All antibodies used for flow cytometry are listed in Supplementary Table S1.
scRNA-seq
Cell sorting for scRNA-seq
Immune cells were isolated from tumors and stained as previously described. Then, CD45+ (Cat. #103106, BioLegend) and live (7AAD−, Cat. #420404, BioLegend) cells were sorted using a BD FACSMelody cell sorter.
scRNA-seq library generation
About 1 × 104 sorted cells for each sample were loaded onto the 10× Chromium Controller (10× Genomics). The scRNA-seq libraries were generated by the Chromium Single Cell 3′ v2 Reagent Kit (Cat. #PN-120267, 10× Genomics) and sequenced using a NextSeq 500/550 High Output Kit v2 (150 cycles; Illumina) according to the manufacturer's protocol. The single-cell suspensions of live CD45+ cells were diluted in nuclease-free water and then loaded to Chromium Controller. RNA transcripts were uniquely barcoded and reverse-transcribed within nanoliter-scale droplets. Barcoded (10×) full-length cDNAs were then pooled and amplified via PCR. Then, the amplified cDNA went through an end repair process, A-tailing, adaptor ligation, and sample index PCR. The libraries were sequenced using a NextSeq 500/550 High Output Kit v2 (150 cycles; Illumina, 20024904).
scRNA-seq analysis
Raw sequencing data were demultiplexed and converted to gene-barcode matrices using the Cell Ranger (version 2.2.0) mkfastq and count functions, respectively (10x Genomics). The mouse reference genome mm10 was used for alignment. Data were further analyzed in R (version 3.6.1) using Seurat (version 3.1.0; ref. 15). A total of 3,214 cells from the nonresponder and 1,864 cells from the responder mice were recovered and merged into one Seurat object. The number of genes detected per cell, number of unique molecular identifiers (UMI), and percentage of transcripts derived from the mitochondria were plotted; cells that expressed less than 200 or more than 2,500 genes and cells with percent mitochondrial genes over 5% were removed to filter out doublets and cells with low read quality. Differences in cell library sizes (number of UMIs) and percentage of reads derived from the mitochondria were regressed out to prevent these technical variables from influencing cell clustering. Raw UMI counts were normalized and log-transformed. Principal component analysis was performed using variable genes, and the top 10 principal components were used for t-distributed stochastic neighbor embedding (t-SNE) analysis. These first 10 principal components were used to cluster the cells with Seurat's implementation of a shared nearest neighbor modularity optimization–based clustering algorithm (Louvain's original algorithm described in ref. 16). To identify marker genes, the FindAllMarkers function was used with Wilcoxon test for single-cell gene expression. For each cluster, only genes that were expressed in more than 25% of cells with at least 0.25-fold difference were considered. To characterize cell types, we performed annotation in SingleR for each single cell independently based on the ImmGen, GeneQuery, and Enrichr databases. Normalized data were used in feature plots or violin plots. Mean expression of markers inside each cluster was used to perform gene set enrichment analysis (GSEA) using the fgsea R package (17). The volcano plot was generated by R package EnhancedVolcano.
Bulk RNA sequencing
For each biological replicate, 0.5 × 106 live (7AAD−) bone marrow–derived MDSCs (BM-MDSC) were FACS-sorted. Total RNA was extracted with the RNeasy Plus Micro Kit per the manufacturer's protocol (Cat. #74034, QIAGEN). Library preparation was performed according to the Smart-seq2 protocol (18). Sequencing was performed on a NextSeq 500/550 High Output Kit v2 (75 cycles; Illumina) in a 37 bp paired-end mode. Sequenced reads were mapped to the mouse reference genome sequence (mm10) using TopHat v2.1.1 in combination with Bowtie2 v2.2.8 and Samtools v1.3. Fragments per kilobase of exon per million mapped fragments were calculated, and differential expression analysis was performed using Cufflinks v2.2.2. Heatmaps were generated with the R package pheatmap. GSEA was performed using the fgsea R package. Ingenuity pathway analysis (IPA; QIAGEN Bioinformatics) was used for bioinformative analysis of the data.
Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request. RNA sequencing (RNA-seq) data were deposited in the NCBI Gene Expression Omnibus database under accession codes GSE166308 and GSE166309.
Clinical data analysis
To further test if myeloid cells with gene signatures of high FAO metabolism correlated with ICB resistance, we analyzed scRNA-seq data from nivolumab-treated patients with melanoma (accession: GSE120575). Data from all 10 ICB responder and all 19 ICB nonresponder patients were used for analysis. Cells were clustered using Seurat in a manner similar to the murine data above. One cluster was identified as myeloid cells based on high expression of Integrin alpha M (ITGAM), also known as CD11b, and CD33. GSEA was performed on these cells after separating the cells by patient status (responder vs. nonresponder) using the clusterProfiler package (v 3.16.1) in R and gene sets from the Kyoto Encyclopedia of Genes and Genomes.
To correlate the expression of PIM1 and survival, the Web-based tool Gene Expression Profiling Interactive Analysis (GEPIA) was used to analyze datasets from The Cancer Genome Atlas (TCGA). Overall survival analysis was performed for colon adenocarcinoma (n = 270), skin cutaneous melanoma (n = 458), and combined 27 types of solid tumors (n = 3,834) datasets, and the significance was determined using the log-rank test in GEPIA.
Cleavage Under Targets and Tagmentation chromatin profiling
A total of 100,000 cells from Control and AZD groups were used for library construction using the Cleavage Under Targets and Tagmentation (CUT&Tag) chromatin immunoprecipitation sequencing approach (version 2; ref. 19; stepwise protocol can be found at https://www.protocols.io/view/bench-top-cut-amp-tag-z6hf9b6). Thirty-seven cycles of paired-end sequencing were performed on an Illumina Nextseq 500, and 5 to 10 million reads were generated for each sample. Each dataset was downsampled to equal read depths. For standardization between experiments, Escherichia coli DNA derived from transposase protein production was used to normalize sample read counts based on the recommendation of CUT&Tag protocol. Reads were aligned to the Mus musculus mm10 genome and E. coli (strain K12) using Bowtie2 (version 2.2.5; ref. 20) with options: –local –very- sensitive-local –no-unal –no-mixed –no-discordant –phred33 -I 10 -X 700. Peaks were called using Sparse Enrichment Analysis for CUT&RUN (SEACR; version 1.1; ref. 21), with options: 0.01 nonstringent. Peaks were visualized by Integrative Genomics Viewer (IGV; version 2.8.2).
Generation of BM-MDSCs and human MDSCs
Femurs were collected from C57BL/6 or Pim1−/− mice and flushed with 10% FBS RPMI medium to harvest bone marrow (BM) precursors. These cells were cultured at 37°C with 5% CO2 in the presence of GM-CSF (40 ng/mL, Cat. #200-15, Shenandoah) and IL-6 (40 ng/mL, 200-02, Shenandoah) for 7 days to generate BM-MDSCs. To generate human MDSCs, the mononuclear cell suspension was prepared from peripheral blood by differential density gradient separation (Ficoll-Hypaque; Cat. #10771, Sigma-Aldrich). These cells were cultured in complete RPMI 1640 media with 30% v/v tumor supernatants, 10% FBS, 1% pen/strep, GM-CSF (Cat. #100-08, Shenandoah), IL-6 (Cat. #100-10, Shenandoah), and G-CSF (Cat. #100-72, Shenandoah) for 4 days.
In vitro suppression assay
For the in vitro suppression assay, activated CD8+ T cells were cocultured with BM-MDSCs. The CD8+ T cells were isolated from naïve C57BL/6 mice (Cat. #19853, Stemcell) and labeled with CellTrace Violet (CTV, Cat. #C34557, Thermo Fisher) followed by activation with anti-CD3/28 mAb. Then BM-MDSCs were added to the culture at different ratios. The proliferation of CTV-labeled CD8+ T cells was analyzed by flow cytometry.
Adoptive transfer of BM-MDSCs
A total of 2 × 106 BM-MDSCs were adoptively transferred into each C57BL/6 CD45.1+ congenic mouse bearing an MC38 tumor. Two weeks later, tumors from recipient mice were harvested and processed for flow cytometry analysis as described above.
Fatty acid uptake assay
BM-MDSCs were incubated with BODIPY FL C16 (1 μg/mL, Cat. #D382, Invitrogen) in glucose-free RPMI 1640 medium (Cat. #11879020, Thermo Fisher Scientific) for 30 minutes at 37°C, and flow cytometry was used to measure its uptake.
Seahorse assay
The Seahorse XF Cell Mito Stress Test Kit (Cat. #103015-100, Seahorse Bioscience) was used according to the manufacturer's protocol to determine the mitochondrial function of BM-MDSCs by measuring the OCR with an XF96 analyzer (Seahorse Bioscience). About 2 × 105 cells were seeded per well in a XF96 cell culture microplate, and compounds were injected during the assay at the following final concentrations: 1.0 μmol/L oligomycin, 1.5 μmol/L carbonyl cyanide p-(trifluoromethoxy) phenylhydrazone (FCCP), 0.5 μmol/L rotenone, 0.5 μmol/L antimycin A, and 100 μmol/L etomoxir. The XF Palmitate-BSA FAO substrate (Cat. #102720-100, Seahorse Bioscience) was used to evaluate the oxidation of exogenously added fatty acids and the proportion of respiration supported by exogenous fatty acids according to the manufacturer's protocol (Seahorse Bioscience, Cat. #102720-100). The BM-MDSCs were inoculated with either palmitate–BSA (200 μmol/L palmitate conjugated with 34 μmol/L BSA) or BSA (34 μmol/L; Seahorse Bioscience), and the OCR was analyzed as described above.
Retroviral transduction
The coding sequences of mouse Pparg were obtained from Addgene (Cat. #17442) and cloned into MSCV-IRES-Thy1.1 (MIT). To produce retroviral supernatant to express Pparg, 293T cells were transfected with either MIT-empty or MIT-PPARg vector along with the pcLEco ecotropic packaging plasmid. After 48 hours, the retrovirus was collected from the culture supernatant and added to Pim1−/− BM-MDSCs followed by spinning transduction. After an additional 3 days of culturing in GM-CSF and IL-6, the positively transduced cells, defined by expression Thy1.1, were sorted for transfer on a FACS Aria III flow cytometer (BD Biosciences).
Statistical analysis
All experiments were performed using randomly assigned mice without investigator blinding. All data points and P values reflect biological replicates from at least three independent experiments. Statistical analysis was performed using GraphPad PRISM 7. Unpaired, two-tailed Student t tests and one-way ANOVA tests with post hoc Tukey−Kramer corrections were used to assess statistical significance. For examining the frequencies of PIM1- and GPR84-expressing myeloid cells in patients, the Fisher exact test was used to evaluate the statistical significance of the ratios of two groups, each with two conditions (e.g., ratio of PIM1+: PIM1− cells in responders vs. nonresponders). P values < 0.05 were considered to indicate a significant difference. The classification and regression trees (CART) is a supervised machine learning method that builds a tree-like predictive model from data and is used to estimate the percentage of MDSCs for predicting responsiveness. The tree is built based on splitting the samples into two subgroups. The splitting rule is to search for a best variable that achieves either the best separation in different groups (classification) or the minimum regression error (regression). The tree will grow as this process is repeated in each subgroup recursively, and will stop when further split does not make meaningful improvement.
Results
Single-cell transcriptomics reveal an association between the immune landscape and ICB resistance
To dissect the transcriptional and metabolic mechanisms of ICB resistance, we employed scRNA-seq to analyze the tumor-infiltrating immune cells in a murine MC38 colorectal cancer model. Consistent with previous work, we found that approximately 50% to 60% of MC38-bearing mice treated with anti–PD-L1 experienced complete tumor regression (22). To obtain immune cells for scRNA-seq analysis before tumors disappeared in responders, we established a bilateral tumor model by inoculating MC38 on both flanks of the mice and demonstrated the response to PD-L1 blockade in a symmetric manner, as previously published (ref. 23; Supplementary Fig. S1A). This model enabled us to examine the immune cells from the left flank tumor via surgical removal after the last dose of anti–PD-L1 administration, when the treatment outcomes were still indistinguishable, yet continuously monitor tumor regression or progression on the right flank to determine responder or nonresponder status, respectively (Fig. 1A).
Single-cell transcriptomics reveal the distinct immune-cell landscapes from an ICB responder and nonresponder. A, Experimental design for establishing a bilateral tumor model to study the immune-cell landscape associated with resistance to ICB immunotherapy in the MC38 colorectal cancer model. B, t-SNE plot showing 12 immune-cell clusters pooled from a responder and a nonresponder. C, Violin plots showing the gene expression distribution of various cell markers in all immune-cell clusters. D, t-SNE plots (top) demonstrating the different tumor-infiltrating immune cells between the nonresponder (left) and responder (right). Proportions of immune-cell subsets found in the tumors are shown using pie charts (bottom). ILC, innate lymphoid cell; NK, natural killer. E, The frequencies of intratumoral CD11b+Gr1+ MDSCs detected by flow cytometry (top) were plotted against the tumor volumes (bottom) for correlation (n = 38, Spearman correlation test, P < 0.05, r = 0.34).
After confirming that the bilateral tumors shared similar immune-cell compositions and phenotypes (Supplementary Fig. S1B and S1C), we next examined how the immune landscape differed between responders and nonresponders. Two days following the end of anti–PD-L1 administration, single-cell transcriptomic datasets from CD45+ immune cells from a responder and a nonresponder were successfully obtained. We recovered 1,864 cells from the responder and 3,214 cells from the nonresponder with an average recovery of approximately 1,500 genes/cell. These separated into 12 clusters using unbiased dimensionality reduction and clustering algorithms in the Seurat R package (Fig. 1B; refs. 24, 25). The expression profile of various myeloid (Adgre1, which encodes F4/80, Itgam, which encodes CD11b, and S100a9) and lymphoid (Cd8a, Cd4, and Ncr1) lineage markers allowed us to assign biological identities to each cluster (Fig. 1C). We identified four distinct CD8+ T-cell clusters (clusters 1, 2, 7, and 8), two distinct CD4+ T-cell clusters (clusters 4 and 11), three myeloid-cell clusters (clusters 0, 5, and 6), one neutrophil cluster (cluster 10), and two natural killer–cell clusters (clusters 3 and 9; Fig. 1B).
To delineate the relationship between the immune landscape and responsiveness to PD-L1 blockade, we performed comparative analysis between the responder and nonresponder (Fig. 1D). We observed that myeloid cells dominated in the nonresponder tumor (61.22%), whereas the percentage of CD8+ T cells (19.26%) was low (Fig. 1D). In contrast, the responder's tumor was predominantly infiltrated by CD8+ T cells (50.75%), whereas myeloid cells comprised only 18.03% of all CD45+ immune cells (Fig. 1D). We also found that tumor progression following PD-L1 blockade positively correlated with the frequency of CD11b+Gr1+ tumor-infiltrating myeloid cells (Fig. 1E).
Given that ICB resistance is associated with tumor-reactive CD8+ T-cell dysfunction (26), we compared the differentially expressed genes in CD8+ T cells from the responder and nonresponder (Supplementary Fig. S1D–S1F). This GSEA revealed that CD8+ T cells from the responder expressed higher levels of genes associated with T-cell costimulation and effector function, such as Cd28, Tnfsf4, Icos, Gzmb, Fasl, and Prf1 (Supplementary Fig. S1E), and had an enriched transcriptional signature of effector CD8+ T cells (Supplementary Fig. S1F). In contrast, CD8+ T cells from the nonresponder expressed higher levels of genes associated with exhaustion such as Pdcd1 (Supplementary Fig. S1E) and positively correlated with an exhausted T-cell signature (Supplementary Fig. S1F).
Myeloid cells with high FAO metabolism correlate with ICB resistance
To gain a more in-depth understanding of how tumor-infiltrating myeloid cells are metabolically regulated in the context of ICB resistance, we further analyzed the differentially expressed genes in the myeloid clusters between responder and nonresponder tumors. Tumor-infiltrating myeloid cells from the nonresponder exhibited higher expression of immunosuppressive molecules such as Arg1, Tgfb1, and Cd274 (Fig. 2A). Furthermore, GSEA revealed that transcriptional signatures related to intratumoral MDSCs were significantly enriched in the nonresponder (Fig. 2B). In contrast, many proinflammatory genes, such as Ccl9, Socs3, and Tnf, were upregulated in the responder along with enriched TNFα signaling, consistent with potential antitumor function.
Heightened FAO metabolism in myeloid cells correlates with ICB resistance. A, Violin plots showing differentially expressed genes in myeloid cells between responder (R) and nonresponder (NR) samples. B, GSEA plots evaluating expression in nonresponders and responders of published datasets of genes downregulated in spleen versus tumor monocytes (top) and genes involved in TNFα signaling (bottom). NES, normalized enrichment score. C, Scatter plot showing the relationship between immunosuppressive activity and FAT activity. The y-axis represents relative expression of immunosuppressive gene signature, and the x-axis represents the relative expression of FAT-related gene signature. Each dot represents a single cell from the myeloid populations pooled from responder (black) and nonresponder (red). D, Volcano plot showing differentially expressed genes between myeloid cells with high or low FAT activity. Each dot denotes one gene. E, The correlation of Pim1 expression and genes related to FAT. F, GSEA of CD33+ ITGAM+ (encodes CD11b) myeloid cells from scRNA-seq data derived from patients with melanoma treated with checkpoint inhibitor therapy reveals the enrichment of fatty acid metabolism and oxidative phosphorylation gene signatures in nonresponder myeloid cells. G, Box plots showing the frequencies of myeloid cells expressing PIM1, GPR84, and PIM1 and GRP84 from nonresponders (n = 19) and responders (n = 10). Statistical significance was determined using the Fisher exact test with a confidence level of 95%. ***, P < 0.001; ****, P < 0.0001.
To better understand mechanisms by which tumor-infiltrating myeloid cells adapt to the nutrient-imbalanced TME, we sought to identify which metabolic processes support their survival and immunosuppressive function. GSEA showed that the signature of oxidative phosphorylation (OXPHOS) was significantly enriched in myeloid cells from the nonresponder (Supplementary Fig. S2A). As OXPHOS in MDSCs is mainly fueled by exogenous fatty acids (27), we explored the association between the fatty acid transport (FAT) pathway and immunosuppressive function using the R package gene set variation analysis to quantify the activity of these biological pathways at the single-cell level (28). The activity score of immunosuppressive activity was calculated for each cell based on the expression level of the genes encoding inhibitory modules, and fatty acid metabolic activity was determined by the level of genes in the Gene Ontology FAT signature. As expected, we found a significant positive correlation between these two pathways (Fig. 2C), which corroborated previous findings showing an essential role for lipid metabolism in the inhibitory functions of MDSCs (29).
By comparing the differentially expressed genes between total myeloid cells with high and low FAT activity, we found that Cd36 and Pim1 were highly upregulated in the high FAT group (Fig. 2D). Previous research has shown that PIM kinases, predominantly PIM1, regulate various aspects of cellular metabolism in myocardiocytes, adipocytes, and T cells (10, 12–14). Thus, we hypothesized that PIM1 regulates FAT activity in tumor-infiltrating myeloid cells for their metabolic adaptation and function in the TME. In support of this idea, Pim1 expression significantly correlated with the expression of FAT-related genes such as Cd36 and Gpr84, which encodes GPR84, a G protein–coupled receptor for medium- and long-chain fatty acids (Fig. 2E; refs. 30, 31). Furthermore, we stratified myeloid cells into Pim1high and Pim1low cells based on their Pim1 expression and identified differentially regulated genes between the two. GSEA demonstrated that Pim1high cells were highly enriched for FAT genes and IL-6/STAT3 signaling pathway genes (Supplementary Fig. S2B). Taken together with the fact that PIM1 is a well-known downstream target of STAT3 signaling (32, 33) and that IL-6 is abundant in the TME (34), we hypothesize that the IL-6/STAT3 pathway is a major driver of PIM1 upregulation in myeloid cells in the TME.
To further test if myeloid cells with gene signatures of high FAO metabolism correlated with ICB resistance, we analyzed scRNA-seq data from nivolumab-treated patients with melanoma (ref. 35; accession: GSE120575). Similar to our murine colorectal carcinoma data, tumor-infiltrating myeloid cells in nonresponder patients exhibited significant enrichment of fatty acid metabolism and oxidative phosphorylation genes (Fig. 2F). In addition, greater frequencies of PIM1- and/or GPR84-expressing myeloid cells in patients significantly and positively correlated with ICB resistance (Fig. 2G). Expression of PIM1 and GPR84 significantly and positively correlated with each other (Supplementary Fig. S2C), suggesting that these could be used as biomarkers either alone or in combination to predict the outcome of ICB therapy. To further evaluate the correlation between PIM1 expression and patient survival, we used the GEPIA Web server to analyze RNA-seq expression data from tumor samples from the TCGA and the GTEx projects (36). In the 27 types of solid tumors cohorts analyzed, high levels of PIM1 were associated with reduced disease-free survival rates, whereas lower expression of PIM1 was associated with significantly longer survival (P = 0.002 by log-rank test; Supplementary Fig. S2D). Collectively, our results suggest that PIM1 expression and associated metabolic reprogramming in tumor-infiltrating myeloid cells might be responsible for ICB resistance.
PIM1 is required for the immunosuppressive properties of myeloid cells
To test the necessity of PIM1 for the function of immunosuppressive myeloid cells, we analyzed BM-MDSCs (37). In this in vitro model, expression of PIM1 increased with the gradual enrichment of MDSCs by day 7 (Supplementary Fig. S3A). Lack of PIM1 led to a reduced formation of total CD11b+Gr1+ MDSCs and frequency of Ly6G+CD11b+Gr1+ polymorphonuclear MDSCs (PMN-MDSC), but it did not affect Ly6C+CD11b+Gr1+ monocytic MDSCs (M-MDSC; Supplementary Fig. S3A). To further assess the role of PIM1 in MDSCs, we sorted CD11b+Gr1+ cells and performed bulk RNA-seq analysis to compare gene expression profiles between wild-type (WT) and Pim1−/− BM-MDSCs. The expressions of genes encoding inhibitory molecules, such as indolamine-2, 3-dioxygenase (Ido1), arginase 1 (Arg1), TGFβ1 (Tgfb1), and PD-L1 (Cd274), were substantially reduced in the absence of PIM1 (Fig. 3A). In agreement with these observations, flow cytometry analysis revealed that protein levels of Arg1 were significantly decreased in Pim1−/− BM-MDSCs (Supplementary Fig. S3B). In addition, the ability of Pim1−/− MDSCs to inhibit T-cell proliferation in vitro was significantly compromised (Fig. 3B). Likewise, pharmacologic inhibition of PIM kinases by AZD1208, a potent and selective pan-PIM kinase inhibitor (38, 39), also led to similar functional defects of BM-MDSCs (Supplementary Fig. S3C and S3D). Together, these results uncover an essential role of PIM1 in regulating the suppressive function of MDSCs in vitro.
PIM1 is required for the immunosuppressive phenotype of MDSCs. A, Heatmap showing expression of genes associated with immunosuppressive functions in BM-MDSCs generated from WT (n = 3) and Pim1−/− (n = 3) mice. B, WT (n = 3) and Pim1−/− (n = 3) MDSC-suppressive activities were evaluated by the ability to inhibit CD8+ T-cell proliferation. Bar graphs were plotted from three independent experiments. C–K, WT (n = 3) and Pim1−/− (n = 5) BM-MDSCs were purified and transferred into CD45.1+ congenic recipient mice bearing an MC38 tumor. C, Experimental design. D, Levels of Arg1 protein expression in intratumoral myeloid cells derived from WT and Pim1−/− donor cells. gMFI, geometric mean fluorescence intensity. E, Representative plot showing the frequencies of M2 macrophages evaluated by CD206 and Relmα expression from WT and Pim1−/− donor cells. F, Quantification of E. G and H, Protein expression levels of MHCII and CD86. I and J, Intratumoral CD8 T cells from WT and Pim1−/− MDSC recipients were stimulated and evaluated for IFNγ production by flow cytometry as shown in representative (I) and quantitative (J) plots. K, Tumor growth was measured using calipers for 2 weeks and plotted. Data are expressed as mean ± SEM from two experiments with three mice per group, and significance was determined by t test. *, P < 0.05; **, P < 0.01; ***, P < 0.001.
Next, we evaluated whether the loss of PIM1 in myeloid cells would affect the immunosuppressive TME in vivo. The necessary role of PIM1 in hematopoietic stem cell homing and engraftment precludes the possibility of generating BM chimeric mice with specific deletion of PIM1 in myeloid cells (40, 41). To overcome this challenge, we adoptively transferred WT or Pim1−/− (CD45.2+) BM-MDSCs into congenic (CD45.1+) WT mice (Fig. 3C). We then inoculated these two groups of recipient mice with MC38 tumor cells. Two weeks later, we found comparable frequencies of CD11b+ cells between adoptively transferred CD45.2+ WT and Pim1−/− cells in tumors (Supplementary Fig. S3E). The majority of donor cells isolated from tumors had differentiated into F4/80+CD11b+ tumor-associated macrophages (TAM; Supplementary Fig. S3F and S3G). We then assessed whether PIM1 deficiency affected the ability of transferred myeloid cells to inhibit antitumor immunity. First, we found that protein levels of Arg1 were significantly reduced in CD11b+ intratumoral myeloid cells from Pim1−/− compared with WT donor (Fig. 3D). Furthermore, Pim1−/− CD11b+ cells exhibited significantly lower expression of the M2 macrophage markers Relmα and CD206 (Fig. 3E and F), but expressed significantly higher levels of CD86 and MHC II than their WT counterparts (Fig. 3G and H). This is consistent with our in vitro observations that Pim1−/− BM-derived macrophages (BMDM) exhibited higher expression of M1 markers than WT BMDMs (Supplementary Fig. S3H). Furthermore, another important M1 gene signature, antigen presentation, was also enriched in Pim1−/− BMDMs (Supplementary Fig. S3I). These data suggest that PIM1-deficient myeloid cells fail to maintain their immunosuppressive phenotype and instead differentiate toward a more immune-stimulating state. Supporting this notion, we observed that tumor-infiltrating CD8+ T cells from mice that received Pim1−/− MDSCs produced more effector cytokines, such as IFNγ, than those from mice that received WT MDSCs (Fig. 3I and J). As a result, tumor progression in mice that received Pim1−/− MDSCs was significantly slower than that in mice that received WT MDSCs (Fig. 3K). Collectively, these data revealed a previously unappreciated role of PIM1 in regulating immunosuppressive properties of MDSCs in vivo.
PIM1 regulates FAO metabolism in MDSCs
To evaluate if PIM1-mediated immunosuppression relates to cellular metabolism of MDSCs, we compared the metabolic gene expression profiles of WT and Pim1−/− BM-MDSCs. Genetic ablation of Pim1 resulted in reduced gene expression related to multiple aspects of fatty acid metabolism, including fatty acid degradation (e.g., Acsl1, Acadm, Acca2, and Hadh), fatty acid transportation (e.g., Cpt1a, Cpt1c, and Cpt2), and mitochondrial beta-oxidation (e.g., Acox1, Acads, and Acat1; Fig. 4A). GSEA confirmed that fatty acid metabolism and oxidative phosphorylation signatures were suppressed in Pim1−/− BM-MDSCs (Fig. 4B). To determine whether these transcriptional alternations led to any metabolic changes, we first measured fatty acid uptake upon the addition of fluorescently labeled palmitate (Bodipy FL C16) and observed that Pim1−/− MDSCs exhibited diminished fatty acid uptake (Fig. 4C and D). We then measured OCR using a Seahorse extracellular flux analyzer and found that WT MDSCs had significantly higher basal respiratory capacity and nonsignificantly higher spare respiratory capacity than Pim1−/− MDSCs (Fig. 4E and F). Exposure to exogenous BSA-conjugated XF Palmitate substrate resulted in a significant increase in the basal OCR in WT MDSCs, whereas such an increase was largely diminished in PIM1-deficient MDSCs (Fig. 4G and H). Conversely, glycolytic pathways, as measured by extracellular acidification rate (ECAR) and glucose uptake, were not affected by PIM1 deficiency (Supplementary Fig. S4A–S4C). Consistent with these results, pharmacologic inhibition of PIM1 by AZD1208 in BM-MDSCs resulted in similar metabolic phenotypes (Supplementary Fig. S4D–S4G). Overall, these data suggest that PIM1 regulates FAO metabolism in MDSCs, at least in part, by regulating fatty acid uptake.
Pim1−/− BM-MDSCs exhibit impaired FAO metabolism. A, Heatmap depicting genes related to fatty acid metabolism differentially expressed between WT (n = 3) and Pim1−/− (n = 3) BM-MDSCs. B, GSEA demonstrating that Pim1−/− BM-MDSCs negatively correlate with gene signatures related to FAO and OXPHOS metabolism. NES, normalized enrichment score. C, Histograms showing representative fatty acid uptake from WT (n = 3) and Pim1−/− (n = 3) BM-MDSCs as assessed by flow cytometry. D, Quantification of C. gMFI, geometric mean fluorescence intensity. E, Line graphs depicting the OCR of WT and Pim1−/− BM-MDSCs in response to the Mito Stress assay. Oligo, oligomycin; R+A, rotenone + antimycin A. F, Bar graphs quantifying basal respiratory capacity (BRC) and spare respiratory capacity (SRC) in WT and Pim1−/− BM-MDSCs. N.S., not significant. G and H, Graphs depicting OCR from WT and Pim1−/− BM-MDSCs following exposure to exogenous XF Palmitate-BSA. Data are expressed as mean ± SEM from two experiments, and significance was determined by t test. *, P < 0.05; ***, P < 0.001; ****, P < 0.0001.
The PIM1–PPARγ axis regulates FAO metabolism in MDSCs
To explore the possible molecular mechanisms by which PIM1 affects the metabolism and function of MDSCs, we performed Ingenuity Pathway Analysis (IPA) using differentially expressed genes between WT and Pim1−/− BM-MDSCs. This analysis revealed that PPAR signaling was significantly reduced by the ablation of PIM1 (Fig. 5A). This was further confirmed by GSEA (Fig. 5B). The expression of PPARγ targets, which include genes that encode transcription factors related to lipid metabolism, FAO enzymes, and lipid transporters, were substantially downregulated in Pim1−/− MDSCs (Fig. 5C). Furthermore, protein levels of PPARγ and one of its targets, CD36, were significantly downregulated in Pim1−/− MDSCs (Fig. 5D–G). Together, these results suggest that PPARγ may function as a major downstream target of PIM1 to regulate fatty acid metabolism in MDSCs.
The PIM1–PPARγ signaling axis regulates cellular metabolism and function of MDSCs. IPA (A) and GSEA (B) on bulk RNA-seq data from WT (n = 3) and Pim1−/− (n = 3) MDSCs reveal that PPARγ signaling is significantly disturbed in the absence of Pim1. NES, normalized enrichment score. C, Heatmap showing gene expression of Pparg and its targets in WT and Pim1−/− cells. D–G, Flow cytometry analyses of CD36 (D and E) and PPARγ (F and G) expression in WT and Pim1−/− cells. gMFI, geometric mean fluorescence intensity. H–L, Metabolic and functional analysis of Pim1-deficient MDSCs following ectopic expression of Pparg. Expression levels of CD36 (H and I), fatty acid uptake (J and K), and in vitro suppression activity (L) were measured by flow cytometry in Pim1−/− BM-MDSCs expressing empty vector (MIT; n = 3) or Pparg (n = 3). Data are expressed as mean ± SEM from two experiments, and significance was determined by t test. *, P < 0.05; **, P < 0.01; ***, P < 0.001.
To determine whether reduced PPARγ expression was responsible for the metabolic and functional abnormalities in the PIM1-deficient MDSCs, we attempted to rescue these defects by ectopically expressing PPARγ in Pim1−/− BM-MDSCs (Supplementary Fig. S5A). We found that forced expression of PPARγ bypassed the need for PIM1, restored the expression of CD36 (Fig. 5H and I), and subsequently led to increased fatty acid uptake (Fig. 5J and K). PPARγ overexpression also significantly increased Arg1 protein expression (Supplementary Fig. S5B) and improved the ability of Pim1−/− MDSCs to inhibit T-cell proliferation (Fig. 5L). Taken together, these findings support a mechanism by which PIM1 promotes FAO metabolism and immunosuppressive function in MDSCs via the induction of PPARγ expression.
To determine the mechanisms by which PIM1 regulates the transcription of Pparg, the gene locus of which is known to contain STAT3-binding sites (42), we postulated that PIM1 could regulate STAT3 phosphorylation. In support of this idea, AZD1208 treatment inhibited phosphorylation of the serine 727 (S727) site without affecting phosphorylation of tyrosine 705 (Y705) or total STAT3 protein expression levels (Supplementary Fig. S5C). To further test whether this altered phosphorylation would affect the binding of STAT3 to the cis-regulatory elements of the Pparg locus, CUT&Tag chromatin profiling was performed on MDSCs treated with either PBS or AZD1208. Consistent with published work (42), we observed similar binding patterns of STAT3 to the enhancer regions of Pparg and Csf2ra in the PBS-treated control group (Supplementary Fig. S5D). In contrast, AZD1208 treatment significantly diminished STAT3 binding to the enhancer regions of Pparg without affecting binding to the Csf2ra locus (Supplementary Fig. S5D). These observations collectively suggest that PIM1-mediated phosphorylation of STAT3 at S727 enhances STAT3 transcriptional activity and results in heightened PPARγ expression in MDSCs.
PIM kinase inhibition improves the antitumor efficacy of PD-L1 blockade
Previous research has demonstrated that PIM kinases function as an emergency backup of the AKT–mTOR pathway in effector T cells and are therefore dispensable for normal T-cell activation, expansion, and function (10, 13). Thus, we reasoned we could selectively target PIM-dependent metabolic pathways in immunosuppressive myeloid cells while sparing the antitumor activity of T cells within the tumor. To test this idea, we first verified that pharmacologic inhibition of PIM kinase by AZD1208 did not have direct toxicity against MC38 tumor cells in vitro compared with AZD1208-sensitive AML-AF9 leukemia cells (Supplementary Fig. S6A). Next, we treated MC38 tumor–bearing mice with AZD1208 for 9 days, which sufficiently induced complete tumor regression in 30% of treated mice, whereas anti–PD-L1 treatment alone led to a complete response in 50% of treated mice (Fig. 6A). Coadministration of AZD1208 and anti–PD-L1 induced robust tumor eradication in 70% of mice, greater than either single-agent treatment alone (Fig. 6A). These findings suggest that pharmacologic inhibition of PIM kinases leads to an effective antitumor immune response and enhances ICB response.
PIM kinase inhibition enhances PD-L1 blockade for cancer immunotherapy. A, Bilateral MC38 tumor–bearing mice received the following treatments: vehicle control (n = 9), AZD1208 (n = 10), anti–PD-L1 (n = 10), or a combination of AZD1208 and anti–PD-L1 (n = 10). Individual tumor growth curves are shown. Dotted lines indicate the level of detection. B–I, Tumors on the left flank were surgically removed for immune-cell analysis. Representative flow plots (B) and quantification of CD11b+Gr1+F4/80low MDSCs (C) and CD11b+Gr1−F4/80high macrophages (D) are shown. ns, not significant. E and F, Within the MDSCs, frequencies of CD36highPPARγhigh cells are shown in representative flow plots (E) and quantitative scatter plot (F). G, Scatter plot depicting Arg1 levels as determined by flow cytometry. gMFI, geometric mean fluorescence intensity. H–J, Representative flow plots and quantitative scatter plots of IFNγ and Gzmb expression in tumor-infiltrating CD8+ T cells from four treatment groups. K, Kaplan–Meier survival curves show B16-Bl6 melanoma mice treated with vehicle control, anti–PD-L1, AZD1208, or a combination of AZD1208 and anti–PD-L1. Data represent cumulative results from two independent experiments with n = 9 for vehicle control and combination groups and n = 10 for AZD1208 and α–PD-L1 mAb groups. Cellular significance was determined by t test. Survival curve was analyzed by the log-rank (Mantel–Cox) test. *, P < 0.05; **, P < 0.01; ***, P < 0.001; ****, P < 0.0001.
To delineate the cellular mechanisms underlying the observed effects of PIM kinase inhibition, we analyzed the immune cells in the tumors. AZD1208 alone or AZD1208 + anti–PD-L1 treatment resulted in a significant reduction in the abundance of CD11b+Gr1+ MDSCs (Fig. 6B and C) and M-MDSCs (Supplementary Fig. S6B and S6C) and mildly increased PMN-MDSCs (Supplementary Fig. S6B and S6C) and CD11b+F4/80+ macrophages (Fig. 6D). To further evaluate the effects of PIM kinase inhibition on myeloid cells, protein levels of CD36 and PPARγ were measured to quantify FAO metabolic activity. Consistent with our in vitro data, CD36highPPARγhigh myeloid cells expressed higher levels of PIM1 (Supplementary Fig. S6D) and were significantly decreased in abundance in AZD1208-treated groups, both alone and in combination with anti–PD-L1 (Fig. 6E and F). We also observed downregulation of Arg1 in these two groups (Fig. 6G) and reduced ability to suppress CD8+ T cells (Supplementary Fig. S6E). Next, we characterized the phenotypes of TAMs. Although PIM inhibition has minimal impact on the frequency of CD11b+F4/80high TAMs, their expression of immunosuppressive M2 macrophage markers Relmα and CD206 was significantly reduced (Supplementary Fig. S6F and S6G). To evaluate potential off-target effects of AZD1208, Pim1−/− BM-MDSCs were transferred into MC38 tumor–bearing mice and treated with either vehicle or AZD1208. Pharmacologic PIM inhibition did not further retard tumor progression compared with vehicle control in Pim1−/− MDSCs (Supplementary Fig. S6H). Collectively, these findings suggest that AZD1208 treatment subverts the immunosuppressive phenotype of myeloid cells in the TME.
When disrupting the suppressive phenotype of myeloid cells, we also observed significant improvement in CD8+ T-cell function. AZD1208 treatment alone induced a 2-fold expansion of IFNγ-producing CD8+ T cells in comparison with the PBS control group, and PD-L1 blockade alone led to a 2.8-fold increase compared with PBS control group (Fig. 6H and I). Combination therapy with AZD1208 + anti–PD-L1 further augmented the frequency of IFNγ+ CD8+ T cells to 14%, a 3.5-fold increase over the control group (Fig. 6I). A similar additive effect was observed with the expression of another effector molecule, granzyme B (Gzmb; Fig. 6J). Motivated by these results, we further evaluated the therapeutic effect of this combination therapy on a syngeneic B16-BL6 melanoma model. Coadministration of AZD1208 and anti–PD-L1 promoted strong antitumor activity and enhanced the survival of tumor-bearing mice when compared with either treatment alone (Fig. 6K; Supplementary Fig. S7A–S7E). Taken together, our data suggest that pharmacologic inhibition of PIM kinase diminishes the function of immunosuppressive myeloid cells in the TME, which in turn releases the antitumor CD8+ T-cell response from inhibition. This treatment further augments the therapeutic efficacy of PD-L1 blockade.
PIM kinase inhibition overcomes ICB resistance
The scRNA-seq analyses shown in Fig. 2 revealed a strong correlation between FAT signature and ICB resistance, which prompted us to ask whether the presence of myeloid cells with FAT phenotypes could be used to predict tumor progression after treatment. We thus employed PPARγ and CD36 as indicators of FAT and identified a subset of MDSCs coexpressing both these markers in the circulating blood of MC38 tumor–bearing mice (Fig. 7A). Notably, the frequencies of CD36highPPARγhigh myeloid cells positively correlated with tumor volumes following various treatments (Fig. 7B). The CART-supervised machine learning method estimated that mice with low levels (<50%) of PPARγhighCD36high circulating MDSCs were more likely to experience tumor shrinkage compared with mice with high levels (>50%) of PPARγhighCD36high MDSCs (Supplementary Fig. S8A).
PIM kinase inhibition overcomes PD-L1 blockade resistance. A, Proportions of CD36highPPARγhigh myeloid cells are shown in representative flow plots. B, Scatter plot depicting the correlation between the frequencies of CD36highPPARγhigh myeloid cells and tumor progression (n = 38, Spearman correlation test, P = 0.04, r = 0.3). C–L, Bilateral MC38 tumor–bearing mice were treated with two doses of anti–PD-L1 and then were stratified into responders (R) and nonresponders (NR) based on frequencies of CD36highPPARγhigh cells. All mice received one more dose of anti–PD-L1, whereas NRs received additional daily treatments of either vehicle control or AZD1208 for 10 days (see Supplementary Fig. S7B). Upon completion of the full therapeutic regimen, left-sided tumors were harvested for flow cytometry analysis. C–G, CD11b+Gr1−F4/80low MDSCs (C and D), CD11b+Gr1−F4/80high macrophages (C and E), and CD36highPPARγhigh MDSCs (F and G) are shown in representative contour plots and quantitative scatter plots. H, Scatter plot shows the levels of Arg1 expression in MDSCs as measured by flow cytometry. gMFI, geometric mean fluorescence intensity; ns, not significant. I and J, Intratumoral CD8+ T-cell IFNγ secretion in response to anti-CD3/CD28 stimulation is shown in representative contour plots (I) and quantitative scatter plot (J). K, Quantification of Gzmb expression in MDSCs as measured by flow cytometry (R = 10, NR = 6, and NR+AZD1208 = 8). L, Individual tumor growth curves of nonresponder mice treated with vehicle control (black dotted line) or AZD1208 (red solid line). Curves represent pooled data from three independent experiments with n = 6 for NR+vehicle and n = 8 for NR+AZD1208. Significance was determined by t test. *, P < 0.05; **, P < 0.01; ***, P < 0.001.
Next, we sought to determine whether PIM kinase inhibition can reverse tumor progression in nonresponders. To test this idea, mice bearing bilateral MC38 tumors were treated with two doses of anti–PD-L1, and the nonresponder group received daily administration of either vehicle control or AZD1208 for 10 days (Supplementary Fig. S8B). Compared with the responder group, we observed a significant accumulation of Gr-1+CD11b+ MDSCs in nonresponders (Fig. 7C and D). This enrichment was subverted by AZD1208 treatment, with mild effects on F4/80+CD11b+ TAMs (Fig. 7E). With AZD1208 treatment, the frequencies of CD36highPPARγhigh myeloid cells from the nonresponder group were significantly reduced to a level similar to the responder group (Fig. 7F and G). More interestingly, this reduced FAT activity correlated with decreased expression of M2 markers (CD206 and Relmα) in TAMs (Supplementary Fig. S8C and S8D) and Arg1 in MDSCs (Fig. 7H). Consistent with these data, the effector function of CD8+ T cells in AZD1208-treated nonresponders was significantly increased, as evidenced by their enhanced production of IFNγ and Gzmb (Fig. 7I–K). Subsequently, four out of seven AZD1208-treated mice originally deemed nonresponders achieved complete tumor regression, whereas none of the vehicle-treated mice originally deemed nonresponders eradicated their tumors (Fig. 7L). In addition, we found that PIM inhibition induced similar functional and metabolic changes in human MDSCs via reduction of PPARγ (Supplementary Fig. S8E–S8G). In conclusion, these findings suggest that PIM kinase inhibition can selectively target immunosuppressive myeloid cells with high FAO metabolism and in turn restore sensitivity to ICB in nonresponders.
Discussion
To dissect the cellular and molecular mechanisms of ICB resistance, we employed scRNA-seq and identified that distinct features of the immune landscape in the TME, especially enhanced infiltration of immunosuppressive myeloid cells, correlated with nonresponsiveness to ICB. In-depth computational analysis further unveiled a previously unappreciated role of PIM1 in regulating FAO metabolism in myeloid cells. Mechanistically, PIM1 regulates lipid metabolism by inducing PPARg and its downstream targets, such as CD36, which are responsible for lipid uptake in tumor-infiltrating myeloid cells. Consistent with these data, PIM1-deficient myeloid cells exhibited an impaired ability to inhibit the proliferation and antitumor function of effector CD8+ T cells. Targeting this pathway in MDSCs by pharmacologic inhibition of PIM kinases significantly increased the antitumor immune response and sensitized hosts to anti–PD-L1 therapy, overcoming ICB resistance.
A large body of clinical observations reveals that patients with elevated levels of MDSCs have impaired T-cell function and respond poorly to ICB (43, 44). Supporting this notion, selective elimination of MDSCs in various preclinical and clinical studies shows beneficial effects when combined with immunotherapy for treating solid tumors (45–47). As the importance of MDSCs in promoting resistance to immunotherapy is becoming increasingly recognized, understanding the fundamental biology of MDSCs has garnered an incredible amount of interest (47). In recent years, it has become well appreciated that metabolic pathways in these myeloid cells are not only required for fulfilling their energetic demands but also critical for their immunosuppressive functions (48, 49). It is increasingly clear that lipid uptake for fueling mitochondrial OXPHOS is enhanced in M2 macrophages (50–52), and inhibition of lipid uptake significantly diminishes FAO metabolism, thereby impeding the immunosuppressive activity of MDSCs and enhancing tumor control in a T-cell–dependent manner (7, 29). Consistent with this, we demonstrated that an increased FAO signature in tumor-infiltrating myeloid cells positively correlates with tumor progression and ICB resistance. More importantly, our scRNA-seq analyses provided some mechanistic insights into how the distinct metabolic features of MDSCs are regulated. To this end, a previously unappreciated role of PIM1 was identified: that of a key inducible kinase that regulates fatty acid uptake and oxidative metabolism in MDSCs to facilitate their survival and function in the TME. Our data have revealed significant deficiencies in fatty acid metabolism and suppressive function of MDSCs following genetic ablation or pharmacologic inhibition of PIM1. However, the impact of PIM1 on other metabolic pathways, such as glutaminolysis, is still unclear. Future investigations in identifying the alternative energy source or sources employed by Pim1−/− MDSCs are warranted to better ameliorate the suppressive activity of these cells.
The role of PIM1 in energy metabolism has been implicated in various types of cells, such as adipocytes, cardiomyocytes, colorectal cancer cells, and T cells (14, 53, 54). Nevertheless, little is known about how PIM1 regulates multiple aspects of cellular metabolism (55). Genetic deletion of Pim1 enabled us to identify possible metabolic pathways downstream of PIM kinases in MDSCs. To this end, our transcriptomic data revealed that expression levels of PPARγ, a key transcription factor that regulates lipid metabolism, are critically regulated by PIM1 activity. Although the detailed molecular mechanisms by which PIM1 regulates PPARγ expression remain unknown, our data have shown that overexpression of PPARγ in Pim1−/− MDSCs can sufficiently rescue their metabolic and functional defects and promote immunosuppression. Mechanistically, the downregulation of PPARγ in the absence of Pim1 is likely a consequence of decreased PIM1-mediated STAT3 S727 phosphorylation, which is integral for the maximal transcription activity (7, 29). In response to ligand stimulation such as IL-6, activated JAK proteins phosphorylate the Tyr705 residue of STAT3 to induce its homodimerization and nuclear translocation (53). Another regulatory residue on STAT3 is Ser727, whose phosphorylation takes place in the C-terminal transactivation domain to maximize transcriptional activity (54, 55). Our results showed that pharmacologic inhibition of PIM kinases specifically decreases STAT3Ser727, but not STAT3Tyr705, phosphorylation. More importantly, this altered STAT3 phosphorylation impaired STAT3 binding and transcription of Pparg. Taken together, our findings support a novel mechanism by which PIM1 enhances phosphorylation of STAT3 at the Ser727 site to maximize Pparg expression, which in turn promotes fatty acid metabolism and immunosuppressive activity in MDSCs.
Metabolic reprogramming of immune cells, such as CD8+ T cells, has emerged as a promising approach to enhance the effectiveness of immunotherapies (56, 57). However, a few strategies have been developed to metabolically target MDSCs. Given that PIM1 is highly expressed in the myeloid lineage (Supplementary Fig. S2E) but only serves as a backup for the AKT–mTOR pathway in T cells (13, 58), PIM kinase inhibition should preferentially disrupt the metabolism and function of MDSCs and have little effect on tumor-reactive T cells. Supporting this idea, we observed improved antitumor effector function in CD8+ T cells when the pan-PIM kinase inhibitor AZD1208 and anti–PD-L1 were coadministered. This is consistent with previous observations (13, 58) and suggests that PIM kinase inhibition–promoted function of CD8+ T cells is likely due to their release from MDSC-mediated suppression. Interestingly, a recent report has shown that PIM kinase inhibition improves the efficacy of immunotherapy by promoting memory phenotypes of T cells (59). Nonetheless, whether these phenotypic changes are intrinsic to tumor reactive T cells remains unclear.
PIM kinases are broadly expressed in both hematopoietic and nonhematopoietic cells (40, 41). Thus, it is possible that PIM kinase inhibition could cause off-target effects. For example, given that PIM kinases are expressed in cardiomyocytes and exert potent cardiac protective function by protecting mitochondrial integrity (55), it is possible that AZD1208 treatment may lead to cardiovascular-related adverse effects. Further work on selectively targeting PIM kinases specifically in MDSCs would be ideal for future clinical applications.
Taken together, our findings demonstrate that the PIM1–PPARγ axis regulates FAO metabolism in suppressive myeloid cells and promotes ICB resistance. This provides new insights into developing targeted approaches to abolish immunosuppression and in turn unleash the T-cell response to control tumors.
Authors' Disclosures
P.-C. Ho reports grants from Cancer Research Institute during the conduct of the study. M.B. Dwinell reports other from Protein Foundry, LLC (cofounder and owner) and Xlock Biosciences, LLC (cofounder and owner) outside the submitted work. W. Cui reports a patent for Mitochondrial Targeted PIM Kinase Inhibitors and Methods of Use pending. No disclosures were reported by the other authors.
Authors' Contributions
G. Xin: Conceptualization, data curation, formal analysis, investigation, methodology, writing–original draft, project administration, writing–review and editing. Yao Chen: Formal analysis. P. Topchyan: Resources, formal analysis. M.Y. Kasmani: Resources, formal analysis. R. Burns: Formal analysis. P.J. Volberding: Validation. X. Wu: Data curation. A. Cohn: Formal analysis. Yiliang Chen: Validation. C.-W. Lin: Formal analysis, validation, methodology. P.-C. Ho: Resources. R. Silverstein: Resources. M.B. Dwinell: Writing–review and editing. W. Cui: Conceptualization, resources, project administration.
Acknowledgments
This work is supported by NIH grants AI125741 (W. Cui), AI148403 (W. Cui), CA246920 (P. Topchyan), and DK127526 (M.Y. Kasmani); American Cancer Society Research Scholar Grant (W. Cui); and a Healthier Wisconsin Grant (W. Cui). G. Xin is supported by The Elizabeth Elser Doolittle Postdoctoral Fellowship and Pelotonia Institute for Immuno-Oncology-OSU Comprehensive Cancer Center (PIIO-OSUCCC) start-up funding (46050 061052). W. Cui is supported by the Medical College of Wisconsin Cancer Center and the Bartlog Endowment Fund. This research was completed in part with computational resources and technical support provided by the Research Computing Center at Medical College of Wisconsin.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Footnotes
Note: Supplementary data for this article are available at Cancer Immunology Research Online (http://cancerimmunolres.aacrjournals.org/).
Cancer Immunol Res 2021;9:454–69
- Received May 25, 2020.
- Revision received November 18, 2020.
- Accepted February 10, 2021.
- Published first February 12, 2021.
- ©2021 American Association for Cancer Research.