Inhibition of the CCL2 receptor, CCR2, enhances tumor response to immune checkpoint therapy

1Advances in and approvals for the use of immune checkpoint-based therapies for the treatment of cancer has led a new wave of available second-, and more recently, first-line therapies to traditional chemotherapy and radiation regimens. Programmed cell death protein 1 (PD-1) directed therapies have been shown to re-invigorate CD8+ T cells, rescuing them from an exhausted state, leading to a renewed immune response and in many cases, durable clearance of the tumor. Targeting PD-1 both as a monotherapy or in combination with chemotherapy, targeted agents, or other types of immunotherapy– has been or is being tested in over 3000 clinical trials. While many of the patients in these trials have improved responses relative to standard of care, a significant number will not benefit from such therapies. This has led to a focus on identifying more effective therapeutic combinations,.

2The development of new anti-cancer drugs targeting cytokines and/or their cognate receptors has been an area of great interest in cancer treatment, whether for monotherapy or as an adjuvant in combination with other therapeutic agents. The use of IL-2 to activate the immune system of cancer patients was an early milestone in current cancer immunotherapy. IL-2-based immunotherapy proved that the immune system could eradicate tumor cells–. While relatively infrequent, complete remission and long-term disease-free survival has been reported in patients with melanoma and renal cell carcinoma following IL-2 treatment.

3To date, the majority of the novel combinations that have or are being tested clinically have involved empirical combinations of immune checkpoint-based therapies, with each other as well as with other agents in complementary therapeutic classes,. Newer approaches that identify successful combinations are based on an understanding of the mechanisms underlying the failure of patients to respond to monotherapy,, as well as studies employing functional genomic screens in experimental models,. Using an in vivo functional genomics screen, we identified the DDR2 receptor tyrosine kinase as an important determinant of the efficacy of immune checkpoint blockade therapy, where inhibition in combination with PD-1-targeted antibodies leads to enhanced tumor response compared to monotherapy. Similarly, Manguo et al. performed a loss-of-function in vivo genetic screen which identified PTPN2 as a potential target for enhancing immunotherapy. In another screen, Patel et al. identified APLNR as a critical gene for enhanced response to immunotherapy, where the loss of APLNR reduced the efficacy of T cell-based immunotherapies.

4One of the more compelling candidates from our screen was the cytokine, CCL2/MCP-1, monocyte chemotactic and activating factor/monocyte chemoattractant protein-1. CCL2 is a key molecule in macrophage chemotaxis and activation, and is implicated in the pathogenesis of several diseases including psoriasis, rheumatoid arthritis, asthma, and atherosclerosis–. Blockade of CCL2/CCR2 signaling provides protective immunity in murine models of OVA-induced allergic asthma. CCL2, produced by both cancer and stromal cells, preferentially binds to CCR2, which is expressed to varying degrees in a wide range of organs and tissues including blood, brain, heart, kidney, liver, lung, ovary, pancreas, spinal cord, spleen, and thymus. High levels of CCL2 have been identified in patients with lung adenocarcinoma who traditionally have poorer prognosis. Elevated expression of tumor and systemic CCL2 is associated with poor prognosis in breast cancer patients,. CCL2 is also overexpressed in human liver cancers and is prognostic for hepatocellular carcinoma patients. Interactions between CCL2/CCR2 have been shown to recruit immunosuppressive cells such as myeloid-derived suppressor (MDSC) cells and metastasis-promoting monocytes,. There is also compelling evidence for the targeting of CCL2/CCR2 in the treatment of various cancers. Knockout or blockade of CCL2/CCR2 inhibits primary liver tumor and metastatic growth leading to prolonged survival.

5Here we examined the biological significance of our screen findings which indicated blockade of the CCL2 receptor, CCR2, enhances the therapeutic efficacy of PD-1 inhibition in tumors. We confirmed, in multiple cancer types, that the combination anti-PD-1 and CCR2-targeted therapy leads to enhanced efficacy compared to either agent alone. Exploring the mechanisms of action underlying this finding revealed distinct differences in the immune cell populations and cytokine profile of the combination-treated tumors compared to that of monotherapies.

