ORIGINAL RESEARCH article

Front. Immunol., 21 May 2025

Sec. Systems Immunology

Volume 16 - 2025 | https://doi.org/10.3389/fimmu.2025.1579205

This article is part of the Research TopicAdvances in Management of Aggressive Thyroid Cancer: Medullary and Advanced Thyroid CancerView all 6 articles

Spatial immune profiling reveals distinct microenvironments in medullary thyroid carcinoma

Maria Eduarda de CastroMaria Eduarda de Castro1Gustavo Forlin de SiqueiraGustavo Forlin de Siqueira1Jean Ferrante MarianoJean Ferrante Mariano1Marina Malta Letro KizysMarina Malta Letro Kizys1Lucieli CeolinLucieli Ceolin1Fernando Augusto SoaresFernando Augusto Soares2Rodrigo NatalRodrigo Natal2Humberto CarneiroHumberto Carneiro2Rodrigo Nalio Ramos,Rodrigo Nalio Ramos2,3Laura Sterian WardLaura Sterian Ward4Niels Olsen Saraiva CamaraNiels Olsen Saraiva Camara3Cleber Pinto Camacho,Cleber Pinto Camacho1,5Flavia de Oliveira Facuri ValenteFlavia de Oliveira Facuri Valente1Susan Chow LindseySusan Chow Lindsey1Diego Dias dos SantosDiego Dias dos Santos6Cristiane Damas GilCristiane Damas Gil6Joo Roberto Maciel MartinsJoão Roberto Maciel Martins1Rui Monteiro de Barros Maciel*Rui Monteiro de Barros Maciel1*Lucas Leite Cunha,Lucas Leite Cunha1,7
  • 1Laboratory of Molecular and Translational Endocrinology, Department of Medicine, Escola Paulista de Medicina, Universidade Federal de São Paulo, São Paulo, Brazil
  • 2Pathology Division, D’Or Research Institute, Rede D’Or Hospitals Network, São Paulo, Brazil
  • 3Department of Immunology, Institute of Biomedical Sciences, University of São Paulo, (USP), São Paulo, Brazil
  • 4Laboratory of Cancer Molecular Genetics, School of Medical Sciences, University of Campinas (UNICAMP), Campinas, São Paulo, Brazil
  • 5Molecular Innovation and Biotechnology Laboratory, Postgraduate Medicine Program, Universidade Nove de Julho (UNINOVE), São Paulo, Brazil
  • 6Department of Morphology and Genetics, Escola Paulista de Medicina (EPM), Universidade Federal de São Paulo (UNIFESP), SP, São Paulo, Brazil
  • 7Emergency Medicine and Evidence Based Medicine Unit, Department of Medicine, Escola Paulista de Medicina, Universidade Federal de São Paulo, São Paulo, Brazil

Introduction: Medullary thyroid carcinoma (MTC) is a rare and aggressive thyroid cancer with a challenging prognosis. While the immune microenvironment plays a crucial role in cancer progression, its role in MTC remains underexplored compared to more common thyroid cancers.

Methods: this study investigates the immune landscape of MTC by systematically evaluating immune cell infiltration and expression of immune markers across various tissue topographies. We utilized advanced immunohistochemical techniques to analyze tissue samples from 24 MTC patients, focusing on the tumor core, interface with healthy tissue, adjacent normal thyroid tissue, and lymph node metastases.

Results: our findings reveal a distinct immune profile with increased CD3+, CD4+, CD8+ and CD20+ lymphocytes in normal tissues adjacent to tumors and a notable presence of granzyme B+ cells in the tumor interface, particularly in patients with structural disease. Additionally, we observed a significant enrichment of mast cells in metastatic tissues.

Discussion: these results highlight the complex and spatially dependent immune landscape of MTC, suggesting implications for targeted immunotherapy. This study provides novel insights into the immune microenvironment of MTC and emphasizes the need for further research to elucidate its impact on disease progression and therapeutic response.

Introduction

Medullary thyroid carcinoma (MTC) is a rare yet aggressive form of thyroid cancer arising from parafollicular C cells (1, 2). Clinically, it is characterized by thyroid nodule and elevated calcitonin levels (3). Despite advances in diagnosis and treatment, the prognosis of MTC remains challenging, with variable survival rates and frequent recurrence, even after surgery and adjuvant therapies (4). The heterogeneous prognosis among patients underscores the urgent need for more precise predictive tools to stratify risks and guide personalized therapeutic decisions.

The immunological microenvironment plays a crucial role in the progression of thyroid cancer, particularly evident in papillary thyroid carcinoma (PTC), the most common type (5). Studies have revealed a complex interplay between tumor cells and components of the immune system, including tumor-infiltrating lymphocytes (TILs) and tumor-associated macrophages (TAM) (6). These interactions influence not only the initial treatment response but also the propensity for metastasis and therapeutic resistance (7).

Nonetheless, studies on the immunological microenvironment of MTC are limited. Previous research suggests a significant presence of TILs and regulatory T cells (Tregs) in the tumor microenvironment of MTC, potentially affecting disease progression and treatment response (8). Furthermore, the expression of immune molecules, such as PD-L1, has been associated with immune evasion in MTC, highlighting its potential as a therapeutic target (9). A comprehensive understanding of these immunological mechanisms is crucial for a better prognostic definition and the development of novel therapeutic strategies for MTC.

Several prognostic markers and treatment protocols have been established and are currently being employed (9, 10). However, comprehensive data integrating histopathological, immunohistochemical, and molecular analyses to understand the tumor microenvironment in MTC, especially with regard to immune cell infiltration and its clinical implications, remain limited. We aimed to fill this gap by systematically evaluating the immune landscape of the MTC tissues. Our objectives were to quantify immune cell infiltration in different tissue regions, evaluate the expression of immune checkpoints, and correlate these findings with the clinical outcomes.

