Front. Oncol. Frontiers in Oncology Front. Oncol. 2234-943X Frontiers Media S.A. 10.3389/fonc.2021.771749 Oncology Original Research Identification of Immune-Cell-Related Prognostic Biomarkers of Esophageal Squamous Cell Carcinoma Based on Tumor Microenvironment Cui Yiyao 1 Hou Ruiqin 2 Lv Xiaoshuo 1 Wang Feng 1 Yu Zhaoyan 3 * Cui Yong 1 * 1 Department of Thoracic Surgery, Beijing Friendship Hospital, Affiliated to the Capital University of Medical Sciences, Beijing, China 2 Department of Blood Transfusion, Peking University People’s Hospital, Beijing, China 3 Department of Otorhinolaryngology, Shandong Public Health Clinical Center, Jinan, China

Edited by: Meng Zhou, Wenzhou Medical University, China

Reviewed by: Chao Song, Harbin Medical University, China; Dechao Bu, Institute of Computing Technology (CAS), China

*Correspondence: Zhaoyan Yu, yuzhaoyan@163.com; Yong Cui, cuiyong@ccmu.edu.cn

This article was submitted to Cancer Genetics, a section of the journal Frontiers in Oncology

 25 10 2021 2021 11 771749 07 09 2021 04 10 2021 Copyright © 2021 Cui, Hou, Lv, Wang, Yu and Cui 2021 Cui, Hou, Lv, Wang, Yu and Cui

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.

 Background

Esophageal squamous cell carcinoma (ESCC) is one of the most fatal cancers in the world. The 5-year survival rate of ESCC is <30%. However, few biomarkers can accurately predict the prognosis of patients with ESCC. We aimed to identify potential survival-associated biomarkers for ESCC to improve its poor prognosis.

 Methods

ImmuneAI analysis was first used to access the immune cell abundance of ESCC. Then, ESTIMATE analysis was performed to explore the tumor microenvironment (TME), and differential analysis was used for the selection of immune-related differentially expressed genes (DEGs). Weighted gene coexpression network analysis (WGCNA) was used for selecting the candidate DEGs. Least absolute shrinkage and selection operator (LASSO) Cox regression was used to build the immune-cell-associated prognostic model (ICPM). Kaplan–Meier curve of survival analysis was performed to evaluate the efficacy of the ICPM.

 Results

Based on the ESTIMATE and ImmuneAI analysis, we obtained 24 immune cells’ abundance. Next, we identified six coexpression module that was associated with the abundance. Then, LASSO regression models were constructed by selecting the genes in the module that is most relevant to immune cells. Two test dataset was used to testify the model, and we finally, obtained a seven-genes survival model that performed an excellent prognostic efficacy.

 Conclusion

In the current study, we filtered seven key genes that may be potential prognostic biomarkers of ESCC, and they may be used as new factors to improve the prognosis of cancer.

 esophageal squamous cell carcinoma LASSO regression WGCNA tumor microenvironment prognostic biomarkers Introduction

There are more than half a million new esophageal cancer cases diagnosed each year, leading it to be one of the most universal cancers (1). Esophageal cancer has two main subtypes: esophageal squamous cell carcinoma (ESCC) and esophageal adenocarcinoma (EAC) (2). In Asia, most esophageal cancer is ESCC and had higher mortality rates than other regions (3). Unfortunately, distal metastases have occurred in a great mass of ESCC patients when diagnosed, and the 5-year survival rate is <30% (1). Hence, it is necessary to identify novel biomarkers as improving the prognosis and providing new targets for the therapy of ESCC.

In recent years, the tumor microenvironment (TME) has attracted public attention as a novel therapeutic strategy (4). TME contains numerous cells and acts a vital role in the development and invasion of cancers (5). With the development of tumor immunology, an in-depth understanding of TME is essential to improve immunotherapy (6, 7). More recently, an algorithm called ESTIMATE was able to calculate the abundance of various immune cells in a tissue based on high-throughput second-generation sequencing data (8, 9). This provides a powerful aid for the mining and analysis of the existing numerous cancer data. There were also many research focus on the TME in cancers (10–12).