6We recently reported a functional genomics study evaluating selected U.S. Food and Drug Administration (FDA)–approved drugs, whose corresponding target genes, when inhibited could potentiate the response to anti–PD-1 immunotherapy. In this study, a murine bladder tumor cell line was established from N-butyl-N-(4-hydroxybutyl)nitrosamine (BBN)-induced tumors and adapted to in vitro cell culture to be used for screening purposes. The cells were utilized in an in vivo shRNA-based screen identified genes, that when knocked-down, showed enhanced tumor cell death mediated by the immune-activating anti-PD-1. To determine this, tumor samples were sequenced to quantify shRNA constructs which were absent and present, comparing between the isotype control and anti-PD-1-treated groups, with a goal of identifying genes that are preferentially lost in the anti-PD-1-treated. The study discovered and then validated DDR2 kinase as a promising target leading to the enhancement of response to anti-PD-1 immunotherapy.

7In addition to DDR2, several of the other genes, when depleted by shRNA seemed to enhance the effect of anti–PD-1 immunotherapy and thus may also be viable targets to consider. Hence, we began a multistep process to select additional targets from this screen for further investigation.

8First, we looked for genes whose expression pattern is consistently associated with the expression pattern of DDR2 in human tumors since these may conceivably be members of signaling pathways that a tumor could use to escape DDR2 depletion. Using two human bladder cancer datasets, CNUH and MSKCC, we identified 18 genes in each dataset whose expression pattern grouped with that of DDR2 (Fig. 1a). Interestingly, 16 of the 18 genes (80%) grouped with DDR2 in both datasets (p < 0.001 by Fisher’s Exact Test) (Fig. 1b) and from here on will be called the “DDR2 cluster”.

9To further determine the clinical relevance of this gene set, we developed an expression score from these 16 genes in the DDR2 cluster, and found this to be statistically significantly different in patients with less aggressive non-muscle invasive (NMI) compared to aggressive and advanced muscle-invasive (MI) tumors across 12 different bladder cancer patient cohorts (Fig. 1c). An area under the curve (AUC) < 0.50 indicates higher expression in NMI tumors, while an AUC > 0.50 indicates higher expression in MI tumors, and an AUC = 0.5 is what is expected by chance. High expression score of the 16 gene DDR2 cluster is associated with the more advanced MI disease (Fig. 1c). In contrast, a score generated from the genes outside the DDR2 cluster has an AUC < 0.5 in 10 of the 12 datasets indicating that these genes are not associated with the generally non-fatal NMI disease (Fig. 1c).

10Since PD-1 treatment alone does not affect tumor growth, we reasoned that any genes among the 16 that are higher in the PD-1-treated tumors compared to IgG controls, but unchanged or lower in responder tumors harboring shDDR2 and PD-1 treated, could conceivably be ideal additional targets for therapy. This gene selection strategy implies the expression of these genes may contribute to resistance to anti-PD-1 treatment. To determine if such genes exists, we grew control and shRNA depleted NA13 tumors, treated them with IgG or anti-PD-1 (Supplementary Fig. 1) and carried out RNA-seq as described. Interestingly, three genes survived these strict selection criteria: fibroblast activation protein alpha (FAP), CCL2 and CXCR4 (Table 1). Interestingly, while CCR2 is the primary receptor for CCL2, this ligand can also bind and activate CXCR4. Given the implication of macrophages via CCL2 as contributors to metastatic bladder cancer and potential for immunotherapy, this was of significant interest.

11To evaluate this further in human tumors, we used the BC-BET tool, and found that CCL2 expression is higher in patients with MI compared to NMI tumors in all 12 datasets examined (Fig. 1d). High CCL2 expression is also associated with worse disease-specific or overall survival in multiple datasets (Fig. 1e). Interestingly, it has been shown that CCL2 expression is associated with resistance to anti-PD-1 therapy. Taken together, these data suggest a strong rationale for investigating the efficacy of combination CCL2 inhibition with anti-PD-1 therapy in models of cancer.

12We looked to test the efficacy of targeting CCR2 in combination with PD-1 in multiple tumor models. The tumor cell lines, NA13, B16F10 and E0771, express CCL2 and PD-L1 in vitro (Supplementary Fig. 2) and PD-L1 in vivo (Supplementary Fig. 3) to varying levels. To determine the impact of reduced CCL2 activity on tumor growth when combined with anti-PD-1, we depleted CCL2 expression in NA13 cells using either of two different shRNAs (Supplementary Fig. 4a). Knockdown of CCL2 did not affect growth of NA13 in vitro (data not shown). NA13 cells expressing CCL2 shRNA (shCCL2) were inoculated into mice and on day 14 received treatment with IgG control or anti-PD-1. Tumor growth is hindered in the shCCL2 tumors that received anti-PD-1, but not in the shCCL2 tumors receiving IgG (Supplementary Fig. 4b). To determine the translational potential of this approach we targeted CCR2 using RS504393, a small molecule inhibitor of the CCL2 cognate receptor,. RS504393 is a highly selective CCR2 chemokine receptor antagonist used frequently in the prevention of CCL2/CCR2 interactions in vitro and in in vivo murine models–. When RS504393 is utilized in combination with anti-PD-1 immune checkpoint blockade, this limited NA13 tumor growth synergistically, and in some cases, this resulted in complete tumor clearance (Fig. 2).

