Journal of Immunology Research
Volume 2021 (2021), Article ID 9975423, 7 pages
https://doi.org/10.1155/2021/9975423

Research Article

Hyperion Image Analysis Depicts a Preliminary Landscape of Tumor Immune Microenvironment in OSCC with Lymph Node Metastasis

Shang Xie,¹ Xin-Yuan Zhang,¹ Xiao-Feng Shan,¹ Vicky Yau,² Jian-Yun Zhang,³ Wei Wang,¹ Yong-Pan Yan,⁴ and Zhi-Gang Cai¹

¹Department of Oral and Maxillofacial Surgery, Peking University School and Hospital of Stomatology & National Clinical Research Center for Oral Diseases & National Engineering Laboratory for Digital and Material Technology of Stomatology & Beijing Key Laboratory of Digital Stomatology, Beijing 100081, China
²Stony Brook University School of Dental Medicine, South Drive, Stony Brook, NY 11794, USA
³Department of Oral Pathology, Peking University School and Hospital of Stomatology, Beijing 100081, China
⁴Gencode Diagnostics Inc., #7 Liye Road, 5th Floor, Changping District, Beijing 102206, China

Correspondence should be addressed to Zhi-Gang Cai; c2013xs@163.com

Received 29 March 2021; Revised 20 May 2021; Accepted 4 June 2021; Published 21 June 2021

Academic Editor: Xiao-Jie Lu

Copyright © 2021 Shang Xie et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. Oral squamous cell carcinoma (OSCC) constitutes the most common types of oral cancer. Because its prognosis varies significantly, identification of a tumor immune microenvironment could be a critical tool for treatment planning and predicting a more accurate prognosis. This study is aimed at utilizing the Hyperion imaging system to depict a preliminary landscape of the tumor immune microenvironment in OSCC with lymph node metastasis. Methods. We collected neoplasm samples from OSCC patients. Their formalin-fixed, paraffin-embedded (FFPE) tissue sections were obtained and stained utilizing a panel of 26 clinically relevant metal-conjugated antibodies. Detection and analysis were performed for these stained cells with the Hyperion imaging system. Results. Four patients met our inclusion criteria. We depicted a preliminary landscape of their tumor immune microenvironment and identified 25 distinct immune cell subsets from these OSCC patients based on phenotypic similarity. All these patients had decreased expression of CD8⁺ T cells in tumor specimens. Variety in cell subsets was seen, and more immune activated cells were found in patient A and patient B than those in patient C and patient D. Such differences in tumor immune microenvironments can contribute to forecasting of individual prognoses. Conclusion. The Hyperion imaging system helped to delineate a preliminary and multidimensional landscape of the tumor immune microenvironment in OSCC with lymph node metastasis and provided insights into the influence of the immune microenvironment in determination of prognoses. These results reveal possible contributory factors behind different prognoses of OSCC patients with lymph node metastasis and provide reference for individual treatment planning.

1. Introduction

Oral squamous cell carcinoma (OSCC) contributes the most to neoplasms in the oral and maxillofacial region. Its high mortality leads to a heavy burden to the global health system [1]. Despite rapid advancement in treatments, assessment of accurate prognoses for OSCC has entered a bottleneck period and the five-year overall rate has not significantly ameliorated in the past two decades [2, 3]. One potential reason could be the inability of current treatments to change the existing tumor immune microenvironment; the TNM staging system, on which treatment planning is mostly based, remains inaccurate. In clinical practice, it is indeed common for a small part of OSCC patients with lymph node metastasis to have good prognoses. One possible reason might be that the intrinsic heterogeneity of the tumor immune microenvironment has been overlooked by the TNM system for the mismatch between TMN stage classification and clinical advancement. Thus, there is an important clinical need to identify tools for accurate prognosis evaluation and new treatment strategy formulation, especially for patients with lymph node metastasis.

