Next Article in Journal
In Vivo Application of Carboranes for Boron Neutron Capture Therapy (BNCT): Structure, Formulation and Analytical Methods for Detection
Next Article in Special Issue
Recent Advances in Therapeutic Strategies to Improve Colorectal Cancer Treatment
Previous Article in Journal
Selective Internal Radiotherapy (SIRT) and Chemosaturation Percutaneous Hepatic Perfusion (CS-PHP) for Metastasized Uveal Melanoma: A Retrospective Comparative Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

TIM-3 Expression and M2 Polarization of Macrophages in the TGFβ-Activated Tumor Microenvironment in Colorectal Cancer

1
Department of Gastrointestinal Tract Surgery, Fukushima Medical University School of Medicine, Fukushima 960-1295, Japan
2
Department of Multidisciplinary Treatment of Cancer and Regional Medical Support, Fukushima Medical University School of Medicine, Fukushima 960-1295, Japan
3
Department of Blood Transfusion and Transplantation Immunology, Fukushima Medical University School of Medicine, Fukushima 960-1295, Japan
*
Author to whom correspondence should be addressed.
Cancers 2023, 15(20), 4943; https://doi.org/10.3390/cancers15204943
Submission received: 16 September 2023 / Revised: 29 September 2023 / Accepted: 9 October 2023 / Published: 11 October 2023
(This article belongs to the Special Issue Novel Strategies in the Prevention/Treatment of Colorectal Cancer)

Abstract

:

Simple Summary

TGFβ and TIM-3 have been implicated in driving immune evasion and poor response to anti-PD-1/PD-L1 immune checkpoint inhibitors (ICIs). Using large-scale transcriptomic data analysis, followed by immunohistochemistry and in vitro experiments, we demonstrated a strong link between TIM-3 and M2-like polarization of macrophages in the TGFβ-rich tumor microenvironment in colorectal cancer (CRC). We found that the transcriptional signatures of TIM-3, M2 macrophage and TGFβ-responsive stromal activation were specifically derived from the stromal components and were correlated with each other, and that they were robustly upregulated in CMS4, CMS1 and MSI CRCs. Our study suggests that the blockade of TGFβ or TIM-3, combined with ICIs, may be a promising strategy to circumvent the resistance to current ICIs for CMS4, CMS1 and MSI CRCs, where TIM-3 is highly expressed along with increased TGFβ.

Abstract

TGFβ signaling in the tumor microenvironment (TME) drives immune evasion and is a negative predictor of immune checkpoint inhibitor (ICI) efficacy in colorectal cancer (CRC). TIM-3, an inhibitory receptor implicated in anti-tumor immune responses and ICI resistance, has emerged as an immunotherapeutic target. This study investigated TIM-3, M2 macrophages and the TGFβ-activated TME, in association with microsatellite instability (MSI) status and consensus molecular subtypes (CMSs). Transcriptomic cohorts of CRC tissues, organoids and xenografts were examined (n = 2240). TIM-3 and a TGFβ-inducible stromal protein, VCAN, were evaluated in CRC specimens using immunohistochemistry (n = 45). TIM-3 expression on monocytes and generated M2 macrophages was examined by flow cytometry. We found that the expression of HAVCR2 (TIM-3) significantly correlated with the transcriptional signatures of TGFβ, TGFβ-dependent stromal activation and M2 macrophage, each of which were co-upregulated in CMS4, CMS1 and MSI CRCs across all datasets. Tumor-infiltrating TIM-3+ immune cells accumulated in TGFβ-responsive cancer stroma. TIM-3 was increased on M2-polarized macrophages, and on monocytes in response to TGFβ treatment. In conclusion, we identified a close association between TIM-3 and M2-like polarization of macrophages in the TGFβ-rich TME. Our findings provide new insights into personalized immunotherapeutic strategies based on the TME for CRCs.

1. Introduction

Colorectal cancer (CRC) is the most common gastrointestinal cancer and the third leading cause of cancer-related deaths worldwide [1]. CRC is a heterogeneous disease, with different sets of genetic events associated with diverse clinical characteristics and drug sensitivity. During the past decade, the genomic landscape of CRC has been well-established, and treatment strategies have been personalized according to genomic alterations, including RAS and BRAF mutations and microsatellite instability (MSI); however, specific therapies are only available for a small proportion of patients [2,3]. Numerous studies have recently shown that the tumor microenvironment (TME), in which cancer cells are surrounded by cancer-associated fibroblasts (CAFs), T-cells, tumor-associated macrophages (TAMs), many other cells and their extracellular matrix (ECM) proteins, may be the key driver of clinical and therapeutic outcomes in CRC [4,5,6]. In 2015, the international consortium identified four consensus molecular subtypes (CMSs) based on the gene expression profiles of bulk tumors that encompass mixed signals derived from both cancer cells and stromal cells [7]. CMS1 is closely associated with MSI, mutant BRAF and immune infiltration. CMS2 shows canonical WNT and MYC signaling activation. CMS3 is enriched for RAS mutations and metabolic dysregulation. CMS4 is defined as the mesenchymal subtype that displays prominent transforming growth factor-β (TGFβ) activation, stromal infiltration, angiogenesis and an immunosuppressive TME [2,7]. The CMS classification not only provides valuable insight into tumor biology and therapeutic strategies, but also stratifies clinical outcomes. Patients with CMS4 tumors exhibited the poorest survival rates. In this context, studies have focused more on the role of TGFβ in shaping the tumor-promoting TME in CRC, supporting the potential of TGFβ inhibition in cancer immunotherapy [4,6,8,9,10].
Dysfunctional exhausted T-cell subsets have been shown to express multiple co-inhibitory or checkpoint receptors, such as programmed cell death-1 (PD-1), cytotoxic T-lymphocyte-associated protein 4 (CTLA-4), the T-cell immunoglobulin domain and mucin-domain 3 (TIM-3), lymphocyte-activation gene 3 (LAG-3) and T-cell immunoglobulin and the ITIM domain (TIGIT) [11]. Cancer immunotherapy with immune checkpoint inhibitors (ICIs), particularly antibodies against PD-1 and its ligand programmed death ligand 1 (PD-L1), as well as CTLA-4, has revolutionized the therapeutic management of several solid cancers, such as lung cancer and melanoma. These ICIs have also been applied to CRC treatment; however, they are only effective for 4–5% of CRC exhibiting MSI [2,12,13]. MSI tumors are likely to be immunogenic; mismatch-repair deficiency (dMMR) causes hypermutation and tumor neoantigens, contributing to the creation of the “immune hot” TME. In several clinical trials, anti-PD-1 antibodies, with or without CTLA-4 blockade, demonstrated dramatic and durable responses with prolonged survival in patients with metastatic MSI CRC, although approximately 50% of them had primary resistance to ICIs [12,13]. On the other hand, microsatellite stable (MSS) CRCs showed limited or no clinical responses to ICIs, mostly because of this type being an “immune cold” tumor [13,14]. The distinct responses to ICIs among MSI and MSS CRCs suggest that ICI resistance may be due to an immunosuppressive TME driven by several molecular programs. To this end, Tauriello et al. used a mouse model of MSS CRC displaying the CMS4 phenotype and TGFβ-activated stroma, demonstrating that increased TGFβ in the TME represented a primary mechanism of immune evasion and that TGFβ inhibitors rendered tumors susceptible to PD-1/PD-L1 blockade [9]. In our previous work, we analyzed multiple patient cohorts and identified a CMS4-like subset within MSI CRC, demonstrating the TGFβ-rich TME along with a high ICI resistance signature, suggesting that strategies of TGFβ inhibition combined with ICIs may have potential therapeutic value for some MSI CRCs [15]. Indeed, TGFβ signaling has emerged as a central mediator of immune evasion and ICI resistance across a wide range of cancers, and many clinical trials are currently testing the combination of TGFβ and PD-1/PD-L1 immune checkpoint inhibition in patients with solid cancers [10,16,17,18,19,20].
TIM-3, encoded by the HAVCR2 gene, has recently emerged as an immunotherapy target, as it is often co-expressed alongside PD-1 on the most dysfunctional and exhausted population of tumor-infiltrating T-cells [21,22]. Moreover, TIM-3 may have a role in adaptive resistance mechanisms to ICIs, as it is upregulated in response to anti-PD-1 therapy, and the co-blockade of TIM-3 and PD-1 can overcome resistance to PD-1 therapy in preclinical models [22,23]. Furthermore, the expression of TIM-3 can serve as a marker of poor prognosis in some cancer types, including CRC [21,24,25]. Although TIM-3 was initially thought to be a sign of dysfunctional T-cells, it has been demonstrated that TIM-3 is highly expressed on several types of immune cells, particularly myeloid cells, and that myeloid cells might play a predominant role in TIM-3-mediated anti-tumor immunity [21,26,27,28,29]. However, the mechanisms underlying TIM-3 regulation in non-T-cell populations remain largely undetermined [30]. We recently demonstrated that MSI CRCs exhibiting the TGFβ-rich TME were not only characterized by ICI resistance and poor survival but also the upregulation of the M2 macrophage signature along with HAVCR2 (TIM-3) expression [15]. In this study, we analyzed the relationship between HAVCR2, M2 macrophages, TGFβ-activated stroma and molecular subtypes in multiple transcriptomic and immunohistochemistry cohorts of CRCs and investigated the expression of TIM-3 in M2-polarized macrophages and TGFβ-stimulated monocytes in vitro.