13To assess the efficacy of CCR2 antagonist and anti-PD-1 combination treatment in models beyond bladder cancer as well as in metastatic disease, we examined the efficacy of this combination in the treatment of B16F10 melanoma pulmonary metastases. B16F10 is known to be poorly controlled by PD-1-based monotherapy–. In support of our hypothesis, we found that combination therapy is effective in reducing the number of pulmonary metastases compared to the control or mono therapy cohorts (Fig. 3a–c). These differences were confirmed through lung mass quantification wherein the control and single treatment groups exhibited higher tumor mass burden (Fig. 3b).

14To better understand the underlying mechanism accounting for the differences in tumor burden, we looked at the lung tumor microenvironment using flow cytometry (Supplementary Fig. 5). Analysis of the lungs from the four treatment cohorts indicate a difference in the Siglec-F+ alveolar macrophage (AM) population, while no statistically significant differences were seen in the monocytes (Mono), interstitial macrophages (IM) or dendritic cell (DC) populations (Fig. 3d, e). Siglec-F, used in conjunction with CD11b and CD11c, is a marker specific for murine lung-resident AMs that is not expressed by interstitial or inflammatory macrophages,. The statistically significant higher ratio of AMs in the combination-treated group (Fig. 3d, e), which had the lowest tumor burden, is in agreement with previous observations in which higher AM presence is associated with lower lung tumor burden. Siglec-F+ AMs preferentially remain localized in the healthy alveolar space outside of the tumor nodules and are progressively reduced as tumor burden increases. Thus, their greater presence in the combination-treated group provides further support to our observations of lower tumor burden in this group.

15Since the pro-tumorigenic role of the CCR2/CCL2 axis has been widely reported in breast cancer,,,,, we looked to test the combination therapy in a breast cancer model utilizing the E0771 cell line, which was originally derived from a spontaneous mammary tumor in a female C57BL/6 mouse. Once again, combination therapy with CCR2 antagonist and anti-PD-1 is more efficacious than monotherapy in reducing tumor growth (Fig. 4a–c). Flow cytometry-based analysis of the tumors reveal differences in the expression level of the T cell exhaustion markers LAG3 and PD-1 (Fig. 4d–f, Supplementary Fig. 6). A higher proportion of PD-1+ LAG3− on CD8+ T cells were identified in the tumors from mice treated with CCR2 antagonist and anti-PD-1 (Fig. 4e). In contrast, no statistically significant differences were seen in the PD-1+ LAG3+, PD-1− LAG3−, and PD-1− LAG3+ populations between untreated control and combination therapy (Fig. 4f, Supplementary Fig. 7). Overall expression levels of LAG3 and PD-1 are reduced following treatment with CCR2 antagonist and/or anti-PD-1 compared to control (Fig. 5a, b). The PD-1-treated and combination-treated tumors have lower levels of PD-1+ LAG3+ TIGIT+ CD8+ T cells (Supplementary Fig. 8a), though no differences were seen based on the expression of TIGIT alone (Supplementary Fig. 8b). Furthermore, the frequency of inhibitory FoxP3+ CD4+ regulatory T cells are reduced following single agent and combination treatment compared to control, though more statistically significant in the combination-treated group (Fig. 5c). The frequency of CD8+ T cells varies amongst the tumors from the four different treatment cohorts, with the control group possessing the lowest frequency of CD8+ T cells (Fig. 5d). The combination-treated group possesses a higher ratio of PD-1+LAG3− and IFNγ+ CD8+ T cells relative to Tregs (Supplementary Fig. 9). All three of the tumor cell lines, NA13, B16F10, E0771, express class I MHC at various levels indicating that they can be direct targets of CD8+ T cell recognition (Supplementary Fig. 10a). E0771-OVA cells expressing the ovalbumin peptide SIINFEKL were utilized to determine whether the CD8+ T cell response is tumor-specific. Tumor-infiltrating cells were stained with a tetramer to identify CD8+ T cells that recognize the tumor antigen. There is an increased average number of tetramer+ CD8+ T cells in the tumors from the combination-treated group, though not statistically significant. Of note, three of the five mice show a very strong response with increased recruitment of antigen-specific T cells with combination therapy (Supplementary Fig. 10b). No differences are seen in the frequency of monocytes and macrophages between control and combination-treated tumors. (Supplementary Fig. 11).

16Looking at bulk tumors through a cytokine panel analysis, statistically significant differences were seen in three cytokines across the treatment groups (Fig. 5e, Supplementary Fig. 12). Tumors from mice treated with combination therapy exhibit lower levels of IL-1rα and sICAM-1, but higher levels of IFNγ (Fig. 5e, Supplementary Fig. 12). Furthermore, while no differences are seen in the frequency of IFNγ+ CD8+ T cells, the raw number of cells per gram of tumor reveals a higher number of IFNγ+ CD8+ T cells in the combination-treated group when comparing on an equal tumor weight basis (Fig. 5f). Next, we performed CIBERSORT analysis to estimate relative immune cell populations of human bladder cancer RNA sequencing data from TCGA (Supplementary Fig. 13). In comparing immune population differences between the bottom and top quartile of CCL2 tumor expression, a higher presence of activated DCs and natural killer cells, both producers of IFNγ, are observed in the CCL2 low tumors (Supplementary Fig. 13). While CD8+ T cells are inferred to be higher in the CCL2 high tumors, analysis with CIBERSORT does not account for the activation or exhaustion status of these T cells (Supplementary Fig. 13). Though a higher presence of Tregs are observed in the CCL2 low tumors, inferences between the experimental data and CIBERSORT analysis possess a limitation in the conclusions that can be drawn. A key difference that should be noted is that CIBERSORT is looking at intrinsic, baseline CCL2 levels versus the mouse experiments which are treatment blockade-based. These human data will have considerably more variability/heterogeneity due to the fact that these are inferred estimates of immune populations that come from clinical sequencing samples from a tumor type, MIBC, that is known to have a high amount of inter- and intra-tumoral heterogeneity.

17Overall, we observed that tumors from mice treated with combination therapy possessed significant increases in cytotoxic T cell recruitment and activation, with a concomitant decrease in suppressive regulatory T cells. We believe this change in the T cell population to be correlative with the observed reduction in tumor and metastatic burden in combination-treated mice.

18Inhibition of CCL2/CCR2 has been explored in the clinical setting as a possible cancer therapeutic. While the available data are still sparse so far, no clear benefit for this therapy has emerged. A phase 1b trial in non-metastatic pancreatic cancer patients suggests that CCR2 inhibition decreases tumor-infiltrating macrophages and regulatory T cells, while also increasing effector T cells. Combination therapy of CCX872, a CCR2 specific antagonist, with FOLFIRINOX (fluorouracil [5-FU], leucovorin, irinotecan, oxaliplatin) suggests better overall survival with the combination compared to monotherapy (29% vs 19% at 18 months weeks) with no safety concerns. In contrast, a phase 1 trial (NCT00537368) of a human IgG1k monoclonal antibody which binds CCL2, carlumab, also known as CNTO888, in patients with advanced solid tumors refractory to conventional treatments showed no objective anti-tumor response in any of the 44 patients enrolled. A phase 2 study (NTC00992186) to assess the efficacy of carlumab in patients with metastatic prostate cancer reported that no patients responded. These human studies are in line with results that we report here, where we observed modest or no responses to monotherapy with CCL2/CCR2 blockade. An important finding common to all these studies is the general tolerability of the treatment with few mild-to-moderate adverse events. This makes CCL2/CCR2 inhibition ideal for combining with other therapies, such as with immune checkpoint blockade.

19The improved efficacy of the combination treatment when CCR2 antagonist and anti-PD-1 are used concurrently to target the tumor (Fig. 2) compared to CCL2 knockdown (Supplementary Fig. 4b) could be attributed to several reasons: better reduction in CCR2 receptor signaling that is achieved by reduction in CCL2; blockade of CCL2 being expressed by other cell types, such as CAFs, which are known to secrete CCL2 and promote macrophages in the tumor; and finally, non-specific effects of the inhibitor that serendipitously affect anti-PD-1 response. Since tumor intrinsic CCL2 may not be the only relevant source of CCL2, in our study we chose to target the CCL2 receptor, CCR2. Furthermore, CCR2 antagonists provide greater real-world relevance in which the findings from our study could more readily be applied due to the use of CCR2 antagonists in clinical trials. BMS-741672, BMS-813160, PF-04634817 are three CCR2 antagonists which have been or are currently being tested in clinical trials (ClinicalTrials.gov Identifier: NCT04123379, NCT00699790, NCT01712061, NCT00699790). Indeed, based on our findings, therapy targeting CCL2/CCR2 in combination with blockade of PD-1 may produce synergistic responses in some patients. In addition to blocking CCL2/CCR2 and PD-1 separately, development of bispecific antibodies would be a novel and elegant solution towards inhibiting these molecules in the cancer cell.