The tumor microenvironment (TME) is comprised of the tumor cells, surrounded by immune cells, fibroblasts, extracellular matrix, signaling molecules, etc. [4, 5]. Tumor cells and the components of TME, including CD8⁺ T cells [6, 7], CD4⁺ T cells [8, 9], NK cells [10], dendritic cells [11, 12], B cells [13, 14], macrophages [15], FOXP3 [16], and Collagen I [17] have close correlations to tumor treatment and prognosis. In past decades, due to the limitations of tools and clinical methodologies, it was difficult to analyze multiple components and assess tumor heterogeneity of the TME. Yet, with the newly developed Hyperion imaging system, we are now able to include such components to determine more accurate prognoses [18].

In this study, we depict a preliminary and multidimensional landscape of the tumor immune microenvironment of OSCC patients with lymph node metastasis and provide insights into its influence in prognosis by utilizing the Hyperion imaging system.

2. Materials and Methods

2.1. Sample Selection

Medical records of patients admitted for operation to the Peking University Hospital of Stomatology between October, 2018 and April, 2019 were reviewed. In order to focus on specific stages of OSCC (T1-T2, lymph node metastasis) and decrease the clinical heterogeneity, the following criteria were established: (1) 25 to 45 years old, (2) Asian, (3) tumor size less than 2 cm, (4) tumors were located at the tongue, (5) pathological diagnosis were SCC, and (6) tumors have well-moderately differentiated SCC and with lymph node metastasis. No patient has undergone radiotherapy, chemotherapy, or immunotherapy before operation. Tumor tissues were obtained from the Department of Pathology, Peking University Hospital of Stomatology. The included patients were followed up on April 19th, 2021.

2.2. Panel and Antibody Conjugation

The antibody panel was designed with 26 protein markers, as described in our previous publication in which descriptions of these antibodies, metal labels, clones, catalog, and vendor can be found [19]. The panel information can be also found in Supplementary S1. Conjugation of metal-conjugated antibodies was performed utilizing a Maxpar labeling kit (Fluidigm, USA).

2.3. Tissue Antibody Labeling

Tissue samples were obtained and fixed in formalin and embedded in paraffin. Protocols of tissue antibody labeling were similar to those described in the published articles [19, 20].

2.4. Hyperion Tissue Imager Scan

The largest square area of an imager scan is . Images were acquired and analyzed utilizing a Hyperion imaging system (Fluidigm) [19, 21, 22].

2.5. Data Procedure and Normalization

Data were divided into single cells using CellProfiler v3.1.8 [23]. Detailed data procedures were similar to those in the published article [19]. For t-SNE and PhenoGraph, the data were normalized by Harmony [24].

2.6. Clustering and Barnes–Hut t-SNE

Single cells were clustered into groups based on their phenotypical similarity and molecular markers utilizing PhenoGraph [25]. For visualization, high-dimensional data were decreased to two dimensions utilizing the algorithm t-SNE [26–28]. We used the Barnes–Hut achievement of t-SNE to Harmony-normalized data with default parameters ().

3. Results

3.1. Characteristics of OSCC Patients

After reviewing, four OSCC patients met our inclusion criteria. Their clinical characteristics, treatment strategies, and follow-up information are shown in Table 1. All four patients underwent extensive resection of tongue neoplasm and neck lymph node dissection. During follow-up, one patient (patient C, female) died of tumor recurrence at 21 months after the operation and this patient had no history of smoking and drinking.

Table 1: Clinicopathological characteristics of the 4 included patients.

Study NO
Demographic information
Tumor site
Tumor size
TNM stage
Smoking
Alcohol
Pathological features
Treatment strategy
Follow-up (months)

Patient A
Male, 36 years, Asian
Tongue
Ulceration, diameter 2.0 cm
T2N1M0
Yes
Yes
Well-moderately differentiated SCC, lymph node metastasis, 1.0 cm depth of invasion, invasion of striated muscles
Operation: extensive resection of tongue cancer and cervical lymph node dissection; radiotherapy but without chemotherapy
No recurrence, 25 months

