Abstract
The human cutaneous metastatic melanoma is the deadliest skin cancer. Partial, or less frequently complete spontaneous regressions could be observed, mainly mediated by T cells. Nevertheless, the underlying mechanisms are not fully unraveled. We investigated the first events of the immune response related to cancer regression in Melanoma-bearing Libechov Minipigs (MeLiM), a unique swine model of cutaneous melanoma that regresses spontaneously. Using a multiparameter flow cytometry strategy and integrating new clinical and histological criteria of the regression, we show that T cells and B cells are present only in the late stages, arguing against their role in the initial destruction of malignant cells. NK cells infiltrate the tumors before T cells and therefore might be involved in the induction of the regression process. Myeloid cells represent the main immune population within the tumor microenvironment regardless of the regression stage. Among those, MHCII+ CD163− macrophages that differ phenotypically and functionally compared to other tumor-associated macrophages, increase in number together with the first signs of regression suggesting their crucial contribution to initiating the regression process. Our study supports the importance of macrophage reprogramming in humans to improve current immunotherapy for metastatic melanoma.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
Introduction
The cutaneous melanoma is the deadliest skin cancer. Its incidence has steadily increased since the 1970s with a predominance in fair-skinned populations from Australia, New Zealand, North America and northern Europe [1]. The high mortality rate caused by the disease is due to its high rate of metastasis. Nevertheless, spontaneous regressions, partial or less frequently complete, are observed macroscopically and histologically with an overall incidence ranging from 10 to 35% and even reaching 58% in primary melanoma with a Breslow thickness below 1 mm [2, 3]. The progressive disappearance of malignant cells is associated with fibrosis, presence of lymphocytes and melanophages and varying degrees of neovascularization, depigmentation around the tumors, and inflammation. Although the mechanisms of spontaneous regression are not fully unraveled, the host immune system plays a key role, attested to by the early activation of immunity against primary cutaneous melanoma leading to lymphocytic and histiocytic infiltration in regressing lesions [4, 5]. The immune reactions toward “self” cancer antigens may lead to the destruction of normal melanocytes and subsequently depigmentation occasionally occurring in patients. Such auto-immunity has even been associated with a favorable melanoma prognosis [3, 6, 7].
In-depth understanding of new immune parameters could contribute to guide the development of novel cancer treatments [8, 9]. In particular, the identification of immune cells related to the spontaneous tumor regression is of great interest but remains difficult in humans due to its low prevalence and the limited access to lesions with active regression. In addition, the extent of cancer regression is not objectively defined by a universal scheme, leading to an unreliable reporting of regression that may explain the conflicting data regarding on prognostic importance of regression [10]. Melanoma-bearing Libechov Minipigs (MeLiM) spontaneously develop multifocal primary cutaneous melanomas around birth, with clinical and histopathological features comparable to human counterparts [11, 12]. 80% of animals exhibit multiple lesions from benign to highly invasive, eventually leading to lymph nodes and visceral metastasis. A regression occurs without any treatment 2 to 4 months after birth, corresponding to a complete disappearance of primary tumors and metastases. Overall, only a 4% mortality rate is observed, likely due to metastatic complications appearing before the regression onset. Our team has already described clinically and histologically the regression process often accompanied by a local or systemic depigmentation of hair, skin and eyes [12, 13] and reported its association with the modulation of immune-related genes with humoral, monocyte/macrophage-like, and T/NK signatures [14]. In the MeLiM model, we have recently evidenced changes in the circulating immune cell composition due to melanoma occurrence: young melanoma-bearing piglets harboring higher proportions of NK cells, CD4+ and CD4+ CD8α+ T cells, and CD21− B cells among B cells consistent with the immune-mediated spontaneous regression [15]. Others have shown the involvement of fibronectin and tenascin C in forming fibrous tissue during the spontaneous tumor regression [16].
Here, we performed a longitudinal analysis of melanoma lesions from a large cohort of MeLiM pigs, from the initiation to latest stages of tumor regression. Using our recent multiparametric flow cytometry strategy [15] and integrating new clinical and histological criteria of the regression, we studied the tumor microenvironment and identified the different subsets of T, B, NK cells and myeloid cells infiltrating melanoma lesions. We also evidenced a macrophage subset exhibiting very different functions compared to other tumor-associated macrophages (TAMs) and that increases with the first signs of regression, supporting their contribution to the initiation of the regression process.
Methods
Experimental pigs
MeLiM pigs (14 males and 20 females) from 13 different litters were examined for the presence of cutaneous melanoma every 1 to 2 weeks until 3 months of age and then every 2 to 3 weeks until 5 months. All animals were genotyped for a specific haplotype segregating in MeLiM, responsible for epitope deficiency to the anti-CD4 antibody [13]. Clinical examination consisted in identifying new pigmented melanocytic lesions, monitoring the evolution of previously identified lesions and examining the presence of palpable lymphadenopathies (supplementary Tables 1 and 2). Some tumors were excised surgically under general anesthesia and either stored in liquid nitrogen for protein extraction, in 10% buffered formalin for histological analyses, or in PBS containing 5% of Vetedine solution (Vetoquinol), 200 units/mL of penicillin, 200 µg/mL of streptomycin and 0.2% of Fungizone (all from Gibco) for further processing of enzymatic digestion.
Histopathology
Hematoxylin–eosin-saffron—stained paraffin-embedded sections were evaluated histologically according to the human classification. Regression was characterized histologically by the presence of a zone of dermal fibrosis and if so, other criteria were examined in the regressive zone (supplementary Table 2).
Immunohistochemical staining and image analysis
Immunohistochemical stains for CD3 (mouse IgG1 anti-pig CD3, clone PPT3, 2.5 µg/mL, SouthernBiotech) were performed on 5-μm thick paraffin-embedded tissue sections using standard techniques after unmasking at 90 °C for 1h, ImmPRESS® HRP Anti-Mouse IgG (Peroxidase) kit (LSBio) and HRP green and PAS hematoxylin counterstain. Images were acquired with the Pannoramic SCAN digital slide scanner (objective magnification X40). Zones were defined manually and CD3 expression was quantified in each zone using with the CaseViewer and QuantCenter softwares (all from 3DHISTECH) and expressed as HRP green colored pixels relatively to 1000 total pixels. Partition of connective tissue vs tumoral cells was performed on Trichrome- or HES-stained sections using the same softwares.
Tumoral protein extraction and cytokine assays
Proteins were extracted from frozen tumors in protein lysis buffer (20mM Tris HCl pH7.5, 150 mM NaCl, 1mM EDTA, 1% NP-40, 0.1% SDS, 0.5% sodium deoxycholate) containing protease and phosphatase inhibitor (Pierce). Protein concentrations were measured with the BCA protein assay kit (Pierce) and cytokines were assessed on 12.5 µg of total protein. IFNα, IFNγ, IL-2, IL-4, IL-6, IL-8, IL-10, IL-17, IL-1β and TNFα were measured by a Cytokine Bead Assay as previously described [17]. CCL2 and TGFβ were assessed using Swine CCL2 VetSet ELISA Development Kit (KingFisher Biotech) and Mouse/Rat/Porcine/Canine TGF-beta 1 Quantikine ELISA Kit (R&D systems) following the manufacturers’ instructions (supplementary Table 3).
Enzymatic digestion of lesions and identification of cells
Epidermis and hypodermis were carefully removed and tumor tissues were weighted. Single cell suspensions were obtained by an enzymatic treatment for 1 h at 37 °C under agitation with collagenase B at 4 mg/mL and DNase I at 0.1 mg/mL (Roche) in DMEM containing 100 units/mL penicillin, 100 µg/mL streptomycin, 2 mM l-Glutamine, 0.5 mM EDTA and 2% FBS (all from Gibco). Extensive washings were performed in the same medium without enzymes to remove released melanin.
Absolute number and viability of cells were determined with the ViaCount Assay performed on easyCyte 6HT-2L Guava flow cytometer (Millipore) following manufacturer’s instructions. Cells were stained as previously described for blood cell analysis [15]. Briefly, three combinations of antibodies, listed in supplementary Table 4, with the aqua LIVE/DEAD® Fixable Dead Cell Stain Kit (Thermo Fisher Scientific) for viability, were used to identify lymphoid cells (combinations A: CD45, CD3, CD8α, CD4, γδTCR and CD16 and B: CD45, MHC II, CD21 and CD79a) and myeloid cells (combination C: CD45, MHC II, PG68A, CD163, CD172a, CD14 and CADM1). Five million cells were processed in each combination and finally fixed in BD CellFIX solution before analysis on a BD LSR Fortessa cytometer (BD Biosciences). Data were analyzed using FlowJo V10 software (supplementary Table 5).
T and NK cell functionality
After enzymatic digestion of tumoral lesions, five million cells were incubated in 1 mL of RPMI medium containing 5% FBS, 100 units/mL penicillin, 100 µg/mL streptomycin, 2 mM l-Glutamine, 10 mM Hepes (all from Gibco) with or without PMA (50 ng/mL, from Sigma-Aldrich) and ionomycin (1 µg/mL, from Sigma-Aldrich) for 3 h at 37 °C in the presence of monensin (1/1000, eBioscience) and brefeldin A (1/1000, Invitrogen). Cells were then stained for viability with the Near IR LIVE/DEAD® Fixable Dead Cell Stain Kit (Thermo Fisher Scientific) and stained with a combination of antibodies targeting CD45, CD3, CD8α and CD16 (combination D in supplementary Table 4). Detection of IFNγ was done by intracellular staining using Foxp3/Transcription Factor Staining Buffer Set (eBioscience) as previously described [15]. IFNγ intracellular staining specificity was determined by incubating cells with labeled isotype control. Cells were then fixed in BD CellFIX solution before analysis on a BD LSR Fortessa cytometer (BD Biosciences). The data analysis was performed using FlowJo V10 software.
Macrophage subpopulation cell sorting and qRT-PCR
After enzymatic digestion of tumoral lesions, cells were stained with the combination for myeloid cells (combination C in supplementary Table 4: CD45, MHC II, PG68A, CD163, CD172a, CD14 and CADM1) but not fixed. Macrophages subpopulations were sorted using a FACS Aria III cell sorter (BD Biosciences; 100 µm nozzle). RNA from sorted cells were extracted using the Arcturus PicoPure RNA Isolation kit (ThermoFisher Scientific) according to manufacturer’s instructions. Their concentrations were measured with a NanoDrop 2000 spectrophotometer (16.9 ± 13.8 µg were obtained per sample) and RNA integrity was assessed by an Agilent 2100 Bioanalyzer using Eukaryote total RNA 6000 Nano Kit (RIN obtained were 8.7 ± 0.6, ranging from 7.2 to 9.7). RT-PCR was performed from 10 to 100 ng of RNA using TaqMan Reverse Transcription Reagents (Applied Biosystems, ThermoFisher Scientific). Relative mRNA expression was evaluated by qPCR using either the Fast SYBR™ Green Master Mix or TaqMan Fast Advanced Master Mix on the QuantStudio 12K Flex system (all from Applied Biosystems) with the primers described in supplementary Table 6. A Ct value equal to 40 was set for non-detected genes. Gene expressions were normalized using the ΔCt method by subtracting the Ct obtained for the gene of interest by the geometric mean of the Ct obtained for Ribosomal protein S24 (RPS24) and Ribosomal protein L32 (RPL32) reference genes. Relative gene expression (RQ) were calculated by subtracting the ΔCt(sample) from the ΔCt(mean) across all samples and the 2−ΔΔCt method.
Statistics
Data analyses were performed using R (v3.6.1). For statistical inference, different transformations have been applied on data: a square root transformation for percentages of connective tissue, a log base 10 + 1 transformation for cytokine concentration and a log base 2 transformation for qPCR relative gene expression. Factorial Analysis of Mixed Data (FAMD) and Principal Component Analysis (PCA) were performed using FactoMineR (v2.1) package and visualized with the factoextra R package (v1.0.6). Statistical description of groups by the quantitative and qualitative variables were performed using v test obtained using the catdes function. Eventually, the imputePCA function of missMDA (v1.16) package was used to handle missing values. Heatmaps were built using ComplexHeatmap (v2.0.0) package.
Effect of sex and regression stages was evaluated on the different proportions of cell populations among immune cells in the tumors using a linear mixed model using the lm function with sex (two levels) and regression stages (5 levels) as fixed effects. p values for fixed effects were obtained using the lmerTest package (v3.1-1) by type III ANOVA tables with Satterthwaite’s approximation to degrees of freedom using the ANOVA function. Sex effect was not significant. Pairwise comparisons with Tukey’s adjustment were performed to assess the differences within regression stages using the emmeans package (v1.5.2-1). A significance level of 0.05 was applied. All significant F and p-values are reported in the supplementary Tables 7–9 as mentioned in figure legends. Results in the text are expressed as mean ± SD. Other statistical tests used are mentioned in the appropriate figure legend.
Results
Clinical and histological features of spontaneous melanoma regression in MeLiM pigs
A total of 92 cutaneous pigmented lesions from 34 pigs aged 14 to 154 days were analyzed. Supplementary Tables 1 and 2 present the pig’s individual clinical data with the histological and clinical parameters, respectively. Lesions corresponded to invasive nodular melanoma equivalent to Clark’s level IV or V. The tumor area consisted of highly pigmented cells. Regression was characterized histologically by the presence of a zone of dermal fibrosis. Melanoma without any microscopic regression, collected from animals aged 14 to 48 days, were classified as R0. Melanoma from age-matched animals with a partial zone of regression were classified as R1. The stage of regression was then determined kinetically: R2 from 51 to 69 days, R3 from 72 to 94 days and R4 from 100 to 154 days. Based on the melanoma clinical data and on the histological features of tumoral and regressive zones (Table 1), FAMD allowed to identify many characteristics of progressive R0 lesions and R1 lesions in partial regression (Fig. 1a). Indeed, in animals belonging to the group R1, lesions size levels off, adenomegalies are more frequent, and lesions display more connective tissue and clustered melanophages corresponding to dense nests of heavily pigmented melanin-laden histiocytes. Both R0 and R1 lesions are frequently ulcerated but R0 lesions more intensely.
Separation of lesions in regression from R1 to R4 groups is more continuous (Fig. 1b and Table 1). While the regression occurs, lesion size stabilizes and then decreases, Clark’s level and ulceration reduce, and tumors become grayish. Histologically, progressively with age, cutaneous melanoma becomes flat and the tumoral zone decreases as the lymphoid infiltration and proportion of connective tissue increase. A grenz zone could be observed at the last time-point. Regression first occurs partially in the dermis and then extends to the epidermis. Lymphoid infiltration is more frequent in R3 and R4 stages of regression.
Evolution of cytokine environment in melanoma lesions along the regression process
Thirteen cytokines have been measured in tumoral lesions along the regression process (supplementary Table 3). IL-17 and IFNα, weakly expressed (52.3 and 214.4 pg/mg of tumor, respectively), have not been included in further analyses. IL-4, TGFβ and IL-8 were assessed in low amounts (471 to 820 pg/mg of tumor); CCL2, IL-12, IL-6 and IL-10 in moderate amounts (2825 to 3440 pg/mg of tumor); and IFNγ, TNFα, IL-1β and IL-2 in high amounts (6293 to 13,492 pg/mg of tumor). Cytokine environment in tumors has then been studied by a multiple factor analysis. The first two principal components after a PCA explain 90.6% of the variation within the dataset. The plot in Fig. 2a illustrates the relative contribution of each cytokine to the two principal components and highlights two clusters of cytokines. IL-2, TNFα, IFNγ, IL-6, IL-10, IL-12 and IL-4 correlate together and are main contributors to the first dimension. IL-1β, CCL2, IL-8 and TGFβ are more spread along the second-dimension axis. The evolution profile of these clusters along the regression is represented Fig. 2b. Tumors with no regression (R0 group) express more IL-8 and tumors initiating regression (R1) express more TGFβ. On the opposite, less IL-2, TNFα, IL-6, IL-1β, IL-10, CCL2, IL-8 and TGFβ are found at the latest stage of regression (R4). Thus, the cytokine profile in tumors evolves along the regression process, with a significant decrease for nearly all the cytokines measured.
Phenotypic characterization of melanoma-infiltrating immune cells during the regression process
Immune cells infiltrating melanoma lesions have been analyzed by multicolor flow cytometry using the gating strategy we recently published [15]. All individual data are reported in supplementary Table 5. The number of non-immune cells (CD45−) per gram of tumor tends to decrease along the regression, from 27.2 million on average at R1 to 15.8 (p = 0.023) and 19.1 million at R3 and R4, respectively (Fig. 3a). The number of immune cells (CD45+) is quite stable with a mean of 14.5 million of CD45+ cells extracted per gram of tumor (Fig. 3b). They represent 38.4% of total cells (ranging from 18.2 to 70.2%) and slightly increase between R1 and R4 stages (Fig. 3c). Interestingly, the proportions of immune cell subtypes differ during the regression process (Fig. 3d). Myeloid cells account for 93.3% of immune cells associated with the tumor in R0 lesions, but only 54.4% in R4 lesions. In contrast, lymphoid cells represent only 3.4% of cells in R0 lesions, but reach respectively 11.4% and even 29.1% of them at R3 and R4 stages. Immune cell subtypes that have evolved during the regression process are detailed below.
Lymphoid cells infiltrate the tumor at the latest stage of regression
B cells poorly infiltrate tumors, and reach 0.47% of immune cells at stage 4 of regression (Fig. 4a). CD3+ T cells are present in few amounts at the beginning of regression (R1 to R3) and only in a high proportion (27.9% of CD45+ cells) at stage 4 (Fig. 4b). CD3+ T cell functionality of R2 stage lesions has been evaluated by measuring their secretion of IFN-γ by intracellular staining following stimulation of single cell suspension of tumor cells with PMA/Ionomycin (Fig. 4c) and 10.4% of them (range 5–18%) produce IFN-γ after stimulation vs 2.3% when non-stimulated. Immunohistochemistry reveals that CD3+ T cells are located in regressive and intermediate zones from the beginning of the regression and particularly highly concentrated at stage 4. They were also present in remaining tumor zones only at the latest stage of regression (Fig. 4d–f). All T cell subtypes (CD4− CD8α+, CD4+ CD8α−, CD4+ CD8α+, CD4− CD8α−, γδ and NKT T cells) exhibit the same evolution profile along regression, with a statistically significant presence in R4 lesions (Fig. 4g). The majority of them at the latest stage of regression are CD4− CD8α+ T cells representing on average 11.0% of immune cells, followed by γδ T cells (6.4% of immune cells), CD4+ CD8α− T cells (4.5% of immune cells), CD4− CD8α− T cells (3.1% of immune cells), CD4+ CD8α+ T cells (2.1% of immune cells) and NKT cells (0.82% of immune cells). Indeed, in the MeLiM model, B and T cells are detected within melanoma lesions only at the latest stage of regression and may thus be the consequence of this process.
NK cells (CD45+ CD3− CD8α+ CD16+ cells) represent 0.17 and 0.27% of immune cells in tumor lesions of pigs at R0 and R1 stages, respectively. This proportion increases regularly along the regression process reaching 0.59 and 0.75% at R3 and R4 stages (Fig. 4h). NK cell functionality of R2 stage lesions has been evaluated by measuring their secretion of IFN-γ by intracellular staining following stimulation of single cell suspension of tumor cells with PMA/Ionomycin (Fig. 4i). Interestingly, 30.0% of NK cells (range 11.4–42.3%) produce IFN-γ after stimulation, supporting their putative involvement in the tumor regression in vivo.
Slight modulation of PMN and DC during the course of regression
Despite the decrease in proportion of myeloid cells among tumor-infiltrating immune cells in the course of the regression process, myeloid cells represent the main population within the tumor microenvironment whatever the regression stage (Fig. 3d). Among those, polymorphonuclear (PMN) cells represent 20.0% of cells before regression, and decrease all along the regression process reaching 8.2 and 3.9% at R3 and R4 stages, respectively (Fig. 5a). DC remained at low levels in the tumors (ranging from 2.2 to 3.9% of immune cells) without statistically significant evolution during regression (Fig. 5a). All subtypes of DC we have investigated evolve similarly (supplementary Fig. 1).
MHCII+ TAMs accumulate from the beginning of the regression process
TAMs proportion remains stable from R0 to R3 (ranging from 66.5 to 71.0% on average among CD45+ cells) and decreased at the R4 stage where they represent 55.9% of immune cells (Fig. 5a). Four major subsets of TAMs were identified considering their expression of MHCII, CD163 and CD14: MHCII− CD14−, MHCII− CD14+, MHCII+ CD163− and MHCII+ CD163+ (Fig. 5b). MHCII− TAMs (CD14− or CD14+) express low levels of CD163 and MHCII+ TAMs (CD163− or CD163+) low levels of CD14. Interestingly, as illustrated in Fig. 5c, the proportion of the different subtypes among immune cells evolves along the regression. Indeed, MHCII− CD14− and MHCII− CD14+ TAMs significantly decrease with regression, particularly at R3 stage and at R4 stage, respectively. On the opposite, the proportion of MHCII+ TAMs among immune cells increases with regression. Interestingly, that of MHCII+ CD163− TAMs doubles between R0 and R1 lesions derived from age-matched animals, at the earliest step of the regression, and significantly increases at R2 stage (Fig. 5C). When representing the four TAMs subsets in percentages of all TAMs (Fig. 5d), their phenotypes evolve in huge proportion: MHCII− TAMs decrease by half, representing 43 and 34% of total macrophages at R0 for CD14− and CD14+, respectively, and 27 and 14% of them at the end of regression. Conversely, MHCII+ TAMs increase, from 16 to 34% for CD163+ and from 6 to 25% for CD163− when comparing R0 and R4 stages.
CD163− MHCII+ TAMs are functionally very different from other TAMs subsets
To investigate their functionality, we sorted the four different subsets of TAMs from pigs belonging to R1 and R2 stages of regression and analyzed their gene expression profiles by RT-qPCR. Globally, MHCII+ CD163− TAMs clustered and differed significantly from the other subsets (Fig. 6a). The detailed gene expression is presented in Fig. 6b for the genes for which a statistically significant difference could be revealed, others are in supplementary Fig. 2. MHCII+ CD163− TAMs express lower levels of melanoma genes (LYZ, S100A8), cytokines and chemokines genes (IL1B, IL8, IL10, CCL2 and CCL5), of antigen presentation related genes (CD80 and CD274), and of genes implicated in angiogenesis, invasion and epithelial-mesenchymal transition (CSF1, MMP14 and VEGFA). They show also lower expression of ARG1 and IDO1 (even though not significantly) and higher levels of IL4I1. MHCII+ CD163+ TAMs also differ from other TAMs subpopulations, with higher expression of AIF, MPEG1, FCGR2B, VEGFA and lower expression of NOS2, MMP9 and CSF1. The increase of MHCII+ CD163− TAMs, occurring concomitantly with the first signs of regression, and functionally very different from other TAMs subpopulations supports their contribution in the initiation of the regression process.
Discussion
The activation of immune system plays a crucial role in the control of malignant cell growth in cutaneous melanoma. Nevertheless, the immune actors are in most cases inefficient with the exception of infrequent clinical situations characterized by partial or complete spontaneous tumor regression. The regression process is associated with an important tumor lymphoid infiltration [18]. Thus, CD8+ cytotoxic T cells are considered as the main effectors recognizing specifically tumor-associated antigens and their significant presence in melanoma lesions in regression supported their direct involvement in the immune-mediated regression [19, 20]. However, other analyses revealed no differences in the CD8+ T cell infiltration, when comparing either melanoma with or without regression, or regressed and not regressed areas from the same lesions [21, 22]. Such a result discrepancy may be due either to various parameters required for the regression classification or to the lesion analysis at different stages of the regression. Anyway, whether the CD8+ T cell infiltration is a cause or a consequence of tumor regression remains unclear. Beside CD8+ T cells, the contribution of CD56+ NK cells as cytotoxic effectors during early regression phase is supported by the activation of NK cell receptors in the tumor bed [23, 24] and the increased NK cell number in regressing melanocytic lesions [21]. Alternatively, NK cells might eliminate melanoma cells through the secretion of TNFα and IFNγ promoting Th1 differentiation and the recruitment of CD8+ T cells into regressing tumors [25, 26]. The number of global CD4+ T cells was similar in human melanoma lesions with or without any regression [22]. Among CD4+ T cells, FOXP3+ regulatory T cells are rare in areas with regression, and much more frequent in areas without regression [22, 27]. A recent study in the MeLiM pig model showed that CD4−CD8+ T cells accounted for the majority of tumor-infiltrating lymphocytes (TILs) during late phases of melanoma regression and that only CD4+CD8hi TILs expanded during latest stages [28]. Here we deciphered more finely the kinetics of the initiation of regression. We found that neither CD4+CD8α−, CD4+CD8α+ and CD4−CD8α+ T cells, nor γδ T cells, were present in the active phase of the regression and accumulated only in late regression. This data supports the notion of a regression process initiated independently of these T cell subsets presence. We further reported that NK cells infiltrated tumors earlier than T cells consistent with their contribution in melanoma regression and their potential role in the recruitment of CD8+ T cells within the tumor microenvironment.
B cells are unique producers of antibodies, but several B-cell subsets can differently shape the responses to melanoma through their capacity to produce pro-inflammatory or anti-inflammatory cytokines and to present tumor-associated antigens [29]. Data on B cells infiltrating human melanoma regression are still limited, and their involvement in this process remains elusive. The number of global CD20+ B cells was similar in human melanoma lesions with or without any evidence of regression [22]. Even if B lymphocyte infiltration was very limited in MeLiM and restrained to the late phases of regression, we have previously reported an immunoglobulin signature [14] and that tumor regression is frequently accompanied by a local or systemic depigmentation of the skin, hair, and eyes that is associated with a CD4 haplotype affecting the concentration of seric Ig [13].
Regarding the role of myeloid cells, human melanoma regression areas were characterized by a high number of DC that sustains their involvement in the spontaneous tumor cell destruction [27, 30]. In all studied MeLiM tumors, the proportion of different subtypes of DC remained low and quite stable arguing against crucial role for DC in melanoma regression in the pig model. In contrast, PMN cells decreased all along the regression process. Whether such a reduction contributes to the tumor control requires further investigations.
TAMs account for a substantial but variable proportion of the immune cell infiltrate in human cutaneous melanoma [31]. Even though TAMs promote tumor growth in many solid tumor types [32,33,34], they may also exhibit protective roles in specific disease stages or cancer types [35, 36]. Their significance for clinical outcome remains unclear probably due to the large spectrum of their phenotypes ranging from classically activated (or M1-like) and alternatively activated (or M2-like) associated with anti- and pro-tumoral functions, respectively, and their heterogeneity within the tumor microenvironment [37,38,39]. Thus, Etzerodt et al. recently identified in a melanoma mouse model four different TAM subsets depending on their expression of the scavenger receptor CD163 and MHCII (CD163− MHCII−, CD163− MHCII+, CD163+ MHCII−, CD163+ MHCII+) [40]. CD163+ TAMs represent the most mature and immunosuppressive, characterized by an increased expression of genes associated with a M2-like signature (including Il4ra, Mrc1, Stab1 and Slco2b1) or T cell suppression (Il10, Ido1 and Lgals1), and their specific depletion reduced significantly melanoma growth by driving the recruitment of inflammatory monocytes and subsequently CD4+ and CD8+ T cell recruitment and activation. Notably, CD163+ TAMs are associated with a poor prognosis in human melanoma [41].
Compared to the abundant knowledge of TAMs in tumorigenesis [34, 39], the role of the TAMs in tumor regression is only poorly investigated. For instance, TAMs can switch from pro-tumoral to anti-tumoral status and cooperate with TILs leading to reduced tumor burden after STING/Type 1 IFN activation or inhibition of Class IIa HDAC or Clever-1 [42,43,44]. Our data in the MeLiM model highlight for the first time the role of TAMs in the spontaneous melanoma regression. Indeed, they represented about 70% among tumor-infiltrating CD45+ cells from R0 to R3 stages. In particular, MHCII+ CD163− TAMs were the earliest tumor-infiltrating immune population during the regression process that doubled between R0 and R1 lesions and significantly increased at R2 stage. More interestingly, this subset differed significantly from MHCII+ CD163+, MHCII− CD14− and MHCII− CD14+ TAMs in R1 and R2 stages. It exhibited lower levels of IL1B, IL8 and IL10 and CCL5 genes and of genes implicated in angiogenesis (i.e. VEGFa). Their low expression of antigen presentation related genes (CD80 and CD274) and the absence of T cells at the earliest stage of melanoma regression support a T-cell independent action of macrophages on tumor cells. Nonetheless, MHCII + CD163− TAMs accumulate within the regressing lesions (R3, R4) and could also promote the activation of tumor-specific T cells. Compared to other TAMs subsets, they display very low levels of Arginase 1 and IDO1 and a high level of IL4I1 gene expression.
IL4I1 is a secreted l-amino acid oxidase that mainly catabolizes l-phenylalanine, known to be expressed by TAMs of most solid tumors [45, 46]. This enzyme promotes tumor growth by shaping the immune microenvironment in transplanted and spontaneous melanoma murine models [47, 48]. Our recent data also showed that IL4I1 expression is related to a microenvironment enriched in regulatory T cells and poor in granzyme B-positive CD8+ T cells in human primary cutaneous melanoma and is associated with a higher risk of poor outcome in patients [49]. Altogether, these data support a key role for IL4I1 in melanoma aggressiveness. The high transcriptional IL4I1 expression by MHCII+ CD163− TAMs in MeLiM pigs with regressing melanoma was unexpected, but may result from the enriched IFNγ, TNFα, IL12 and IL6 tumor microenvironment between R0 and R1 stages. Interestingly, programmed cell death 1 receptor is a marker of both activation and exhaustion of melanoma-infiltrating T cells. Whether IL4I1 expression impacts the properties of TAMs similarly requires further investigations. The conflicting effects of the IL4I1 enzyme may depend on the tumor context, either deleterious in the context of a chronic inflammation (tumor progression) or anti-tumoral in the context of tumor regression. Our present data shows that IL4I1 produced by MHCII+ CD163− TAMs does not prevent tumor control in our model. This enzyme, known for its primarily immunosuppressive effect on T cells, may even contribute to the regression of melanoma, at the earliest stages of the process when T cells are absent from the melanoma microenvironment. Finally, it is crucial to unravel whether IL4I1 derived from MHCII+ CD163− TAMs initiates the tumor control or is merely a consequence of the tumor regression to develop an accurate therapy targeting this enzyme.
In summary, spontaneously regressing melanoma represents a fascinating model to improve our knowledge in the natural immune process of tumor clearance. Our study clearly identified distinct waves of immune cells from the earliest stage of regression to the latest one and evidenced the presence of a particular TAMs subset in very early regressing lesions. Although current immunotherapies mainly aim to reinvigorating the TILs activity, harnessing both the adaptive and innate arms of the immune system emerged to improve treatment. In particular, the appropriate activation of MHCII+ CD163− TAMs might lead to more efficient treatment for melanoma.
Data availability
The data generated in this study are available within the article and its supplementary data files.
Abbreviations
- FAMD:
-
Factorial analysis of mixed data
- MeLiM:
-
Melanoma-bearing Libechov Minipigs
- PCA:
-
Principal component analysis
- TAM:
-
Tumor-associated macrophage
- TIL:
-
Tumor-infiltrating lymphocyte
References
Schadendorf D, van Akkooi ACJ, Berking C et al (2018) Melanoma. Lancet 392:971–984. https://doi.org/10.1016/S0140-6736(18)31559-9
Blessing K, Mclaren KM (1992) Histological regression in primary cutaneous melanoma: recognition, prevalence and significance. Histopathology 20:315–322
Ribero S, Moscarella E, Ferrara G et al (2016) Regression in cutaneous melanoma: a comprehensive review from diagnosis to prognosis. J Eur Acad Dermatol Venereol 30:2030–2037. https://doi.org/10.1111/jdv.13815
McGovern VJ, Shaw HM, Milton GW (1983) Prognosis in patients with thin malignant melanoma: influence of regression. Histopathology 7:673–680. https://doi.org/10.1111/j.1365-2559.1983.tb02279.x
Ronan SG (1987) Thin malignant melanomas with regression and metastases. Arch Dermatol 123:1326. https://doi.org/10.1001/archderm.1987.01660340088026
Ribero S, Gualano MR, Osella-Abate S et al (2015) Association of histologic regression in primary melanoma with sentinel lymph node status: a systematic review and meta-analysis. JAMA Dermatol 151:1301–1307. https://doi.org/10.1001/jamadermatol.2015.2235
Zugna D, Senetta R, Osella-Abate S et al (2018) Favourable prognostic role of histological regression in stage III positive sentinel lymph node melanoma patients. Br J Cancer 118:398–404. https://doi.org/10.1038/bjc.2017.397
Zhang Y, Zhang Z (2020) The history and advances in cancer immunotherapy: understanding the characteristics of tumor-infiltrating immune cells and their therapeutic implications. Cell Mol Immunol 17:807–821. https://doi.org/10.1038/s41423-020-0488-6
Attrill GH, Ferguson PM, Palendira U et al (2021) The tumour immune landscape and its implications in cutaneous melanoma. Pigment Cell Melanoma Res 34:529–549. https://doi.org/10.1111/pcmr.12926
Gardner LJ, Strunck JL, Wu YP, Grossman D (2019) Current controversies in early-stage melanoma: Questions on incidence, screening, and histologic regression. J Am Acad Dermatol 80:1–12. https://doi.org/10.1016/j.jaad.2018.03.053
Hruban V, Horák V, Fortýn K et al (2004) Inheritance of malignant melanoma in the MeLiM strain of miniature pigs. Vet Med (Praha) 49:453–459. https://doi.org/10.17221/5739-VETMED
Vincent-Naulleau S, Le Chalony C, Leplat J-J et al (2004) Clinical and histopathological characterization of cutaneous melanomas in the melanoblastoma-bearing libechov minipig model. Pigment Cell Res 17:24–35. https://doi.org/10.1046/j.1600-0749.2003.00101.x
Blanc F, Créchet F, Bruneau N et al (2018) Impact of a CD4 gene haplotype on the immune response in minipigs. Immunogenetics 70:209–222. https://doi.org/10.1007/s00251-017-1037-z
Rambow F, Piton G, Bouet S et al (2008) Gene expression signature for spontaneous cancer regression in melanoma pigs. Neoplasia 10:714–726. https://doi.org/10.1593/neo.08344
Blanc F, Prévost-blondel A, Piton G et al (2020) The Composition of circulating leukocytes varies with age and melanoma onset in the MeLiM pig biomedical model. Front Immunol 11:1–14. https://doi.org/10.3389/fimmu.2020.00291
Planska D, Kovalska J, Cizkova J, Horak V (2018) Tissue rebuilding during spontaneous regression of melanoma in the melanoma-bearing libechov minipig. Anticancer Res 38:4629–4636. https://doi.org/10.21873/anticanres.12767
Bordet E, Blanc F, Tiret M et al (2018) Porcine reproductive and respiratory syndrome virus type 1.3 Lena triggers conventional dendritic cells 1 activation and T helper 1 immune response without infecting dendritic cells. Front Immunol 9:1–14. https://doi.org/10.3389/fimmu.2018.02299
Nathanson (1976) Spontaneous regression of malignant melanoma: a review of the literature on incidence, clinical features, and possible mechanisms. Natl Cancer Inst Monogr 44:67–76
Aung PP, Nagarajan P, Prieto VG (2017) Regression in primary cutaneous melanoma: etiopathogenesis and clinical significance. Lab Investig 97:657–668. https://doi.org/10.1038/labinvest.2017.8
Mackensen A, Carcelain G, Viel S et al (1994) Direct evidence to support the immunosurveillance concept in a human regressive melanoma. J Clin Invest 93:1397–1482. https://doi.org/10.1172/jci117116
McKay K, Moore PC, Smoller BR, Hiatt KM (2011) Association between natural killer cells and regression in melanocytic lesions. Hum Pathol 42:1960–1964. https://doi.org/10.1016/j.humpath.2011.02.019
Osella-Abate S, Conti L, Annaratone L et al (2019) Phenotypic characterisation of immune cells associated with histological regression in cutaneous melanoma. Pathology 51:487–493. https://doi.org/10.1016/j.pathol.2019.04.001
Moretta L, Bottino C, Pende D et al (2006) Surface NK receptors and their ligands on tumor cells. Semin Immunol 18:151–158. https://doi.org/10.1016/j.smim.2006.03.002
Gauthier L, Morel A, Anceriz N et al (2019) Multifunctional natural killer cell engagers targeting NKp46 trigger protective tumor immunity. Cell 177:1701-1713.e16. https://doi.org/10.1016/j.cell.2019.04.041
Habif G, Crinier A, André P et al (2019) Targeting natural killer cells in solid tumors. Cell Mol Immunol 16:415–422. https://doi.org/10.1038/s41423-019-0224-2
Cheng M, Chen Y, Xiao W et al (2013) NK cell-based immunotherapy for malignant diseases. Cell Mol Immunol 10:230–252. https://doi.org/10.1038/cmi.2013.10
Cioplea M, Nichita L, Georgescu D et al (2020) FOXP3 in melanoma with regression: between tumoral expression and regulatory T cell upregulation. J Immunol Res. https://doi.org/10.1155/2020/5416843
Cizkova J, Sinkorova Z, Strnadova K et al (2019) The role of αβ T-cells in spontaneous regression of melanoma tumors in swine. Dev Comp Immunol 92:60–68. https://doi.org/10.1016/j.dci.2018.10.001
Prevost-Blondel A, Richard Y (2019) Interleukin 4–induced gene 1 as an emerging regulator of B-cell biology and its role in cutaneous melanoma. Crit Rev Immunol 39:39–57. https://doi.org/10.1615/CritRevImmunol.2019030020
Nedelcu RI, Ion DA, Holeab CA et al (2015) Dendritic cells in melanoma—immunohistochemical study and research trends. Rom J Morphol Embryol 56:997–1002
Erdag G, Schaefer JT, Smolkin ME et al (2012) Immunotype and immunohistologic characteristics of tumor-infiltrating immune cells are associated with clinical outcome in metastatic melanoma. Cancer Res 72:1070–1080. https://doi.org/10.1158/0008-5472.CAN-11-3218
DeNardo DG, Ruffell B (2019) Macrophages as regulators of tumour immunity and immunotherapy. Nat Rev Immunol 19:369–382. https://doi.org/10.1038/s41577-019-0127-6
Bercovici N, Guerin MV, Trautmann A, Donnadieu E (2019) The remarkable plasticity of macrophages: a chance to fight cancer. Front Immunol. https://doi.org/10.3389/fimmu.2019.01563
Duhan V, Smyth MJ (2021) Innate myeloid cells in the tumor microenvironment. Curr Opin Immunol 69:18–28. https://doi.org/10.1016/j.coi.2021.01.001
Ruffell B, Coussens LM (2015) Macrophages and therapeutic resistance in cancer. Cancer Cell 27:462–472. https://doi.org/10.1016/j.ccell.2015.02.015
Nalio Ramos R, Missolo-Koussou Y, Gerber-Ferder Y et al (2022) Tissue-resident FOLR2+ macrophages associate with CD8+ T cell infiltration in human breast cancer. Cell 185:1189-1207.e25. https://doi.org/10.1016/j.cell.2022.02.021
Falleni M, Savi F, Tosi D et al (2017) M1 and M2 macrophages’ clinicopathological significance in cutaneous melanoma. Melanoma Res 27:200–210. https://doi.org/10.1097/CMR.0000000000000352
Tremble LF, McCabe M, Walker SP et al (2020) Differential association of CD68+ and CD163+ macrophages with macrophage enzymes, whole tumour gene expression and overall survival in advanced melanoma. Br J Cancer 123:1553–1561. https://doi.org/10.1038/s41416-020-01037-7
Mantovani A, Marchesi F, Jaillon S et al (2021) Tumor-associated myeloid cells: diversity and therapeutic targeting. Cell Mol Immunol 18:566–578. https://doi.org/10.1038/s41423-020-00613-4
Etzerodt A, Tsalkitzi K, Maniecki M et al (2019) Specific targeting of CD163+ TAMs mobilizes inflammatory monocytes and promotes T cell-mediated tumor regression. J Exp Med 216:2394–2411. https://doi.org/10.1084/jem.20182124
Lee WJ, Lee MH, Kim HT et al (2019) Prognostic significance of CD163 expression and its correlation with cyclooxygenase-2 and vascular endothelial growth factor expression in cutaneous melanoma. Melanoma Res 29:501–509. https://doi.org/10.1097/CMR.0000000000000549
Guerriero JL, Sotayo A, Ponichtera HE et al (2017) Class IIa HDAC inhibition reduces breast tumours and metastases through anti-tumour macrophages. Nature 543:428–432. https://doi.org/10.1038/nature21409
Weiss JM, Guerin MV, Regnier F et al (2017) The STING agonist DMXAA triggers a cooperation between T lymphocytes and myeloid cells that leads to tumor regression. Oncoimmunology 6:e1346765. https://doi.org/10.1080/2162402X.2017.1346765
Viitala M, Virtakoivu R, Tadayon S et al (2019) Immunotherapeutic blockade of macrophage clever-1 reactivates the CD8+ T-cell response against immunosuppressive tumors. Clin Cancer Res 25:3289–3303. https://doi.org/10.1158/1078-0432.CCR-18-3016
Carbonnelle-Puscian A, Copie-Bergman C, Baia M et al (2009) The novel immunosuppressive enzyme IL4I1 is expressed by neoplastic cells of several B-cell lymphomas and by tumor-associated macrophages. Leukemia 23:952–960. https://doi.org/10.1038/leu.2008.380
Mulder K, Patel AA, Kong WT et al (2021) Cross-tissue single-cell landscape of human monocytes and macrophages in health and disease. Immunity 54:1883-1900.e5. https://doi.org/10.1016/j.immuni.2021.07.007
Lasoudris F, Cousin C, Prevost-Blondel A et al (2011) IL4I1: an inhibitor of the CD8+ antitumor T-cell response in vivo. Eur J Immunol 41:1629–1638. https://doi.org/10.1002/eji.201041119
Bod L, Lengagne R, Wrobel L et al (2017) IL4-induced gene 1 promotes tumor growth by shaping the immune microenvironment in melanoma. Oncoimmunology 6:e1278331. https://doi.org/10.1080/2162402X.2016.1278331
Ramspott JP, Bekkat F, Bod L et al (2018) Emerging role of IL-4–induced gene 1 as a prognostic biomarker affecting the local T-cell response in human cutaneous melanoma. J Invest Dermatol 138:2625–2634. https://doi.org/10.1016/j.jid.2018.06.178
Acknowledgements
The authors wish to acknowledge the technical support provided by Edwige Bouguyon (Université Paris‐Saclay, INRAE, VIM, Jouy-en-Josas, France), Jean-Luc Coville and Fabrice Andréoletti (Université Paris‐Saclay, INRAE, AgroParisTech, GABI, Jouy‐en‐Josas, France).
The authors greatly acknowledge the Cytometry and Immunobiology Facility (Cybio) of the Cochin Institute, the MetaFun platform of MetaGénoPolis for access to the flow cytometers apparatus. This work has benefited from the facilities and expertise of @BRIDGe (GABI, INRAE, AgroParisTech, Université Paris‐Saclay) for histology and virtual slide scanner and of the porcine experimental unit (IE1313, GABI, INRAE, AgroParisTech, Université Paris‐Saclay) housing the MeLiM pigs.
Funding
This work was supported by Institut National du Cancer, in the frame of the project INCa-ITMO Cancer inside the plan cancer 2009–2013—project MeLiMun. F. Blanc was supported by a post-doctoral fellowship from the “Institut National du Cancer”.
Author information
Authors and Affiliations
Contributions
FB, NB, SV-N and AP-B designed experiments. FB, GP, SV-N and EB did the experiments. J-JL was in charge of the animal experimentations. FB analyzed the flow cytometry data. FB, GE and SV-N analyzed in situ data. FB made the figures. FB and AP-B wrote the paper. All authors read and approved the submitted version.
Corresponding authors
Ethics declarations
Conflict of interest
The authors have no relevant financial or non-financial interests to disclose.
Ethics approval
Animal care and use were carried out under licenses issued by the Direction Départementale de la Protection des Populations des Yvelines (DDPPY, agreement number C78-719) and experiments were ethically approved by the Committee on the Ethics of Animal Experiments of AgroParisTech and INRAE Jouy-en-Josas (COMETHEA, authorization number 12/091).
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Blanc, F., Bertho, N., Piton, G. et al. Deciphering the immune reaction leading to spontaneous melanoma regression: initial role of MHCII+ CD163− macrophages. Cancer Immunol Immunother 72, 3507–3521 (2023). https://doi.org/10.1007/s00262-023-03503-6
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00262-023-03503-6