20In our analysis of the tumor-infiltrating immune cells, higher proportions of PD-1+ LAG3+ CD8+ T cells were present in the tumors of control mice. Like the immune checkpoint PD-1, LAG3 is a co-inhibitory receptor which suppresses T cell activation and cytokine secretion. LAG3 and PD-1 are commonly described as markers of exhausted T cells with combination targeting of the two being effective in eliciting a strong T cell response leading to tumor clearance. LAG3 has been reported to synergize with PD-1, potentially amplifying the PD-1 inhibitory effect. Thus, it is possible that downregulation of LAG3 by this novel combination treatment might alter responsiveness to immune checkpoint blockade. Recently, progenitor, transitional, and terminally exhausted subsets of CD8+ T cells were identified in response to chronic stimulation, with differential responsiveness to checkpoint blockade. In that study, LAG3 was most highly expressed in the terminally exhausted subset. The appearance of PD-1+ LAG3− cells correlates with treatment effectiveness, but whether these cells represent a specific type of effector cell, are direct progenitors of effector cells, or are another transitional population of cells will require more extensive analysis of their transcriptome in the future.

21While we observed few changes in the myeloid compartment after blocking the CCL2/CCR2 axis combined with or without anti-PD-1, an increase in the Siglec-F+ AM population was seen (Fig. 3d, e). Siglec-F is a marker specific for murine lung-resident AMs that is not expressed by interstitial or inflammatory macrophages,. AMs play an important role in maintaining lung homeostasis. AMs remove clear apoptotic cells, environmental particles, and pathogens by interacting with cells in the alveolar epithelium via cell-surface receptors and cytokines.

22Siglecs belong to the transmembrane lectins family that are expressed on immune cells and binds to sialic acids. Cancer cells express high levels of sialic acids that can interact with Siglecs on immune cells. Siglecs can contribute to both facilitation as well as attenuation of anti-tumor activity. Their anti-tumor effect is due to the inhibition of tumor-promoting inflammation. For example, macrophages that express Siglec-1 (CD169) have been shown to present tumor antigens to cytotoxic T cells, and Siglec-1 deficient macrophages result in inhibition of anti-tumor immunity. The function of Siglec-E, the murine equivalent of human Siglec-9, have been tested in subcutaneous syngeneic mouse model. Once tumors are established, Siglec-E deficiency results in more aggressive tumor growth compared to wildtype.

23In a study by Poczobutt et al., multiple subsets of myeloid cells were identified: MacA cells, which bear markers of AMs (Siglec-F+/CD11c+), MacB1 cells (CD11b+/CD64lo/CD11c+), MacB2 cells (CD11b+/CD64int/CD11c2), and MacB3 cells (CD11b+/CD64hi/CD11c+). Using data from PRECOG, they reported that MacA represents a resident AM population that is enriched in genes predicting good clinical outcomes.

24We observed a higher number of IFNγ+ CD8+ T cells in the combination-treated group. Concurrently, we also observed increased levels of IFNγ, as measured through a multi-cytokine panel of bulk tumor. In the clinical setting, high IFNγ expression or associated signatures is associated with improved outcomes to immune checkpoint blockade,. Metastatic non-small cell lung cancer (NSCLC) and urothelial cancer patients who received PD-L1 inhibitor and subsequently exhibited an increased IFNγ gene signature (IFNγ, CD274, LAG3, and CXCL9) had better overall response rates and longer median progression‐free survival, independent of PD-L1 expression. Moreover, in melanoma and NSCLC patients treated with anti-PD-1 therapies, higher IFNγ protein expression is associated with longer progression-free survival. Importantly, this suggests that the higher observed proportion of IFNγ and IFNγ-expressing cells in combination-treated tumors from our study are indicative of not only reduced tumor growth in mice but is more likely to be clinically translatable to patient outcomes. Taken together, this finding, the known tolerability in patients of CCR2-targeted therapies, and evidence that the combination is more effective across multiple tumor models, we believe clinical trials to be strongly justified.