Patient B
Male, 39 years, Asian
Tongue
Ulceration, diameter 1.3 cm
T1N1M0
Yes
No
Well differentiated SCC, lymph node metastasis, 0.2 cm depth of invasion
Operation: extensive resection of tongue cancer and cervical lymph node dissection; without radiotherapy and chemotherapy
No recurrence, 24 months

Patient C
Female, 41 years, Asian
Tongue
Ulceration, diameter 1.3 cm
T2N1M0
No
No
Well differentiated SCC, lymph node metastasis, 0.6 cm depth of invasion, invasion of striated muscles and nerves
Operation: extensive resection of tongue cancer and cervical lymph node dissection; without radiotherapy and chemotherapy
Death, 21 months

Patient D
Female, 29 years, Asian
Tongue
Ulceration, diameter 1.8 cm
T2N1M0
No
No
Well differentiated SCC, lymph node metastasis, 1.0 cm depth of invasion
Operation: extensive resection of tongue cancer and cervical lymph node dissection; without radiotherapy and chemotherapy
No recurrence, 26 months

TNM stage: T = tumor; N = node; M = metastasis. SCC: squamous cell carcinoma.

3.2. Sample Procedure and Data Acquisition

All sample slices were marked by a panel of 26 metal-labeled antibodies. Labeling antibodies were detected, and the data were analyzed by the Hyperion imaging system. The scheme of the overall study design is shown in Figure 1.

Figure 1: Overall study design. Hyperion imaging system was utilized to detect OSCC tissue samples, and the output data were used for cluster analyses. TSCC: tongue squamous cell carcinoma; FFPE: formalin-fixed, paraffin-embedded.

3.3. Analysis of Immune Phenotypes of OSCC Samples by Hyperion Imaging System

These metal antibody-labeled slices were scanned to identify different markers of tumor cells, lymphocyte types, vasculature, and structures of tissue cells by the Hyperion imaging system. To visualize the characteristics of TME among the four OSCC patients (patients A, B, C, and D), a t-SNE map was firstly utilized to illustrate single cells from each PhenoGraph cluster according to samples (Figure 2(a)). This panel showed that the components in TME were obviously distinctive among these patients despite similar clinicopathological characteristics, indicating differences in their future prognoses.

(a)

(b)

Figure 2: Cluster analysis of OSCC samples. (a) Single cells from each PhenoGraph cluster in four samples of OSCC with lymph node metastasis displayed by t-SNE descending dimension map were colored-coded in heat map images. The colors represent different samples. (b) t-SNE descending dimension map displays single cells from each PhenoGraph cluster. They are identified in heat map images colored according to individual cluster. Colors represent different clusters of immune cells. OSCC: oral squamous cell carcinoma.

To visualize the panorama of immune cell subsets and extracellular matrix among OSCC patients with lymph node metastasis, we also showed a t-SNE map according to different clusters of cell subsets. In this panorama, we identified 25 unequal cell clusters in TME based on phenotypic similarity (Figure 2(b)).

To further explore these cell clusters, we further visualized these markers via heat map. Clusters and tumor markers were grouped based on expression profiles, and the result is shown in Figure 3(a). From the heat map, we inferred 25 different cell clusters representing the CD8⁺ T cell subsets, CD4⁺ T cell subsets, Collagen I, CD19⁺ B cell subsets, CD56⁺ NK cell subsets, and FOXP3 cell subsets.

(a)

(b)

Figure 3: Heat map of OSCC samples. (a) The heat map shows the differential expression of the panel markers of 25 subsets. Clusters and markers are grouped according to expression profiles. Certain cell clusters are identified as known immune cell types according to typically expressed markers. (b) Differential cluster analyses of the 4 OSCC samples with lymph node metastasis. OSCC: oral squamous cell carcinoma.

Based on differences in cell clusters, cluster and patient analyses were performed (Figure 3(b)). From Figure 3(b), we found that cell clusters had a significant disparity among different patients, but all CD8⁺ cell populations were found to be underexpressed in all included patients, implying that all of them had an incomplete tumor immune microenvironment.