2. Materials and Methods

2.1. Microarray and RNA-Seq Data Analysis, Gene Signatures

The transcriptomic datasets used in this study are summarized in Table S1. We used the preprocessed values obtained from each dataset. If a gene was represented by multiple probes, only the probe with the highest mean expression was used. We used the following six datasets of CRC patient samples, for which CMS classification was obtained from the CRCSC synapse data portal accessed on 3 February 2022 (https://www.synapse.org/#!Synapse:syn2623706/files/) [7]. The Cancer Genome Atlas (TCGA) RNA-seq dataset (COADREAD) was obtained through cBioPortal (https://www.cbioportal.org/) and LinkedOmics (http://linkedomics.org/cptac-colon/), as previously described [3,15,31,32]. Four microarray datasets, GSE39582, GSE33113, GSE14333 and GSE37892, were downloaded from the Gene Expression Omnibus (GEO) repository (http://www.ncbi.nlm.nih.gov/geo), and a microarray dataset KFSYCC was downloaded from the CRCSC synapse data portal (Table S1). In addition, we used an RNA-seq dataset (GSE171682) of CRC organoids and matched primary CRC tissues [33], and a microarray dataset (GSE100550) that profiled CRC patient tumors, patient-derived xenografts, patient-derived organoids and cell lines (Table S1) [34].
The gene expression-based signatures used in this study were previously described [15]. Briefly, we built a TGFβ-stromal signature, designated TBSS, on the basis of a set of TGFβ-responsive stromal genes, including VCAN, SERPINE1, CALD1, FAP, POSTN and IGFBP7 [35]. Tumors were dichotomized into TBSSLow or TBSSHigh based on the median. We also estimated the overall TGFβ levels based on the expression of TGFB1 and TGFB3, designated TGFB [4]. The M2 macrophage signature was built based on a set of genes encoding M2 surface markers, including CD163, MSR1 (CD204), MRC1 (CD206), CD200, FCER2 (CD23), CD36, CD209 and SLAMF1 (CD150) [15]. We also generated a set of genes assembled from published signatures, including Stromal ESTIMATE, which infers the faction of stromal cells, the C-ECM-UP signature, correlating with the dysregulation of ECM, and mesenchymal genes associated with epithelial-to-mesenchymal transition (EMT), designated stromal/ECM/mesenchymal genes [16,36,37].

2.2. Patient Samples

We enrolled forty-five consecutive patients with stage I–IV colorectal adenocarcinoma who underwent surgery at Fukushima Medical University Hospital between 2019 and 2021 without preoperative treatment, such as self-expanding metal stents, chemotherapy or radiotherapy. We used formalin-fixed paraffin-embedded (FFPE) whole tissue sections for immunohistochemistry. Tumors were classified according to the Japanese Classification of Colorectal, Appendiceal and Anal Carcinoma. The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of Fukushima Medical University.

2.3. Immunohistochemistry

Four-μm thick FFPE sections were deparaffinized in xylene and rehydrated in ethanol. Endogenous peroxidases were blocked with 0.3% hydrogen peroxide in methanol. Antigens were retrieved by autoclave in 10 mM citrate buffer solution (105 °C, pH 6.0) for 5 min for VCAN, or by autoclave in Target Retrieval Solution (Agilent Technologies, Santa Clara, CA, USA) (100 °C, pH 9.0) for 10 min for TIM-3. Slides were incubated with primary rabbit polyclonal anti-VCAN antibody (HPA004726, Prestige Antibodies® Powered by Atlas Antibodies, Sigma-Aldrich, St Louis, MO, USA; 1:500) or primary rabbit monoclonal anti-TIM-3 antibody (D5D5R XP®; #45208; Cell Signaling Technology, Danvers, MA, USA; 1:400) at 4 °C overnight, and then detected using a horseradish peroxidase-coupled anti-rabbit polymer (K4003; Envision+ system, Agilent Technologies). Peroxidase was visualized with diaminobenzidine (Dojindo, Kumamoto, Japan), and nuclei were counterstained with hematoxylin.

2.4. Assessment of Staining

Digital images were acquired using a digital slide scanner (NanoZoomer-SQ; Hamamatsu Photonics, Hamamatsu, Japan). The scanned images were independently evaluated by two observers (M. Katagata and H. Okayama), who had no prior knowledge of the clinical data. The assessment of VCAN staining was previously described [35]. Briefly, the percentage positivity of stromal VCAN was scored as 0 (0–5%), 1 (5–25%) or 2 (>25%) in the tumor stromal area, and the intensity of staining was scored as 0 (negative), 1 (moderate) or 2 (strong). The scores were then combined to generate a VCAN score (min 0, max 4). For TIM-3 staining, we evaluated five independent high-power microscopic fields at the hotspot areas of each tumor. The proportion of TIM-3-stained tumor-infiltrating immune cells (TIICs) was graded as 0 (<1%), 1 (1–5%), 2 (5–20%), 3 (20–50%) or 4 (>50%), and the intensity of staining was graded as 0 (negative), 1 (weak), 2 (moderate) or 3 (strong). We then multiplied the scores to obtain the total staining score for one field (min 0, max 12), and the final TIM-3+ TIICs score for each tumor was the mean of the five fields. We further selected tumors with a VCAN score of 3 or 4 that showed a heterogeneous VCAN staining pattern, and then evaluated the TIM-3+ TIICs in stromal areas with strong VCAN staining (VCANHigh areas) and those with negative or moderate VCAN staining (VCANLow areas).

2.5. Immunofluorescence

The FFPE sections were deparaffinized and rehydrated. Antigens were retrieved by autoclave for 20 min in Target Retrieval Solution (Agilent Technologies) (100 °C, pH 9.0). Blocking was performed with 2% donkey serum in TBS/T for 30 min. Thereafter, the sections were stained with anti-human CD163 antibody (NCL-L-CD163; Leica biosystems, Wetzlar, Germany; 1:200) and anti-human TIM-3 antibody (D5D5R XP®; #45208; Cell Signaling Technology; 1:200) at 4 °C overnight. The sections were incubated for 1 h with Alexa Fluor 488-conjugated anti-mouse (A-21202; Thermo Fisher Scientific, Waltham, MA, USA; 1:200) and Alexa Fluor 555-conjugated anti-rabbit (A-31572; Thermo Fisher Scientific; 1:200) secondary antibodies. Nuclei were stained with DAPI (D9542; Sigma-Aldrich; 40 ng/mL) for 10 min. Finally, the slides were mounted with ProLong Glass Antifade Mountant (P36984; Thermo Fisher Scientific), and images were obtained using an Olympus FV1000-D confocal microscope.

2.6. Flow Cytometry Analysis

Cells were stained with antibodies, including APC/Cyanine7-conjugated anti-human CD14 mAb (325620, BioLegend, San Diego, CA, USA), PE/Cyanine7-conjugated anti-human CD163 mAb (333613, BioLegend), Pacific Blue-conjugated anti-human CD86 mAb (305423, BioLegend), 7-AAD Viability Staining Solution (420404, BioLegend), PE-conjugated anti-human TIM-3 mAb (345005, BioLegend) and PE-conjugated isotype control mAb (400111, BioLegend). Unstained samples were used as the negative controls. The stained cells were measured using a BD FACS Canto II flow cytometer (BD Biosciences, San Jose, CA, USA) and the data were analyzed using FlowJo software (version 10.8.1, FlowJo, Ashland, OR, USA).

2.7. Generation of M2 and M1 Macrophages

M2 and M1 macrophages were generated by the generally established methods [38,39]. Human peripheral blood mononuclear cells (PBMCs) were collected from the peripheral blood of healthy volunteers using a BD Vacutainer® CPT™ (BD Bioscience), and monocytes were isolated using a MidiMACS™ Separator (130-042-302, Miltenybiotec, Bergisch Gladbach, Germany), LS column (130-042-401, Miltenybiotec) and CD14 MicroBeads human (130-050-201, Miltenybiotec). CD14+ monocytes were cultured in a 12-well plate at a density of 5 × 105 per well at 37 °C in 5% CO2 atmosphere in RPMI1640 (Sigma-Aldrich) supplemented with 20 ng/mL of M-CSF (216-MC, R&D systems, Minneapolis, MN, USA), penicillin/streptomycin and 10% heat-inactivated FBS for 6 days. Subsequently, we continued to culture the cells in the same medium for 3 days to generate M0 macrophages, while we added 20 ng/mL each of IL-4 (204-IL, R&D Systems), IL-10 (217-IL, R&D systems) and IL-13 (213-ILB, R&D systems) into the culture medium and continued to culture the cells for 3 days for the generation of M2 macrophages. To differentiate M1 macrophages, we utilized 20 ng/mL of GM-CSF (7954-GM, R&D Systems) and 20 ng/mL each of LPS (L2630, Sigma-Aldrich), IFN-γ (285-IF, R&D Systems) and IL-6 (206-IL, R&D Systems).

2.8. TGFβ Treatment

PBMCs were cultured in a 12-well plate at a density of 5 × 105 per well at 37 °C in 5% CO2 atmosphere in RPMI1640 (Sigma-Aldrich) supplemented with 10% heat-inactivated FBS and penicillin/streptomycin, and were exposed to 0, 0.1, 1 or 10 ng/mL of TGFβ1 (7754-BH, R&D systems). Cells were harvested 48 h after treatment initiation and used in each experiment. To generate IL-2-stimulated T-cells, PBMCs were stimulated with complete medium alone or with human IL-2 (589102, BioLgend) at 25 ng/mL for 6 days, and the cells were then treated with 10 ng/mL of TGFβ1 in the presence of IL-2 for 48 h.

2.9. Statistical Analysis

Correlations were assessed using Spearman’s coefficient. The Mann–Whitney U test or paired Wilcoxon signed-rank test was used to determine differences between two variables. The Kruskal-Wallis test with Dunn’s post hoc test was used to compare multiple groups. All statistical analyses were conducted using GraphPad Prism v9.5.1 (GraphPad Software Inc., San Diego, CA, USA). All p-values were two-sided, and we considered p-values less than 0.05 to be statistically significant.

3. Results

3.1. Correlation between the Expression of HAVCR2, M2 Macrophage and TBSS in MSI and MSS CRC

In our recent study, which focused on MSI CRC, we developed the TGFβ-stromal signature (TBSS), by which to infer the levels of TGFβ-dependent stromal activation, based on six TGFβ-responsive, stroma-specific genes, consisting of VCAN, POSTN, FAP, CALD1, SERPINE1 and IGFBP7, and identified a CMS4-like, TGFβ-activated stromal subset within MSI CRCs characterized by high levels of M2 macrophage signature, ICI resistance signature and poor prognosis [15]. Notably, this stromal MSI subtype also consistently displayed the upregulation of HAVCR2 (encoding TIM-3), but not of the other immune checkpoint genes.
In the present study, we first examined whether the association of the TGFβ-rich TME with HAVCR2 expression and M2 macrophage signature is characteristic of MSI tumors or irrespective of MSI status. Using the TCGA RNA-seq dataset, although MSI tumors showed higher levels of HAVCR2 expression and M2 macrophage signature than those of MSS, we found that the levels of HAVCR2 expression, M2 macrophage signature and TBSS were significantly correlated with TGFB levels in the analysis of both MSI and MSS tumors (Figure 1A–C and Figure S1A–F). Therefore, we decided to analyze MSI and MSS CRCs together. We differentially examined expressed genes between TBSSLow and TBSSHigh tumors and found the striking upregulation of HAVCR2 (log2 fold change = 1.80, p = 2.76 × 10−46), along with marked enrichment of TGFB genes, stromal/mesenchymal/ECM genes and M2 macrophage marker genes, such as MSR1 (CD204), CD163 (CD163) and MRC1 (CD206), in TBSSHigh tumors (Figure 1D). It is worth noting that the expression of HAVCR2 had a much closer relationship with high TBSS, among the other co-inhibitory receptor genes, including TIGIT, CTLA4, PDCD1 (PD-1) and LAG3, prompting us to further explore the role of HAVCR2 in the TGFβ-rich TME in CRC.

3.2. The Expression of HAVCR2 and M2 Macrophage Signature in the TGFβ-Activated Stroma in CMS1 and CMS4 CRC

To validate the strong association of TGFB with HAVCR2, M2 macrophage and TBSS, which was found to be independent of MSI status, we utilized multiple independent cohorts (n = 1975, in total), including TCGA RNA-seq and five microarray datasets in conjunction with CMS classification and MSI status (Table S1). In all six cohorts, we confirmed significant positive correlations of TGFB levels with HAVCR2, M2 macrophage signature and TBSS (Figure 1E). These variables were also significantly correlated with one another in each of the cohorts (Figure S2). These data consistently indicated that the expression of HAVCR2 along with M2 macrophage signature, was strongly upregulated in the TGFβ-activated TME. When focusing on the CMS subtypes and MSI status, TGFB, HAVCR2, M2 macrophage and TBSS were together found to be significantly higher in both CMS1 (mainly MSI) and CMS4 tumors, compared to those of CMS2 and CMS3 in all cohorts (Figure 2). Collectively, although TGFB, HAVCR2, M2 macrophage and TBSS were robustly correlated with each other regardless of MSI status or CMS subtype, they were found to be strongly co-upregulated in MSI, CMS1 and CMS4 CRCs.
The transcriptomic data from bulk tumor samples (Figure 1 and Figure 2) consisted of mixed signals originating from both epithelial tumor cells and many types of stromal cells in the TME, in which most cells have the potential to respond to TGFβ, and the responses largely depend on cell types [8]. Therefore, we sought to determine whether the TGFB-dependency of HAVCR2, M2 macrophage and TBSS found in patient cohorts were stroma-specific and unrelated to the signals derived from the epithelial component. To this end, we further analyzed transcriptomic datasets of CRC bulk tissues, patient-derived organoids, patient-derived xenografts and cell lines. Organoids and cell lines lack a stromal component, whereas, in xenografts, cancer stroma is replaced by murine counterparts and is not or only partially detected by human-specific expression platforms. Based on pairwise analysis of patient-derived CRC organoids and their matched CRC bulk tissues (GSE171862), we found that the levels of TGFB were significantly correlated with HAVCR2, M2 macrophage and TBSS in bulk tumor tissues (Spearman’s r: 0.6716–0.8463, p < 0.0001) but not in organoids (Spearman’s r: −0.01699–0.3112) (Figure 3A–C). Correspondingly, in GSE100550, a significant positive correlation of TGFB with HAVCR2, M2 macrophage and TBSS was specifically found in bulk CRC tissues (Spearman’s r: 0.5792–0.7909, p < 0.0001), but not in CRC organoids, xenografts or cell lines (Spearman’s r: −0.2239–0.3333) (Figure 3D–F). These data clearly support the notion that the TGFB-dependent increase in HAVCR2, M2 macrophage and TBSS in CRC exclusively arises from the cancer stroma.

3.3. TIM-3+ Tumor-Infiltrating Immune Cells in TGFβ-Activated Stroma

We next conducted immunohistochemistry for TIM-3 and VCAN, a TGFβ-responsive stromal protein also known as versican, using surgically resected CRC tissues. Stromal VCAN immunostaining can indicate TGFβ-activated stroma, as it was exclusively stained in the cancer stroma, mainly expressed by CAFs in association with phospho-SMAD3 staining [15,35]. VCAN IHC scores and TIM-3+ tumor-infiltrating immune cells (TIICs) were evaluated using 45 CRC specimens. Representative images of VCAN and TIM-3 are shown in Figure 4A–F and Figure S3. Overall, TIM-3+ TIICs were often found in the tumor center, tumor invasive front and lymphoid structures in the TME, and TIM-3 was mainly expressed by macrophage-like TIICs (Figure S3). Moreover, immunofluorescent staining for TIM-3 and an M2 macrophage marker, CD163, in CRC tissues confirmed that TIM-3 exhibited a high degree of overlap with CD163+ macrophages (Figure 4G). We found that tumors with higher stromal VCAN scores were significantly associated with higher levels of TIM-3+ TIICs in the TME (Figure 4H, Spearman’s r = 0.584, p < 0.0001). Considering the heterogeneous patterns of VCAN and TIM-3 staining within tumors, we used 24 tumors showing heterogeneous VCAN staining, and the stromal areas of each tumor were further divided into those with strong VCAN staining (VCANHigh areas) and those with negative or moderate VCAN staining (VCANLow areas) (Figure S4). We found significantly higher levels of TIM-3+ TIICs in VCANHigh stromal areas than in VCANLow areas (Figure 4I, p < 0.0001).

3.4. M2 Macrophages Express TIM-3 and TGFβ Induces TIM-3 on Monocytes

We demonstrated that the close relationship between M2 macrophage signature and TGFB is consistent with the fact that TGFβ can contribute to inducing the M2-like macrophage phenotype, and that M2-like macrophages can be a major source of TGFβ in the TME [8,10,40,41]. Given the strong correlation between HAVCR2 expression and M2 macrophage signature at the transcriptional level, we investigated the surface expression of TIM-3 in M2-polarized macrophages. Human CD14+ monocytes, isolated from PBMCs, were differentiated into the M2 phenotype in vitro, showing high levels of cell surface TIM-3 together with an M2 marker, CD163, compared to M0 macrophages by flow cytometry (Figure 5A,B and Figure S5A,B). Although M1 macrophages were also generated from PBMCs, cell surface TIM-3 expression did not significantly differ between generated M0 and M1 macrophages (Figure 5C,D and Figure S5C,D).
Since we found the TGFB-dependent increase in HAVCR2 and M2 macrophage signature in transcriptomic datasets and dense infiltration of TIM-3+ TIICs in VCAN+ TGFβ-activated stroma by immunohistochemistry, we hypothesized that TGFβ may induce TIM-3 expression in immune cells, particularly in myeloid cells. To this end, we analyzed the expression of TIM-3 in monocytes, CD4+ T-cells and CD8+ T-cells obtained from PBMCs treated with varying concentrations of TGFβ in culture (Figure S6A). We found that TGFβ induced a significant dose-dependent increase in TIM-3 expression in monocytes (Figure 5E,F), but not in CD4+ T-cells or CD8+ T-cells (Figure S6B,C). Although we further utilized IL-2-stimulated T-cells expressing TIM-3, treatment with TGFβ did not significantly alter the levels of TIM-3 on the surface of CD4+ or CD8+ T-cells regardless of IL-2 stimulation (Figure S6D,E).

4. Discussion

In the present study, we found a strong correlation between the TGFβ-dependent active TME and HAVCR2 expression in CRC, and showed consistent results across multiple datasets, which consisted of approximately 2200 CRC tissues, organoids and xenografts. This study revealed that HAVCR2 was strikingly upregulated in a considerable fraction of CMS4 and CMS1/MSI CRCs, both of which were characterized by the TGFβ-activated TME. These correlations at transcriptional levels were further recapitulated by immunohistochemistry. We found that the levels of tumor-infiltrating TIM-3+ immune cells, mainly M2-like macrophages, were correlated with tumors exhibiting TGFβ-activated stroma, as detected by VCAN staining, and that TIM-3+ cells particularly accumulated in VCANHigh stromal areas. We also demonstrated a significant correlation of TGFB with the expression of HAVCR2 and M2 macrophage signature, and the upregulation of TIM-3 was observed on the surface of M2-polarized macrophages and TGFβ-stimulated monocytes. Therefore, in the TGFβ-rich TME, TIM-3 may be mainly expressed by M2-like macrophages at least in part induced by TGFβ. To our knowledge, this is the first study to address TIM-3 and M2 macrophages in the TGFβ-rich TME, in association with MSI status and CMS subtypes in CRC. Our results are in good agreement with a previous study on hepatocellular carcinoma by Yan et al., which reported that TGFβ not only induced TIM-3 on macrophages but also facilitated M2-like polarization in cultured macrophages [27]. Yan et al. also identified a SMAD-binding site in the 5′ region of TIM-3, and demonstrated that TGFβ enhanced the TIM-3 promotor activity examined by luciferase assays [27]. TIM-3 can promote M2 macrophage polarization through several mechanisms in CRC and HCC [21,27,29]. Collectively, several lines of evidence indicate a mechanistic link between the M2-like polarization of macrophages and increased TIM-3 expression in the TGFβ-rich TME. It is also worth noting that TIM-3, M2-like macrophages and the TGFβ-activated TME are known to have tumor-promoting and immunosuppressive functions, as well as poor prognostic impact in solid cancers, particularly in CRC [4,5,6,7,8,24,30,42].
TGFβ is a pleiotropic cytokine with complex and context-dependent functions in cancer. TGFβ is known to be a potent growth inhibitor of cells of epithelial origin, thus, it acts as a tumor-suppressor in pre-malignant cells in the early stages of tumorigenesis [10]. During CRC progression, TGFβ levels increase in the TME [6]. Although it is well documented that TGFβ promotes cancer cell-intrinsic effects on invasion and metastasis associated with EMT, cancer cells often acquire loss-of-function mutations in TGFβ signaling pathway genes to avoid the growth inhibitory effect of TGFβ [6,10]. Consequently, in advanced-stage tumors, TGFβ plays a key role in shaping the tumor-promoting TME and driving immune evasion, where CAFs are activated and TAMs are polarized toward the M2 phenotype with both cell types being the major sources of TGFβ [8,10,40,41]. TGFβ has also emerged as a negative predictor of ICI response. In clinical trials and animal models, TGFβ-activated CAFs and TGFβ-induced ECM genes are associated with a poor response to anti-PD-1/PD-L1 therapies in several solid cancers [8,16,17,18,19]. Moreover, several strategies for the dual inhibition of TGFβ and PD-1/PD-L1 signaling have been extensively studied in preclinical models, showing uniformly positive results across multiple types of solid cancers [17,20,43,44,45]. In mouse models of MSS CRC with CMS4 phenotype, TGFβ inhibition exerted synergic anti-tumor effects with anti-PD-L1 therapy [9], while we previously identified a CMS4-like subgroup within MSI CRCs displaying TGFβ-active stroma as well as an ICI resistance signature [15]. Therefore, the combined blockade of TGFβ and PD-1/PD-L1 may provide clinically meaningful responses in some patients with MSI and MSS CRC, possibly by making the immunosuppressive TME sensitive to ICIs. However, despite the success in preclinical models, clinical trials of anti-TGFβ combined with ICIs in several solid cancers, including CRC, have suggested that this approach may be limited, as the addition of TGFβ inhibitors often failed to show meaningful clinical benefits beyond ICIs alone in unselected patients [20]. In addition, the systemic adverse effects of inhibiting TGFβ signaling and its narrow therapeutic window require careful consideration [8,20]. Although some trials have already been terminated and many others are currently underway, it may be suggested that patient selection and additional therapeutic combinations are required to maximize the benefit of these strategies. For example, CMS4 CRCs (mostly MSS), as well as a CMS4-like subset of MSI and/or CMS1 CRCs, might be a candidate population that may be effectively treated with anti-TGFβ combined with anti-PD-1/PD-L1 therapies, as they have a worse prognosis and exhibit much higher levels of TGFβ-dependent stromal activation than those of CMS2 and CMS3. More importantly, TIM-3 can be an actionable target in combination with a PD-1/PD-L1 blockade for these patients. The dual inhibition of PD-1 and TIM-3 showed synergic effects in animal models [22,46], and adaptive resistance to the PD-1 blockade was associated with the upregulation of TIM-3 [23]. Anti-TIM-3 antibodies, mostly in combination with an anti-PD-1/PD-L1 blockade, are currently being investigated in clinical trials for the treatment of advanced cancers, and have shown encouraging early efficacy results with no serious adverse events or dose-limiting toxicities [21,30,47]. Concerning the adverse events observed in clinical trials of TGFβ inhibition [8,20], the blockade of TIM-3 may be a safe and promising alternative, at least in selected patients with CRC. Therefore, we suggest that, in future trials, dual TIM-3 and PD-1/PD-L1 inhibition should be considered for patients with CMS4 and MSI CRC, where TIM-3 is highly expressed in the TGFβ-rich TME.

5. Conclusions

In conclusion, using large-scale transcriptomic data analysis, followed by immunohistochemistry and in vitro experiments, we demonstrated a mechanistic link between TIM-3 expression and M2-like macrophages in the TGFβ-activated TME. We found that TGFB, HAVCR2, M2 macrophage and TBSS signatures were specifically derived from the stromal components and robustly correlated with each other, and that they were co-upregulated in CMS4, CMS1 and MSI CRCs. Since TGFβ and TIM-3 are implicated in immune evasion and poor response to PD-1/PD-L1 blockade, immunotherapeutic strategies targeting TIM-3 and/or TGFβ combined with PD-1/PD-L1 signaling may overcome the resistance to current ICIs, particularly for CMS4 and MSI CRCs characterized by high levels of TIM-3 and TGFβ.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15204943/s1, Figure S1: Levels of M2 macrophage signature and TGFβ-stromal signature in association with TGFB in MSI and MSS CRCs in TCGA dataset; Figure S2: Pairwise correlations between the expression of TGFB, TGFβ-stromal signature, M2 macrophage signature and HAVCR2, in six independent cohorts of CRC; Figure S3: Representative images of immunohistochemistry for TIM-3 in CRC tissues; Figure S4: Heterogenous VCAN and TIM-3 staining; Figure S5: Generation of M2 and M1 macrophages; Figure S6: CD14+ monocytes, CD3+CD4+ T-cells and CD3+CD8+ T-cells from PBMCs; Table S1: Datasets used in this study. References [3,31,32,33,34,48,49,50,51,52,53] are cited in the Supplementary Materials.

Author Contributions

Conceptualization, H.O., E.E. and K.K.; methodology, H.O. and S.N.; software, H.O.; validation, H.O., M.K. and S.N.; formal analysis, H.O. and S.N.; investigation, H.O., M.K., S.N., K.S., T.S., M.S. (Mei Sakuma) and S.F.; resources, W.S., M.S. (Motonobu Saito), Z.S. and T.M.; data curation, H.O.; writing—original draft preparation, H.O., M.K. and S.N.; writing—review and editing, H.O., K.M. and K.K.; visualization, H.O. and M.K.; supervision, H.O. and K.K.; project administration, H.O., S.N., K.M. and K.K.; funding acquisition, H.O., E.E. and W.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Japan Society for the Promotion of Science KAKENHI grant numbers, 22K08780, 21K16454 and 21K08802.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki, and approved by the Institutional Review Board of Fukushima Medical University (Reference No. 2289 and 2847).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The public datasets used in this study are available as described earlier. The data generated in this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Siegel, R.L.; Miller, K.D.; Fuchs, H.E.; Jemal, A. Cancer statistics, 2022. CA Cancer J. Clin. 2022, 72, 7–33. [Google Scholar] [CrossRef] [PubMed]
  2. Sveen, A.; Kopetz, S.; Lothe, R.A. Biomarker-guided therapy for colorectal cancer: Strength in complexity. Nat. Rev. Clin. Oncol. 2020, 17, 11–32. [Google Scholar] [CrossRef] [PubMed]
  3. Cancer Genome Atlas Network. Comprehensive molecular characterization of human colon and rectal cancer. Nature 2012, 487, 330–337. [Google Scholar] [CrossRef] [PubMed]
  4. Calon, A.; Lonardo, E.; Berenguer-Llergo, A.; Espinet, E.; Hernando-Momblona, X.; Iglesias, M.; Sevillano, M.; Palomo-Ponce, S.; Tauriello, D.V.; Byrom, D.; et al. Stromal gene expression defines poor-prognosis subtypes in colorectal cancer. Nat. Genet. 2015, 47, 320–329. [Google Scholar] [CrossRef]
  5. Isella, C.; Terrasi, A.; Bellomo, S.E.; Petti, C.; Galatola, G.; Muratore, A.; Mellano, A.; Senetta, R.; Cassenti, A.; Sonetto, C.; et al. Stromal contribution to the colorectal cancer transcriptome. Nat. Genet. 2015, 47, 312–319. [Google Scholar] [CrossRef]
  6. Calon, A.; Espinet, E.; Palomo-Ponce, S.; Tauriello, D.V.; Iglesias, M.; Céspedes, M.V.; Sevillano, M.; Nadal, C.; Jung, P.; Zhang, X.H.; et al. Dependency of colorectal cancer on a TGF-β-driven program in stromal cells for metastasis initiation. Cancer Cell 2012, 22, 571–584. [Google Scholar] [CrossRef]
  7. Guinney, J.; Dienstmann, R.; Wang, X.; de Reynies, A.; Schlicker, A.; Soneson, C.; Marisa, L.; Roepman, P.; Nyamundanda, G.; Angelino, P.; et al. The consensus molecular subtypes of colorectal cancer. Nat. Med. 2015, 21, 1350–1356. [Google Scholar] [CrossRef]
  8. Tauriello, D.V.F.; Sancho, E.; Batlle, E. Overcoming TGFbeta-mediated immune evasion in cancer. Nat. Rev. Cancer 2022, 22, 25–44. [Google Scholar] [CrossRef]
  9. Tauriello, D.V.F.; Palomo-Ponce, S.; Stork, D.; Berenguer-Llergo, A.; Badia-Ramentol, J.; Iglesias, M.; Sevillano, M.; Ibiza, S.; Canellas, A.; Hernando-Momblona, X.; et al. TGFbeta drives immune evasion in genetically reconstituted colon cancer metastasis. Nature 2018, 554, 538–543. [Google Scholar] [CrossRef]
  10. Batlle, E.; Massague, J. Transforming Growth Factor-beta Signaling in Immunity and Cancer. Immunity 2019, 50, 924–940. [Google Scholar] [CrossRef]
  11. Anderson, A.C.; Joller, N.; Kuchroo, V.K. Lag-3, Tim-3, and TIGIT: Co-inhibitory Receptors with Specialized Functions in Immune Regulation. Immunity 2016, 44, 989–1004. [Google Scholar] [CrossRef]
  12. Jin, Z.; Sinicrope, F.A. Mismatch Repair-Deficient Colorectal Cancer: Building on Checkpoint Blockade. J. Clin. Oncol. 2022, 40, 2735–2750. [Google Scholar] [CrossRef]
  13. Ganesh, K.; Stadler, Z.K.; Cercek, A.; Mendelsohn, R.B.; Shia, J.; Segal, N.H.; Diaz, L.A., Jr. Immunotherapy in colorectal cancer: Rationale, challenges and potential. Nat. Rev. Gastroenterol. Hepatol. 2019, 16, 361–375. [Google Scholar] [CrossRef] [PubMed]
  14. Le, D.T.; Uram, J.N.; Wang, H.; Bartlett, B.R.; Kemberling, H.; Eyring, A.D.; Skora, A.D.; Luber, B.S.; Azad, N.S.; Laheru, D.; et al. PD-1 Blockade in Tumors with Mismatch-Repair Deficiency. N. Engl. J. Med. 2015, 372, 2509–2520. [Google Scholar] [CrossRef] [PubMed]
  15. Endo, E.; Okayama, H.; Saito, K.; Nakajima, S.; Yamada, L.; Ujiie, D.; Kase, K.; Fujita, S.; Endo, H.; Sakamoto, W.; et al. A TGFbeta-Dependent Stromal Subset Underlies Immune Checkpoint Inhibitor Efficacy in DNA Mismatch Repair-Deficient/Microsatellite Instability-High Colorectal Cancer. Mol. Cancer Res. 2020, 18, 1402–1413. [Google Scholar] [CrossRef] [PubMed]
  16. Chakravarthy, A.; Khan, L.; Bensler, N.P.; Bose, P.; De Carvalho, D.D. TGF-beta-associated extracellular matrix genes link cancer-associated fibroblasts to immune evasion and immunotherapy failure. Nat. Commun. 2018, 9, 4692. [Google Scholar] [CrossRef] [PubMed]
  17. Mariathasan, S.; Turley, S.J.; Nickles, D.; Castiglioni, A.; Yuen, K.; Wang, Y.; Kadel Iii, E.E.; Koeppen, H.; Astarita, J.L.; Cubas, R.; et al. TGFbeta attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells. Nature 2018, 554, 544–548. [Google Scholar] [CrossRef] [PubMed]
  18. Grauel, A.L.; Nguyen, B.; Ruddy, D.; Laszewski, T.; Schwartz, S.; Chang, J.; Chen, J.; Piquet, M.; Pelletier, M.; Yan, Z.; et al. TGFbeta-blockade uncovers stromal plasticity in tumors by revealing the existence of a subset of interferon-licensed fibroblasts. Nat. Commun. 2020, 11, 6315. [Google Scholar] [CrossRef]
  19. Dominguez, C.X.; Muller, S.; Keerthivasan, S.; Koeppen, H.; Hung, J.; Gierke, S.; Breart, B.; Foreman, O.; Bainbridge, T.W.; Castiglioni, A.; et al. Single-Cell RNA Sequencing Reveals Stromal Evolution into LRRC15+ Myofibroblasts as a Determinant of Patient Response to Cancer Immunotherapy. Cancer Discov. 2020, 10, 232–253. [Google Scholar] [CrossRef]
  20. Metropulos, A.E.; Munshi, H.G.; Principe, D.R. The difficulty in translating the preclinical success of combined TGFbeta and immune checkpoint inhibition to clinical trial. EBioMedicine 2022, 86, 104380. [Google Scholar] [CrossRef]
  21. Wolf, Y.; Anderson, A.C.; Kuchroo, V.K. TIM3 comes of age as an inhibitory receptor. Nat. Rev. Immunol. 2020, 20, 173–185. [Google Scholar] [CrossRef] [PubMed]
  22. Sakuishi, K.; Apetoh, L.; Sullivan, J.M.; Blazar, B.R.; Kuchroo, V.K.; Anderson, A.C. Targeting Tim-3 and PD-1 pathways to reverse T cell exhaustion and restore anti-tumor immunity. J. Exp. Med. 2010, 207, 2187–2194. [Google Scholar] [CrossRef] [PubMed]
  23. Koyama, S.; Akbay, E.A.; Li, Y.Y.; Herter-Sprie, G.S.; Buczkowski, K.A.; Richards, W.G.; Gandhi, L.; Redig, A.J.; Rodig, S.J.; Asahina, H.; et al. Adaptive resistance to therapeutic PD-1 blockade is associated with upregulation of alternative immune checkpoints. Nat. Commun. 2016, 7, 10501. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, Y.; Cai, P.; Liang, T.; Wang, L.; Hu, L. TIM-3 is a potential prognostic marker for patients with solid tumors: A systematic review and meta-analysis. Oncotarget 2017, 8, 31705–31713. [Google Scholar] [CrossRef]
  25. Chen, K.; Gu, Y.; Cao, Y.; Fang, H.; Lv, K.; Liu, X.; He, X.; Wang, J.; Lin, C.; Liu, H.; et al. TIM3+ cells in gastric cancer: Clinical correlates and association with immune context. Br. J. Cancer 2022, 126, 100–108. [Google Scholar] [CrossRef]
  26. Dixon, K.O.; Tabaka, M.; Schramm, M.A.; Xiao, S.; Tang, R.; Dionne, D.; Anderson, A.C.; Rozenblatt-Rosen, O.; Regev, A.; Kuchroo, V.K. TIM-3 restrains anti-tumour immunity by regulating inflammasome activation. Nature 2021, 595, 101–106. [Google Scholar] [CrossRef]
  27. Yan, W.; Liu, X.; Ma, H.; Zhang, H.; Song, X.; Gao, L.; Liang, X.; Ma, C. Tim-3 fosters HCC development by enhancing TGF-beta-mediated alternative activation of macrophages. Gut 2015, 64, 1593–1604. [Google Scholar] [CrossRef]
  28. Komohara, Y.; Morita, T.; Annan, D.A.; Horlad, H.; Ohnishi, K.; Yamada, S.; Nakayama, T.; Kitada, S.; Suzu, S.; Kinoshita, I.; et al. The Coordinated Actions of TIM-3 on Cancer and Myeloid Cells in the Regulation of Tumorigenicity and Clinical Prognosis in Clear Cell Renal Cell Carcinomas. Cancer Immunol. Res. 2015, 3, 999–1007. [Google Scholar] [CrossRef]
  29. Jiang, X.; Zhou, T.; Xiao, Y.; Yu, J.; Dou, S.; Chen, G.; Wang, R.; Xiao, H.; Hou, C.; Wang, W.; et al. Tim-3 promotes tumor-promoting M2 macrophage polarization by binding to STAT1 and suppressing the STAT1-miR-155 signaling axis. Oncoimmunology 2016, 5, e1211219. [Google Scholar] [CrossRef]
  30. Zeidan, A.M.; Komrokji, R.S.; Brunner, A.M. TIM-3 pathway dysregulation and targeting in cancer. Expert Rev. Anticancer Ther. 2021, 21, 523–534. [Google Scholar] [CrossRef]
  31. Gao, J.; Aksoy, B.A.; Dogrusoz, U.; Dresdner, G.; Gross, B.; Sumer, S.O.; Sun, Y.; Jacobsen, A.; Sinha, R.; Larsson, E.; et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci. Signal 2013, 6, pl1. [Google Scholar] [CrossRef]
  32. Vasaikar, S.V.; Straub, P.; Wang, J.; Zhang, B. LinkedOmics: Analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018, 46, D956–D963. [Google Scholar] [CrossRef] [PubMed]
  33. Cho, E.J.; Kim, M.; Jo, D.; Kim, J.; Oh, J.H.; Chung, H.C.; Lee, S.H.; Kim, D.; Chun, S.M.; Kim, J.; et al. Immuno-genomic classification of colorectal cancer organoids reveals cancer cells with intrinsic immunogenic properties associated with patient survival. J. Exp. Clin. Cancer Res. 2021, 40, 230. [Google Scholar] [CrossRef]
  34. Linnekamp, J.F.; Hooff, S.R.V.; Prasetyanti, P.R.; Kandimalla, R.; Buikhuisen, J.Y.; Fessler, E.; Ramesh, P.; Lee, K.; Bochove, G.G.W.; de Jong, J.H.; et al. Consensus molecular subtypes of colorectal cancer are recapitulated in in vitro and in vivo models. Cell Death Differ. 2018, 25, 616–633. [Google Scholar] [CrossRef]
  35. Chida, S.; Okayama, H.; Noda, M.; Saito, K.; Nakajima, T.; Aoto, K.; Hayase, S.; Momma, T.; Ohki, S.; Kono, K.; et al. Stromal VCAN expression as a potential prognostic biomarker for disease recurrence in stage II-III colon cancer. Carcinogenesis 2016, 37, 878–887. [Google Scholar] [CrossRef]
  36. Yoshihara, K.; Shahmoradgoli, M.; Martinez, E.; Vegesna, R.; Kim, H.; Torres-Garcia, W.; Trevino, V.; Shen, H.; Laird, P.W.; Levine, D.A.; et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat. Commun. 2013, 4, 2612. [Google Scholar] [CrossRef]
  37. Mak, M.P.; Tong, P.; Diao, L.; Cardnell, R.J.; Gibbons, D.L.; William, W.N.; Skoulidis, F.; Parra, E.R.; Rodriguez-Canales, J.; Wistuba, I.I.; et al. A Patient-Derived, Pan-Cancer EMT Signature Identifies Global Molecular Alterations and Immune Target Enrichment Following Epithelial-to-Mesenchymal Transition. Clin. Cancer Res. 2016, 22, 609–620. [Google Scholar] [CrossRef] [PubMed]
  38. Zarif, J.C.; Hernandez, J.R.; Verdone, J.E.; Campbell, S.P.; Drake, C.G.; Pienta, K.J. A phased strategy to differentiate human CD14+monocytes into classically and alternatively activated macrophages and dendritic cells. Biotechniques 2016, 61, 33–41. [Google Scholar] [CrossRef]
  39. Min, A.K.T.; Mimura, K.; Nakajima, S.; Okayama, H.; Saito, K.; Sakamoto, W.; Fujita, S.; Endo, H.; Saito, M.; Saze, Z.; et al. Therapeutic potential of anti-VEGF receptor 2 therapy targeting for M2-tumor-associated macrophages in colorectal cancer. Cancer Immunol. Immunother. 2021, 70, 289–298. [Google Scholar] [CrossRef] [PubMed]
  40. Aras, S.; Zaidi, M.R. TAMeless traitors: Macrophages in cancer progression and metastasis. Br. J. Cancer 2017, 117, 1583–1591. [Google Scholar] [CrossRef]
  41. Cassetta, L.; Pollard, J.W. Targeting macrophages: Therapeutic approaches in cancer. Nat. Rev. Drug Discov. 2018, 17, 887–904. [Google Scholar] [CrossRef]
  42. Vayrynen, J.P.; Haruki, K.; Lau, M.C.; Vayrynen, S.A.; Zhong, R.; Dias Costa, A.; Borowsky, J.; Zhao, M.; Fujiyoshi, K.; Arima, K.; et al. The Prognostic Role of Macrophage Polarization in the Colorectal Cancer Microenvironment. Cancer Immunol. Res. 2021, 9, 8–19. [Google Scholar] [CrossRef]
  43. Martin, C.J.; Datta, A.; Littlefield, C.; Kalra, A.; Chapron, C.; Wawersik, S.; Dagbay, K.B.; Brueckner, C.T.; Nikiforov, A.; Danehy, F.T., Jr.; et al. Selective inhibition of TGFbeta1 activation overcomes primary resistance to checkpoint blockade therapy by altering tumor immune landscape. Sci. Transl. Med. 2020, 12, eaay8456. [Google Scholar] [CrossRef] [PubMed]
  44. Ravi, R.; Noonan, K.A.; Pham, V.; Bedi, R.; Zhavoronkov, A.; Ozerov, I.V.; Makarev, E.; Artemov, V.A.; Wysocki, P.T.; Mehra, R.; et al. Bifunctional immune checkpoint-targeted antibody-ligand traps that simultaneously disable TGFbeta enhance the efficacy of cancer immunotherapy. Nat. Commun. 2018, 9, 741. [Google Scholar] [CrossRef]
  45. Lan, Y.; Zhang, D.; Xu, C.; Hance, K.W.; Marelli, B.; Qi, J.; Yu, H.; Qin, G.; Sircar, A.; Hernandez, V.M.; et al. Enhanced preclinical antitumor activity of M7824, a bifunctional fusion protein simultaneously targeting PD-L1 and TGF-beta. Sci. Transl. Med. 2018, 10, eaan5488. [Google Scholar] [CrossRef] [PubMed]
  46. Ngiow, S.F.; von Scheidt, B.; Akiba, H.; Yagita, H.; Teng, M.W.; Smyth, M.J. Anti-TIM3 antibody promotes T cell IFN-gamma-mediated antitumor immunity and suppresses established tumors. Cancer Res. 2011, 71, 3540–3551. [Google Scholar] [CrossRef]
  47. Curigliano, G.; Gelderblom, H.; Mach, N.; Doi, T.; Tai, D.; Forde, P.M.; Sarantopoulos, J.; Bedard, P.L.; Lin, C.C.; Hodi, F.S.; et al. Phase I/Ib Clinical Trial of Sabatolimab, an Anti-TIM-3 Antibody, Alone and in Combination with Spartalizumab, an Anti-PD-1 Antibody, in Advanced Solid Tumors. Clin. Cancer Res. 2021, 27, 3620–3629. [Google Scholar] [CrossRef] [PubMed]
  48. Marisa, L.; de Reyniès, A.; Duval, A.; Selves, J.; Gaub, M.P.; Vescovo, L. Gene expression classification of colon cancer into molecular subtypes: Characterization, validation, and prognostic value. PLoS Med 2013, 10, e1001453. [Google Scholar] [CrossRef] [PubMed]
  49. de Sousa EMelo, F.; Colak, S.; Buikhuisen, J.; Koster, J.; Cameron, K.; de Jong, J.H. Methylation of cancer-stem-cell-associated Wnt target genes predicts poor prognosis in colorectal cancer patients. Cell Stem Cell 2011, 9, 476–485. [Google Scholar] [CrossRef]
  50. Jorissen, R.N.; Gibbs, P.; Christie, M.; Prakash, S.; Lipton, L.; Desai, J. Metastasis-Associated Gene Expression Changes Predict Poor Outcomes in Patients with Dukes Stage B and C Colorectal Cancer. Clin. Cancer Res. 2009, 15, 7642–7651. [Google Scholar] [CrossRef]
  51. Laibe, S.; Lagarde, A.; Ferrari, A.; Monges, G.; Birnbaum, D.; Olschwang, S. A seven-gene signature aggregates a subgroup of stage II colon cancers with stage III. OMICS 2012, 16, 560–565. [Google Scholar] [CrossRef] [PubMed]
  52. Guinney, J.; Ferte, C.; Dry, J.; McEwen, R.; Manceau, G.; Kao, K.J. Modeling RAS phenotype in colorectal cancer uncovers novel molecular traits of RAS dependency and improves prediction of response to targeted agents in patients. Clin. Cancer Res. 2014, 20, 265–272. [Google Scholar] [CrossRef] [PubMed]
  53. Lal, N.; White, B.S.; Goussous, G.; Pickles, O.; Mason, M.J.; Beggs, A.D. KRAS Mutation and Consensus Molecular Subtypes 2 and 3 Are Independently Associated with Reduced Immune Infiltration and Reactivity in Colorectal Cancer. Clin. Cancer Res. 2018, 24, 224–233. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Association of TGFB with HAVCR2, M2 macrophage signature and TGFβ-stromal signature (TBSS) in multiple RNA-seq or microarray cohorts of CRC in conjunction with MSI status and CMS subtypes. (A,B) Correlations between the expression of TGFB and HAVCR2 in MSI (A) and MSS (B) CRC. Spearman’s coefficients are indicated (p < 0.0001). (C) HAVCR2 expression in MSI and MSS CRC. (D) Volcano plot displaying differentially expressed genes between TBSSLow and TBSSHigh tumors. The x-axis shows the log2-fold change, and the y-axis displays the -log10 of the p-values. (E) Correlations of TGFB with HAVCR2 (red), M2 macrophage (yellow) and TBSS (green) in six independent cohorts of CRC. Spearman’s coefficients are shown as colors corresponding to HAVCR2 (red), M2 macrophage (yellow) and TBSS (green) for each dataset. p-values for all correlations < 0.0001.
Figure 1. Association of TGFB with HAVCR2, M2 macrophage signature and TGFβ-stromal signature (TBSS) in multiple RNA-seq or microarray cohorts of CRC in conjunction with MSI status and CMS subtypes. (A,B) Correlations between the expression of TGFB and HAVCR2 in MSI (A) and MSS (B) CRC. Spearman’s coefficients are indicated (p < 0.0001). (C) HAVCR2 expression in MSI and MSS CRC. (D) Volcano plot displaying differentially expressed genes between TBSSLow and TBSSHigh tumors. The x-axis shows the log2-fold change, and the y-axis displays the -log10 of the p-values. (E) Correlations of TGFB with HAVCR2 (red), M2 macrophage (yellow) and TBSS (green) in six independent cohorts of CRC. Spearman’s coefficients are shown as colors corresponding to HAVCR2 (red), M2 macrophage (yellow) and TBSS (green) for each dataset. p-values for all correlations < 0.0001.
Cancers 15 04943 g001
Figure 2. Association between CMS subtypes and the levels of TGFB, TGFβ-stromal signature (TBSS), M2 macrophage signature and HAVCR2 expression in six independent cohorts of CRC, including TCGA, GSE39582, GSE33113, GSE14333, GSE37892 and KFSYSCC. MSI status was available for the former three cohorts; red dots represent MSI tumors. Kruskal-Wallis test, followed by Dunn’s post hoc test to compare each CMS subtype with CMS2. **** p < 0.0001, *** p < 0.001, ** p < 0.01, * p < 0.05, ns p > 0.05.
Figure 2. Association between CMS subtypes and the levels of TGFB, TGFβ-stromal signature (TBSS), M2 macrophage signature and HAVCR2 expression in six independent cohorts of CRC, including TCGA, GSE39582, GSE33113, GSE14333, GSE37892 and KFSYSCC. MSI status was available for the former three cohorts; red dots represent MSI tumors. Kruskal-Wallis test, followed by Dunn’s post hoc test to compare each CMS subtype with CMS2. **** p < 0.0001, *** p < 0.001, ** p < 0.01, * p < 0.05, ns p > 0.05.
Cancers 15 04943 g002
Figure 3. Stroma-specific association of TGFB with TGFβ-stromal signature (TBSS), M2 macrophage signature and HAVCR2. (AC) Primary CRC tissues (orange) and matched organoids (purple) in an RNA-seq dataset (GSE171682). (DF) CRC tissues (orange), organoids (purple), xenografts (blue) and cell lines (black) in a microarray dataset (GSE100550). Spearman’s coefficients are indicated as colored numbers corresponding to tissues (orange), organoids (purple), xenografts (blue) and cell lines (black) for each dataset.
Figure 3. Stroma-specific association of TGFB with TGFβ-stromal signature (TBSS), M2 macrophage signature and HAVCR2. (AC) Primary CRC tissues (orange) and matched organoids (purple) in an RNA-seq dataset (GSE171682). (DF) CRC tissues (orange), organoids (purple), xenografts (blue) and cell lines (black) in a microarray dataset (GSE100550). Spearman’s coefficients are indicated as colored numbers corresponding to tissues (orange), organoids (purple), xenografts (blue) and cell lines (black) for each dataset.
Cancers 15 04943 g003
Figure 4. Immunohistochemistry for a TGFβ-responsive stromal protein, VCAN (A,B), and TIM-3+ tumor-infiltrating immune cells (TIICs) (CF) in CRC tissues. (A,C,E) Serial sections from a representative case of CRC showing negative stromal VCAN staining (A) with low infiltration of TIM-3+ TIICs (C,E). (B,D,F) Serial sections from a representative case of CRC showing strong stromal VCAN staining (B) with dense TIM-3+ TIICs (D,F). (G) Representative immunofluorescence of CD163 (green), TIM-3 (red) and DAPI (blue) in CRC tissues. Scale bar, 100μm. (H) Association between the levels of TIM-3+ TIICs and stromal VCAN in CRC (n = 45) by Spearman’s correlation. (I) Association between TIM-3+ TIICs and heterogenous VCAN expression patterns (VCANLow areas and VCANHigh areas) within the tumor (n = 24). Pairwise significance was tested using the Wilcoxon signed rank test.
Figure 4. Immunohistochemistry for a TGFβ-responsive stromal protein, VCAN (A,B), and TIM-3+ tumor-infiltrating immune cells (TIICs) (CF) in CRC tissues. (A,C,E) Serial sections from a representative case of CRC showing negative stromal VCAN staining (A) with low infiltration of TIM-3+ TIICs (C,E). (B,D,F) Serial sections from a representative case of CRC showing strong stromal VCAN staining (B) with dense TIM-3+ TIICs (D,F). (G) Representative immunofluorescence of CD163 (green), TIM-3 (red) and DAPI (blue) in CRC tissues. Scale bar, 100μm. (H) Association between the levels of TIM-3+ TIICs and stromal VCAN in CRC (n = 45) by Spearman’s correlation. (I) Association between TIM-3+ TIICs and heterogenous VCAN expression patterns (VCANLow areas and VCANHigh areas) within the tumor (n = 24). Pairwise significance was tested using the Wilcoxon signed rank test.
Cancers 15 04943 g004
Figure 5. TIM-3 expression on M2 and M1-polarized macrophages and TGFβ-treated monocytes. (A,B) TIM-3 expression on M0 and M2-polarized macrophages. Representative histograms for TIM-3 are shown (A), and the expression levels of TIM-3 on M2 macrophages present as the mean ± s.e.m. relative to those of M0 macrophages (n = 4) (B). Mann–Whitney U test. * p < 0.05. (C,D) TIM-3 expression on M0 and M1-polarized macrophages. Representative histograms for TIM-3 (C), and the levels of TIM-3 on M1 macrophages as the mean ± s.e.m. relative to those of M0 macrophages (n = 3) (D). n.s. p > 0.05. (E,F) Representative histograms displaying TIM-3 expression on monocytes treated with various concentrations of TGFβ (E), and TIM-3 expression levels on TGFβ-treated monocytes present as the mean ± s.e.m. relative to those of controls (n = 10) (F). Kruskal-Wallis test with Dunn’s post hoc test. *** p < 0.001.
Figure 5. TIM-3 expression on M2 and M1-polarized macrophages and TGFβ-treated monocytes. (A,B) TIM-3 expression on M0 and M2-polarized macrophages. Representative histograms for TIM-3 are shown (A), and the expression levels of TIM-3 on M2 macrophages present as the mean ± s.e.m. relative to those of M0 macrophages (n = 4) (B). Mann–Whitney U test. * p < 0.05. (C,D) TIM-3 expression on M0 and M1-polarized macrophages. Representative histograms for TIM-3 (C), and the levels of TIM-3 on M1 macrophages as the mean ± s.e.m. relative to those of M0 macrophages (n = 3) (D). n.s. p > 0.05. (E,F) Representative histograms displaying TIM-3 expression on monocytes treated with various concentrations of TGFβ (E), and TIM-3 expression levels on TGFβ-treated monocytes present as the mean ± s.e.m. relative to those of controls (n = 10) (F). Kruskal-Wallis test with Dunn’s post hoc test. *** p < 0.001.
Cancers 15 04943 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Katagata, M.; Okayama, H.; Nakajima, S.; Saito, K.; Sato, T.; Sakuma, M.; Fukai, S.; Endo, E.; Sakamoto, W.; Saito, M.; et al. TIM-3 Expression and M2 Polarization of Macrophages in the TGFβ-Activated Tumor Microenvironment in Colorectal Cancer. Cancers 2023, 15, 4943. https://doi.org/10.3390/cancers15204943

AMA Style

Katagata M, Okayama H, Nakajima S, Saito K, Sato T, Sakuma M, Fukai S, Endo E, Sakamoto W, Saito M, et al. TIM-3 Expression and M2 Polarization of Macrophages in the TGFβ-Activated Tumor Microenvironment in Colorectal Cancer. Cancers. 2023; 15(20):4943. https://doi.org/10.3390/cancers15204943

Chicago/Turabian Style

Katagata, Masanori, Hirokazu Okayama, Shotaro Nakajima, Katsuharu Saito, Takahiro Sato, Mei Sakuma, Satoshi Fukai, Eisei Endo, Wataru Sakamoto, Motonobu Saito, and et al. 2023. "TIM-3 Expression and M2 Polarization of Macrophages in the TGFβ-Activated Tumor Microenvironment in Colorectal Cancer" Cancers 15, no. 20: 4943. https://doi.org/10.3390/cancers15204943

APA Style

Katagata, M., Okayama, H., Nakajima, S., Saito, K., Sato, T., Sakuma, M., Fukai, S., Endo, E., Sakamoto, W., Saito, M., Saze, Z., Momma, T., Mimura, K., & Kono, K. (2023). TIM-3 Expression and M2 Polarization of Macrophages in the TGFβ-Activated Tumor Microenvironment in Colorectal Cancer. Cancers, 15(20), 4943. https://doi.org/10.3390/cancers15204943

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop