Skip to main content

ORIGINAL RESEARCH article

Front. Immunol., 28 June 2023
Sec. Cancer Immunity and Immunotherapy
This article is part of the Research Topic Community Series in Biomarkers in the Era of Cancer Immunotherapy: Zooming in from Periphery to Tumor Microenvironment, Volume II View all 14 articles

Targeting immunosuppressive Ly6C+ classical monocytes reverses anti-PD-1/CTLA-4 immunotherapy resistance

B. Leticia Rodriguez&#x;B. Leticia Rodriguez1†Limo Chen&#x;Limo Chen1†Yanli Li,Yanli Li1,2Shucheng Miao,Shucheng Miao1,3David H. PengDavid H. Peng1Jared J. FradetteJared J. Fradette1Lixia DiaoLixia Diao4Jessica M. KonenJessica M. Konen1Frank R. Rojas AlvarezFrank R. Rojas Alvarez5Luisa M. SolisLuisa M. Solis5Xiaohui Yi,Xiaohui Yi6,7Aparna Padhye,Aparna Padhye1,3Laura A. GibsonLaura A. Gibson1Joshua K. OchiengJoshua K. Ochieng1Xiaofei ZhouXiaofei Zhou7Jing WangJing Wang4Don L. Gibbons,*Don L. Gibbons1,8*
  • 1Department of Thoracic/Head and Neck Medical Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, United States
  • 2Department of Obstetrics and Gynecology, Shanghai General Hospital, Shanghai Jiao Tong University, School of Medicine, Shanghai, China
  • 3United of Texas (UT) Health Graduate School of Biomedical Sciences, The University of Texas MD Anderson Cancer Center, Houston, Texas, United States
  • 4Department Bioinformatics and Computational Biology, The University of Texas MD Anderson Cancer Center, Houston, TX, United States
  • 5Department of Translational Molecular Pathology, The University of Texas MD Anderson Cancer Center, Houston, TX, United States
  • 6Bellicum Pharmaceuticals, Inc., Houston, TX, United States
  • 7Department of Immunology, The University of Texas MD Anderson Cancer Center, Houston, TX, United States
  • 8Department of Molecular and Cellular Oncology, The University of Texas MD Anderson Cancer Center, Houston, TX, United States

Introduction: Despite significant clinical advancement with the use of immune checkpoint blockade (ICB) in non-small cell lung cancer (NSCLC) there are still a major subset of patients that develop adaptive/acquired resistance. Understanding resistance mechanisms to ICB is critical to developing new therapeutic strategies and improving patient survival. The dynamic nature of the tumor microenvironment and the mutational load driving tumor immunogenicity limit the efficacy to ICB. Recent studies indicate that myeloid cells are drivers of ICB resistance. In this study we sought to understand which immune cells were contributing to resistance and if we could modify them in a way to improve response to ICB therapy.

Results: Our results show that combination anti-PD-1/CTLA-4 produces an initial antitumor effect with evidence of an activated immune response. Upon extended treatment with anti-PD-1/CTLA-4 acquired resistance developed with an increase of the immunosuppressive populations, including T-regulatory cells, neutrophils and monocytes. Addition of anti-Ly6C blocking antibody to anti-PD-1/CTLA-4 was capable of completely reversing treatment resistance and restoring CD8 T cell activity in multiple KP lung cancer models and in the autochthonous lung cancer KrasLSL-G12D/p53fl/fl model. We found that there were higher classical Ly6C+ monocytes in anti-PD-1/CTLA-4 combination resistant tumors. B7 blockade illustrated the importance of dendritic cells for treatment efficacy of anti-Ly6C/PD-1/CTLA-4. We further determined that classical Ly6C+ monocytes in anti-PD-1/CTLA-4 resistant tumors are trafficked into the tumor via IFN-γ and the CCL2-CCR2 axis. Mechanistically we found that classical monocytes from ICB resistant tumors were unable to differentiate into antigen presenting cells and instead differentiated into immunosuppressive M2 macrophages or myeloid-derived suppressor cells (MDSC). Classical Ly6C+ monocytes from ICB resistant tumors had a decrease in both Flt3 and PU.1 expression that prevented differentiation into dendritic cells/macrophages.

Conclusions: Therapeutically we found that addition of anti-Ly6C to the combination of anti-PD-1/CTLA-4 was capable of complete tumor eradication. Classical Ly6C+ monocytes differentiate into immunosuppressive cells, while blockade of classical monocytes drives dendritic cell differentiation/maturation to reinvigorate the anti-tumor T cell response. These findings support that immunotherapy resistance is associated with infiltrating monocytes and that controlling the differentiation process of monocytes can enhance the therapeutic potential of ICB.

GRAPHICAL ABSTRACT
www.frontiersin.org

Graphical Abstract

Introduction

Immunotherapy has considerably changed the treatment landscape for patients with lung cancer over the past decade. Multiple checkpoint inhibitors have been FDA approved alone, in combination (anti-PD-(L)1/CTLA-4) or with chemotherapy for NSCLC in the locally-advanced or metastatic setting including nivolumab, pembrolizumab, durvalumab, atezolizumab, ipilimumab, and tremelimumab (17). Checkpoint inhibitors have recently been approved for early-stage disease and are expected to make remarkable improvements in NSCLC patient survival (8, 9). Efforts are now geared towards expanding the use of ICB to benefit the largest number of patients, because despite the remarkable clinical improvements, a major subset of patients do not respond to therapy due to adaptive or acquired resistance (10). The dynamic nature of the tumor microenvironment and the mutational load driving tumor immunogenicity limit the efficacy to ICB (11). Since the majority of NSCLC patients will eventually progress on ICB therapy, the exploration of novel concurrent or sequential combinatorial treatment approaches is essential.

There has been increasing evidence showing that myeloid cells (MC) are a potential driver of ICB resistance (12, 13). The MC population includes granulocytes, monocytes, macrophages and dendritic cells, which are far more abundant than T cells in tumors (14). They are considered critical modulators of the immune tumor microenvironment due to their dynamic differentiation capacity, complex specialized function and polarization (15). Here we describe a novel strategy of modulating the differentiation capacity of a progenitor myeloid cell population, namely monocytes, into dendritic cells and attaining durable therapeutic benefit to ICB resistance in lung cancer models.

We found that adaptive resistance to combination anti-PD-1/CTLA-4 was due in part to immunosuppressive cells including regulatory T cells, neutrophils, and monocytes. Addition of Ly6C blockade to the anti-PD-1/CTLA-4 combination was critical to prevent resistance from occurring and restored T cell activity. Importantly blockade of Ly6C is sufficient to change the differentiation capacity of the monocyte subsets into an anti-tumor dendritic cell population. The classical/nonclassical monocyte subsets identified have diverse differentiation potential. Our study identified that classical Ly6C+ monocytes in ICB resistant tumors have loss of Flt3/PU.1, which prevents their differentiation into dendritic cells. We show for the first time that blocking Ly6C on an intratumoral progenitor monocyte can enhance immunotherapy and ultimately influence the ability of the immune system to durably control cancer cells.

Materials and methods

Reagents and cancer cell lines

Anti-mouse PD-1 (clone RMP1-14), anti-mouse CTLA-4 (clone 9D9), anti-mouse CD80 (clone 16-10A1), anti-mouse CD86 (clone GL-1), anti-mouse Ly6C (clone Monts 1), anti-mouse CCL2 (clone 2H5), and the isotype-matched IgG controls were purchased from BioXCell. Recombinant murine CCL2 (catalog #250-10), IFN-γ (catalog #315-05), IL4 (catalog #214-14), GM-CSF (catalog #315-03), M-CSF (catalog #315-02) were purchased from PeproTech. Murine lung cancer cell lines 393P, 412P, 344P, 344SQ, and 307P were derived from K-rasLA1/+p53R172HΔg/+ mice as previously described in multiple prior lung cancer studies (1623). The LLC-JSP Lewis lung tumor cell line as well as the B16 and B16-OVA melanoma cell lines were maintained in our laboratories (24, 25).

Lung cancer patient samples and mRNA expression profiling

Experimental details regarding TCGA datasets including RNA extraction, mRNA library preparation, sequencing (Illumina HiSeq platform), quality control, data processing and quantification of gene expression are previously published (26). For the PROSPECT and BATTLE-2 samples, the mRNA was extracted from frozen tumor tissue corresponding to the same specimen from which the formalin fixed, paraffin embedded blocks were made. Array based expression profiling of PROSPECT tumors was performed using the Illumina Human WG-6 v3 BeadChip, according to the manufacturer’s protocol. Gene expression data for the PROSPECT dataset have been previously deposited in the GEO repository (GSE42127) (27, 28). The raw data files of transcriptomes were analyzed using Bioconductor R packages. Scatter plots were generated through the cBioPortal (29, 30).

Animal tumor models

Animal studies were reviewed and approved by the University of Texas MD Anderson Cancer Center Institutional Animal Care and Use Committee (IACUC). Six- to 8-week-old 129/Sv, C57BL/6, OT-I TCR transgenic mice were obtained from Charles River Laboratories and maintained in our laboratories. CD11c-DTR mice were from Jackson Laboratories.

For studying the tumor growth and animal survival in syngeneic immunocompetent models, cancer cells with indicated numbers in 100 µl of culture medium were injected subcutaneously into the mouse flank. Tumor growth was measured with digital calipers once a week, and tumor sizes were calculated using the formula: ½ (length × width × width) at indicated time points. At end points, mice were sacrificed to examine tumor weight, metastatic lung nodules as described previously (17). The genetically engineered mouse models (GEMM) (KrasLSL-G12D/p53fl/fl) were generated as previously described (19, 31). Mice were infected at 4 months of age with a 2.5x107 titer pfu adenovirus-cre by intratracheal intubation. Treatment with antibodies were performed at 2-3 months post-infections when tumors were visualized by micro-CT imaging. All mice were immunocompetent and assessed for health daily by the department of veterinary medicine and surgery at The University of Texas MD Anderson Cancer Center.

In vivo treatments

Mice were treated with antibodies (200 µg of anti-PD-1 per mouse; 200 µg of anti-CTLA-4 per mouse; 150 µg of anti-CD25 per mouse; 300 µg of anti-Gr-1 per mouse; 250 µg of anti-Ly6G per mouse; 200 µg of anti-Ly6C per mouse; 100 µg of anti-CCL2 per mouse; or combination) or their Ig G control via IP injection once a week for indicated weeks beginning on the indicated day after tumor cells were subcutaneously implanted. Mice were humanely euthanized at week 4 or week 8 (Figure 1). For the depletion experiments in Figure 2 mice were treated until the IgG CTL group reached maximum allowable size per IACUC protocol. Mice were euthanized at week 9-14 depending on tumor growth from controls. For the anti-PD-1_CTLA4_Ly6C (Figure 3A) studies mice were euthanized at 6 weeks post implantation in order to evaluate immune cell infiltration. GEMM experiment (Figures 3E–I) were sacrificed at 8 weeks post infection when a few of the mice began to show signs of dyspnea/respiratory distress. B16 bone marrow chimeras experiment mice were scarified at week 4. The CD80/86 blockade experiment (Figure 3K) mice were euthanized at week 8.

FIGURE 1
www.frontiersin.org

Figure 1 PD-1 and CTLA-4 blockade produces an initial antitumor effect with development of acquired resistance upon extended treatment (A) anti-PD-1 (200 µg per mouse), anti-CTLA-4 (200 µg per mouse), combination (200 µg of anti-PD-1 plus 200 µg of anti-CTLA-4), or their IgG control mixture was intraperitoneally injected into 129/Sv mice once a week for 3 weeks beginning on day 7 after 344SQ lung cancer cells were subcutaneously implanted (0.1 x 106 cells per mouse). Tumors were measured once a week for 4 weeks. The tumor growth curves are shown in the left panel. Tumor weight at week 4 shown on right panel. (B) Tumors from (A) were harvested at the end point. The indicated immune markers were determined by FACS analysis and summarized data is shown as percentage of CD8, CD8+/Ki-67, Lag3+/Tim3+/CD8+, total CD4 T-cells, and ICOS+/CD4 T-cells. (C) 129/Sv mice were intraperitoneally injected with 400 µg per mouse of anti-CD4 (α -CD4), anti-CD8 (α -CD8), or anti-CD4 plus anti-CD8 (α -CD4/CD8) antibody when 344SQ cancer cells were injected. 200 µg of T cell depletion antibody per mouse was injected into the mice once weekly for 3 weeks to maintain T cell depletion. Mice were treated with combined anti-PD-1 and anti-CTLA-4 (Combination: 200 µg of anti-PD-1 plus 200 µg of anti-CTLA-4 per mouse) beginning on day 7 after a subcutaneous 344SQ cancer cell injection (0.1 x 106 cells per mouse; n = 5) for 3 weeks. The tumor growth was monitored once a week for 4 weeks. The tumor growth curves are shown on the left and the tumor weights are shown on the right. (D) 129/Sv mice were treated weekly with combined anti-PD-1 and anti-CTLA-4 (Combination: 200 µg of anti-PD-1 plus 200 µg of anti-CTLA-4 per mouse) or their IgG control mixture (IgG control) beginning on day 7 after a subcutaneous 344SQ cancer cell injection (0.1 x 106 cells per mouse; n = 5) for 7 weeks. The tumor growth was monitored once a week for 8 weeks. The tumor growth curves are shown. (E) tSNE CD45 plots shown from 344SQ tumors obtained from week 8 from IgG control or combo treated 129/Sv mice from (D). (F) The enrichment of CD45+CD4+CD25+FoxP3+ T-Regulatory cells, CD45+CD3-CD11b+Ly6G+ neutrophils, and CD45+CD3-CD14+Ly6C+ monocyte cells in tumors were analyzed by FACS. p values were calculated with t-test. n.s., no significant difference; *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001.

FIGURE 2
www.frontiersin.org

Figure 2 Anti-PD-1_CTLA-4 treatment resistance is eradicated by anti-Ly6C antibody treatment. (A-H) 129/Sv mice (n = 5) were treated weekly with indicated antibodies beginning on day 7 after a subcutaneous cancer cell injection (0.1 x 106 344SQ cells, 1 x 106 344P cells, 2 x 106 393P cells, or 1 x 106 412P cells per mouse). Dosage of all animal studies: 200 µg of anti-PD-1, 200 µg of anti-CTLA-4, 150 µg of anti-CD25, 300 µg of anti-Gr-1, 250 µg of anti-Ly6G, 200 µg of anti-Ly6C per mouse per injection. p values were calculated with t-test. n.s., no significant difference; *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001.

FIGURE 3
www.frontiersin.org

Figure 3 Anti-PD-1/CTLA-4/Ly6C treatment reduces tumor growth and increases CD8 T cells/dendritic cells. (A) 344SQ tumor bearing 129/Sv mice were weekly treated with anti-PD-1_CTLA-4 (200 µg of anti-PD-1 plus 200 µg of anti-CTLA-4 per mouse), anti-Ly6C alone (200 µg of anti-Ly6C per mouse), anti-PD-1_CTLA-4_Ly6C (200 µg of anti-PD-1 plus 200 µg of anti-CTLA-4 plus 200 µg of anti-Ly6C per mouse), or their IgG control mixture (IgG control) beginning on week 2 after a subcutaneous cancer cell injection (1 x 106 cells per mouse; n = 5) for 4 weeks. Mice received total of 4 treatments starting at week 2 post tumor cell implantation. (B) tSNE CD45 plots from 344SQ tumors treated with IgG control, anti-PD-1_CTLA-4, or anti-PD-1_CTLA-4_Ly6C from (A). Tumors from (A) were harvested to prepare single cell suspensions for FACS analysis. (C) Percentages of CD8 T cells (left), Effector/memory CD8 T cells (middle), and dendritic cells (right). (D) Percentage of Ly6C- (left) and Ly6C+ (right) CD14+CD115+ monocytes. KrasLSL-G12D/p53fl/fl mice generated through intratracheal administration of adenovirus expressing Cre recombinase were treated with either anti-PD-1/CTLA-4 or anti-PD-1/CTLA-4/Ly6C for 8 weeks (E) Micro-CT images shown at week 0 (baseline) and week 8 (endpoint) for anti-PD-1/CTLA-4 or anti-PD-1/CTLA-4/Ly6C. Dashed yellow circles indicate lung tumors. H indicates of heart. (F) Percentage change of tumor area was calculated taking into account prior time point and normalized to the baseline measurement. (G) H&E stained lung sections at week 8 from anti-PD-1/CTLA-4 or anti-PD-1/CTLA-4/Ly6C treated mice from (E). bar= 5 mm. (H) Number of lung tumors on H&E sections from (G) treated with anti-PD-1/CTLA-4 or anti-PD-1/CTLA-4/Ly6C weekly for 8 weeks. (I) IHC of CD8 stained lung tumors from (F) week 8. Percentage of CD8 T cells found in the lung tumors from KrasLSL-G12D/p53fl/fl mice treated for 8 weeks with anti-PD-1/CTLA-4 or anti-PD-1/CTLA-4/Ly6C. (J) The dendritic cells from CD11c-DTR mice were transferred into C57BL/6 mice to generate the chimerical mice, after 8-week stabilization, the mice were treated with diphtheria toxin twice a week to maintain the depletion of dendritic cells. The B16 melanoma-bearing mice were treated weekly with anti-PD-1_CTLA-4_Ly6C for 3 weeks staring on week 1 after tumor cells inoculation. (K) 344SQ tumor bearing 129/Sv mice were treated weekly with blockade of PD-1, CTLA-4, and Ly6C (200 µg of anti-PD-1 plus 200 µg of anti-CTLA-4 plus 200 µg of anti-Ly6C per mouse) or their IgG control mixture (IgG control) beginning on day 7 after a subcutaneous cancer cell injection (0.1 x 106 cells per mouse; n = 5) for 4 weeks. For blocking B7 signal, antibodies (anti-B7: 300 μg of anti-CD80 and 300 μg of anti-CD86 per mouse) were intraperitoneally administered 1 day before the first dose of therapy, and then once a week to maintain the blockade. ANOVA test was used to analyze the data. n.s., not significance; *p < 0.05; **p < 0.01; ***p< 0.001; ****p < 0.0001.

Chemotaxis assay

As modified as previously reported (32), the sorted Ly6Chigh or Ly6Clow monocytes were performed using 24-well transwell plates (Cell Biolabs) with a 5 µm pore size polycarbonate filter. Monocytes were added to the upper chamber of the transwell (1 × 104 cells/transwell) in the presence of CCL2 with indicated concentrations. Transwell plates were incubated at 37°C for 90 minutes. The transmigrated cells were collected from the lower chamber and counted on a flow cytometer. To determine the absolute number of cells in each sample, a standard number of 20 µM size fluorescent microspheres was added to each tube and counted along with the cells. The total number of transmigrated cells equals the number of counted monocytes × the total number of beads/the number of beads counted. The results are expressed as the mean ± SD of the chemotaxis index. The chemotaxis index represents the fold increase in the number of migrated cells in response to CCL2 over the spontaneous cell migration (to control medium).

qRT-PCR assay

Total RNA was isolated from cultured cells or sorted cells from tumor tissues with TRIzol (Invitrogen). cDNA was synthesized by using the SuperScript III kit (Invitrogen), and qPCR was performed with the SYBR Green PCR Master Mix (Applied Biosystems). Primers were designed with National Center for Biotechnology Information (NCBI) primer design software. Relative expression levels were normalized by L32 and calculated by the 2-ΔΔCt method (33). The primers used for amplification are found in the Supplementary Table 1.

Enzyme-linked immunosorbent assay

To measure the concentration of IFN-α, IFN-γ, TNF-α, IL2, IL12 p40, IL23, IL27, and CCL2, tumor lysates or co-culture supernatants were used to perform ELISA assays with the kits for IFN-α (R&D Systems, 42120-1), IFN-γ (eBioscience, 88-7314-88), TNF-α (eBioscience, 88-7324-86), IL-2 (Thermo Fisher Scientific, BMS601), IL-12 p40 (abcam, ab171179), IL23 (R&D Systems, M2300), IL-27 (Thermo Fisher Scientific, BMS6024), and CCL2 (R&D Systems, DY479). The experiments were conducted in triplicates.

Fluorescence-activated cell sorting analysis

Single-cell suspensions were prepared and stained according to standard protocols for flow cytometry including RBC lysis (Biolegend) following manufacturer recommendations. Cells were incubated with anti-CD16/32 to block nonspecific binding and then stained for 1 hr at room-temperature with appropriate dilutions of various combinations of fluorochrome conjugated antibodies. For intracellular staining, cells were fixed and permeabilized with BD Cytofix/Cytoperm (BD Biosciences). Cells were stained with intracellular mixture of fluorochrome conjugated antibodies for 30 minutes. The stained cells were acquired on a BD fortessa, Canto (BD Biosciences), or Cytek Aurora and the data was analyzed using FlowJo software (version 7.6; Tree Star). Doublets were excluded based upon forward and side scatter and dead cells were excluded by using BV510 Ghost dye. All flow cytometry plots shown and data quantitated from flow cytometry were gated upon doublet exclusion (FSC SSC) and viable (BV510 ghost -) CD45+ cells. Gating schemes found in (Figure S10). All antibodies used including dilutions are found in Supplementary Table 1.

Monocyte differentiation assay

The assays were performed as previously described (32, 3437). Briefly, isolated monocytes at 4 × 104/mL were cultured in the 6-well plates with 10% fetal calf serum Iscoves modified Dulbecco medium containing designated cytokines. Cytokines were used at the following concentrations: 10 ng/mL of IL-4, 20 ng/mL of M-CSF, 100 ng/mL of GM-CSF, and 50 ng/mL of IFN-γ. After 6 days, cells were harvested to analyze the populations of myeloid-derived suppressor cells, M1 macrophages, and M2 macrophages by FACS. Dendritic cell maturation was induced by adding TNF-α (10 ng/mL) to the cultures on day 6. Mature dendritic cells were measured by FACS 2 days later.

Bone marrow chimeras and diphtheria toxin treatment

The experimental protocol was modified according to a previous study (38). Briefly, femurs and tibias were taken from CD11c-DTR mice. The bone marrow was flushed out with a syringe and passed through a 70 μm nylon mesh to generate the single cell suspension. After red blood cells were lysed, the cells were washed twice and re-suspended (20 × 106 cells/ml) in HBSS Supplementaryed with HEPES, L-glutamine, penicillin, streptomycin, and gentamycin sulfate. The recipient mice were irradiated (1000 rad) before 4 × 106 bone marrow cells were intravenously transferred. The mice were allowed to rest for 8 weeks before B16 melanoma cell inoculation. Diphtheria toxin was intraperitoneally administered to mice at 4 ng/g body weight. Diphtheria toxin were administered 3-4 days apart to bone marrow chimeras for 6 times.

Antibody-mediated T cell depletion and B7 blockade

For T cell depletion, mice were intraperitoneally pretreated with anti-CD4, anti-CD8, or combination antibodies (400 μg per mouse) one week before tumor cell injection. 200 μg of antibody per mouse was injected into the mice once weekly for indicated weeks beginning on day 1 after a subcutaneous cancer cell injection as previously described (18). For blockade of B7 signal, the protocol was modified according to the previous report (39). 300 μg of anti-CD80 and 300 μg of anti-CD86 per mouse were intraperitoneally administered 1 day before the first dose of combined anti-PD-1_CTLA-4, and then once a week to maintain the blockade.

Statistical analysis

All statistical analyses were performed using GraphPad Prism version 6.00 (GraphPad Software, La Jolla). The un-paired t-test was used for comparisons between two-group means. The ANOVA test was used for comparisons between more than two-group means. All p values are two-tailed and for all analyses, p < 0.05 is considered statistically significant, unless otherwise specified.

Results

PD-1 and CTLA-4 blockade produce an initial antitumor effect with development of acquired resistance upon extended treatment

Due to the distinct signaling pathways of PD-1 and CTLA-4, we wanted to evaluate whether an antibody combination in NSCLC models would improve antitumor immunity as compared to single agent treatment. We evaluated a syngeneic KP tumor model (344SQ) derived from KrasLA1/+; TP53R172HΔG (KP) mice in the 129/Sv background (17, 19), a highly metastatic lung cancer model that has been previously shown to partially respond to PD-(L)1 blockade (17, 18, 40). We observed tumor growth inhibition with anti–PD-1 treatment, but not anti-CTLA-4, and enhanced tumor growth inhibition with combination blockade (Figure 1A). We further observed that both anti-PD-1 and combination blockade strongly increased T cell infiltration into and proliferation within the tumors, whereas CTLA-4 blockade alone did so to a lesser extent (Figure 1B). Combination PD-1/CTLA-4 blockade increased the proliferation and infiltration of CD4+ and CD8+ T-cells more than PD-1 blockade alone (Figure 1B), and significantly decreased the percentage of tumor infiltrating exhausted LAG3+ TIM3+ CD8+ T cells compared to either alone (Figure 1B). Single agent anti-PD-1 or anti-CTLA-4 increased the infiltrating ICOS+ CD4+ T cells to a similar degree, but combination blockade interestingly had a greater effect on the percentage of tumor-infiltrating lymphocytes (TILs) composed of ICOS+ CD4+T cells (Figure 1B). Evaluation of cytokine changes in the tumor microenvironment by ELISA-based analysis of tumor lysates demonstrated that combined blockade elevated the levels of cytokines IFN-γ, TNF-α, IL4, M-CSF, GM-CSF, and CCL2, while it suppressed the immune suppressive cytokine TGF-β (Figure S1A). The effect of anti-PD-1 on IFN-γ and TNF-α production was significantly stronger than that of anti-CTLA-4, while the effect on TGF-β, IL4, M-CSF, GM-CSF, and CCL2 was similar (Figure S1A).

To evaluate the contribution of either CD4+ or CD8+ T cells in the anti–PD-1/CTLA-4 combination, we tested the in vivo effects of CD4+ and/or CD8+ T cell depletion on syngeneic tumor growth. Response to combination anti-PD-1/CTLA-4 diminished when CD4+ T cells, CD8+ T cells, or both T cell subsets were depleted by antibody treatment (Figure 1C). Although dual depletion of CD4+ and CD8+ T cells trended toward a larger tumor size, there was no significant difference from the single CD4+ or CD8+ depletion (Figure 1C). Anti-PD-1/CTLA-4 treatment inhibits tumor progression in a CD4+ and CD8+ T cell-dependent manner, with CD4+ and CD8+ T cells working in an orchestrated fashion to control tumors.

Despite the enhanced response to anti–PD-1/CTLA-4 combination treatment, the effect was transient when observed over an extended treatment period (Figure 1D). Tumors in the treatment group ultimately developed adaptive resistance to the combination therapy and by ~7-8 weeks had primary tumor size similar to the controls. We analyzed harvested tumor tissue at the point of treatment resistance (week 8) using flow cytometry and found that the immunosuppressive populations of T-regulatory cells, neutrophils, and monocytes were highly enriched in the treatment resistant tumors (Figures 1E, F).

Anti-PD-1/CTLA-4 treatment resistance is reversed by addition of anti-Ly6C antibody

In tumors with resistance to the anti-PD-1/CTLA-4 treatment, there was a significant enrichment of CD45+CD4+CD25+FoxP3+ T-Regulatory cells, CD45+CD3-CD11b+Ly6G+ neutrophils, and CD45+CD3-CD14+Ly6C+ monocytes (Figure 1F). These populations have been well documented as immune suppressive in cancer immunotherapy (4143). We evaluated their role in the resistance to combination PD-1/CTLA-4 blockade by adding: anti-CD25 (clone PC61 for depleting CD25+ T-Regulatory cells (44)), anti-Gr1 (clone RB6-8C5 for depleting CD45+CD3-CD11b+Gr1+ myeloid cells (45)), or anti-Ly6G (clone 1A8 for depleting CD45+CD3-CD11b+Ly6G+ neutrophils (46)). In each case the resistance to anti-PD-1/CTLA-4 treatment was substantially inhibited (Figures 2A-D), demonstrating a role for multiple immune cell populations in the antitumor effect, but in each case the tumors demonstrated resistance and the ability to grow out over time.

Given that Ly6C was highly up-regulated on infiltrating monocytes in tumors with resistance to the anti-PD-1/CTLA-4 treatment, we tested whether Ly6C blockade could reverse the resistance. Strikingly, when we used the anti-Ly6C antibody (clone Monts 1) to block Ly6C, the resistance to anti-PD-1/CTLA-4 was completely suppressed in the aggressive KP 344SQ tumor model (Figures 2E, S2A). We tested the concurrent anti-PD-1/CTLA-4/Ly6C therapy in multiple other KP NSCLC animal models and consistently observed that the resistance to anti-PD-1/CTLA-4 treatment was reversed by anti-Ly6C antibody treatment in 393P, 344P, and 412P KP syngeneic models (Figures 2F-H, S2A-D).

To further understand the efficacy of this combination, we tested whether dual blockade of Ly6C/PD-1 or Ly6C/CTLA-4 could produce the same complete tumor regression as triple blockade. Although anti-Ly6C/CTLA-4 treatment achieved a significantly better outcome in terms of tumor growth and immune infiltration than that of anti-Ly6C/PD-1, dual blockade of Ly6C/CTLA-4 did not produce a curative effect (Figures S3A, B). Importantly, the anti-Ly6C antibody Monts 1 in combination with clodronate liposomes (a commonly used macrophage depleting agent), was reportedly used for in vivo macrophage depletion (47) (https://bxcell.com/product/ly6c/). However, although we observed that Monts 1 clearly binds to Ly6C (Figure S4A), the antibody treatment alone did not demonstrate depleting effects (Figure S4B).

Taken together, multiple factors contribute to the immunosuppressive tumor microenvironment and resistance to anti-PD-1/CTLA-4 treatment. However, CD14+ Ly6C+ monocytes play a dominant role in the acquired resistance to combination ICB therapy.

Tumor-infiltrating monocyte-derived dendritic cells play a key role in maximizing the antitumor immunity of anti-PD-1/CTLA-4 treatment

Triple blockade of Ly6C, PD-1, and CTLA-4 resulted in complete tumor responses, thus precluding further analysis of tumor tissues to explore the underlying mechanism of response (Figures 2E-H). To address this, we increased the number of implanted tumor cells and started treatment at two weeks post implantation (Figure 3A). The triple blockade of Ly6C, PD-1, and CTLA-4 caused a significant tumor regression by week 6 (Figure 3A).

Multi-parameter spectral flow cytometry from 344SQ tumors treated with anti-PD-1/CTLA-4 or anti-PD-1/CTLA-4/Ly6C and harvested at week 6 revealed drastic changes to immune cell populations. We observed an increase in classical monocytes (CD45+CD3-CD14+Ly6G-Ly6C+) and neutrophils (CD45+CD3-CD11b+Ly6C-Ly6G+) in anti-PD-1/CTLA-4 treated tumors (Figure 3B). Anti-PD-1/CTLA-4/Ly6C tumors had a significant reduction of monocytes, neutrophils, and an increase in dendritic cells (CD45+CD11b-CD11c+MHCII+F4/80-CD209a+), as compared to anti-PD-1/CTLA-4 treated tumors. Other immune cell populations such as M1/M2 macrophages and MDSC showed no significant difference between treatment groups.

We found that anti-PD-1/CTLA-4/Ly6C treated tumors had increased percentages of CD8 T cells and effector/memory CD8 T cells (Figure 3C) as compared to anti-PD-1/CTLA-4. We found no difference in M1 or M2 macrophages as compared to anti-PD-1/CTLA-4 or IgG CTL (Figure S5A), but interestingly we observed increased dendritic cells in the anti-PD-1/CTLA-4/Ly6C treated tumors (Figure 3C). Intratumoral monocytes from anti-PD-1/CTLA-4/Ly6C treated tumors were composed of nonclassical ~30% Ly6C- and ~10% classical Ly6C+ monocytes, but with the anti-PD-1/CTLA-4 combination we observe a decrease in nonclassical monocytes Ly6C- monocytes and increase in classical Ly6C+ monocytes (Figure 3D). Intriguingly, addition of anti-Ly6C to the combination was capable of restoring the levels of classical and non-classical monocytes to that of the IgG treated controls.

We further tested the efficacy of the triple blockade of Ly6C/PD-1/CTLA-4 in a clinically relevant autochthonous lung cancer GEMM (KrasLSL-G12D/p53fl/fl) where tumors are induced through intratracheal administration of adenoviral-Cre recombinase (23, 48). Primary lung tumors developed ~8 weeks after infection and were confirmed by micro-CT imaging (Figure 3E, left panel). Once tumors were confirmed, mice were randomized and treated with anti-PD-1/CTLA-4 or anti-PD-1/CTLA-4/Ly6C for 8 weeks (Figure 3E, right panel). Micro-CT image analysis revealed that mice treated with anti-PD-1/CTLA-4/Ly6C had a smaller percentage change of lung adenocarcinomas compared to mice treated with anti-PD-1/CTLA-4 by week 8 (Figure 3F). H&E stained lung sections at week 8 revealed fewer lung tumors in the mice treated with anti-PD1-1/CTLA-4/Ly6C as compared to anti-PD-1/CTLA-4 (Figures 3G, H, S5B). To assess immune response in the autochthonous lung tumors from KrasLSL-G12D/p53fl/fl mice treated with anti-PD-1/CTLA-4 or anti-PD-1/CTLA-4/Ly6C we probed for CD8 and quantified the percentage of CD8 cells within each tumor. We found a higher percentage of CD8 T cells in lung tumors from KrasLSL-G12D/p53fl/fl mice treated with anti-PD-1/CTLA-4/Ly6C compared to anti-PD-1/CTLA-4 (Figures 3I, S5C). Overall, these data demonstrate the efficacy of combining anti-PD-1/CTLA-4/Ly6C immunotherapy to control KP lung tumor growth.

Lastly, we sought to test the importance of dendritic cells within anti-PD-1/CTLA-4/Ly6C treated tumors. Using the B16 melanoma tumor model we observed similar response/resistance to combination anti-PD-1/CTLA-4 treatment in the KP lung model, as well as infiltration of Ly6C+ monocytes (Figure S6A). After inoculation of mice with B16-OVA tumor cells, we injected OVA-pulsed ex vivo generated tumor monocyte-derived dendritic cells (TMDC) into the tumor-bearing mice at weeks 1 and 2. Compared with the anti-PD-1/CTLA-4 treatment alone, the anti-PD-1/CTLA-4 treatment plus TMDCs was significantly more potent in suppressing the tumor growth and inducing cytokines from tumor-infiltrating effector T cells and dendritic cells (Figures S6B, C). Second, to confirm the role of dendritic cells in the anti-PD-1/CTLA-4/Ly6C antitumor response, we performed an in vivo dendritic cell depletion assay. The dendritic cells from CD11c-DTR mice were transferred into C57BL/6 mice to generate chimeric mice. After 8-weeks to allow for reconstitution, mice were implanted with B16 tumor cells and treated with combination anti-PD-1/CTLA-4/Ly6C with or without addition of diphtheria toxin (DT) twice a week to deplete the dendritic cells. The results demonstrated that the B16 tumors from mice treated with diphtheria toxin lacked dendritic cells and the efficacy of the anti-PD-1/CTLA-4/Ly6C was lost (Figure 3J).

Dendritic cell-mediated B7 signaling plays a critical role in the antitumor immune response (39), and we further tested the role of antigen expressing cells such as dendritic cells and macrophages in the tumor regression with anti-PD-1/CTLA-4/Ly6C treatment. When B7 signaling (CD80/CD86) was blocked, the treatment effect on KP tumors was reversed and demonstrated no efficacy (Figures 3K, S6D). B7 signaling (CD80/CD86) are expressed on dendritic cells and macrophages. APC-mediated B7 signaling is critical for anti-PD-1/CTLA-4/Ly6C treatment mediated tumor regression. Taken together the data reveal that addition of Ly6C blockade to anti-PD-1/CTLA-4 was capable of suppressing Ly6C+ monocyte levels, restoring T cell responses and delineated a critical role for dendritic cells within the anti-PD-1/CTLA-4/Ly6C treated tumors.

Classical Ly6C+ monocytes in anti-PD-1/CTLA-4 resistant tumors are mediated by IFN-γ and the CCL2-CCR2 axis

Our data reveal that anti-Ly6C addition to the combination treatment can restore immune cell responses. How the monocytes are trafficked and which type of monocytes (classical vs. non-classical) is unknown. First, we evaluated Ly6C expression on CD115+ CD14+ monocytes within the IgG control or anti-PD-1/CTLA-4 treated tumors and found increased Ly6C+ expression and a decrease in Ly6C- on monocytes within combo treated tumors (Figure 4A). Monocytes are heterogenous and consist of a classical/intermediate Ly6C+CD11b+CD115+ population or a nonclassical Ly6C-CD11b+CD115+ population. Published single cell transcriptomic analysis of tumor-infiltrating myeloid cells in human and mouse lung cancers revealed that monocytes are conserved across the species with evidence of a one-to-one equivalence (15), with the murine Ly6C+ classical monocyte being equivalent to the human CD14++ classical monocyte. We used the human gene signatures from Zilionis et al. that are equivalent to the Ly6C+ classical and the Ly6C- non-classical monocytes and interrogated a dataset from melanoma patients treated with nivolumab (49). We found that classical monocytes were higher in patients showing progressive disease/stable disease as compared to patients showing partial response/complete response to immunotherapy and the reverse trend was found with the non-classical monocytes (Figure 4B). A similar trend was found in NSCLC patients treated with nivolumab or pembrolizumab, although not statistically different due to low sample numbers in the responder group (50) (Figure S7A). Overall, these data demonstrate that tumor-infiltrating Ly6C+ classical monocytes are associated with immunotherapy resistance.

FIGURE 4
www.frontiersin.org

Figure 4 Classical Ly6C+ monocytes in anti-PD-1/CTLA-4 resistant tumors is mediated by IFN-γ and the CCL2-CCR2 axis (A) 129/Sv mice were treated weekly with combined anti-PD-1 and anti-CTLA-4 (Combination: 200 µg of anti-PD-1 plus 200 µg of anti-CTLA-4 per mouse) or their IgG control mixture (IgG control) beginning on day 7 after a subcutaneous 344SQ cancer cell injection (0.1 x 106 cells per mouse; n = 5) for 7 weeks. FACS analysis from intertumoral monocytes CD115+CD14+ monocytes and Ly6C+ (left) or Ly6C- (right). (B) Classical monocytes and non-classical monocytes from melanoma nivolumab treated patients (GSE91061). (C) 129/Sv mice were treated weekly with combined anti-PD-1 and anti-CTLA-4 (Combination: 200 µg of anti-PD-1 plus 200 µg of anti-CTLA-4 per mouse) or their IgG control mixture (IgG control) beginning on day 7 after a subcutaneous 344SQ cancer cell injection (0.1 x 106 cells per mouse; n = 5) for 7 weeks. The IFN-γ concentrations in tumors were measured by ELISA assay. (D) The monocytes were sorted (CD45+CD11b+CD115+Ly6G- for sorting) from untreated 344SQ tumors. 10,000 monocytes were plated in the 12-well plate. After 12 hrs of IFN-γ (20 ng/ml) stimulation, the expression of Ly6C on monocytes was analyzed by FACS. The experiments were repeated at least three times. The representative histograms are shown. (E) The CCL2 concentrations in tumors treated with IgG CTL or combo treatment were measured by ELISA assay. (F) Monocytes from IgG CTL or combo treated 344SQ tumors were harvested at week 7 and FACS for analyzing CCR2 expression on total monocytes. (G) The sorted Ly6Chigh or Ly6Clow monocytes from 344SQ tumors were added to the upper chamber of the transwell (1 × 104 cells/transwell) in the presence of CCL2 with indicated concentrations. Transwell plates were incubated at 37 °C for 90 minutes. The transmigrated cells were collected from the lower chamber and analyzed by a flow cytometer. The chemotaxis index represents the fold increase in the number of migrated cells in response to CCL2 over the spontaneous cell migration. (H) Heat map of association between the mRNA levels of human monocyte gene signature and IFNγ pathway in the TCGA lung adenocarcinoma (n = 1008). Spearman correlation with adjusted p value < 0.05 was used as the criteria to select the most significant IFN-γ-related markers for generating the heat map. (I) Spearman’s rank correlation (Rho) was used to assess the association between CCL2/CCR2 and human monocyte marker expression in lung cancer patients’ samples from TCGA (LUAD). The results are shown with t-test. n.s., no significant difference; *p < 0.05; **p < 0.01; ****p < 0.0001.

It was previously reported that Ly6C is regulated by IFN-α, IFN-γ and IL-27 (5153). We observed an increase of IFN-γ in tumors with both short-term and long-term anti-PD-1/CTLA-4 treatment but did not observe a significant change in IFN-α or IL-27 (Figures 4C, S1A). We therefore sought to determine whether anti-PD-1/CTLA-4-induced IFN-γ could regulate Ly6C on monocytes. Monocytes from the KP 344SQ tumors were sorted (CD45+CD11b+CD115+Ly6G-) and incubated with IFN-γ overnight. Ly6C expression on monocytes showed a significant up-regulation after IFN-γ stimulation (Figure 4D). Furthermore, CCL2 in tumors and CCR2 on tumor-infiltrating monocytes were increased with combination anti-PD-1/CTLA-4 treatment (Figures 4E, F, S1A). The CCL2-CCR2 axis plays an important role in recruiting monocytes into tumors (54) and we set up a chemotaxis assay to evaluate the CCL2-CCR2 axis in our models. CCL2-containing media was placed in the bottom chamber and sorted Ly6Clow or Ly6Chigh monocytes were placed in the upper chamber. More Ly6Chigh monocytes migrated in a manner that was CCL2 concentration dependent than Ly6Clow monocytes (Figure 4G). For in vivo testing we used an anti-CCL2 antibody, as previously reported (54, 55). The CCL2 blocking antibody treatment partially reversed the efficacy of the anti-PD-1/CTLA-4 treatment (Figure S7B), while significantly reducing intratumoral monocytes and mature dendritic cells (Figures S7C, D). Further analysis demonstrated that anti-CCL2 treatment dampened the T cell immune response to anti-PD-1/CTLA-4 therapy, with a decrease in total CD4 and CD8 T cell numbers and inhibition of T-cell proliferation (Figure S7E). Our data shows that IFN-γ stimulation leads to upregulation of Ly6C+ expression on CD115+ monocytes and that the Ly6C+ monocytes are recruited via the CCL2-CCR2 axis.

To investigate whether these observations occur in human lung cancers, we analyzed three independent lung cancer patient tumor datasets; The Cancer Genome Atlas (TCGA), PROSPECT and BATTLE-2 trials from MD Anderson (40). We utilized a monocyte gene signature as previously reported (56) and consistently found a strong correlation between the monocyte gene signature and the IFN-γ signature (5759), (p < 0.05) (Figures 4H, S8A). A positive correlation was also found between monocytes (CD14+/CD16a+/CD16b+) and the CCL2-CCR2 axis (Figures 4I, S8B). Taken together, both murine and human data suggest that intertumoral monocyte infiltration is caused by signaling through the CCL2-CCR2 axis, which is potentiated by monocyte Ly6C levels regulated in an IFN-γ dependent manner.

Anti-PD-1/CTLA-4 resistant monocytes have decreased differentiation capacity due to loss of Flt3/PU.1

Monocytes are heterogeneous and function to differentiate into dendritic cells, macrophages, or MDSC based upon the cytokine milieu (54, 6062). In the present study, we found that the tumor-infiltrating monocytes (CD45+CD11b+CD115+Ly6G-) could be classified into two subpopulations, Ly6C+ and Ly6C- monocytes, with combination immunotherapy treated tumors showing higher percentage of Ly6C+ (Figure 5A). Our aim was to determine if Ly6C plays a role to determine monocyte differentiation or if Ly6C simply serves as a marker. To do so, we set up ex vivo monocyte differentiation assays according to previous reports (32, 3437). Monocytes from KP tumors treated with IgG or anti-PD-1/CTLA-4 antibodies until week 7 were sorted using the CD45+CD11b+CD115+Ly6G- markers and then further sorted into Ly6C+ and Ly6C- subsets. The sorted monocytes with different expression levels of Ly6C were cultured ex vivo in 12-well plates containing various cytokine treatments (Figure 5B). M-CSF+IFNγ was added to the sorted monocytes to induce differentiation into M1 macrophages. IL-4 was added to induce differentiation into M2 macrophages. Dendritic cells were incubated with GM-CSF+IL-4 for 6 days then TNF-α was added for maturation. GM-CSF was added to induce differentiation into MDSC. Differentiated cells were then quantified using flow cytometry.

FIGURE 5
www.frontiersin.org

Figure 5 Immunotherapy resistant Ly6C+ monocytes have decreased antigen presenting cell differentiation capacity due to loss of Flt3/PU.1. (A) Tumor-infiltrating monocytes sorted from 344SQ tumors treated with IgG control or combo (anti-PD-1/CTLA-4) are classified into two populations: Ly6C+monocytes and Ly6C-monocytes. (B) The trans-differentiation assay scheme. (C) The differentiation assays were performed under different conditions. The percent of M1 macrophage, M2 macrophage, MDSC, or dendritic cell in total monocytes were analyzed with their markers by FACS at the end points. (D) heterogeneity of intertumoral Ly6C+ monocytes from 344SQ tumors treated with IgG CTL or combo (anti-PD-1/CTLA-4) for 7 weeks. Flt3 and CD11c flow dot plots were used to characterize Ly6C+ monocytes. R1, R2, and R3 gates are found in box gates. (E) Percentage of R1 monocytes CD11c-Flt3-, R2 monocytes CD11c-Flt3+, and R3 monocytes CD11c+Flt3+ from Ly6C+ intertumoral monocytes treated with IgG CTL or combo until week 7. (F) qPCR of Sfpi1 from Ly6C+ FACS sorted monocytes from IgG CTL or combo (anti-PD-1/CTLA-4) treated tumors. (G) qPCR of Flt3 from Ly6C+ FACS sorted monocytes from IgG CTL or combo (anti-PD-1/CTLA-4) treated tumors. (H) Intracellular percentage of PU.1 levels on R1 monocytes from Ly6C+ FACS sorted monocytes from IgG CTL or combo (anti-PD-1/CTLA-4) treated tumors. (I) schematic showing differentiation results from (C) and loss of Flt3/PU.1 in the combo (anti-PD-1/CTLA-4) resistant Ly6C+ monocytes. ANOVA test is shown. n.s., no significant difference; *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001.

We observed that there was a higher percentage of M1 macrophage differentiation upon addition of M-CSF + IFNγ from monocytes sorted from IgG CTL treated tumors, as compared to monocytes treated with PBS (Figure 5Ci). M1 macrophages differentiated from Ly6C- monocytes but not from the Ly6C+ monocytes sorted from combination treated tumors (Figure 5Ci). M2 macrophage differentiation was not different across monocytes from IgG CTL tumors (Figure 5Cii). However, we observed changes in differentiation within the monocytes sorted from the combo treated tumors in which IL-4 had been added, with a higher percentage of M2 macrophages from the Ly6C+ vs the Ly6C- monocytes from the combo treated tumors (Figure 5Cii). Conversely, mature dendritic cells were obtained upon addition of cytokines to both Ly6C+ and Ly6C- monocytes from the IgG CTL treated tumor, but we did not observe dendritic cell differentiation from the Ly6C+ monocyte population from combo treated tumors, even though the Ly6C- monocyte population produced dendritic cell levels equivalent to the IgG CTL (Figure 5Ciii). Finally, MDSC differentiation was performed with addition of GM-CSF. We found low levels of differentiation in the Ly6C- monocyte population from either IgG CTL or combo treated tumors (Figure 5Civ). Ly6C+ monocytes were the predominant contributor of MDSCs and the Ly6C+ monocytes from combo treated tumors had higher levels than the IgG CTL treated tumors (Figure 5Civ). These data suggest that Ly6C+ and Ly6C- monocytes from IgG CTL treated tumors differentiate into M1 macrophages and dendritic cells (Figure 5C). Ly6C- monocytes from combo treated tumors also differentiate into M1 macrophage and dendritic cells (Figure 5C), but the predominate Ly6C+ monocyte population found in combo treated tumors was unable to differentiate into M1 macrophages or dendritic cells and instead differentiated into M2 macrophages and MDSC (Figure 5C). The differences in differentiation potential are largely due to the Ly6C+ classical monocyte population found in the combination immunotherapy treated tumors.

We evaluated Ly6C+ monocyte heterogeneity to understand how polarization of classical Ly6C+ monocytes into macrophages or dendritic cells may be regulated (63). SIRPα+/Flt3+ staining on Ly6C+ monocytes can be subclassified into R1, R2, and R3 subpopulations that can differentiate into macrophages, dendritic cells, and pre-dendritic cells, respectively (63). We utilized a staining/gating scheme that included Flt3 and SIRPα to establish that the 344SQ tumor infiltrating Ly6C+ monocytes are heterogenous (Figures 5D, S9A). We found that IgG treated tumor-derived Ly6C+ monocytes had similar levels of R1, R2, and R3 monocytes as compared to matched spleens (Figure 5D). By contrast, the percentage of R1 tumor-derived monocytes was higher in combo treated as compared to IgG treated controls or matched spleens (Figures 5D, E), while the infiltrating monocytes from anti-PD-1/CTLA-4 treated tumors had decreased levels of R2 and R3 monocytes (Figures 5D, E). R2 and R3 monocytes give rise to monocyte-derived dendritic cells, whereas R1 monocytes generate macrophages. This data is consistent with the differentiation assay results (Figure 5C) demonstrating that Ly6C+ monocytes from the combo group have decreased levels of dendritic cells as compared to Ly6C+ monocytes from the IgG CTL group.

We next sought to understand whether the transcription factor Sfpi1 (PU.1) was regulating the formation of R2 monocytes within the Ly6C+ subset. PU.1 is known to promote Flt3 expression (64). We evaluated levels of Spfi1 and Flt3 in Ly6C+ sorted monocytes and found lower expression of both in the monocytes from combination treated tumors as compared to the IgG CTL (Figures 5F, G). There was also a decrease in the percentage of intracellular PU.1 in Ly6C+ monocytes from the combination treated resistant tumors (Figure 5H). It is known that Ly6C+ monocytes give rise to Ly6C- monocytes via Nr4a1 (Nur77) (65), with increased Ly6C+ monocytes in Nr4a1-/- mice. To assess for a role of Nr4a1 in tumor infiltrating monocytes we performed flow cytometry using intracellular Nur77 and qPCR from Ly6C sorted monocytes from combo resistant and IgG CTL treated tumors. We observed higher levels of Nur77 expression on the Ly6C- monocytes as compared to Ly6C+ (Figure S9B), while the Ly6C+ monocytes from combo treated tumors had similar Nur77 levels compared to the IgG CTL. Similarly, levels of Nr4a1 were elevated in the Ly6C- sorted monocytes from treated tumors and we found higher levels of Nr4a1 on combo treated Ly6C+ monocytes compared to IgG CTL tumors (Figure S9B). These data suggest that the increased infiltration of Ly6C+ monocytes in combo treated tumors was not due to a decrease in Nr4a1 (Nur77) expression. Taken together, these data suggest that in the combination therapy resistant tumors classical Ly6C+ monocytes have decreased Flt3/PU.1 expression (Figure 5I), which leads to their decreased ability to differentiate into antigen presenting cells and instead differentiate into immunosuppressive cells.

The data from multiple tumor models under treatment conditions that produced complete tumor eradication highlight a model in which the importance of tumor infiltrating monocyte-derived dendritic cells is critical to maximizing the therapeutic outcome of combined PD-1 and CTLA-4 blockade in lung cancer. Ly6C+ monocytes differentiate into immunosuppressive cell types, while their blockade drives dendritic cell differentiation/maturation to boost an anti-tumor T cell response.

Discussion

Combination treatment strategies with PD-(L)1 and CTLA-4 ICB antibodies have shown a higher clinical efficacy than either agent alone (5, 6669) and are currently approved in melanoma and NSCLC. CTLA-4 blockade has an impact on the immune priming phase by inducing the activation and proliferation of effector T cells and by reducing T-reg mediated suppression of T cell responses (69, 70). PD-(L)1 blockade functions during the effector phase to restore immune function in T cells that have high levels of antigen exposure (39, 69, 71). Due to their nonredundant effects on immune cell biology, PD-1 and CTLA-4 blockade have the potential for synergistic effects in the treatment of NSCLC. Despite the therapeutic benefit observed with combination immunotherapies, the majority of patients do not benefit from treatment or initial responders may relapse due to acquired resistance (10, 67, 68, 72).

Here we show that blockade of PD-1 and CTLA-4 resulted in robust tumor suppression with functional augmentation of tumor-infiltrating CD4+/CD8+ T cells and an antitumor cytokine milieu. But the therapeutic effects of combination ICB were observed during the early phase of tumor growth, with resistance to therapy developing over time. FACS analysis revealed an enrichment of immunosuppressive intratumoral regulatory T cells, neutrophils and Ly6C+ monocytes. In addition, the initial functionally active T cell compartment reverted to an exhausted, non-proliferative state.

To explore the optimal strategy for treating the resistant tumors, we evaluated antibody-based depletion of each of the immunosuppressive immune cell populations. We found that addition of a depleting antibody for CD25+ T cells, GR1 expressing myeloid cells, and Ly6G+ cells significantly inhibited the resistance to anti-PD-1/CTLA-4, but did not produce complete tumor regressions. Strikingly, addition of anti-Ly6C to the anti-PD-1/CTLA-4 combination herein proved to not only prevent resistance from occurring but also produced complete tumor rejection. Further analysis of tumors treated with anti-Ly6C/PD-1/CTLA-4 revealed an accumulation of mature dendritic cells. Dendritic cells are antigen-presenting cells with potent T cell stimulatory capacity that direct the activation and differentiation of T cells by providing co-stimulatory signals (73). Dendritic cells are critical regulators of both innate and adaptive immune responses (38, 74). The combination of anti-PD-1/CTLA-4 was shown to enhance T cell function, however without a sufficiently large number of activated dendritic cells to cross-present and activate the T cells, a sub-optimal anti-tumor effect was observed.

We observed a highly enriched classical monocyte population in the tumors resistant to anti-PD-1/CTLA-4 treatment. Similar observations were made in NSCLC and melanoma datasets with patients treated with ICB in that non-responders/patients with progressive disease have higher level of classical monocytes as compared to responders or patients with response. The monocytes were found to be recruited to the tumor via the CCL2-CCR2 axis in an IFN-γ dependent manner. High infiltration of monocytes with immunosuppressive features have been shown to lead to T cell dysfunction and failure to respond to immunotherapy (7578). Interestingly, the depletion of the suppressive monocytes with anti-CCL2 did not enhance antitumor response to anti-PD-1/CTLA-4 therapy, but instead blunted the effect.

Addition of anti-Ly6C to the anti-PD-1/CTLA-4 combination prevented resistance from forming and importantly facilitated complete tumor elimination. Our data suggests that anti-Ly6C in the anti-PD-1/CTLA4 combination has an impact predominantly on classical Ly6C+ monocytes. The heterogenous Ly6C monocyte populations have distinct differentiation potential. Classical Ly6C+ monocytes in the immunotherapy resistant tumors have lost their differentiation potential into M1 and dendritic cells, via Flt3/PU.1. The current study delineates the importance of classical monocytes within immunotherapy resistant tumors and potentially how their differentiation can be directed into tumor antigen recognizing dendritic cells. Tumor-infiltrating monocyte-derived dendritic cells have potent antigen-presenting capacity, thereby enhancing the antitumor activity of anti-PD-1/CTLA-4 therapy.

Murine monocytes are heterogenous and consist of a classical/intermediate population Ly6C+CD11b+CD115+ or a nonclassical Ly6C-CD11b+CD115+. Intriguingly, a human homolog for Ly6C has not been reported and the ligands for Ly6C are unknown. Human monocytes are also heterogenous and are made up of classical CD14++CD16-, intermediate CD14+CD16+, or nonclassical CD14+CD16++ monocytes (79, 80). Single cell transcriptomics analysis of tumor-infiltrating myeloid cells in humans and mouse lung cancers revealed that monocytes are conserved across the species with evidence of a one-to-one equivalence (15). Nonetheless the identification of a Ly6C human homolog or orthologue, or a better understanding of the monocyte population(s) defined by Ly6C expression for which there are human counterparts, is warranted to better define how these cell types can be directed as possible treatment options to mitigate ICB therapy resistance.

In addition, a better understanding of the blocking antibody anti-Ly6C Monts1 and how it binds to Ly6C on the surface of classical monocytes and changes differentiation potential will be necessary. Blocking antibodies are known to bind to ligand or receptor on the surface of cells and block the target signaling pathway. For example, PD-1 pathway blockade restores activity of antitumor T cells that have become quiescent by preventing the binding of PD-1 to PD-L1. Blockade with anti-Ly6C could be potentially preventing binding of a ligand to Ly6C on the cell surface of monocytes. It is also possible that the blocking Ly6C antibody may also block Ly6C on other immune cell types including neutrophils, macrophages, and MDSC. Better understanding of the Ly6C interactome will help uncover interactions that are critical for monocyte development and are currently under investigation.

Myeloid cell subsets can circumvent T cell directed ICB therapies by reducing antigen presentation and attenuation of T cell activation, thereby leading to impaired antitumor immune response (1214, 81). ICB aimed at restoring T cell function only partially addresses the immune dysfunction in the tumor microenvironment. It is also essential to modulate the myeloid populations in order to produce a durable anti-tumor immune response. Our study demonstrates that monocytes, a myeloid cell subset, are recruited into inflamed tumor tissues and proliferate after anti-PD-1/CTLA-4 treatment, where they can mediate suppressor roles or rapidly differentiate to dendritic cells upon Ly6C blockade. These TMDCs can drive T-cell activation and proliferation, thereby enhancing the antitumor immunity. This study not only delineates a previously unrecognized mechanism of monocyte-to-dendritic cell differentiation in the immunotherapy resistant tumor microenvironment but also provides a potential strategy to optimize immunotherapy for cancer.

Data availability statement

The authors confirm that the data supporting the findings in the results are available within the article and the supplementary material. Raw data from this study are available from the corresponding author upon request.

Ethics statement

The animal studies were reviewed and approved by The University of Texas MD Anderson Cancer Center Animal Care and Use Committee.

Author contributions

BR, LC, and DG conceived and designed the study. LC, YL, BR, DP, SM, and XY performed and interpreted the experiments. LD developed and performed bioinformatics and statistical analysis. JF, LG, JK, JO, AP, FA, LS, and XZ helped with experiments. LC, JW, and DG supervised the participants. BR, LC, and DG wrote/edited the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported in part by Elsa U. Pardee Foundation research grant, LCD-615867 Lung Cancer Discovery award from American Lung Association, and Mike Hogg Fund research grant to LC; 1F32CA239292 to JK; CPRIT RP150405, CPRIT RP200235, CPRIT-MIRA RP160652-P3, NCI-R37CA214609, Uniting Against Lung Cancer/Lung Cancer Research Foundation award, Rexanna’s Foundation for Fighting Lung Cancer, the Connie Rasor Endowment for Cancer Research, NCI P50CA070907, and grant research support from Janssen Research & Development to DG. The work was also supported by the generous philanthropic contributions to The University of Texas MD Anderson Lung Cancer Moon Shots Program and Cancer Center Support Grant (CCSG CA016672).

Acknowledgments

We thank the members of our laboratories for critical discussions. Special thanks to Dr. Karen Clise-Dwyer and Dr. Kathryn Cox at the MDACC Advanced Cytometry & Sorting Facility for help in designing and validating of our flow cytometry panels.

Conflict of interest

DG serves on scientific advisory committees for AstraZeneca, GlaxoSmithKline, Sanofi, Menarini Ricerche, 4D Pharma, Onconova, Eli Lilly and Janssen and has received research support from Janssen, Takeda, Ribon Therapeutics, Astellas, NGM Biopharmaceuticals, Boehringer Ingelheim and AstraZeneca. Author XY was employed by Bellicum Pharmaceuticals.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1161869/full#supplementary-material

Supplementary Figure 1 | (A) The short-term treatment of anti-PD-1_CTLA-4 changes the profile of cytokines/chemokines in tumors. Tumors from Figure 1A, anti-PD-1 (200 µg per mouse), anti-CTLA-4 (200 µg per mouse), combination (200 µg of anti-PD-1 plus 200 µg of anti-CTLA-4), or their IgG control mixture was intraperitoneally injected into Sv129 mice once a week for 3 weeks beginning on day 7 after 344SQ lung cancer cells were subcutaneously implanted (0.1 x 106 cells per mouse). Tumors were measured once a week for 4 weeks, were harvested at the end point. 1/3 of each individual tumor was prepared for tumor lysate. The concentration of IFN-α, IFN-γ, TNF-α, TGF-β, IL4, IL27, M-CSF, GM-CSF, and CCL2 in tumors were measured by ELISA assay. The data were summarized and shown. The ANOVA test was used to analyze the data. ns, no significant difference; *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001.

Supplementary Figure 2 | Kaplan-Meier survival curves from 344SQ tumors from Figure 2E (A) survival curve from 393P tumors from Figure 2F (B) survival curve from 344P tumors from Figure 2G (C) survival curve from 412P tumors from Figure 2H (D). p values were calculated with t-test. n.s., no significant difference; *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001.

Supplementary Figure 3 | Dual treatment of anti-Ly6C_PD-1 or anti-Ly6C_CTLA-4. (A) The dual blockade of Ly6C_PD-1 or Ly6C_CTLA-4 does not result in a curative therapeutic effect. 344SQ tumor bearing 129/Sv mice were weekly treated with anti-Ly6C alone (200 µg of anti-Ly6C per mouse), anti-Ly6C_PD-1 (200 µg of anti-Ly6C plus 200 µg of anti-PD-1 per mouse), anti- Ly6C_CTLA-4 (200 µg of anti-Ly6C plus 200 µg of anti-CTLA-4 per mouse), or their IgG control mixture (IgG control) beginning on day 7 after a subcutaneous cancer cell injection (1 x 106 cells per mouse; n = 5) for 4 weeks. The tumor growth was weekly monitored. The tumor growth curves are shown in the left panel. At week 6, mice were euthanized to measure the primary tumor weight (middle) and metastatic lung nodules (right). At the end point, tumors were analyzed by FACS for immune profiling (B). The summarized data with ANOVA test is shown in. ns, no significant difference; *p < 0.05; **p < 0.01.

Supplementary Figure 4 | Monts-1 binds to Ly6C but does not deplete the Ly6C+ population (A) Ly6C or empty control vector expression in HEK293T cells were evaluated for EGFP expression (top) or Monts-1 Ly6C AF647 and HK1.4 Ly6C APC-Cy7 (bottom). Ly6C or empty control vector expression in HEK293T cells were evaluated for EGFP expression (top) or Monts-1 Ly6C AF647 and AL-21 Ly6C PE (bottom). (B) 129/Sv mice were treated with 200 µg of anti-Ly6C (Monts 1) per mouse or the IgG control via IP injection once a week for 4 weeks. The Ly6C+CD11b+ population in spleens was analyzed by FACS. The staining antibody for Ly6C is HK1.4. The representative plots are shown.

Supplementary Figure 5 | 344SQ tumors from Figure 3A were FACS analyzed for M1 macrophage or M2 macrophage markers. (B) H&E stained lung sections at week 8 from anti-PD-1/CTLA-4 or anti-PD-1/CTLA-4/Ly6C treated mice from Figure 3E. Scale bar= 5 mm. (C) IHC of representative CD8 staining in lung tumor from mice treated in Figure 3E, quantitation shown in Figure 3I.

Supplementary Figure 6 | (A) C57BL/6 mice were treated weekly with combined anti-PD-1 and anti-CTLA-4 (Combination: 200 mg of anti-PD-1 plus 200 mg of anti-CTLA-4 per mouse) or their IgG control mixture (IgG control) beginning on day 7 after a subcutaneous melanoma B16 cell injection (0.1 x 106 cells per mouse; n = 5) for 8 weeks. FACS analysis of percentage of CD14+CD45+ monocytes quantitated from B16 at week 9. (B) B16-OVA (B16-expressing OVA) melanoma cells (2 × 105 per mouse) were subcutaneously injected into 6-8 week old C57BL/6 mice. On day 7, the mice were treated with IgG control or anti-PD-1_CTLA-4 (200 mg of anti-PD-1 plus 200 mg of anti-CTLA-4 per mouse) once a week for three times. In the treatment of dendritic cell transfer plus anti-PD-1_CTLA-4 (DC + anti-PD-1_CTLA-4) group, the TMDCs were pulsed with 10 µg/ml OVA overnight and injected via tail vein (1 × 105 per mouse) on day 7 and 14. BSA-treated DCs were used as the control. The tumor growth curves and tumor weights at the end point are shown. The mature dendritic cells in tumors were analyzed by FACS and are shown on right. (C) The tumor lysates were subjected to ELISA for measuring cytokine concentrations in tumors. (D) 393P tumor bearing 129/Sv mice were treated weekly with triple blockade of PD-1, CTLA-4, and Ly6C (Tritherapy: 200 µg of anti-PD-1 plus 200 µg of anti-CTLA-4 plus 200 µg of anti-Ly6C per mouse) or their IgG control mixture (IgG control) beginning on day 7 after a subcutaneous cancer cell injection (0.1 x 106 cells per mouse; n = 5) for 4 weeks. For blocking B7 signal, antibodies (anti-B7: 300 μg of anti-CD80 and 300 μg of anti-CD86 per mouse) were intraperitoneally administered 1 day before the first dose of therapy, and then once a week to maintain the blockade.

Supplementary Figure 7 | (A) Classical monocytes and non-classical monocytes from NSCLC nivolumab treated patients (GSE126044). (B) 344SQ tumor bearing 129/Sv mice were weekly treated with the anti-CCL2 alone (100 mg of anti-CCL2 per mouse), the combined anti-PD-1 and anti-CTLA-4 (anti-PD-1_CTLA-4: 200 mg of anti-PD-1 plus 200 mg of anti-CTLA-4 per mouse), the combined anti-PD-1, anti-CTLA-4, and anti-CCL2 (anti-PD-1_CTLA-4_CCL2: 200 mg of anti-PD-1 plus 200 mg of anti-CTLA-4 plus 100 mg of anti-CCL2 per mouse), or their IgG control mixture beginning on day 7 after cancer cell subcutaneous injection (0.2 x 106 cells per mouse; n = 5) for 5 weeks. Tumors were measured once a week. The tumor growth curves, and primary tumor masses are shown. At the end point, tumors were harvested to prepare single cell suspensions. The single cell suspensions were subjected to FACS analysis. The tumor-infiltrating monocyte populations (C) and mature dendritic cells (mDCs) (D) and the CD4 and CD8 T cells (E) from different groups are shown. The results are shown with t-test. ns, no significant difference; *p < 0.05; **p < 0.01; ***p < 0.001; ****p < 0.0001. In (F, J, and K), the results are shown with ANOVA test. ns, no significant difference; *p < 0.05; **p < 0.01; ****p < 0.0001.

Supplementary Figure 8 | (A) Heat map of association between the mRNA levels of human monocyte gene signature and IFNγ pathway in the BATTLE-2 and PROSPECT datasets. Spearman correlation with adjusted p value < 0.05 was used as the criteria to select the most significant IFN-γ-related markers for generating the heat map. (B) Spearman’s rank correlation (Rho) was used to assess the association between CCL2/CCR2 and human monocyte marker expression in lung cancer patients’ samples from BATTLE-2 and PROSPECT datasets.

Supplementary Figure 9 | (A) Gating scheme used to characterize R1, R2, and R3 Ly6C+ monocytes shown in Figure 5E. (B) (Left) Percentage of Nur77 on monocytes from 344SQ tumors treated with IgG or anti-PD-1/CTLA-4 combo until week 7. (Right) qPCR of Nr4a from Ly6C+ or Ly6C- FACS sorted monocytes from IgG CTL or combo (anti-PD-1/CTLA-4) treated tumors.

Supplementary Figure 10 | FACS analysis gating scheme.

Supplementary Table 1 | Table of antibodies used for flow cytometry and primers used for qPCR’s.

References

1. Reck M, Rodríguez-Abreu D, Robinson AG, Hui R, Csőszi T, Fülöp A, et al. Pembrolizumab versus chemotherapy for PD-L1–positive non–Small-Cell lung cancer. New Engl J Med (2016) 375(19):1823–33. doi: 10.1056/NEJMoa1606774

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Mok TSK, Wu YL, Kudaba I, Kowalski DM, Cho BC, Turna HZ, et al. Pembrolizumab versus chemotherapy for previously untreated, PD-L1-expressing, locally advanced or metastatic non-small-cell lung cancer (KEYNOTE-042): a randomised, open-label, controlled, phase 3 trial. Lancet (2019) 393(10183):1819–30. doi: 10.1016/S0140-6736(18)32409-7

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Gandhi L, Rodríguez-Abreu D, Gadgeel S, Esteban E, Felip E, De Angelis F, et al. Pembrolizumab plus chemotherapy in metastatic non-Small-Cell lung cancer. N Engl J Med (2018) 378(22):2078–92. doi: 10.1056/NEJMoa1801005

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Carbone DP, Reck M, Paz-Ares L, Creelan B, Horn L, Steins M, et al. First-line nivolumab in stage IV or recurrent non-Small-Cell lung cancer. N Engl J Med (2017) 376(25):2415–26. doi: 10.1056/NEJMoa1613493

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Hellmann MD, Paz-Ares L, Bernabe Caro R, Zurawski B, Kim SW, Carcereny Costa E, et al. Nivolumab plus ipilimumab in advanced non-Small-Cell lung cancer. N Engl J Med (2019) 381(21):2020–31. doi: 10.1056/NEJMoa1910231

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Herbst RS, Giaccone G, de Marinis F, Reinmuth N, Vergnenegre A, Barrios CH, et al. Atezolizumab for first-line treatment of PD-L1-Selected patients with NSCLC. N Engl J Med (2020) 383(14):1328–39. doi: 10.1056/NEJMoa1917346

PubMed Abstract | CrossRef Full Text | Google Scholar

7. West H, McCleod M, Hussein M, Morabito A, Rittmeyer A, Conter HJ, et al. Atezolizumab in combination with carboplatin plus nab-paclitaxel chemotherapy compared with chemotherapy alone as first-line treatment for metastatic non-squamous non-small-cell lung cancer (IMpower130): a multicentre, randomised, open-label, phase 3 trial. Lancet Oncol (2019) 20(7):924–37. doi: 10.1016/S1470-2045(19)30167-6

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Forde PM, Spicer J, Lu S, Provencio M, Mitsudomi T, Awad MM, et al. Neoadjuvant nivolumab plus chemotherapy in resectable lung cancer. N Engl J Med (2022) 386(21):1973–85. doi: 10.1056/NEJMoa2202170

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Garassino MC, Mazieres J, Reck M, Chouaid C, Bischoff H, Reinmuth N, et al. Durvalumab after sequential chemoradiotherapy in stage III, unresectable NSCLC: the phase 2 PACIFIC-6 trial. J Thorac Oncol (2022) 17:1415–27. doi: 10.1016/j.jtho.2022.07.1148

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Sharma P, Hu-Lieskovan S, Wargo JA, Ribas A. Primary, adaptive, and acquired resistance to cancer immunotherapy. Cell (2017) 168(4):707–23. doi: 10.1016/j.cell.2017.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Horvath L, Thienpont B, Zhao L, Wolf D, Pircher A. Overcoming immunotherapy resistance in non-small cell lung cancer (NSCLC) - novel approaches and future outlook. Mol Cancer (2020) 19(1):141. doi: 10.1186/s12943-020-01260-z

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Medina-Echeverz J, Aranda F, Berraondo P. Myeloid-derived cells are key targets of tumor immunotherapy. Oncoimmunology (2014) 3:e28398. doi: 10.4161/onci.28398

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Kim IS, Gao Y, Welte T, Wang H, Liu J, Janghorban M, et al. Immuno-subtyping of breast cancer reveals distinct myeloid cell profiles and immunotherapy resistance mechanisms. Nat Cell Biol (2019) 21(9):1113–26. doi: 10.1038/s41556-019-0373-7

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Takahashi H, Varner J. Rel-ating myeloid cells to cancer therapy. Nat Cancer (2020) 1(5):480–1. doi: 10.1038/s43018-020-0069-8

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Zilionis R, Engblom C, Pfirschke C, Savova V, Zemmour D, Saatcioglu HD, et al. Single-cell transcriptomics of human and mouse lung cancers reveals conserved myeloid populations across individuals and species. Immunity (2019) 50(5):1317–34.e10. doi: 10.1016/j.immuni.2019.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Ahn YH, Gibbons DL, Chakravarti D, Creighton CJ, Rizvi ZH, Adams HP, et al. ZEB1 drives prometastatic actin cytoskeletal remodeling by downregulating miR-34a expression. J Clin Invest (2012) 122(9):3170–83. doi: 10.1172/JCI63608

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Chen L, Gibbons DL, Goswami S, Cortez MA, Ahn YH, Byers LA, et al. Metastasis is regulated via microRNA-200/ZEB1 axis control of tumour cell PD-L1 expression and intratumoral immunosuppression. Nat Commun (2014) 5:5241. doi: 10.1038/ncomms6241

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Chen L, Yi X, Goswami S, Ahn YH, Roybal JD, Yang Y, et al. Growth and metastasis of lung adenocarcinoma is potentiated by BMP4-mediated immunosuppression. Oncoimmunology (2016) 5(11):e1234570. doi: 10.1080/2162402X.2016.1234570

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Gibbons DL, Lin W, Creighton CJ, Rizvi ZH, Gregory PA, Goodall GJ, et al. Contextual extracellular cues promote tumor cell EMT and metastasis by regulating miR-200 family expression. Genes Dev (2009) 23(18):2140–51. doi: 10.1101/gad.1820209

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Peng DH, Ungewiss C, Tong P, Byers LA, Wang J, Canales JR, et al. ZEB1 induces LOXL2-mediated collagen stabilization and deposition in the extracellular matrix to drive lung cancer invasion and metastasis. Oncogene (2017) 36(14):1925–38. doi: 10.1038/onc.2016.358

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Ungewiss C, Rizvi ZH, Roybal JD, Peng DH, Gold KA, Shin DH, et al. The microRNA-200/Zeb1 axis regulates ECM-dependent beta1-integrin/FAK signaling, cancer cell invasion and metastasis through CRKL. Sci Rep (2016) 6:18652. doi: 10.1038/srep18652

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Yang Y, Ahn YH, Chen Y, Tan X, Guo L, Gibbons DL, et al. ZEB1 sensitizes lung adenocarcinoma to metastasis suppression by PI3K antagonism. J Clin Invest (2014) 124(6):2696–708. doi: 10.1172/JCI72171

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Peng DH, Kundu ST, Fradette JJ, Diao L, Tong P, Byers LA, et al. ZEB1 suppression sensitizes KRAS mutant cancers to MEK inhibition by an IL17RD-dependent mechanism. Sci Transl Med (2019) 11(483). doi: 10.1126/scitranslmed.aaq1238

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Palumbo JS, Kombrinck KW, Drew AF, Grimes TS, Kiser JH, Degen JL, et al. Fibrinogen is an important determinant of the metastatic potential of circulating tumor cells. Blood (2000) 96(10):3302–9. doi: 10.1182/blood.V96.10.3302

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Palumbo JS, Potter JM, Kaplan LS, Talmage K, Jackson DG, Degen JL. Spontaneous hematogenous and lymphatic metastasis, but not primary tumor growth or angiogenesis, is diminished in fibrinogen-deficient mice. Cancer Res (2002) 62(23):6966–72.

PubMed Abstract | Google Scholar

26. Cancer Genome Atlas Research N. Comprehensive molecular profiling of lung adenocarcinoma. Nature (2014) 511(7511):543–50. doi: 10.1038/nature13385

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Tang H, Xiao G, Behrens C, Schiller J, Allen J, Chow CW, et al. A 12-gene set predicts survival benefits from adjuvant chemotherapy in non-small cell lung cancer patients. Clin Cancer Res (2013) 19(6):1577–86. doi: 10.1158/1078-0432.CCR-12-2321

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Lou Y, Diao L, Cuentas ER, Denning WL, Chen L, Fan YH, et al. Epithelial-mesenchymal transition is associated with a distinct tumor microenvironment including elevation of inflammatory signals and multiple immune checkpoints in lung adenocarcinoma. Clin Cancer Res (2016) 22(14):3630–42. doi: 10.1158/1078-0432.CCR-15-1434

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal (2013) 6(269):pl1. doi: 10.1126/scisignal.2004088

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discovery (2012) 2(5):401–4. doi: 10.1158/2159-8290.CD-12-0095

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Zheng S, El-Naggar AK, Kim ES, Kurie JM, Lozano G. A genetic mouse model for metastatic lung cancer with gender differences in survival. Oncogene (2007) 26(48):6896–904. doi: 10.1038/sj.onc.1210493

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Partida-Sanchez S, Goodrich S, Kusser K, Oppenheimer N, Randall TD, Lund FE. Regulation of dendritic cell trafficking by the ADP-ribosyl cyclase CD38: impact on the development of humoral immunity. Immunity (2004) 20(3):279–91. doi: 10.1016/S1074-7613(04)00048-2

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res (2001) 29(9):e45. doi: 10.1093/nar/29.9.e45

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Ohishi K, Varnum-Finney B, Serda RE, Anasetti C, Bernstein ID. The notch ligand, delta-1, inhibits the differentiation of monocytes into macrophages but permits their differentiation into dendritic cells. Blood (2001) 98(5):1402–7. doi: 10.1182/blood.V98.5.1402

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Morales JK, Kmieciak M, Knutson KL, Bear HD, Manjili MH. GM-CSF is one of the main breast tumor-derived soluble factors involved in the differentiation of CD11b-Gr1- bone marrow progenitor cells into myeloid-derived suppressor cells. Breast Cancer Res Treat (2010) 123(1):39–49. doi: 10.1007/s10549-009-0622-8

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Gordon S. Alternative activation of macrophages. Nat Rev Immunol (2003) 3(1):23–35. doi: 10.1038/nri978

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Murray PJ, Wynn TA. Protective and pathogenic functions of macrophage subsets. Nat Rev Immunol (2011) 11(11):723–37. doi: 10.1038/nri3073

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Zammit DJ, Cauley LS, Pham QM, Lefrancois L. Dendritic cells maximize the memory CD8 T cell response to infection. Immunity (2005) 22(5):561–70. doi: 10.1016/j.immuni.2005.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Kamphorst AO, Wieland A, Nasti T, Yang S, Zhang R, Barber DL, et al. Rescue of exhausted CD8 T cells by PD-1-targeted therapies is CD28-dependent. Science (2017) 355(6332):1423–7. doi: 10.1126/science.aaf0683

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Chen L, Diao L, Yang Y, Yi X, Rodriguez BL, Li Y, et al. CD38-mediated immunosuppression as a mechanism of tumor cell escape from PD-1/PD-L1 blockade. Cancer Discovery (2018) 8:1156–75. doi: 10.1158/2159-8290.CD-17-1033

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Dodagatta-Marri E, Meyer DS, Reeves MQ, Paniagua R, To MD, Binnewies M, et al. Alpha-PD-1 therapy elevates Treg/Th balance and increases tumor cell pSmad3 that are both targeted by alpha-TGFbeta antibody to promote durable rejection and immunity in squamous cell carcinomas. J Immunother Cancer (2019) 7(1):62. doi: 10.1186/s40425-018-0493-9

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Liu Y, Wei G, Cheng WA, Dong Z, Sun H, Lee VY, et al. Targeting myeloid-derived suppressor cells for cancer immunotherapy. Cancer Immunol Immunother (2018) 67(8):1181–95. doi: 10.1007/s00262-018-2175-3

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Giese MA, Hind LE, Huttenlocher A. Neutrophil plasticity in the tumor microenvironment. Blood (2019) 133(20):2159–67. doi: 10.1182/blood-2018-11-844548

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Setiady YY, Coccia JA, Park PU. In vivo depletion of CD4+FOXP3+ treg cells by the PC61 anti-CD25 monoclonal antibody is mediated by FcgammaRIII+ phagocytes. Eur J Immunol (2010) 40(3):780–6. doi: 10.1002/eji.200939613

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Li Z, Pang Y, Gara SK, Achyut BR, Heger C, Goldsmith PK, et al. Gr-1+CD11b+ cells are responsible for tumor promoting effect of TGF-beta in breast cancer progression. Int J Cancer (2012) 131(11):2584–95. doi: 10.1002/ijc.27572

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Daley JM, Thomay AA, Connolly MD, Reichner JS, Albina JE. Use of Ly6G-specific monoclonal antibody to deplete neutrophils in mice. J Leukoc Biol (2008) 83(1):64–70. doi: 10.1189/jlb.0407247

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Rowe AM, Yun H, Treat BR, Kinchington PR, Hendricks RL. Subclinical herpes simplex virus type 1 infections provide site-specific resistance to an unrelated pathogen. J Immunol (2017) 198(4):1706–17. doi: 10.4049/jimmunol.1601310

PubMed Abstract | CrossRef Full Text | Google Scholar

48. DuPage M, Dooley AL, Jacks T. Conditional mouse lung cancer models using adenoviral or lentiviral delivery of cre recombinase. Nat Protoc (2009) 4(7):1064–72. doi: 10.1038/nprot.2009.95

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Riaz N, Havel JJ, Makarov V, Desrichard A, Urba WJ, Sims JS, et al. Tumor and microenvironment evolution during immunotherapy with nivolumab. Cell (2017) 171(4):934–49.e16. doi: 10.1016/j.cell.2017.09.028

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Cho J-W, Hong MH, Ha S-J, Kim Y-J, Cho BC, Lee I, et al. Genome-wide identification of differentially methylated promoters and enhancers associated with response to anti-PD-1 therapy in non-small cell lung cancer. Exp Mol Med (2020) 52(9):1550–63. doi: 10.1038/s12276-020-00493-8

PubMed Abstract | CrossRef Full Text | Google Scholar

51. DeLong JH, Hall AO, Konradt C, Coppock GM, Park J, Harms Pritchard G, et al. Cytokine- and TCR-mediated regulation of T cell expression of Ly6C and sca-1. J Immunol (2018) 200(5):1761–70. doi: 10.4049/jimmunol.1701154

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Seo SU, Kwon HJ, Ko HJ, Byun YH, Seong BL, Uematsu S, et al. Type I interferon signaling regulates Ly6C(hi) monocytes and neutrophils during acute viral pneumonia in mice. PloS Pathog (2011) 7(2):e1001304. doi: 10.1371/journal.ppat.1001304

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Jutila MA, Kroese FG, Jutila KL, Stall AM, Fiering S, Herzenberg LA, et al. Ly-6C is a monocyte/macrophage and endothelial cell differentiation antigen regulated by interferon-gamma. Eur J Immunol (1988) 18(11):1819–26. doi: 10.1002/eji.1830181125

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Qian BZ, Li J, Zhang H, Kitamura T, Zhang J, Campion LR, et al. CCL2 recruits inflammatory monocytes to facilitate breast-tumour metastasis. Nature (2011) 475(7355):222–5. doi: 10.1038/nature10138

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Getts DR, Terry RL, Getts MT, Muller M, Rana S, Shrestha B, et al. Ly6c+ "inflammatory monocytes" are microglial precursors recruited in a pathogenic manner in West Nile virus encephalitis. J Exp Med (2008) 205(10):2319–37. doi: 10.1084/jem.20080421

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Nirmal AJ, Regan T, Shih BB, Hume DA, Sims AH, Freeman TC. Immune cell gene signatures for profiling the microenvironment of solid tumors. Cancer Immunol Res (2018) 6(11):1388–400. doi: 10.1158/2326-6066.CIR-18-0342

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Gao J, Shi LZ, Zhao H, Chen J, Xiong L, He Q, et al. Loss of IFN-gamma pathway genes in tumor cells as a mechanism of resistance to anti-CTLA-4 therapy. Cell (2016) 167(2):397–404 e9. doi: 10.1016/j.cell.2016.08.069

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Roh W, Chen PL, Reuben A, Spencer CN, Prieto PA, Miller JP, et al. Integrated molecular analysis of tumor biopsies on sequential CTLA-4 and PD-1 blockade reveals markers of response and resistance. Sci Transl Med (2017) 9(379). doi: 10.1126/scitranslmed.aah3560

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Gao J, Ward JF, Pettaway CA, Shi LZ, Subudhi SK, Vence LM, et al. VISTA is an inhibitory immune checkpoint that is increased after ipilimumab therapy in patients with prostate cancer. Nat Med (2017) 23(5):551–5. doi: 10.1038/nm.4308

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Geissmann F, Auffray C, Palframan R, Wirrig C, Ciocca A, Campisi L, et al. Blood monocytes: distinct subsets, how they relate to dendritic cells, and their possible roles in the regulation of T-cell responses. Immunol Cell Biol (2008) 86(5):398–408. doi: 10.1038/icb.2008.19

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Geissmann F, Manz MG, Jung S, Sieweke MH, Merad M, Ley K. Development of monocytes, macrophages, and dendritic cells. Science (2010) 327(5966):656–61. doi: 10.1126/science.1178331

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Delneste Y, Charbonnier P, Herbault N, Magistrelli G, Caron G, Bonnefoy JY, et al. Interferon-gamma switches monocyte differentiation from dendritic cells to macrophages. Blood (2003) 101(1):143–50. doi: 10.1182/blood-2002-04-1164

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Menezes S, Melandri D, Anselmi G, Perchet T, Loschko J, Dubrot J, et al. The heterogeneity of Ly6C(hi) monocytes controls their differentiation into iNOS(+) macrophages or monocyte-derived dendritic cells. Immunity (2016) 45(6):1205–18. doi: 10.1016/j.immuni.2016.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Carotta S, Dakic A, D'Amico A, Pang SH, Greig KT, Nutt SL, et al. The transcription factor PU.1 controls dendritic cell development and Flt3 cytokine receptor expression in a dose-dependent manner. Immunity (2010) 32(5):628–41. doi: 10.1016/j.immuni.2010.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Hanna RN, Carlin LM, Hubbeling HG, Nackiewicz D, Green AM, Punt JA, et al. The transcription factor NR4A1 (Nur77) controls bone marrow differentiation and the survival of Ly6C- monocytes. Nat Immunol (2011) 12(8):778–85. doi: 10.1038/ni.2063

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Curran MA, Montalvo W, Yagita H, Allison JP. PD-1 and CTLA-4 combination blockade expands infiltrating T cells and reduces regulatory T and myeloid cells within B16 melanoma tumors. Proc Natl Acad Sci United States America (2010) 107(9):4275–80. doi: 10.1073/pnas.0915174107

CrossRef Full Text | Google Scholar

67. Chae YK, Arya A, Iams W, Cruz MR, Chandra S, Choi J, et al. Current landscape and future of dual anti-CTLA4 and PD-1/PD-L1 blockade immunotherapy in cancer; lessons learned from clinical trials with melanoma and non-small cell lung cancer (NSCLC). J Immunother Cancer (2018) 6(1):39. doi: 10.1186/s40425-018-0349-3

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Hodi FS, Chiarion-Sileni V, Gonzalez R, Grob JJ, Rutkowski P, Cowey CL, et al. Nivolumab plus ipilimumab or nivolumab alone versus ipilimumab alone in advanced melanoma (CheckMate 067): 4-year outcomes of a multicentre, randomised, phase 3 trial. Lancet Oncol (2018) 19(11):1480–92. doi: 10.1016/S1470-2045(18)30700-9

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Wei SC, Levine JH, Cogdill AP, Zhao Y, Anang NAS, Andrews MC, et al. Distinct cellular mechanisms underlie anti-CTLA-4 and anti-PD-1 checkpoint blockade. Cell (2017) 170(6):1120–33 e17. doi: 10.1016/j.cell.2017.07.024

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Wolchok JD, Saenger Y. The mechanism of anti-CTLA-4 activity and the negative regulation of T-cell activation. oncologist (2008) 13 Suppl 4:2–9. doi: 10.1634/theoncologist.13-S4-2

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Topalian SL, Hodi FS, Brahmer JR, Gettinger SN, Smith DC, McDermott DF, et al. Safety, activity, and immune correlates of anti-PD-1 antibody in cancer. N Engl J Med (2012) 366(26):2443–54. doi: 10.1056/NEJMoa1200690

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Ascierto PA, Butterfield LH, Demaria S, Ferris RL, Freeman GJ, Lo RS, et al. The great debate at "Immunotherapy bridge 2018", Naples, November 29th, 2018. J Immunother Cancer (2019) 7(1):221. doi: 10.1186/s40425-019-0683-0

PubMed Abstract | CrossRef Full Text | Google Scholar

73. Sallusto F, Lanzavecchia A. The instructive role of dendritic cells on T-cell responses. Arthritis Res (2002) 4 Suppl 3:S127–32. doi: 10.1186/ar567

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Liu YJ. Dendritic cell subsets and lineages, and their functions in innate and adaptive immunity. Cell (2001) 106(3):259–62. doi: 10.1016/S0092-8674(01)00456-1

PubMed Abstract | CrossRef Full Text | Google Scholar

75. Gebhardt C, Sevko A, Jiang H, Lichtenberger R, Reith M, Tarnanidis K, et al. Myeloid cells and related chronic inflammatory factors as novel predictive markers in melanoma treatment with ipilimumab. Clin Cancer Res (2015) 21(24):5453–9. doi: 10.1158/1078-0432.CCR-15-0676

PubMed Abstract | CrossRef Full Text | Google Scholar

76. Sade-Feldman M, Kanterman J, Klieger Y, Ish-Shalom E, Olga M, Saragovi A, et al. Clinical significance of circulating CD33+CD11b+HLA-DR- myeloid cells in patients with stage IV melanoma treated with ipilimumab. Clin Cancer Res (2016) 22(23):5661–72. doi: 10.1158/1078-0432.CCR-15-3104

PubMed Abstract | CrossRef Full Text | Google Scholar

77. Krieg C, Nowicka M, Guglietta S, Schindler S, Hartmann FJ, Weber LM, et al. High-dimensional single-cell analysis predicts response to anti-PD-1 immunotherapy. Nat Med (2018) 24(2):144–53. doi: 10.1038/nm.4466

PubMed Abstract | CrossRef Full Text | Google Scholar

78. Singhal S, Stadanlick J, Annunziata MJ, Rao AS, Bhojnagarwala PS, O'Brien S, et al. Human tumor-associated monocytes/macrophages and their regulation of T cell responses in early-stage lung cancer. Sci Transl Med (2019) 11(479). doi: 10.1126/scitranslmed.aat1500

PubMed Abstract | CrossRef Full Text | Google Scholar

79. Sprangers S, de Vries TJ, Everts V. Monocyte heterogeneity: consequences for monocyte-derived immune cells. J Immunol Res (2016) 2016:1475435. doi: 10.1155/2016/1475435

PubMed Abstract | CrossRef Full Text | Google Scholar

80. Anbazhagan K, Duroux-Richard I, Jorgensen C, Apparailly F. Transcriptomic network support distinct roles of classical and non-classical monocytes in human. Int Rev Immunol (2014) 33(6):470–89. doi: 10.3109/08830185.2014.902453

PubMed Abstract | CrossRef Full Text | Google Scholar

81. Youn JI, Nagaraj S, Collazo M, Gabrilovich DI. Subsets of myeloid-derived suppressor cells in tumor-bearing mice. J Immunol (2008) 181(8):5791–802. doi: 10.4049/jimmunol.181.8.5791

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: monocytes, immunotherapy, resistance, myeloid cells, PD-1/CTLA-4 immunotherapy

Citation: Rodriguez BL, Chen L, Li Y, Miao S, Peng DH, Fradette JJ, Diao L, Konen JM, Alvarez FRR, Solis LM, Yi X, Padhye A, Gibson LA, Ochieng JK, Zhou X, Wang J and Gibbons DL (2023) Targeting immunosuppressive Ly6C+ classical monocytes reverses anti-PD-1/CTLA-4 immunotherapy resistance. Front. Immunol. 14:1161869. doi: 10.3389/fimmu.2023.1161869

Received: 08 February 2023; Accepted: 08 June 2023;
Published: 28 June 2023.

Edited by:

Michael P. Gustafson, Mayo Clinic Arizona, United States

Reviewed by:

Trupti Vardam-Kaur, Omeros Corporation, United States
Marion Curtis, Mayo Clinic Arizona, United States

Copyright © 2023 Rodriguez, Chen, Li, Miao, Peng, Fradette, Diao, Konen, Alvarez, Solis, Yi, Padhye, Gibson, Ochieng, Zhou, Wang and Gibbons. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Don L. Gibbons, dlgibbon@mdanderson.org

These authors have contributed equally to this work and share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.