Notably, CD56⁺ NK cells were highly expressed in patient A and CD4⁺ in patient B; FOXP3 and CD19 were close to completely absent in both. Collagen I was highly expressed with less expression of CD4⁺ cells in patient C, and FOXP3 and CD19 were highly expressed while CD4⁺ cells were absent in patient D. It suggested the existence of more immunosuppressive cells and less immune activation cells in patient C and patient D.

Taken together, it was hinted that patient A and patient B might have better prognoses than patient C and patient D, reflecting that a thorough classification of OSCC patients might be beneficial for accurate prognostic assessment of OSCC and implementation of precisely targeted therapies.

4. Discussion

Previous studies demonstrated that tumor cells, immune cells, signaling molecules, and extracellular matrix have mutual interactions in TME [29, 30]. Therefore, understanding differences of tumor immune microenvironment components among patients would provide pivotal information for evaluating prognosis and identifying individualized treatments. This study is first of a new focus to depict a landscape of TME in OSCC patients with lymph node metastasis using the Hyperion imaging system.

The Hyperion imaging system is a new scan and analytical tool for analyses of heterogeneity in TME, and it allows simultaneous detection of dozens of metal-labeled antibodies in the same tissue slices. Currently, it has been successfully used to show the components of many tumor microenvironments [22, 31, 32] and explore the immune response of COVID-19 patients [20]. However, to date, there has been no publication on the use of the Hyperion imaging system to scan and analyze the heterogeneity of OSCC with lymph node metastasis. Therefore, our study is aimed at filling in the gap by exploring and analyzing the intratumoral heterogeneity of OSCC with lymph node metastasis, and thereby suggesting strategies for prognosis and precision treatment planning.

In this study, we depicted a landscape of TME according to a panel of 26 antibodies. In this landscape, we identified several critical immune cell phenotypes including CD8⁺ T cells, CD4⁺ T cells, CD56⁺ NK cells, CD19⁺ B cells, FOXP3⁺ cells, and Collagen I cell subsets.

Previous publications had demonstrated that cytotoxic CD8⁺ T cells targeted and killed tumor cells [6, 7, 33]. In our study, we found low expression of CD8⁺ T cells in all included patients, suggesting their poorly activated immune systems. The prediction of poor prognoses in general for these patients is consistent with the publication of Fang et al. [34].

As we all know, CD4⁺ T cells are important for the activation and expansion of CD8⁺ effectors, and also vital for the maintenance and generation of CD8⁺ T cell memory [8]. In the TME, they might also participate to amplify the CD8⁺ T cell response and have the ability to control tumor growth in different ways [8]. In our study, we found that CD4⁺ T cells were highly expressed in patient B but poorly expressed in patient C and patient D, implying that patient C and patient D had lower immunity.

Natural killer cells (NK cells) are also important immune cells, and many authors proposed that the presence of NK cells was correlated to favorable prognosis of OSCC. For example, Ikeda et al. reported that even a small quantity of NK cells were still helpful to determine prognosis of advanced OSCC [35]. In our study, high expression of CD56 NK cells in patient A compared to that in patient C and patient D suggested a better relative prognosis in the patient.

Classically, B cells are deemed to positively regulate immune responses and boost T cell activation and proliferation via antigen presentation [36]. However, numerous studies recently reported that B cells were relevant to unfavorable prognosis and advancement of lymph node metastasis [14, 37]. In this study, B cells were more highly expressed in patient D but poorly expressed in both patient A and patient B. Besides, the presence of high numbers of FOXP3⁺ T cells in cancer tissues predicts worse relapse-free survival and decreased overall survival [16]. Also highly expressed in patient D, FOXP3⁺ might contribute to worse prognosis in this patient comparing to those of patient A and patient B.

Previous publications have reported that Collagen I increased tumorigenesis, promoted cancer metastasis by upregulating N-cadherin expression and increasing motility [38], and promoted epithelial-to-mesenchymal transition by TGF-β signaling [39]. In our study, Collagen I was highly expressed in patient C but no obvious expressions were seen in other patients, suggesting that patient C had the worse immunosuppressive TME. Notably, patient C died of tumor metastasis at 21 months after the surgical operation, potentially caused by her poor tumor microenvironment.

Based on these data, we could deduce that although all patients demonstrated similar clinical characteristics, patient A and B might have potentially better prognoses due to the presence of more immune activated cells and less immunosuppressive cells in their tumor microenvironment.

The current study is a very early study to analyze heterogeneity and depict the landscape of the tumor immune microenvironment of OSCC using the Hyperion imaging system. Hence, further studies would definitely be required to address the following limitations. Firstly, since current data were based on limited numbers of patients, the expression tendency for the abovementioned markers and cell subsets would need further investigation to look for consistency. Secondly, the phenotypes and immune cell clusters were displayed according to bioinformatics analysis and have not been identified by cell and animal experiments. Besides, the fee for utilizing the Hyperion imaging system could be very challenging to investigate individual cases. Based on these limitations, further studies are required to increase the sample size and identify the roles of these phenotypes with a low-cost approach in the future.

5. Conclusion

The Hyperion imaging system depicted a preliminary and multidimensional immune microenvironment landscape for OSCC with lymph node metastasis and provided insights into the influence of the heterogeneity and variable immune cell subsets on determination of prognoses. These results reveal possible factors which contribute to different prognoses in spite of the similar lymph node involvement and clinical characteristics. Furthermore, long term follow-up will be needed in the future, and further study for targeted therapy will greatly improve individualization of OSCC treatments.

Data Availability

Data are available from corresponding author Zhi-Gang Cai (c2013xs@163.com) or first author Shang Xie (xs2013@hsc.pku.edu.cn).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

All authors contributed to this study. Zhi-Gang Cai designed, edited, and revised the article. Shang Xie, Xin-Yuan Zhang, and Xiao-Feng Shan selected the patients. Xin-Yuan Zhang collected the patients’ clinical data and followed up all included patients. Shang Xie performed research, analyzed data, and drafted and revised the manuscript. Vicky Yau provided good comments and edited and polished the manuscript. Jian-Yun Zhang and Wei Wang made FFPE sections. Yong-Pan Yan provided good comments and analyzed data. All the authors approved the final manuscript. Shang Xie and Xin-Yuan Zhang contributed equally to this manuscript.

Acknowledgments

This research was partly funded by the Peking University Hospital of Stomatology Young Foundation (YS020219) and the National Natural Science Foundation of China (82002878 and 81870781).

Supplementary Materials

S1: antibodies targeting tumor immune environment of OSCC. (Supplementary Materials)

References

Z. H. Ren, C. Y. Hu, H. R. He, Y. J. Li, and J. Lyu, “Global and regional burdens of oral cancer from 1990 to 2017: results from the global burden of disease study,” Cancer Communications, vol. 40, no. 2-3, pp. 81–92, 2020.

B. A. van Dijk, M. T. Brands, S. M. Geurts, M. A. Merkx, and J. L. Roodenburg, “Trends in oral cavity cancer incidence, mortality, survival and treatment in the Netherlands,” International Journal of Cancer, vol. 139, no. 3, pp. 574–583, 2016.

J. D. S. Moro, M. C. Maroneze, T. M. Ardenghi, L. M. Barin, and C. C. Danesi, “Oral and oropharyngeal cancer: epidemiology and survival analysis,” Einstein (São Paulo), vol. 16, no. 2, article eAO4248, 2018.

J. A. Joyce and D. T. Fearon, “T cell exclusion, immune privilege, and the tumor microenvironment,” Science (New York, NY), vol. 348, no. 6230, pp. 74–80, 2015.

W. H. Fridman, F. Pages, C. Sautes-Fridman, and J. Galon, “The immune contexture in human tumours: impact on clinical outcome,” Nature Reviews Cancer, vol. 12, no. 4, pp. 298–306, 2012.

D. Ostroumov, N. Fekete-Drimusz, M. Saborowski, F. Kuhnel, and N. Woller, “CD4 and CD8 T lymphocyte interplay in controlling tumor growth,” Cellular and Molecular Life Sciences: CMLS, vol. 75, no. 4, pp. 689–713, 2018.

S. Shimizu, H. Hiratsuka, K. Koike et al., “Tumor-infiltrating CD8(+) T-cell density is an independent prognostic marker for oral squamous cell carcinoma,” Cancer Medicine, vol. 8, no. 1, pp. 80–93, 2019.

J. R. Richardson, A. Schöllhorn, C. Gouttefangeas, and J. Schuhmacher, “CD4+ T cells: multitasking cells in the duty of cancer immunotherapy,” Cancers, vol. 13, no. 4, p. 596, 2021.

D. Y. Oh, S. S. Kwek, S. S. Raju et al., “Intratumoral CD4⁺ T cells mediate anti-tumor cytotoxicity in human bladder cancer,” Cell, vol. 181, no. 7, pp. 1612–1625.e13, 2020.

A. Jewett, J. Kos, Y. Fong et al., “NK cells shape pancreatic and oral tumor microenvironments; role in inhibition of tumor growth and metastasis,” Seminars in Cancer Biology, vol. 53, pp. 178–188, 2018.

N. C. Jung, J. H. Lee, K. H. Chung, Y. S. Kwak, and D. S. Lim, “Dendritic cell-based immunotherapy for solid tumors,” Translational Oncology, vol. 11, no. 3, pp. 686–690, 2018.

R. L. Sabado, S. Balan, and N. Bhardwaj, “Dendritic cell-based immunotherapy,” Cell Research, vol. 27, no. 1, pp. 74–95, 2017.

S. Li, C. Huang, G. Hu et al., “Tumor-educated B cells promote renal cancer metastasis via inducing the IL-1β/HIF-2α/Notch1 signals,” Cell Death & Disease, vol. 11, no. 3, p. 163, 2020.

Y. Gu, Y. Liu, L. Fu et al., “Tumor-educated B cells selectively promote breast cancer lymph node metastasis by HSPA4-targeting IgG,” Nature Medicine, vol. 25, no. 2, pp. 312–322, 2019.

E. Hadler-Olsen and A. M. Wirsing, “Tissue-infiltrating immune cells as prognostic markers in oral squamous cell carcinoma: a systematic review and meta-analysis,” British Journal of Cancer, vol. 120, no. 7, pp. 714–727, 2019.

J. Shou, Z. Zhang, Y. Lai, Z. Chen, and J. Huang, “Worse outcome in breast cancer with higher tumor-infiltrating FOXP3+ Tregs: a systematic review and meta-analysis,” BMC Cancer, vol. 16, no. 1, 2016.

C. E. Barcus, K. A. O’Leary, J. L. Brockman et al., “Elevated collagen-I augments tumor progressive signals, intravasation and metastasis of prolactin-induced estrogen receptor alpha positive mammary tumor cells,” Breast Cancer Research: BCR, vol. 19, no. 1, p. 9, 2017.

H. W. Jackson, J. R. Fischer, V. R. T. Zanotelli et al., “The single-cell pathology landscape of breast cancer,” Nature, vol. 578, no. 7796, pp. 615–620, 2020.

S. Xie, X.-F. Shan, V. Yau et al., “Hyperion imaging system reveals heterogeneous tumor microenvironment of oral squamous cell carcinoma patients at T1N0M0 stage,” Annals of Translational Medicine, vol. 8, article e1513, 2020.

C. Wang, J. Xu, S. Wang et al., “Imaging mass cytometric analysis of postmortem tissues reveals dysregulated immune cell and cytokine responses in multiple organs of COVID-19 patients,” Frontiers in Microbiology, vol. 11, p. 600989, 2020.

Q. Chang, O. I. Ornatsky, I. Siddiqui, A. Loboda, V. I. Baranov, and D. W. Hedley, “Imaging mass cytometry,” Cytometry Part A: The Journal of the International Society for Analytical Cytology, vol. 91, no. 2, pp. 160–169, 2017.

B. Allo, X. Lou, A. Bouzekri, and O. Ornatsky, “Clickable and high-sensitivity metal-containing tags for mass cytometry,” Bioconjugate Chemistry, vol. 29, no. 6, pp. 2028–2038, 2018.

A. E. Carpenter, T. R. Jones, M. R. Lamprecht et al., “CellProfiler: image analysis software for identifying and quantifying cell phenotypes,” Genome Biology, vol. 7, no. 10, p. R100, 2006.

I. Korsunsky, N. Millard, J. Fan et al., “Fast, sensitive and accurate integration of single-cell data with Harmony,” Nature Methods, vol. 16, no. 12, pp. 1289–1296, 2019.

J. H. Levine, E. F. Simonds, S. C. Bendall et al., “Data-driven phenotypic dissection of AML reveals progenitor-like cells that correlate with prognosis,” Cell, vol. 162, no. 1, pp. 184–197, 2015.

L. van der Maaten and G. Hinton, “Visualizing data using t-SNE,” Journal of Machine Learning Research, vol. 9, pp. 2579–2605, 2008.

M. Lvd, “Barnes-Hut-SNE,” Computer Science, vol. 1301, pp. 1–11, 2013.

D. Kobak and P. Berens, “The art of using t-SNE for single-cell transcriptomics,” Nature Communications, vol. 10, no. 1, p. 5416, 2019.

V. Thorsson, D. L. Gibbs, S. D. Brown et al., “The immune landscape of cancer,” Immunity, vol. 48, no. 4, pp. 812–830.e14, 2018.

B. Qiao, J. Huang, Z. Mei, A. K. Lam, J. Zhao, and L. Ying, “Analysis of immune microenvironment by multiplex immunohistochemistry staining in different oral diseases and oral squamous cell carcinoma,” Frontiers in Oncology, vol. 10, p. 555757, 2020.

E. Gerdtsson, M. Pore, J. A. Thiele et al., “Multiplex protein detection on circulating tumor cells from liquid biopsies using imaging mass cytometry,” Convergent Science Physical Oncology, vol. 4, no. 1, p. 015002, 2018.

R. Elaldi, P. Hemon, L. Petti et al., “High dimensional imaging mass cytometry panel to visualize the tumor immune microenvironment contexture,” Frontiers in Immunology, vol. 12, p. 666233, 2021.

W. Wang, M. Green, J. E. Choi et al., “CD8⁺ T cells regulate tumour ferroptosis during cancer immunotherapy,” Nature, vol. 569, no. 7755, pp. 270–274, 2019.

J. Fang, X. Li, D. Ma et al., “Prognostic significance of tumor infiltrating immune cells in oral squamous cell carcinoma,” BMC Cancer, vol. 17, no. 1, p. 375, 2017.

T. Ikeda, S. Seki, M. Fujiwara et al., “Low-risk population among patients with tumor-node-metastasis stage III/IV oral squamous cell carcinoma,” Oncology Letters, vol. 14, no. 3, pp. 3711–3716, 2017.

A. Sarvaria, J. A. Madrigal, and A. Saudemont, “B cell regulation in cancer and anti-tumor immunity,” Cellular & Molecular Immunology, vol. 14, no. 8, pp. 662–674, 2017.

H. P. Dong, M. B. Elstrand, A. Holth et al., “NK- and B-cell infiltration correlates with worse outcome in metastatic ovarian carcinoma,” American Journal of Clinical Pathology, vol. 125, no. 3, pp. 451–458, 2006.

Y. Shintani, M. A. Hollingsworth, M. J. Wheelock, and K. R. Johnson, “Collagen I promotes metastasis in pancreatic cancer by activating c-Jun NH (2)-terminal kinase 1 and up-regulating N-cadherin expression,” Cancer Research, vol. 66, no. 24, pp. 11745–11753, 2006.

Y. Shintani, M. Maeda, N. Chaika, K. R. Johnson, and M. J. Wheelock, “Collagen I promotes epithelial-to-mesenchymal transition in lung cancer cells via transforming growth factor-beta signaling,” American Journal of Respiratory Cell and Molecular Biology, vol. 38, no. 1, pp. 95–104, 2008.