25NA13 cell line was isolated and cultured from N-butyl-N-(4-hydroxybutyl) nitrosamine (BBN) carcinogen-induced bladder tumor of C57BL/6 female mice. E0771 was a gift from Dr. Traci Lyons (University of Colorado). E0771 was maintained in RPMI with 5% fetal bovine serum (FBS). B16F10 was obtained from the American Type Culture Collection (ATCC) through the University of Colorado Tissue Culture Core and maintained in DMEM media with 10% FBS. All lines have been authenticated and tested to be mycoplasma-free. All lines were grown at 37 °C in a humidified atmosphere (5% CO2).

26Cells were homogenized using QIAshredder (Qiagen), followed by RNA extraction using a RNeasy Plus Mini Kit with gDNA Eliminator (Qiagen). cDNA was synthesized using iScript Reverse Transcription Supermix (Bio-Rad). qPCR was then performed (Quant Studio 6 Flex Real-Time PCR system, Applied Biosystems, USA), using iQ SYBR Green Supermix (Bio-Rad). The following primer pairs were used: Mouse CCL2: forward 5′ AGTAGGCTGGAGAGCTACAA 3′; reverse 5′ GTATGTCTGGACCCATTCCTTC 3′ and mouse PD-L1: forward 5′ TCCATCCTGTTGTTCCTCATT 3′ Reverse 5′ TCCACATCTAGC ATTCTCACTTG 3′. To determine the changes in mRNA expression as measured by qRT-PCR, the ΔΔCt method was used. Expression was normalized to internal control β-actin.

27Western Blot analysis was performed using total protein extracted from NA13, B16F10, and E0771 cell lines using RIPA buffer (Sigma, USA) containing protease and phosphatase inhibitors (Roche). Following total protein quantification, equal amounts of protein were separated by 4–20% Mini-PROTEAN® TGX™ Protein Gel (BioRad) and transferred to PVDF membrane. The membrane was probed with antibodies against PD-L1 (MAB90781, R&D Systems, USA) or anti-β-actin (13E5, Cell Signaling, USA) Rabbit monoclonal antibodies diluted at 1:1000 in 5% non-fat milk blocking buffer. The blots were imaged using iBright imaging system according to the manufacturer’s instruction (ThermoFisher Scientific, USA).

28Immunohistochemistry was performed on 5-µm thick sections of formalin-fixed paraffin-embedded (FFPE) NA13 tumors blocks using a Rabbit anti-mouse antibody (D5V3B, Cell Signaling Technology, USA) at a dilution of 1:100. Antigen retrieval was performed according to standardized protocols by heating with 10 mM citrate buffer. The activity of endogenous peroxidase (peroxidase blocking reagent, Dako) were neutralized and non-specific binding were then blocked with goat serum 1:100 in PBS. Tissue sections were incubated with primary antibodies. The slides were then incubated with anti-rabbit HQ secondary antibody followed by anti-HQ HRP linking antibody (Ventana Medical System, USA). Finally, slides were incubated with Discovery ChromoMap DAB reagent. Slides were counterstained with hematoxylin, dehydrated by ethanol and mounted. Slides were imaged on an Olympus BX43 light microscope with ×20 magnification using an Olympus DP26 digital camera.

29The response of mouse tumor cells to rIFNγ (R&D system) was tested after incubation of 200,000 cells/well in 6-well tissue culture plate in the presence of rIFNγ (50 and 250 unites) for 48 h followed by FACS analysis of MHC class I antigen expression using APC anti-mouse H-2Db and APC mouse IgG2a K isotype control (FC) antibodies (Biolegend).

30Female C57BL/6 mice (Charles River) were received at 6-week old and allowed to acclimate for at least one week in sterile micro isolator cages with constant temperature and humidity. Mice had free access to food and water. Mice were housed in specific-pathogen-free conditions and cared for in accordance with US National Institutes of Health guidelines, and all procedures were approved by the University of Colorado Denver Animal Care and Use Committee and carried out according to approved protocols.

31For NA13, mice were injected with 1 × 106 cells in 100 µL sterile PBS subcutaneously in the hind flank. For the E0771 mammary tumor model, mice were injected with 5 × 104 cells in the third thoracic mammary fat pad. For the B16F10-induced pulmonary metastases model, mice were inoculated intravenously with 2 × 105 B16F10 cells in 100 µL sterile PBS. Mice were examined twice weekly. Tumor measurements commenced from when the tumor was first palpable. Tumor size was determined using an electronic caliper to measure the length and width and calculated by (L × W2)/2, where L is the largest diameter measurement of the tumor and W is the shorter perpendicular tumor measurement. Measurements were taken from distinct samples of each tumor-bearing mouse. Animals were randomized into treatment groups ensuring similar average tumor volumes amongst the groups, weighed and identified via ear punch.

32Mouse anti-PD-1 antibody (IgG1-D265A) and isotype control (IgG1, clone 4F7) were produced by Bristol-Myers Squibb laboratories (Redwood City, CA) and was formulated in PBS and administered intraperitoneally at a dose of 50 µg/mouse (NA13, E0771) or 100 µg/mouse (B16F10) for a total of three doses. CCR2 small molecule inhibitor RS504393 (Tocris) was given daily at 2 mg/kg by oral gavage to B16F10 mice and intraperitoneally to E0771 and NA13 mice.

33B16F10 mice were euthanized 4 weeks after tumor cell injection. Lungs were extracted and the number of visible metastases was quantified, then lungs were processed and analyzed. E0771 tumors were mechanically dissociated in Click’s media in the absence of mercaptoethanol or L-glutamine (Irvine Scientific). Cells were digested for 1 h at 37 °C with 500 units/ml collagenase type II and IV and 20 μg/ml DNase (Worthington Biochemical). The digested tissue suspension was then filtered through a 100 μm strainer. Filtered cells were carefully layered into a centrifuge tube containing 5 ml Lympholyte-M (Cedarlane). The cells were centrifuged at 1500 × g for 20 min., then the interface lymphocyte layer was carefully removed. The cells were washed prior to staining. For the CD8 T cell panel, cells were stained with CD8 APC/Cy7 (clone 53-6.7) (1:400), CD3 FITC (clone 17A2) (1:300), CD45 BV510 (clone 30-F11) (1:300), CD44 BV421 (clone IM7) (1:400), PD-1 PE (clone 29F.1A12) (1:200), LAG3 PerCP-Cy5.5 (clone C9B7W) (1:100), and Tigit APC (clone1G9) (1:100). CD8 T cells were gated on live, CD3+/CD8+ double-positive cells. The cells were then further classified based on the expression of PD-1 and Lag-3. For the macrophage panel, cells were stained with CD45 BV510, F4/80 APC/Cy7 (clone BM8) (1:100), CD11b BV421 (clone M1/70) (1:600), CD64 PerCP-Cy5.5 (clone X54-5/7.1) (1:200), MERTK FITC (clone 2B10C42) (1:100), PD-L1 PE (clone 10F.9G2) (1:200), and Ly-6G APC/Cy7 (clone 1A8) (1:100). Live cells were gated for CD45+ cells. Neutrophils were gated by the expression of Ly-6Ghi/CD11b+, and not included in further analysis. Macrophages were gated as F4/80+/CD11bhi population and MERTKhi/CD64hi population. Monocytes were confirmed with two population gates as F4/80lo/CD11b+ and MERTKlo/CD64lo–. For the tetramer panel, cells were stained with CD8 APC/Cy7, CD35 BV410, CD44 BV421, LAG3 PerCP Cy5.5, PD-1 APC (clone 29F.1A12) (1:100), CD4 FITC (clone GK1.5) (1:200), and MHC-I-SIINFEKL tetramer PE (NIH) (1:200). Live CD8+ cells were gated CD44hi/tetramer+. Cells were further analyzed for expression of PD-1 and Lag-3. IMs are identified based MERTKhi CD64hi CD11bhi, and subdivided into three pulmonary IM subtypes: IM1 (CD11clo CD206+ MHCIIlo), IM2 (CD11clo CD206+ MHCIIhi), and IM3 (CD11chi CD206lo MHCIIhi). All flow cytometry antibodies used were purchased from Biolegend unless otherwise indicated. All samples were run on the CyAn ADP flow cytometer, acquired using Summit software (BD Biosciences), and analyzed using FlowJo software (Tree Star). Flow cytometry data were taken from distinct samples and not the same samples repeatedly measured.

34Cells were isolated from the tissue and treated with or without (controls) phorbol 12-myristate 13-acetate (PMA) (20 ng/ml) (Sigma, St. Louis, MO) plus ionomycin (1 µg/ml) (Sigma, St. Louis, MO) for 4–6 h at 37 °C in the presence of 2 µg/ml of brefeldin A (Adipogen, San Diego, CO) in RPMI + 2.5% FBS. Cells were then stained with CD8 APC/Cy7, CD45 BV510, CD44 BV421, and CD4 PerCP-Cy5.5 (clone GK1.5) (1:400). Cells were incubated at 37 °C for 30 min. Following incubation, cells were washed, then fixed with 1% paraformaldehyde and 4% sucrose for 10 min at room temperature in the dark. Cells were then permeabilized with BD Perm Wash (BD Biosciences) and stained for cytokines IFNγ APC/Cy7 (XMG1.2) (1:200). Cells were incubated overnight at 4 °C in the dark, then washed with Perm Wash and resuspended in FACs buffer (0.5% Bovine Serum Albumin and 0.1% Sodium Azide in PBS).

35For FoxP3 analysis, cells were stained for CD8 APC/Cy7, CD45 BV510, CD4 perCP-Cy5.5, B220 FITC (clone RA3-6B2) (1:300), CD25 APC (clone 3C7) (1:200), and PD-1 BV421 (clone 29 F.1A12) (1:100). Cells were incubated in antibody at 37 °C for 30 min. Cells were fixed and permeabilized for 30 min at 4 °C in the dark using the FoxP3 Transcription Factor Staining Buffer Set (eBioscience; 00-5523). The cells were then stained with FoxP3 PE (MF-14) (1:200) and incubated in the dark overnight at 4 °C. Cells were gated for live CD45+, followed by gating for the CD4+/CD8− population. Further classification of activated CD4 T cells as CD25+/FoxP3− and regulatory CD4 T cells as FoxP3+.

36Library preparation and sequencing were performed by Novogene. Sequencing was performed on an Illumina platform paired end 150 bp with 20 million reads per sample.

37Transcript quantification was done using RSEM (v1.2.31) with default parameters and Bowtie2 (v2.1.0) as the read aligner. Reads were mapped directly to mouse transcripts and summarized at the gene level using annotations from Ensembl r91, genome build GRCm38.p5. Quantification of genes as expected counts were compiled. Differential expression was performed using voom function in the limma R package. Genes with an average expected count <5 were removed, normalization factors were calculated, and comparisons between groups were made using the voom function using default parameters.

38Tumors were isolated and flash frozen in liquid nitrogen. Tissue was homogenized in PBS containing complete mini protease inhibitor cocktail (MilliporeSigma). Following homogenization, Triton X-100 was added to a final concentration of 1%, then frozen at –80 °C overnight. Samples were thawed on ice and centrifuged at 10,000 × g for 5 min, and supernatants were collected. A BCA protein assay was performed to determine lysate protein concentration (Pierce, Thermo Scientific). Purified protein lysates were applied to the Proteome Profiler Mouse Cytokine Array Kit, Panel A as per manufacturer’s instructions (R&D Systems, catalog ARY006). The chemiluminescence reaction was measured with Biorad ChemiDoc MP Imaging system.

39Preliminary analysis using the Bladder Cancer Biomarker Evaluation Tool (BC-BET) was used to evaluate the association of CCL2 expression with bladder cancer characteristics. The 12 patient cohorts (n = 1257) included in BC-BET were then downloaded. In brief, processed data were downloaded from the Gene Expression Omnibus (AUH-1 cohort: Accession #GSE3167; AUH-2: GSE5479; CNUH, GSE13507; DFCI: GSE31684; Lindgren: GSE19915; Lindgren-2: GSE32548; MDA-1: GSE48276 and MDA-2: GSE48075; UVA: GSE37317, from Array Express, (Stransky-1 and Stransky-2: E-TABM-147 or as Supplemental Material to Blaveri and MSKCC. Probes for target genes were identified from the microarray platform annotation. When multiple probes for a gene were available, probe with the highest mean expression was used. Hierarchical agglomerative clustering with complete linkage was used, based on the Euclidian distance between samples. The DDR2 cluster score was calculated by normalizing the expression of each of the 16 probes to have a mean of 0 and a standard deviation of 1, and then finding the mean of all signature probes for each sample.

40CIBERSORT using TCGA bladder cancer gene expression was processed according to The Cancer Immune Atlas (TCIA) and the associated data as downloaded directly from the TCIA website.

41The differences between the groups were evaluated by Student’s t test, one-way or two-way ANOVA using GraphPad Prism 7.0 Software. The specific type of statistical analysis utilized is indicated in the corresponding figure legends. Exact p-values are indicated in the figure. A value of p < 0.05 was considered to be statistically significant. Experiments were repeated at least two times to ensure reproducibility.

42Further information on research design is available in the Nature Research Reporting Summary linked to this article.