In recent years, many algorithms based on network analysis have emerged endlessly (13). A widely used algorithm, called weighted gene coexpression network analysis (WGCNA), could calculate the relationship between genes and patients’ information (14). The advantage of WGCNA is that genes with analogous patterns in expression level can be synthesized into coexpression modules and then associated with clinical characteristics so that genes related to specific traits can be screened by dimensionality reduction (15–17). The WGCNA can analyze data from large samples and establish key genes for further validation after identifying expression modules that are relevant to clinical characteristics (18). WGCNA provides a powerful aid in the search for clinically relevant molecular markers. Here, we screened multiple modules and genes closely related to the tumor microenvironment and immune cells and obtained possible biomarkers, which could improve the prognosis of ESCC ( Figure 1A ).

Assessment of immune cell abundance and the tumor microenvironment in ESCC. (A) Flow chart of this study. (B) Immune cell abundance in the training dataset. The p-value is displayed at the top of the bar. We only selected the significantly different cells for further analysis. (C) Immune cell abundance in the test dataset. (D) Differential analysis of immune score based on ESTIMATE analysis. The volcano plot shows the DEGs in high and low immune-score groups. (E) Enrichment analysis of immune-associated DEGs.

 Methods Data Sources

Three cohorts were used in the study. The training dataset was extracted from the GEO database (GSE53625). It contains 179 ESCC patients’ expression data and their clinical trait information. One of the test datasets was downloaded from the TGGA database. We obtained the RNA-seq data of 84 ESCC patients and corresponding clinical information. The gene expression level was normalized by a log2 conversion. Another test dataset containing 125 ESCC patients with clinical information was downloaded and used as another test dataset (GSE121931).

 Differential Expression Analysis

We used the ESTIMATE algorithm to get the immune scores of ESCC patients (19). Then, the patients were divided into high and low immune‐score groups based on the median of the score. Finally, we used the “limma” R package to perform the differential expression analyses between the two groups. The threshold of differentially expressed genes (DEGs) is adjusted p < 0.05 and |log2 fold change| ≥ 0.585.

We also applied the ImmuCellAI database to predict the abundance of 24 immune cell types in ESCC samples (20). The abundance of these immune cells in high and low immune‐score groups was analyzed and used as the trait of ESCC samples for further analysis.

 Weighted Gene Coexpression Network Analysis

We used the “WGCNA” R package to construct the coexpression network of immune-related DEGs based on the automatic network construction function. Then, gene modules with similar expression pattern hierarchical were detected. Finally, patients’ characteristics of immune cell infiltration were associated with these modules, and the key genes from the candidate module were explored for further analysis.

 Enrichment Analysis

We used the “clusterProfiler” R package for functional enrichment analysis of genes of interest. The background functions including the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and GO terms. In the analysis, only functions with a p < 0.05 will be selected.

 Survival Analysis

The R package glmnet was applied to build an immune cell-associated prognostic model (ICPM) of ESCC patients by LASSO analysis. Subsequently, we used the survival R package to depict the survival curve to estimate the efficiency of the ICPM and the genes in the model. We analyzed the overall survival between different clinical subgroups of ESCC patients based on the risk score, too. The subgroups included age (≤59 or >59 years), gender, alcohol, tobacco, and tumor–node–metastasis (TNM) stage. The pROC package was used to assess the prognostic efficiency of the ICPM. Univariate cox regression analysis was performed to filter the prognosis-associated genes in the ICPM, and a Cox regression model was constructed based on its result.

 Statistical Analysis

All statistical analyses were calculated through R software (version 4.0.3). t-test was used to compare the differences between the selected two groups. Adjusted p < 0.05 was regarded as significant.

 Results Evaluation of Immune Cell Abundance in ESCC Patients and Identification of Immune-Associated DEGs

In the training dataset, we split the 179 ESCC patients into high- and low-immune score groups based on the ESTIMATE analysis. Next, we explored the immune cell abundance of the cohort. It showed that most of the 24 immune cells were different in these groups except CD4-naive, cytotoxic T cell, NKT, DC, and B cell ( Figure 1B ). Then, the test data contains 84 ESCC patients who had performed the same analysis. These immune cells had a similar abundance in the test dataset ( Figure 1C ).