Materials and methods

Design

This study was approved by the Research Ethics Committee of Hospital São Paulo and Universidade Federal de São Paulo, São Paulo, Brazil (309–2023). We investigated 24 patients of the BRASMEN project (CAAE: 70952523.8.0000.5505) with histopathological confirmation of MTC, whose clinical data were updated from data collected from their medical records (11). Aggressiveness at diagnosis was determined using the tumor-node-metastasis classification and staging system for MTC (11). All patients followed a treatment protocol established by the American Thyroid Association (4).

The tissues of the patients, fixed in formalin and embedded in paraffin, were reviewed for diagnostic confirmation and selection of the most representative areas for the construction of a Tissue Microarray (TMA) (Beecher Instruments, Silver Springs, MD, USA). From each case, two tissue samples with an estimated area of ​​0.79 mm² were obtained from the following regions: (i) tumor core, (ii) tumor interface with healthy thyroid tissue, (iii) adjacent healthy thyroid tissue, and (iv) lymph node metastasis (when present).

The inclusion criteria were as follows: (i) diagnosis of MTC in any age group, (ii) clinical data available for evaluation, (iii) primary tumor tissues stored in paraffin available for morphological investigation, and (iv) primary tumor tissues stored in paraffin with quality and viability for molecular studies. Patients were excluded if: (i) they did not have clinical data in their medical records; (ii) they did not present tissue samples in a pathological anatomy repository; (iii) the diagnosis of MTC could not be confirmed by histological slide review; (iv) molecular diagnosis of MEN2 could not be performed by sequencing the RET gene and; (v) the paraffin tissue was not of sufficient quality for molecular analysis.

After total thyroidectomy, patients were classified as: i) biochemical disease-free, in other words, no laboratory evidence of residual disease or recurrence when measuring serum calcitonin or serum carcinoembryonic antigen (CEA); ii) structural disease-free, when there was no evidence of visible disease on imaging tests such as PET-CT; iii) disease-free, those patients who have neither biochemical nor structural disease.

Immunohistochemistry

Immunohistochemistry was performed in a hospital clinical routine laboratory, whose immune cell markers analyzed included TAM (CD68) and subsets of TIL, such as CD3, CD4, CD8, and CD20. The activation status of these immune cells was identified using Granzyme-B and PD-L1(22C3) (Figure 1).

Figure 1
www.frontiersin.org

Figure 1. Immunohistochemical. Immunohistochemistry was carried out in the tumor, interface, normal tissue adjacent to the tumor and metastasis regions for the markers CD3, CD4, CD8, CD20, FOXP3, CD68, granzyme B and PDL1. Scale bar, 20 µm, identified in the bottom left corner of each image.

The analyses were carried out by immunohistochemical staining using the ultraview universal DAB detection kit ventana (Roche) and EnVision FLEX Visualization Systen kit AGILENT (DAKO), according to standard operating protocols. The corresponding positive and negative controls were analyzed and found to be adequate. The respective antibodies are: rabbit polyclonal anti-CD3 antibody (DAKO); mouse monoclonal anti-human anti-CD20 (L26 - DAKO); rabbit monoclonal anti-FOX P3 (EP340 - BIOSB); rabbit monoclonal anti-CD4 (SP35 - Ventana); rabbit monoclonal anti-CD8 (SP57 - Ventana); mouse monoclonal anti-CD68 (PG-M1 - DAKO); mouse monoclonal anti-PD-L1 (22C3 - DAKO); rabbit polyclonal anti-granzyme B (Roche).

The independent assessment of immune cell markers was carried out by calculating the count of positive cells per TMA point based on an estimated area of ​​approximately 0.79 mm². The cases were grouped into the following categories for statistical analysis: 0 (complete absence of immune cells), 1 (up to ten immune cells), and 2 (more than 10 immune cells). Visual assessment of PD-L1(22C3) was performed by estimating the percentage of positive cells according to the Combined Positive Score (CPS), which is defined as the number of PD-L1 positive cells divided by the total number of tumor cells × 100. The PD-L1 CPS cutoffs were: <1 (negative) and ≥ 1 (positive). Since only 22C3 clone was positive in the samples, only 22C3 clone was considered for further analysis.

We also assessed the presence of tumor necrosis and quantified mitotic activity and Ki67 expression. Tumor necrosis was evaluated visually on hematoxylin and eosin-stained slides and classified as either present or absent. Mitotic activity was quantified by counting mitotic figures in ten high-power fields (HPFs, 40× objective) within the most mitotically active area of the tumor, and the results expressed as the number of mitosis per 10 HPFs. Ki67 expression was analyzed by immunohistochemistry using a mouse monoclonal anti-Ki67 antibody (MIB-1; Dako, Glostrup, Denmark). The Ki67 labeling index was calculated as the percentage of positively stained tumor nuclei in a total count of at least 1,000 cells in the area of the highest labeling density (“hot spot”) identified at low magnification (10×).

The tumors were classified as low, intermediate, or high based on the mitotic count, ki-67 proliferative index, and presence or absence of necrosis. Low proliferative activity was defined as <3 mitoses/10 HPFs, a Ki-67 proliferative index of <3%, and absence of coagulative necrosis. An intermediate degree of the proliferative index was considered when there were 3 to 20 mitoses/10 HPFs or a Ki-67 index of 3% to 20%, without coagulative necrosis. Cases with <3 mitoses/10 HPFs and a Ki-67 proliferative index of <3%, but with coagulative necrosis were also considered intermediate grade. A high-grade proliferative index was considered when there were 3–20 mitoses/10 HPFs2 or a Ki-67 index of 3–20% and coagulative necrosis. Cases with >20 mitoses/10 HPFs or a Ki-67 index of >20% were also considered high-grade, regardless of whether or not there was coagulative necrosis.

QuPath software made it possible to quantify the number of positively labeled cells over the total number of cells, in the entire region of the section (12). Furthermore, to standardize the area evaluated, we estimated the value of positive cells in the area of ​​1 mm² and this was the number used for subsequent statistical analysis.

We conducted a hierarchical cluster analysis of our immunohistochemical data using Cluster 3.0 software to evaluate the similarities and distances between patients based on the immunological parameters. In accordance with the authors’ guidelines, we applied complete-linkage hierarchical clustering to the following eight parameters: CD3, CD20, FOXP3, CD4, CD8, CD68, PDL1 (22C3 clone), Granzyme B, and mast cells. The resulting dendrogram, which visually represents patient similarities and distances, was generated using Java TreeView software (13).

Assessment of mast cell infiltration

Mast cell infiltration was assessed using a Rapid Panotic technique (Larboclin). Sections obtained from the samples were prepared as described in the immunohistochemistry section. To visualize mast cells, the slide was immersed for 1 min in solution I (0.1% triarylmethane), 1 min in solution II (0.1% xanthenes), and 30 s in solution III (thiazines to 0.1%). Sections were washed in water buffered at pH 7 and dehydrated for slide preparation. Mast cells were identified and characterized by visualization of metachromatic cytoplasmic granules. The cells were quantified using the 40x objective of QuPath software in the total area of ​​the sections, approximately 0.79 mm².

Statistical analysis

All slides were scanned by ScanScope CS² using a 40x objective for subsequent quantification of immunolabeling using the QuPath software. Clinical, anatomopathological, and immunohistochemical data were crossed in search of associations with statistical significance using the chi-square method for quantitative variables and the Kruskal-Wallis test for independent samples to perform non-parametric tests for ordinal and numeric variables. GraphPad software (version 9) was used to present data as the mean ± standard error of the mean (S.E.M.). Non-normally distributed data were used for datasets with more than two groups. Statistical analysis was performed using the Statistical Package for the Social Sciences (SPSS, Chicago, IL, USA)(R) software, version 23.0. Statistical significance was set at P < 0.05. We performed Benjamini-Hochberg correction to reduce the bias caused by multiple testing, adjusting the p-value to reduce the false discovery rate. All adjusted p-values (Q) < 0.05 were considered significant.

Results

Patient characteristics

Table 1 shows the demographic, clinical, and pathological characteristics of the patients. The majority of patients were diagnosed after the age of 45 years and were female. Most included patients had sporadic MTC. Table 2 summarizes the main pathological characteristics of the investigated tumors. The most common histological type was the classic type, with 62.5% of patients presenting with multifocality, surgical margins free in less than 50% of cases, and a minority with tumors larger than 4 cm. Thyroiditis was present in 12.5% of the cases, and TNM staging was well distributed. Central lymph node invasion was observed in 52.6% of patients. The patients were followed-up in our service for 10 years. Most patients had structural or biochemical remission of the disease after surgery; however, only 36.4% were free of concomitant biochemical and structural disease after surgery (Table 1).

Table 1
www.frontiersin.org

Table 1. Descriptive table of demographic and clinical characteristics.

Table 2
www.frontiersin.org

Table 2. Descriptive table of anatomopathological features.

Tumor necrosis was not observed in any of the cases. Mitotic activity varied greatly from 0 to 4 mitoses per 10 high power fields (HPFs), with a median of 0 mitoses/10 HPFs. A high mitotic index (>3 mitoses/10 HPFs) was present in 12.5% of the cases. The Ki67 labeling index ranged from 0% to 5% and the proliferative index was considered low in 95.8% of cases.

Semi-quantitative analysis of immunological markers

Supplementary Table 1 shows the results of the semi-quantitative analysis of the infiltration of different immune cells in various histological topographies. We did not observe any difference in the categories of immune cell infiltration among samples from different topographies. PD-L1 positivity was rare and equally scarce across the different tissues evaluated.

Quantitative analysis of immunological markers

In general, normal thyroid tissue samples adjacent to the tumor were enriched with immune cells of different immunophenotypes. Immune cell infiltration was analyzed using immunohistochemistry and semi-automatic quantification. Normal tissue samples adjacent to the tumors showed greater infiltration of CD3+ lymphocytes than the center of the tumor (p<0.001) and the region of metastasis (p<0.01; Figure 2A). Normal tissues adjacent to the tumor (p<0.01) showed greater infiltration of CD4+ lymphocytes than samples from the tumor center (Figure 2B). Normal thyroid tissues adjacent to the tumor were more enriched with CD8+ lymphocytes than the tumor center (p<0.001) and interface (p<0.05; Figure 2C). Normal thyroid tissues showed greater infiltration of CD20+ lymphocytes than samples from the tumor center (p<0.01; Figure 2D). Metastatic tissue samples were enriched with immune cells of different immunophenotypes, although statistical analysis showed no differences between other topographies (Figures 2B–H). The samples exhibited sparse infiltration of granzyme B+ lymphocytes, with a tendency towards greater infiltration in metastatic tissues (Figure 2G). Immunohistochemical evaluation of PD-L1 was positive for this marker in all tissue topographies analyzed, with clone 22C3 showing a higher number of positive cells per mm² in the normal tissue adjacent to the tumor when compared to the tumor center (p < 0.001) and the interface region (p < 0.05) (Figure 2H).

Figure 2
www.frontiersin.org

Figure 2. Quantitative analysis of immune cell infiltration and expression of immunological markers in different regions of tumor tissue. Quantification of lymphocyte and macrophage markers (A-H). The analysis reveals a general decrease in the infiltration of T lymphocytes (CD3, CD4, CD8) and B cells (CD20) in tumor and metastatic tissue compared to normal tissue, suggesting an immunosuppressive tumor microenvironment. The low expression of PD-L1 in the tumor shows that this is not one of the main mechanisms of immune escape. The Y-axis indicates the marker analyzed, while the X-axis shows the region analyzed. Data are presented as mean ± S.E.M of protein expression in positive cells per mm², and dot plots represent individual values per group (n = 24 patients/group). *P < 0.05; **P < 0.01; ***P < 0.001 vs. Normal (Kruskal-Wallis, Dunn’s post-test).

Analysis of tissues using rapid panoptic staining revealed the presence of a few mast cells across different histological topographies (Figure 3A). Metastatic tissue samples showed greater infiltration of these cells compared to other topographies, especially when compared to the tumor center (p<0,05; Figure 3B).

Figure 3
www.frontiersin.org

Figure 3. Mast cell analysis. Representative image of histological staining, showing mast cells identified by their characteristic purple staining indicated by red arrows (A). Quantitative analysis revealed a significant increase in the number of mast cells in the region of metastasis compared to the other regions, suggesting a possible association between the presence of mast cells and the metastatic process. Quantification of mast cells per 0.79 mm² (B). Data are represented as mean ± standard error of the mean (S.E.M) of the number of cells per 0.79 mm², and dot plots show individual values by group (n = 24 patients/group). *P < 0.05 (Kruskal-Wallis, Dunn’s post-test).

Analysis of clinical status and clinical characteristics

Supplementary Table 2 summarizes the clinical description of the patients included in this study. Patients with T1 tumors (tumor limited to the thyroid of 2 cm or less) were predominantly found in the group without evidence of active disease (62.5%), followed by those with biochemical disease (37.5%), and none in the group with structural disease (P=0.366). The distribution of T2 tumors (tumor limited to the thyroid larger than 2 cm and smaller than 4 cm) was 25% without active disease, 37.5% with biochemical disease, and 37.5% with structural disease. T3 tumors (tumor with extrathyroidal extension or larger than 4 cm limited to the thyroid) accounted for 25% of patients without active disease, 25% with biochemical disease, and 50% with structural disease. T4 tumors (moderately or very advanced disease, invading other tissues) showed 50% with biochemical disease and 50% with structural disease.

In terms of nodal status, N0 patients (no metastases in regional lymph nodes) were mainly in the group with no evidence of active disease (100%), and none with biochemical or structural disease (P=0.004). N1 patients (metastases in lymph nodes at level VI) were distributed as follows: 28.6% without active disease, 57.1% with biochemical disease, and 14.3% with structural disease. N2 patients (lymph node metastases to other levels) showed 0% without active disease, 50% with biochemical disease, and 50% with structural disease.

Patients with M0 status (no distant metastasis) were predominantly without evidence of active disease (66.7%) and some had biochemical disease (33.3%), with no evidence of structural disease (P=0.230). M1 status (distant metastasis) patients were more distributed: 14,3% without active disease, 28,6% with biochemical disease, and 57,1% with structural disease. Patients with Mx status (it was not possible to identify) were more distributed: 37.5% without active disease, 43.8% with biochemical disease, and 18.8% with structural disease.

The absence of multifocality was associated with 66.7% of patients without evidence of active disease, 22.2% with biochemical disease, and 11.1% with structural disease (P=0.046). The presence of multifocality was associated with 15.4% of patients without active disease, 46.2% with biochemical disease, and 38.5% with structural disease.

Patients who underwent total thyroidectomy (TT) without central emptying had 66.7% without active disease and 33.3% with biochemical disease (P=0.092). Those who had TT with central emptying comprised 57.1% without active disease, 42.9% with biochemical disease, and none with structural disease. Patients with TT with central and lateral emptying showed 16.7% without active disease, 33.3% with biochemical disease, and 50% with structural disease.

Analysis of immunological markers versus pathological and clinical characteristics

Mitotic activity showed a positive relationship with CD68+ infiltration in the metastasis region (Spearman = 0.470; P = 0.032) and a positive relationship with CD4+ (Spearman = 0.949; P = 0.051) and CD8+ (Spearman = 0.949; P = 0.051) infiltration as well as in the metastasis region, but with a non-significant marginal p-value. Tumor size was negatively correlated with infiltration and CD3+ cells in the tumor region (Spearman = -0.497; P = 0.022). Ki-67 expression in the interface region showed a negative relationship with CD68+ infiltration in the normal tissue region adjacent to the tumor (Spearman = -0.458; P = 0.042).

Patients were classified according to their clinical status as follows: (i) no evidence of active disease, (ii) biochemical disease, or (iii) structural disease. Following this classification, we analyzed the infiltration of immune cells according to each status. Normal thyroid tissue showed greater infiltration of CD3 cells, especially in individuals with structural disease (Figure 4A). In addition, metastatic tissues showed a sparse infiltration of CD3, CD4 and CD8 cells in all clinical statuses (Figures 4A–C). Patients with biochemical disease presented with normal thyroid tissue with a higher infiltration of CD20 cells than the tumor core (P<0.01). The same trend was observed in patients with structural disease (Figure 4D). The expressions of FOXP3 and CD68 were not significant (Figures 4E, F). The interface tissue of patients with structural disease showed greater infiltration of granzyme B+ cells than the interface area of individuals with biochemical disease (p<0.05; Figure 4G). PDL1 (22C3) showed higher expression in normal tissues adjacent to the tumor, but not significantly (Figure 4H). Mast cell showed no statistical relevance or significant trends (Figure 4I).

Figure 4
www.frontiersin.org

Figure 4. Marker analysis according to clinical status. Cross-analysis between cell expression (Y-axis) and patient clinical status (X-axis), including no evidence of active disease, biochemical disease, and structural disease (A-I). The legend indicates the analyzed tissue regions. Data are shown as means ± S.E.M of protein expression in positive cells per mm², and dot plots represent individual values by group (n = 24 patients/group). *P < 0.05; **P < 0.01 vs Tumor (Kruskal-Wallis, Dunn’s post-test).

Analysis of the mixture of immune microenvironment markers

Cluster analysis showed that each histological topography exhibited a unique pattern of similarities and discrepancies among the markers, as shown in Figure 5. Furthermore, the correlation matrix analysis revealed that each microenvironment topography presents unique patterns of bivariate correlations between the biomarkers, as visually demonstrated in Figure 6.

Figure 5
www.frontiersin.org

Figure 5. Cluster analysis of the immune microenvironment in different topographies. Hierarchical clustering of immune markers in the four topographies. The immune response profiles were grouped into clusters based on the similarities and differences between the immune markers. (A) center of the tumor. (B) invasive margin. (C) normal thyroid tissue adjacent to the tumor. (D) metastasis. Each column represents an individual sample and each row corresponds to a specific immune marker. GRZ, granzyme B; MC, mast cell. The color scale indicates the normalized expression intensity, with red representing high expression and black representing low expression. The greater the intensity of the color, the greater the relationship between the markers.

Figure 6
www.frontiersin.org

Figure 6. Correlation matrix. The values express Spearman’s correlation coefficients between the density of immunological markers and mast cells (MC) in four different regions of the tumor microenvironment. (A) center of the tumor. (B) invasive margin. (C) normal thyroid tissue adjacent to the tumor. (D) metastasis. Positive values indicate a positive relationship between the markers and negative values indicate an inverse relationship. GRZ, granzyme B; MC, mast cell. The bar next to it indicates the intensity of the relationship between the markers. The greater the intensity of the color, the greater the relationship between the markers.

First, analyzing the microenvironment of the tumor core (Figure 6A), the sparse infiltration of mast cells in the tumor core resulted in this marker showing greater discrepancy compared to other immune markers. Immune function markers such as PD-L1 and Granzyme, demonstrated a stronger correlation of positivity within the tumor core. Mast cells, on the other hand, showed an intense negative relationship with CD3.

In the interface region (Figure 6B), the profile of immune response markers differed from that of the tumor core by exhibiting two major groups of immune responses that diverged in terms of Granzyme B+ cell infiltration. In addition, CD3 and CD8 showed a strong positive correlation.

In normal tissue adjacent to the tumor (Figure 6C), the immune response appeared relatively uniform, with minimal variation in similarities and distinctions among the immune response markers.

The metastasis microenvironment was characterized by strong correlations (Figure 6D). In the metastatic microenvironment, we observed a strong inverse correlation between PD-L1 expression and Granzyme B+ lymphocyte infiltration. In this microenvironment, mast cell infiltration is higher than in other topographies and positively correlates with the enrichment of various other immune cells, such as CD3+, CD4+, CD8 and CD20+ lymphocytes. Conversely, mast cell infiltration was inversely correlated with FOXP3, CD68 and PD-L1 expression levels.

We corrected the p-values using the Benjamini-Hochberg test to obtain the corrected p-values (Q). The correlation of the markers in the intratumoral region revealed a positive and significant correlation coefficient (Spearman = 0.774; Q = 0.001) of Granzyme B with PDL1 and CD4 with CD8 (Spearman = 0.695; Q = 0.006). CD3 and mast cells showed a negative relationship (Spearman = -0.755; Q = 0.001). The interface region showed a positive and significant correlation (Spearman = 0.762; Q = 0.005) between CD3 and CD8. While the normal region adjacent to the tumor showed no significant Q.

Discussion

To our knowledge, this is the first study to demonstrate the spatial dependence of the immune microenvironment on the immune profile of MTC and its relationship with the demographic and pathological characteristics of the patients. Immune profiling indicated that normal thyroid tissues adjacent to tumors were enriched with CD3+, CD4+, CD8+, and CD20+ lymphocytes compared with tumor centers. Metastatic tissues exhibited a tendency to greater infiltration of CD20+ immune cells, although this was not consistently significant. The immune profile highlighted that the immune microenvironment of the interface tissue, where normal thyroid tissue meets the tumor, showed a distinct pattern of immune cell infiltration. This region showed a greater infiltration of CD3+ and CD4+ lymphocytes than the center of the tumor, suggesting a robust immune response at the margins of the tumor. In addition, the interface tissue showed greater infiltration of granzyme B+ cells, particularly in patients with structural disease. Figure 7 illustrates the CMT tumor microenvironment. These results indicate that the interface tissue plays a critical role in the immune response against MTC, with potential implications for prognosis and therapeutic strategies.

Figure 7
www.frontiersin.org

Figure 7. Illustration of the tumor microenvironment. The figure shows the four areas of the tumor analyzed: normal tissue, interface, tumor center and metastasis. Normal tissue shows greater infiltration of T and B lymphocytes. The interface follows the same pattern. In contrast, the tumor center is poor in these cells. Metastasis differs from other tissues, showing infiltration not only of T and B lymphocytes, but also of macrophages and mast cells.

The spatial-dependent immune landscape of the inflammatory immune microenvironment in MTC is consistent with findings from other studies, which highlight a complex interaction between immune cells and tumor progression (14, 15). Studies have shown that normal thyroid tissue adjacent to tumors often exhibits a higher density of immune cells, particularly CD3+, CD4+, CD8+, and CD20+ lymphocytes, indicating a robust immune response in non-tumor tissues (7). Similarly, metastatic tissues frequently demonstrate increased infiltration of mast cells and granzyme B+ lymphocytes, reflecting a localized inflammatory response associated with advanced disease stages (16, 17). The distinct clustering patterns of immune markers, such as the correlation of CD4 and CD8 cells with inflammatory responses in the tumor interface, have also been reported in the literature (18, 19), emphasizing the dynamic and heterogeneous nature of the immune microenvironment in MTC. These findings corroborate the notion of an intricate immune landscape within MTC, characterized by diverse immune cell populations and marker expressions, underscoring the potential for targeted immunotherapeutic strategies, as discussed in recent studies.

While the expression of immune checkpoint markers, such as PD-L1, is generally sparse across different tissue topographies, its elevated levels at metastatic sites suggest an immune evasion mechanism, consistent with observations in other cancer types (20, 21). The knowledge about PD-L1 expression in MTC is limited, and the few existing studies are inconsistent, with varying results (9, 2224). In most of these studies, PD-L1 is minimally expressed in tumor cells, and in one study, its very low expression did not correlate with clinical characteristics (9, 2224). Our study showed that PD-L1 (clone 22C3) is more highly expressed in normal tissue adjacent to the tumor than intratumoral tissue, revealing that this is not the main mechanism of lymphocyte inactivation in MTC and is therefore not a good way to treat patients. Some studies have similar results to ours; analysis with the SP142 and SP263 clones revealed almost no PD-L1 expression in MTC cells and inflammatory cells, and no clinicopathological or prognostic association (22, 24). A Chinese study with 49 patients analyzed PD-L1 expression for clones 22C3 and SP142; neither clone was significantly expressed in MTC and was not related to clinicopathological parameters (25).

Furthermore, we demonstrated for the first time that mast cells are abundant in the metastatic immunological environment of MTC. This enrichment was especially noted among patients with aggressive disease and biochemical manifestations of cancer. In the context of tumor immune biology, mast cells are a negligible subpopulation of immune cells, and few studies have emphasized their role in tumor progression. The spatial accumulation of mast cells in a metastatic microenvironment may indicate that mast cells promote tumor progression (26). Mast cell density is a key regulatory point for the local immune response (27, 28) and it is associated with angiogenesis and poor prognosis in some solid tumors (2932). Although further studies are needed to evaluate the translational impact of inhibiting mast cell activity, especially in patients with an aggressive MTC phenotype, this finding opens interesting perspectives for new therapeutic approaches.

In conclusion, our study presents a comprehensive analysis of the immunological landscape of MTC, revealing critical spatial dependencies in immune cell infiltration across different tissue topographies. Our findings underscore the enriched presence of CD3+, CD4+, CD8+, and CD20+ lymphocytes in normal thyroid tissues adjacent to tumors compared to tumor centers, suggesting a distinct and potentially more robust immune response in non-tumor regions. The identification of higher immune cell infiltration at the tumor interface and in metastatic tissues further emphasizes the dynamic interaction between immune cells and tumor progression in MTC. Notably, our study is the first to document the enrichment of mast cells in the metastatic immune microenvironment of MTC, particularly in patients with aggressive disease, suggesting a novel role for mast cells in tumor progression. These insights into the immune microenvironment of MTC could inform the development of targeted immunotherapies and improve prognostic strategies, ultimately contributing to more personalized and effective management of this challenging cancer.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by Research Ethics Committee of Hospital São Paulo and Universidade Federal de São Paulo, São Paulo, Brazil (309–2023). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

MD: Writing – original draft, Writing – review & editing. GD: Writing – original draft. JFM: Writing – original draft. MK: Writing – review & editing. LC: Writing – review & editing. FS: Writing – original draft. RN: Writing – review & editing. HC: Writing – review & editing. RR: Writing – original draft. LW: Writing – review & editing. NC: Writing – original draft. CC: Writing – original draft. FF: Writing – review & editing. SL: Writing – original draft. DS: Writing – original draft. CG: Writing – original draft. JM: Writing – review & editing. RM: Writing – review & editing. LLC: Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. The São Paulo Research Foundation (FAPESP) to RMBM (21/02752-6), 2024/06697-8 (LLC), MMLK (2022/10804-9), 2024/01435-5 (MEC), 2022/13451-0 (GFS), 2022/11509-0 (LC). National Council for Scientific and Technological Development (CNPq) and the Coordination for higher Education Staff Development (CAPES) to financial support.

Conflict of interest

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

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

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

Supplementary material

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

References

1. Gild LM, Clifton-Bligh JR, Wirth JL, and Robinson GB. Medullary thyroid cancer: updates and challenges. Endocr Rev. (2023) 44:934–46. doi: 10.1210/endrev/bnad013

PubMed Abstract | Crossref Full Text | Google Scholar

2. Mathiesen SJ, Kroustrup PJ, Vestergaard P, Stochholm K, Poulsen LP, Rasmussen KÅ, et al. Survival and long-term biochemical cure in medullary thyroid carcinoma in Denmark 1997–2014: A nationwide study. Thyroid. (2019) 29:368–77. doi: 10.1089/thy.2018.0564

PubMed Abstract | Crossref Full Text | Google Scholar

3. Schmid WK. Histopathology of C cells and medullary thyroid carcinoma. Medullary Thyroid Carcinoma. (2015) 204:41–60. doi: 10.1007/978-3-319-22542-5_2

PubMed Abstract | Crossref Full Text | Google Scholar

4. Wells SA Jr., Asa SL, Dralle H, Elisei R, Evans DB, Gagel RF, et al. Revised American Thyroid Association guidelines for the management of medullary thyroid carcinoma. Thyroid. (2015) 25:567–610. doi: 10.1089/thy.2014.0335

PubMed Abstract | Crossref Full Text | Google Scholar

5. Cunha LL, Marcello AM, and Ward SL. The role of the inflammatory microenvironment in thyroid carcinogenesis. Endocr-Relat Cancer. (2014) 21:R85–R103. doi: 10.1530/ERC-13-0431

PubMed Abstract | Crossref Full Text | Google Scholar

6. Cunha LL, Morari CE, Guihen TCA, Razolli D, Gerhard R, Nonogaki S, et al. Infiltration of a mixture of different immune cells may be related to molecular profile of differentiated thyroid cancer. Endocr-Relat Cancer. (2012) 19:L31–L6. doi: 10.1530/ERC-11-0285

PubMed Abstract | Crossref Full Text | Google Scholar

7. Cunha LL, Morari CE, Guihen TCA, Razolli D, Gerhard R, Nonogaki S, et al. Infiltration of a mixture of immune cells may be related to good prognosis in patients with differentiated thyroid carcinoma. Clin Endocrinol. (2012) 77:918–25. doi: 10.1111/j.1365-2265.2012.04482.x

PubMed Abstract | Crossref Full Text | Google Scholar

8. Müller S, Poehnert D, Müller AJ, Scheumann FWG, Koch M, and Lück R. Regulatory T cells in peripheral blood, lymph node, and thyroid tissue in patients with medullary thyroid carcinoma. World J Surge. (2010) 34:1481–7. doi: 10.1007/s00268-010-0484-6

PubMed Abstract | Crossref Full Text | Google Scholar

9. Wusiman D, Guo L, Li L, Zhang X, Zhao X, An Z, et al. Clinicopathological and prognostic significance of PD-L1 and TIM-3 expression in medullary thyroid carcinoma: a retrospective immunohistochemistry study. J Endocrinol Invest. (2023) 47:91–100. doi: 10.1007/s40618-023-02126-z

PubMed Abstract | Crossref Full Text | Google Scholar

10. Franc B, Rosenberg-Bourgin M, Caillou B, Dutrieux-Berger N, Floquet J, Houcke-Lecomte M, et al. Medullary thyroid carcinoma: Search for histological predictors of survival (109 proband cases analysis). Hum Pathol. (1998) 29:1078–84. doi: 10.1016/s0046-8177(98)90417-8

PubMed Abstract | Crossref Full Text | Google Scholar

11. Maciel BMR, Camacho PC, Assumpção MVL, Bufalo EN, Carvalho LA, Carvalho DAG, et al. Genotype and phenotype landscape of MEN2 in 554 medullary thyroid cancer patients: the BrasMEN study. Endocr Connect. (2019) 8:289–98. doi: 10.1530/EC-18-0506

PubMed Abstract | Crossref Full Text | Google Scholar

12. Bankhead P, Loughrey BM, Fernández AJ, Dombrowski Y, Mcart GD, Dunne DP, et al. QuPath: Open source software for digital pathology image analysis. Sci Rep. (2017) 7:16878. doi: 10.1038/s41598-017-17204-5

PubMed Abstract | Crossref Full Text | Google Scholar

13. Saldanha JA. Java Treeview—extensible visualization of microarray data. Bioinformatics. (2004) 20:3246–8. doi: 10.1093/bioinformatics/bth349

PubMed Abstract | Crossref Full Text | Google Scholar

14. Menicali E, Guzzetti M, Morelli S, Moretti S, and Puxeddu E. Immune landscape of thyroid cancers: new insights. Front Endocrinol (Lausanne). (2020) 11:637826. doi: 10.3389/fendo.2020.637826

PubMed Abstract | Crossref Full Text | Google Scholar

15. Papachristos JA, Serrao-Brown H, Gill JA, Clifton-Bligh R, and Sidhu BS. Medullary thyroid cancer: molecular drivers and immune cellular milieu of the tumour microenvironment—Implications for systemic treatment. Cancers. (2024) 16:2296. doi: 10.3390/cancers16132296

PubMed Abstract | Crossref Full Text | Google Scholar

16. Hou Y, Wang Q, Su L, Zhu Y, Xiao Y, and Feng F. Increased tumor-associated mast cells facilitate thyroid cancer progression by inhibiting CD8+ T cell function through galectin-9. Braz J Med Biol Res. (2023) 56:e12370. doi: 10.1590/1414-431X2023e12370

PubMed Abstract | Crossref Full Text | Google Scholar

17. Cunha LL, Nonogaki S, Soares AF, Vassallo J, and Ward SL. Immune escape mechanism is impaired in the microenvironment of thyroid lymph node metastasis. Endocr Pathol. (2017) 28:369–72. doi: 10.1007/s12022-017-9495-2

PubMed Abstract | Crossref Full Text | Google Scholar

18. Chang Z, Li R, Zhang J, An L, Zhou G, Lei M, et al. Distinct immune and inflammatory response patterns contribute to the identification of poor prognosis and advanced clinical characters in bladder cancer patients. Front Immunol. (2022) 13:1008865. doi: 10.3389/fimmu.2022.1008865

PubMed Abstract | Crossref Full Text | Google Scholar

19. Nirmal JA, Maliga Z, Vallius T, Quattrochi B, Chen AA, Jacobson AC, et al. The spatial landscape of progression and immunoediting in primary melanoma at single-cell resolution. Cancer Discov. (2022) 12:1518–41. doi: 10.1158/2159-8290.CD-21-1357

PubMed Abstract | Crossref Full Text | Google Scholar

20. Cunha LL, Marcello AM, Vassallo J, and Ward SL. Differentiated thyroid carcinomas and their B7H1 shield. Future Oncol. (2013) 9:1417–9. doi: 10.2217/fon.13.89

PubMed Abstract | Crossref Full Text | Google Scholar

21. Cunha LL, Marcello AM, Morari CE, Nonogaki S, Conte FF, Gerhard R, et al. Differentiated thyroid carcinomas may elude the immune system by B7H1 upregulation. Endocr-Relat Cancer. (2013) 20:103–10. doi: 10.1530/ERC-12-0313

PubMed Abstract | Crossref Full Text | Google Scholar

22. Bongiovanni M, Rebecchini C, Saglietti C, Bulliard J-L, Marino L, Leval DL, et al. Very low expression of PD-L1 in medullary thyroid carcinoma. Endocr-Relat Cancer. (2017) 24:L35–L8. doi: 10.1530/ERC-17-0104

PubMed Abstract | Crossref Full Text | Google Scholar

23. Karkheiran B, Jahanbin B, Shahverdi G, and Soleimani V. Evaluation of PD-L1 antigen expression using immunohistochemistry technique in medullary thyroid carcinoma samples. Int J Physiol Pathophysiol Pharmacol. (2023) 23:12–20.

PubMed Abstract | Google Scholar

24. Kemal Y, Caliskan S, Gun S, and Kefeli M. Pd-l1 expression in medullary thyroid carcinoma and its association with clinicopathological findings. Turkish J Pathol. (2021) 38:106–13. doi: 10.5146/tjpath.2021.01558

PubMed Abstract | Crossref Full Text | Google Scholar

25. Bai Y, Guo T, Niu D, Zhu Y, Ren W, Yao Q, et al. Clinical significance and interrelations of PD-L1 expression, Ki-67 index, and molecular alterations in sporadic medullary thyroid carcinoma from a Chinese population. Virchows Archiv. (2022) 481:903–11. doi: 10.1007/s00428-022-03390-9

PubMed Abstract | Crossref Full Text | Google Scholar

26. Shi S, Ye L, Yu X, Jin K, and Wu W. Focus on mast cells in the tumor microenvironment: Current knowledge and future directions. Biochim Biophys Acta (BBA) - Rev Cancer. (2023) 1878:188845. doi: 10.1016/j.bbcan.2022.188845

PubMed Abstract | Crossref Full Text | Google Scholar

27. Galli JS, Grimbaldeston M, and Tsai M. Immunomodulatory mast cells: negative, as well as positive, regulators of immunity. Nat Rev Immunol. (2008) 8:478–86. doi: 10.1038/nri2327

PubMed Abstract | Crossref Full Text | Google Scholar

28. Gregory DG, Robbie-Ryan M, Secor HV, Sabatino JJ, and Brown AM. Mast cells are required for optimal autoreactive T cell responses in a murine model of multiple sclerosis. Eur J Immunol. (2005) 35:3478–86. doi: 10.1002/eji.200535271

PubMed Abstract | Crossref Full Text | Google Scholar

29. Takanami I, Takeuchi K, and Naruke M. Mast cell density is associated with angiogenesis and poor prognosis in pulmonary adenocarcinoma. Cancer. (2000) 88:2686–92. doi: 10.1002/1097-0142(20000615)88:12<2686::AID-CNCR6>3.0.CO;2-6

PubMed Abstract | Crossref Full Text | Google Scholar

30. Yodavudh S, Tangjitgamol S, and Puangsa-art S. Prognostic significance of microvessel density and mast cell density for the survival of Thai patients with primary colorectal cancer. J Med Assoc Thailand. (2008) 91:5:723–32.

PubMed Abstract | Google Scholar

31. Ju MJ, Qiu SJ, Gao Q, Fan J, Cai MY, Li YW, et al. Combination of peritumoral mast cells and T-regulatory cells predicts prognosis of hepatocellular carcinoma. Cancer Sci. (2009) 100:1267–74. doi: 10.1111/j.1349-7006.2009.01182.x

PubMed Abstract | Crossref Full Text | Google Scholar

32. Melillo MR, Guarino V, Avilla E, Galdiero RM, Liotti F, Prevete N, et al. Mast cells have a protumorigenic role in human thyroid cancer. Oncogene. (2010) 29:6203–15. doi: 10.1038/onc.2010.348

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: medullary thyroid carcinoma (MTC), immune microenvironment, PD-L1 expression, tumor-infiltrating lymphocytes (TILs), mast cells

Citation: de Castro ME, de Siqueira GF, Mariano JF, Kizys MML, Ceolin L, Soares FA, Natal R, Carneiro H, Ramos RN, Ward LS, Camara NOS, Pinto Camacho C, Facuri Valente FdO, Lindsey SC, dos Santos DD, Gil CD, Martins JRM, Maciel RMdB and Cunha LL (2025) Spatial immune profiling reveals distinct microenvironments in medullary thyroid carcinoma. Front. Immunol. 16:1579205. doi: 10.3389/fimmu.2025.1579205

Received: 18 February 2025; Accepted: 29 April 2025;
Published: 21 May 2025.

Edited by:

Angeliki Chorti, Aristotle University of Thessaloniki, Greece

Reviewed by:

Zhen Dong, Southwest University, China
Dmitry Aleksandrovich Zinovkin, Gomel State Medical University, Belarus

Copyright © 2025 de Castro, de Siqueira, Mariano, Kizys, Ceolin, Soares, Natal, Carneiro, Ramos, Ward, Camara, Pinto Camacho, Facuri Valente, Lindsey, dos Santos, Gil, Martins, Maciel and Cunha. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Rui Monteiro de Barros Maciel, cnVpLm1hY2llbEB1bmlmZXNwLmJy

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