Subsequently, we identified DEGs between these groups and found 1,489 differentially upregulated genes and 213 downregulated genes in the training dataset ( Figure 1D ). The enrichment analysis also showed that these DEGs were involved in the Th17 cell differentiation, Th1 and Th2 cell differentiation, cell adhesion, and so on ( Figure 1E ). It hints that the DEGs participate in the immune-associated functions in ESCC patients.

 Identification of Immune-Cells-Associated Modules

The immune-related DEGs were then utilized to cluster the coexpression modules. We used the pickSoftThreshold function to calculate the cutoff of soft power. Soft power is a key parameter to cluster the modules. In the study, the power threshold was calculated at 7 ( Figure 2A ). When this power is used, the R2 of the scale-free topology model under the soft threshold is 0.82. It hints that the network we constructed conformed to scale-free characteristics ( Figure 2B ). Next, we constructed the network with default module size and built six coexpression modules (green, yellow, turquoise, blue, brown, and red) ( Figure 2C ). Subsequently, we calculated the correlation between these six modules and immune cell abundance to filter the most important modules and crucial genes. It showed that the blue module was significantly associated with most immune cells, especially iTreg and Tfh cells ( Figure 3A ). The blue module contained 240 genes, and the enrichment analysis showed that they are related to the T-cell receptor signaling pathway, cell adhesion, antigen-receptor-mediated signaling pathway, etc. ( Figure 3B ). The enrichment results hinted that these genes are immune-associated in the tumor microenvironment of ESCC.

WGCNA to construct the coexpression modules based on immune-related DEGs. (A, B) The soft power is set at 7 based on the R2 of the scale-free topology model, which reached 0.82. (C) Six modules were clustered based on under the parameter.

Coexpression modules were related to immune cells. (A) Heatmap shows the relationship between coexpression modules and immune cell abundance. (B) Enrichment analysis of candidate DEGs in the blue module.

 Construct the Immune-Cell-Associated Prognostic Model

To select the potential prognostic genes in the blue module, we used LASSO analysis on the 240 genes to generate an ICPM. This is the best model, which contains 13 genes based on LASSO analysis ( Figures 4A, B ). Figure 4C shows the coefficient of the 13 genes in the LASSO model. Subsequently, survival analysis revealed that high-risk-score patients had less overall survival time ( Figure 4D ). We then analyzed the relationship between the patients’ clinical traits and the model to evaluate the ICPM’s efficacy. We found that this model could distinguish the TNM stage ( Figure 4E ). Thus, a stratified survival analysis was used to access the prognostic efficacy in different TNM stages. Interestingly, high-risk-score patients also had a bad prognosis in TNM stages ( Figure 4F ). The result substantiated that the model could be used as a predictor and helps to improve the prognosis of ESCC.

LASSO regression to construct the ICPM and its stratified survival analysis. (A, B) The best criteria to build the model based on LASSO regression. (C) The coefficients of the candidate genes in the survival model. (D) The ICPM has an excellent prognostic efficacy. (E) The ICPM can separate patients with different TNM stages into different subgroups based on risk score. (F) The ICPM could predict the prognosis in patients with different TNM stages. ns, not significant, **p < 0.01.

 Survival Analysis of the 13 Genes in ICPM

The 13 crucial genes in ICPM were then tested by survival analysis. It indicated that most of these genes are prognosis associated in the training dataset except RASGRF1 ( Figure 5 ). To be specific, 4 out of the crucial 13 genes (FFAR4, HTR2A, RASA3, and SH3KBP1) were risk factors of ESCC patients, while the other 8 genes (CD38, CLEC6A, DPP4, STX11, TRAV26-2, TRBV4-2, TRBV5-3, and XCR1) are protective factors.

Survival analysis of the 13 candidate genes in ICPM.

Then, two test datasets were used to test the 13 genes’ survival ability. In the Cancer Genome Atlas (TCGA) test dataset, eight genes were detected, and six of them were survival associated (DPP4, FFAR4, RASA3, RASGRF1, XCR1, and STX11) ( Supplementary Figure 1 ). In another test dataset with 125 ESCC patients, 10 genes were detected, and 6 of them were survival associated (CLEC6A, DPP4, HTR2A, RASA3, SH3KBP1, and XCR1) ( Supplementary Figure 2 ).

 Construct the Ultimate Survival Model Based on the ICPM

After the filtration by two test datasets, we obtained eight stable prognostic genes that were survival associated in at least two datasets ( Figure 6A ). Next, we used univariate Cox regression analysis to select independent prognostic genes. Seven out of the eight crucial genes were independent prognostic factors of overall survival (OS) in ESCC patients ( Figure 6B ). Thus, we used the seven genes to construct a survival model. It showed that the model performed an excellent efficacy in prognosis ( Figure 6C ), and the receiver operating characteristic (ROC) curves also proved the efficacy of the model ( Figure 6D ).

Validation and optimization of ICPM through multidatasets. (A) Stable-survival-associated candidate genes in three cohorts. (B) Cox regression analysis showed that six candidate genes were independent prognostic factors. (C) The optimized survival model has a good effect on ESCC patients’ prognosis. (D) AUC indicates the well-predicted efficacy of the optimized model.

 Discussion

ESCC is one of the most malignant cancers in the world. However, few biomarkers can accurately predict the prognosis of patients with ESCC (21). In addition, some studies have already demonstrated the effect of immune cells on cancers, and the immune cells in ESCC can be applied to access its therapeutic and prognostic effects (10, 22). In the current study, by performing an integrated analysis of immune microenvironment and gene expression pattern, we investigated potential prognostic biomarkers in ESCC base on the following steps: (1) predict the immune cell abundance in ESCC, (2) assess the tumor microenvironment of ESCC and identify differentially expressed genes, (3) enrich DEGs through Gene Ontology and KEGG pathway analysis, (4) construct the coexpression network of immune-related DEGs through WGCNA, (5) associate immune cell infiltration and coexpression modules, (6) construct the survival model, and (7) validate the model and improve its efficacy.

The tumor microenvironment is important for the study of immune-related target molecules and prognostic markers (23, 24). There have been many studies using the tumor microenvironment to search for prognostic markers (25–27). Furthermore, the integration of multiple bioinformatics analyses has been widely used with ESCC (28). Thus, the integration of immune cells and ESCC is feasible in the discovery of prognostic biomarkers. In recent years, immunotherapy has been recognized as a treatment strategy that performs well in many types of cancer (29). It works against cancer cells by inhibiting immune checkpoints. The development of monoclonal antibodies that inhibit programmed death 1 (PD-1) or programmed death-ligand 1 (PD-L1) has also been shown to produce convincing clinical efficacy in a variety of malignancies, including ESCC (30). However, the efficacy of these drugs is limited. Therefore, there is an urgent need to identify new biomarkers associated with immune cells to select patients sensitive to these drugs and thus improve the prognosis of ESCC.

Here, we obtained seven crucial genes that may act as the prognostic factor of ESCC. CLEC6A is a gene of human innate immunity. It can directly mediate intracellular signaling, recognize a variety of endogenous and exogenous ligands, and drive both innate and adaptive immunity (31). DPP4 is an inherent type II transmembrane glycoprotein and serine peptidase involved in glucose and insulin metabolism and immune regulation (32). DPP4 has been found to participate in thyroid papillary carcinoma cell proliferation by inhibiting the mitogen-activated protein kinase (MAPK) pathway (33). FFAR4 encodes G-protein-coupled receptor (GPR) and is involved in anti-inflammatory responses (34). FFAR4 can promote cell proliferation and migration and has been identified as a potential prognostic biomarker for laryngeal cancer (35). HTR2A encodes one of the serotonin receptors and could activate the PI3K-Akt-MTOR signal (36). SH3KBP1 is involved in several cellular processes, such as apoptosis, cytoskeletal rearrangement, and cell adhesion (37). It has been reported to promote tumor proliferation and invasion in ESCC (38). STX11 is involved in intracellular vesicle transport. It has been found to play a tumor-suppressive role in peripheral blood T-cell lymphoma (39). XCR1 signals by increasing intracellular calcium levels. It has been found to promote the migration of non-small cell lung cancer (40). At present, there are no key targets for the prognosis of ESCC. Although these genes are not specifically expressed in ESCC, they have been able to show the potential ability to predict the prognosis of cancer through the integration of the survival model, which is helpful to determine its prognostic targets. Moreover, subsequent studies can further determine the expression levels of these key genes at low throughput levels and further study their important mechanisms in esophageal cancer. In short, our study provides a solid foundation for subsequent prognostic and mechanism research.

All in all, our findings provide several potential prognostic biomarkers of ESCC and may improve the treatment of this type of 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.

 Author Contributions

YiC, ZY, and YoC designed the whole study. YiC and RH performed most of the analysis. XL and FW collected the high throughput data. ZY search the analysis tools. YiC wrote the manuscript. All authors contributed to the article and approved the submitted version.

 Funding

This work was supported by a grant from Wu Jieping Medical Foundation (320.6750.2020-17-4) and project (RDE2019-18) supported by Peking University People’s Hospital Scientific Research Development Funds.

 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.

 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/fonc.2021.771749/full#supplementary-material

Survival analysis of the thirteen candidate genes of ICPM in TCGA dataset.

Survival analysis of the thirteen candidate genes of ICPM in another test dataset.

 References Enzinger PC Mayer RJ . Esophageal Cancer. N Engl J Med (2003) 349(23):2241–52. doi: 10.1056/NEJMra035010 Pennathur A Gibson MK Jobe BA Luketich JD . Oesophageal Carcinoma. Lancet (2013) 381(9864):400–12. doi: 10.1016/S0140-6736(12)60643-6 Brown LM Devesa SS Chow WH . Incidence of Adenocarcinoma of the Esophagus Among White Americans by Sex, Stage, and Age. J Natl Cancer Inst (2008) 100(16):1184–7. doi: 10.1093/jnci/djn211 Yang L Song X Gong T Jiang K Hou Y Chen T . Development a Hyaluronic Acid Ion-Pairing Liposomal Nanoparticle for Enhancing Anti-Glioma Efficacy by Modulating Glioma Microenvironment. Drug Delivery (2018) 25(1):388–97. doi: 10.1080/10717544.2018.1431979 Cavallo F De Giovanni C Nanni P Forni G Lollini PL . 2011: The Immune Hallmarks of Cancer. Cancer Immunol Immunother (2011) 60(3):319–26. doi: 10.1007/s00262-010-0968-0 Najem H Khasraw M Heimberger AB . Immune Microenvironment Landscape in CNS Tumors and Role in Responses to Immunotherapy. Cells (2021) 10(8). doi: 10.3390/cells10082032 Zhao S Peralta RM Avina-Ochoa N Delgoffe GM Kaech SM . Metabolic Regulation of T Cells in the Tumor Microenvironment by Nutrient Availability and Diet. Semin Immunol (2021), 101485. doi: 10.1016/j.smim.2021.101485 Yoshihara K Shahmoradgoli M Martinez E Vegesna R Kim H Torres-Garcia W . Inferring Tumour Purity and Stromal and Immune Cell Admixture From Expression Data. Nat Commun (2013) 4:2612. doi: 10.1038/ncomms3612 Li B Severson E Pignon JC Zhao H Li T Novak J . Comprehensive Analyses of Tumor Immunity: Implications for Cancer Immunotherapy. Genome Biol (2016) 17(1):174. doi: 10.1186/s13059-016-1028-7 Baba Y Nomoto D Okadome K Ishimoto T Iwatsuki M Miyamoto Y . Tumor Immune Microenvironment and Immune Checkpoint Inhibitors in Esophageal Squamous Cell Carcinoma. Cancer Sci (2020) 111(9):3132–41. doi: 10.1111/cas.14541 DeNardo DG Galkin A Dupont J Zhou L Bendell J . GB1275, a First-in-Class CD11b Modulator: Rationale for Immunotherapeutic Combinations in Solid Tumors. J Immunother Cancer (2021) 9(8). doi: 10.1136/jitc-2021-003005 Gu Z Du Y Zhao X Wang C . Tumor Microenvironment and Metabolic Remodeling in Gemcitabine-Based Chemoresistance of Pancreatic Cancer. Cancer Lett (2021) 521:98–108. doi: 10.1016/j.canlet.2021.08.029 Bu D Xia Y Zhang J Cao W Huo P Wang Z . FangNet: Mining Herb Hidden Knowledge From TCM Clinical Effective Formulas Using Structure Network Algorithm. Comput Struct Biotechnol J (2021) 19:62–71. doi: 10.1016/j.csbj.2020.11.036 Langfelder P Horvath S . WGCNA: An R Package for Weighted Correlation Network Analysis. BMC Bioinf (2008) 9:559. doi: 10.1186/1471-2105-9-559 Giulietti M Occhipinti G Righetti A Bracci M Conti A Ruzzo A . Emerging Biomarkers in Bladder Cancer Identified by Network Analysis of Transcriptomic Data. Front Oncol (2018) 8:450. doi: 10.3389/fonc.2018.00450 Luo Z Wang W Li F Songyang Z Feng X Xin C . Pan-Cancer Analysis Identifies Telomerase-Associated Signatures and Cancer Subtypes. Mol Cancer (2019) 18(1):106. doi: 10.1186/s12943-019-1035-x Yuan Y Chen J Wang J Xu M Zhang Y Sun P . Identification Hub Genes in Colorectal Cancer by Integrating Weighted Gene Co-Expression Network Analysis and Clinical Validation In Vivo and Vitro . Front Oncol (2020) 10:638. doi: 10.3389/fonc.2020.00638 Dai Y Sun X Wang C Li F Zhang S Zhang H . Gene Co-Expression Network Analysis Reveals Key Pathways and Hub Genes in Chinese Cabbage (Brassica Rapa L.) During Vernalization. BMC Genomics (2021) 22(1):236. doi: 10.1186/s12864-021-07510-8 Becht E Giraldo NA Lacroix L Buttard B Elarouci N Petitprez F . Estimating the Population Abundance of Tissue-Infiltrating Immune and Stromal Cell Populations Using Gene Expression. Genome Biol (2016) 17(1):218. doi: 10.1186/s13059-016-1070-5 Miao YR Zhang Q Lei Q Luo M Xie GY Wang H . ImmuCellAI: A Unique Method for Comprehensive T-Cell Subsets Abundance Prediction and Its Application in Cancer Immunotherapy. Adv Sci (Weinh) (2020) 7(7):1902880. doi: 10.1002/advs.201902880 Murphy G McCormack V Abedi-Ardekani B Arnold M Camargo MC Dar NA . International Cancer Seminars: A Focus on Esophageal Squamous Cell Carcinoma. Ann Oncol (2017) 28(9):2086–93. doi: 10.1093/annonc/mdx279 Yasuda T Nishiki K Hiraki Y Kato H Iwama M Shiraishi O . Phase II Adjuvant Cancer-Specific Vaccine Therapy for Esophageal Cancer Patients Curatively Resected After Preoperative Therapy With Pathologically Positive Nodes; Possible Significance of Tumor Immune Microenvironment in Its Clinical Effects. Ann Surg (2020). doi: 10.1097/SLA.0000000000003880 Sun J Zhang Z Bao S Yan C Hou P Wu N . Identification of Tumor Immune Infiltration-Associated lncRNAs for Improving Prognosis and Immunotherapy Response of Patients With Non-Small Cell Lung Cancer. J Immunother Cancer (2020) 8(1). doi: 10.1136/jitc-2019-000110 Zhang Z Bao S Yan C Hou P Zhou M Sun J . Computational Principles and Practice for Decoding Immune Contexture in the Tumor Microenvironment. Brief Bioinform (2021) 22(3). doi: 10.1093/bib/bbaa075 Math MV LaBrooy SJ Jacob E Ti TK . Increase in Bile Secretion After Water Ingestion in Cholecystectomised Patients. Jpn J Physiol (1987) 37(4):729–33. doi: 10.2170/jjphysiol.37.729 Lieberman NAP DeGolier K Kovar HM Davis A Hoglund V Stevens J . Characterization of the Immune Microenvironment of Diffuse Intrinsic Pontine Glioma: Implications for Development of Immunotherapy. Neuro Oncol (2019) 21(1):83–94. doi: 10.1093/neuonc/noy145 Job S Rapoud D Dos Santos A Gonzalez P Desterke C Pascal G . Identification of Four Immune Subtypes Characterized by Distinct Composition and Functions of Tumor Microenvironment in Intrahepatic Cholangiocarcinoma. Hepatology (2020) 72(3):965–81. doi: 10.1002/hep.31092 Chen Y Liao LD Wu ZY Yang Q Guo JC He JZ . Identification of Key Genes by Integrating DNA Methylation and Next-Generation Transcriptome Sequencing for Esophageal Squamous Cell Carcinoma. Aging (Albany NY) (2020) 12(2):1332–65. doi: 10.18632/aging.102686 Ribas A Wolchok JD . Cancer Immunotherapy Using Checkpoint Blockade. Science (2018) 359(6382):1350–5. doi: 10.1126/science.aar4060 Janjigian YY Bendell J Calvo E Kim JW Ascierto PA Sharma P . CheckMate-032 Study: Efficacy and Safety of Nivolumab and Nivolumab Plus Ipilimumab in Patients With Metastatic Esophagogastric Cancer. J Clin Oncol (2018) 36(28):2836–44. doi: 10.1200/JCO.2017.76.6212 Kerscher B Willment JA Brown GD . The Dectin-2 Family of C-Type Lectin-Like Receptors: An Update. Int Immunol (2013) 25(5):271–7. doi: 10.1093/intimm/dxt006 Enz N Vliegen G De Meester I Jungraithmayr W . CD26/DPP4 - a Potential Biomarker and Target for Cancer Therapy. Pharmacol Ther (2019) 198:135–59. doi: 10.1016/j.pharmthera.2019.02.015 Hu X Chen S Xie C Li Z Wu Z You Z . DPP4 Gene Silencing Inhibits Proliferation and Epithelial-Mesenchymal Transition of Papillary Thyroid Carcinoma Cells Through Suppression of the MAPK Pathway. J Endocrinol Invest (2021) 44(8):1609–23. doi: 10.1007/s40618-020-01455-7 Stuttgen GM Sahoo D . FFAR4: A New Player in Cardiometabolic Disease? Endocrinology (2021) 162(8). doi: 10.1210/endocr/bqab111 Zhou LZ Cui YX Wang WZ Wu J Sun Z Ma SY . FFAR4 Promotes Cell Proliferation and Migration and Servers as a Potential Biomarker for Clinicopathological Characteristics and Prognosis in Laryngocarcinoma. Eur Rev Med Pharmacol Sci (2019) 23(17):7438–44. doi: 10.26355/eurrev_201909_18853 Gao W Guo N Zhao S Chen Z Zhang W Yan F . HTR2A Promotes the Development of Cardiac Hypertrophy by Activating PI3K-PDK1-AKT-mTOR Signaling. Cell Stress Chaperones (2020) 25(6):899–908. doi: 10.1007/s12192-020-01124-x Ahmad G Mohapatra BC Schulte NA Nadeau SA Luan H Zutshi N . Cbl-Family Ubiquitin Ligases and Their Recruitment of CIN85 Are Largely Dispensable for Epidermal Growth Factor Receptor Endocytosis. Int J Biochem Cell Biol (2014) 57:123–34. doi: 10.1016/j.biocel.2014.10.019 Hua XY Bie XX Cheng X Zhang SG . High Expression of CIN85 Promotes Proliferation and Invasion of Human Esophageal Squamous Cell Carcinoma. Mol Med Rep (2021) 23(1). doi: 10.3892/mmr.2020.11650 Yoshida N Tsuzuki S Karube K Takahara T Suguro M Miyoshi H . STX11 Functions as a Novel Tumor Suppressor Gene in Peripheral T-Cell Lymphomas. Cancer Sci (2015) 106(10):1455–62. doi: 10.1111/cas.12742 Wang T Han S Wu Z Han Z Yan W Liu T . XCR1 Promotes Cell Growth and Migration and Is Correlated With Bone Metastasis in non-Small Cell Lung Cancer. Biochem Biophys Res Commun (2015) 464(2):635–41. doi: 10.1016/j.bbrc.2015.06.175