# Characterization of novel anoikis-related genes as prognostic biomarkers and key determinants of the immune microenvironment in esophageal cancer

**Authors:** Yani Su, Ming Zhang, Peng Xu, Pengfei Wen, Ke Xu, Jiale Xie, Xianjie Wan, Lin Liu, Zhi Yang, Mingyi Yang

PMC · DOI: 10.3389/fimmu.2025.1599171 · Frontiers in Immunology · 2025-07-11

## TL;DR

This study identifies new genes linked to esophageal cancer survival and immune response, offering potential for better treatment strategies.

## Contribution

A novel prognostic model using anoikis-related genes (ARGs) is developed to predict survival and analyze the immune microenvironment in esophageal cancer.

## Key findings

- A risk model with six ARGs (CDK1, IL17A, FOXC2, OLFM3, PIP5K1C, MAPK1) accurately predicts survival in EC patients.
- High-risk EC patients show reduced immune cell infiltration and suppressed immune functions compared to low-risk patients.
- Six potential drugs (BIRB.0796, Camptothecin, etc.) and eight immune checkpoint genes are identified as therapeutic targets.

## Abstract

Esophageal cancer (EC) ranks among the most prevalent malignancies globally and represents a significant and growing public health burden. This study aimed to construct a prognostic model leveraging anoikis-related genes (ARGs) to predict patient survival and elucidate the immunological microenvironment in EC. The findings are anticipated to enhance prognostic accuracy and inform therapeutic strategies, ultimately improving patient outcomes and treatment efficacy.

A comprehensive analysis was conducted using 11 control samples and 159 EC samples obtained from The Cancer Genome Atlas (TCGA) database, alongside associated clinical features. A total of 794 ARGs were curated from GeneCards database. Functional enrichment analyses of EC-related differentially expressed ARGs were performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Prognostic differential ARGs associated with EC were identified through univariate Cox regression analysis, while LASSO regression was employed to minimize overfitting and construct a robust risk prognostic model. The EC cohort was stratified into training and testing groups for model development and verification. Model performance was evaluated through risk curves, survival curves, time-dependent receiver operating characteristic (ROC) curves, ROC curves for the riskscore and clinical features, and independent prognostic analysis. A nomogram with high predictive accuracy was also developed to estimate the prognosis of EC patients. To assess the impact of the risk prognosis model on the immune microenvironment of EC, analyses included tumor microenvironment analysis, single-sample gene set enrichment analysis (ssGSEA), immune cell infiltration correlation analysis, and differential analysis of immune checkpoint expression. Drug sensitivity profiling was conducted to identify potential therapeutic agents for EC. Finally, the expression of selected ARGs was validated at the mRNA level in EC cell lines using real-time quantitative PCR (RT-qPCR).

The ARG-based risk prognostic model was constructed incorporating four high-risk ARGs (CDK1, IL17A, FOXC2, and OLFM3) and two low-risk ARGs (PIP5K1C and MAPK1). This model demonstrated strong predictive accuracy for the survival outcomes of EC patients. Immune correlation analyses revealed that the high-risk group exhibited significantly lower immunological scores compared to the low-risk group. Notably, immune cells such as macrophages and mast cells were markedly downregulated in the high-risk group. Additionally, key immunological functions, including APC co-inhibition, parainflammation, Type I IFN Response, and Type II IFN Response, were significantly suppressed in the high-risk group. Eight immune checkpoint-related genes (TNFRSF25, TNFRSF14, CD70, TNFSF15, TMIGD2, CD160, TNFSF18, and HHLA2) displayed distinct expression differences between high- and low-risk groups. The nomogram developed from this model demonstrated high efficacy in predicting EC patient prognosis. Furthermore, six potential therapeutic agents for EC were identified: BIRB.0796, Camptothecin, CHIR.99021, Methotrexate, PF.4708671, and Vorinostat. Finally, the mRNA expression levels of ARGs were validated using RT-qPCR in EC cell lines. Compared to normal esophageal epithelial cells (NE-2), CDK1 and MAPK1 were significantly upregulated in two EC cell lines (KYSE-30 and KYSE-180).

This study provides valuable insights into the prognostic outcomes and immune microenvironment of EC through the analysis of ARGs. Furthermore, several potential therapeutic agents for EC were identified, offering promising avenues for treatment. These findings hold significant potential for enhancing the survival outcomes of EC patients and provide meaningful guidance for clinical decision-making in managing this malignancy.

## Linked entities

- **Genes:** CDK1 (cyclin dependent kinase 1) [NCBI Gene 983], IL17A (interleukin 17A) [NCBI Gene 3605], FOXC2 (forkhead box C2) [NCBI Gene 2303], OLFM3 (olfactomedin 3) [NCBI Gene 118427], PIP5K1C (phosphatidylinositol-4-phosphate 5-kinase type 1 gamma) [NCBI Gene 23396], MAPK1 (mitogen-activated protein kinase 1) [NCBI Gene 5594], TNFRSF25 (TNF receptor superfamily member 25) [NCBI Gene 8718], TNFRSF14 (TNF receptor superfamily member 14) [NCBI Gene 8764], CD70 (CD70 molecule) [NCBI Gene 970], TNFSF15 (TNF superfamily member 15) [NCBI Gene 9966], TMIGD2 (transmembrane and immunoglobulin domain containing 2) [NCBI Gene 126259], CD160 (CD160 molecule) [NCBI Gene 11126], TNFSF18 (TNF superfamily member 18) [NCBI Gene 8995], HHLA2 (HHLA2 member of B7 family) [NCBI Gene 11148]
- **Chemicals:** BIRB.0796 (PubChem CID 156422), Camptothecin (PubChem CID 2538), CHIR.99021 (PubChem CID 9956119), Methotrexate (PubChem CID 4112), PF.4708671 (PubChem CID 51371303), Vorinostat (PubChem CID 5311)
- **Diseases:** esophageal cancer (MONDO:0007576)

## Full-text entities

- **Genes:** EGFR (epidermal growth factor receptor) [NCBI Gene 1956] {aka ERBB, ERBB1, ERRP, HER1, NISBD2, NNCIS}, NFKB1 (nuclear factor kappa B subunit 1) [NCBI Gene 4790] {aka CVID12, EBP-1, KBF1, NF-kB, NF-kB1, NF-kappa-B1}, CD163 (CD163 molecule) [NCBI Gene 9332] {aka M130, MM130, SCARI1}, CTNNB1 (catenin beta 1) [NCBI Gene 1499] {aka CTNNB, EVR7, MRD19, NEDSDV, armadillo}, CAV1 (caveolin 1) [NCBI Gene 857] {aka BSCL3, CGL3, LCCNS, MSTP085, PPH3, VIP21}, CD68 (CD68 molecule) [NCBI Gene 968] {aka GP110, LAMP4, SCARD1}, ITGB6 (integrin subunit beta 6) [NCBI Gene 3694] {aka AI1H}, HHLA2 (HHLA2 member of B7 family) [NCBI Gene 11148] {aka B7-H5, B7-H7, B7H7, B7y}, TNF (tumor necrosis factor) [NCBI Gene 7124] {aka DIF, IMD127, TNF-alpha, TNFA, TNFSF2, TNLG1F}, SLC14A1 (solute carrier family 14 member 1 (Kidd blood group)) [NCBI Gene 6563] {aka HUT11, HUT11A, HsT1341, JK, Jk(a), Jk(b)}, GSK3B (glycogen synthase kinase 3 beta) [NCBI Gene 2932], P2RX7 (purinergic receptor P2X 7) [NCBI Gene 5027] {aka P2X7}, FOXC2 (forkhead box C2) [NCBI Gene 2303] {aka FKHL14, LD, MFH-1, MFH1}, LINC00511 (long intergenic non-protein coding RNA 511) [NCBI Gene 400619] {aka LCAL5, onco-lncRNA-12}, HCRP-1 [NCBI Gene 387535], SIGLEC1 (sialic acid binding Ig like lectin 1) [NCBI Gene 6614] {aka CD169, SIGLEC-1, SN}, KRAS (KRAS proto-oncogene, GTPase) [NCBI Gene 3845] {aka 'C-K-RAS, C-K-RAS, CFC2, K-RAS2A, K-RAS2B, K-RAS4A}, CAMKK2 (calcium/calmodulin dependent protein kinase kinase 2) [NCBI Gene 10645] {aka CAMKK, CAMKKB}, STC1 (stanniocalcin 1) [NCBI Gene 6781] {aka STC}, OLFM3 (olfactomedin 3) [NCBI Gene 118427] {aka NOE3, NOELIN3, OPTIMEDIN}, PIP5K1C (phosphatidylinositol-4-phosphate 5-kinase type 1 gamma) [NCBI Gene 23396] {aka LCCS3, PIP5K-GAMMA, PIP5K1-gamma, PIP5Kgamma}, MAPK13 (mitogen-activated protein kinase 13) [NCBI Gene 5603] {aka MAPK 13, MAPK-13, PRKM13, SAPK4, p38delta}, IL1RAP (interleukin 1 receptor accessory protein) [NCBI Gene 3556] {aka C3orf13, IL-1RAcP, IL1R3}, CCNB1 (cyclin B1) [NCBI Gene 891] {aka CCNB}, METTL14 (methyltransferase 14, N6-adenosine-methyltransferase non-catalytic subunit) [NCBI Gene 57721] {aka hMETTL14}, SERPINA2 (serpin family A member 2 (gene/pseudogene)) [NCBI Gene 390502] {aka ARGS, ATR, PIL, SERPINA2P, psiATR}, TNFRSF25 (TNF receptor superfamily member 25) [NCBI Gene 8718] {aka APO-3, DDR3, DR3, LARD, TNFRSF12, TR3}, RPS6KB1 (ribosomal protein S6 kinase B1) [NCBI Gene 6198] {aka PS6K, S6K, S6K-beta-1, S6K1, STK14A, p70 S6KA}, YAP1 (Yes1 associated transcriptional regulator) [NCBI Gene 10413] {aka COB1, YAP, YAP-1, YAP2, YAP65, YKI}, ATG7 (autophagy related 7) [NCBI Gene 10533] {aka APG7-LIKE, APG7L, GSA7, SCAR31}, NCAPG (non-SMC condensin I complex subunit G) [NCBI Gene 64151] {aka CAPG, CHCG, NY-MEL-3, YCG1}, PLAG1 (PLAG1 zinc finger) [NCBI Gene 5324] {aka PSA, SGPA, SRS4, ZNF912}, CD8A (CD8 subunit alpha) [NCBI Gene 925] {aka CD8, CD8alpha, IMD116, Leu2, p32}, MVB12B (multivesicular body subunit 12B) [NCBI Gene 89853] {aka C9orf28, FAM125B}, HNF4A (hepatocyte nuclear factor 4 alpha) [NCBI Gene 3172] {aka FRTS4, HNF4, HNF4a7, HNF4a8, HNF4a9, HNF4alpha}, AKT1 (AKT serine/threonine kinase 1) [NCBI Gene 207] {aka AKT, PKB, PKB-ALPHA, PRKBA, RAC, RAC-ALPHA}, PRKAB1 (protein kinase AMP-activated non-catalytic subunit beta 1) [NCBI Gene 5564] {aka AMPK, HAMPKb}, TNFSF18 (TNF superfamily member 18) [NCBI Gene 8995] {aka AITRL, GITRL, TL6, TNLG2A, hGITRL}, CDK1 (cyclin dependent kinase 1) [NCBI Gene 983] {aka CDC2, CDC28A, P34CDC2}, CD160 (CD160 molecule) [NCBI Gene 11126] {aka BY55, NK1, NK28}, POTEF (POTE ankyrin domain family member F) [NCBI Gene 728378] {aka A26C1B, POTE2alpha, POTEACTIN}, STK11 (serine/threonine kinase 11) [NCBI Gene 6794] {aka LKB1, PJS, hLKB1}, GLUD1 (glutamate dehydrogenase 1) [NCBI Gene 2746] {aka GDH, GDH1, GLUD, hGDH1}, DHFR (dihydrofolate reductase) [NCBI Gene 1719] {aka DHFR1, DYR}, MIR574 (microRNA 574) [NCBI Gene 693159] {aka MIR574-3p, MIRN574, hsa-mir-574, mir-574}, APC (APC regulator of Wnt signaling pathway) [NCBI Gene 324] {aka BTPS2, DESMD, DP2, DP2.5, DP3, GS}, CD4 (CD4 molecule) [NCBI Gene 920] {aka CD4mut, IMD79, Leu-3, OKT4D, T4}, CDH1 (cadherin 1) [NCBI Gene 999] {aka Arc-1, BCDS1, CD324, CDHE, ECAD, LCAM}, IL17A (interleukin 17A) [NCBI Gene 3605] {aka CTLA-8, CTLA8, IL-17, IL-17A, IL17, ILA17}, SOX2 (SRY-box transcription factor 2) [NCBI Gene 6657] {aka ANOP3, MCOPS3}, PIK3CB (phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit beta) [NCBI Gene 5291] {aka P110BETA, PI3K, PI3KBETA, PIK3C1}, MAPK1 (mitogen-activated protein kinase 1) [NCBI Gene 5594] {aka ERK, ERK-2, ERK2, ERT1, MAPK2, NS13}, NFIX (nuclear factor I X) [NCBI Gene 4784] {aka CTF, MALNS, MRSHSS, NF-I/X, NF1-X, NF1A}, BCL2L11 (BCL2 like 11) [NCBI Gene 10018] {aka BAM, BIM, BOD}, FAM3C (FAM3 metabolism regulating signaling molecule C) [NCBI Gene 10447] {aka GS3786, ILEI}, CD70 (CD70 molecule) [NCBI Gene 970] {aka CD27-L, CD27L, CD27LG, LPFS3, TNFSF7, TNLG8A}, TNFSF15 (TNF superfamily member 15) [NCBI Gene 9966] {aka TL1, TL1A, TNLG1B, VEGI, VEGI192A}, IFNG (interferon gamma) [NCBI Gene 3458] {aka IFG, IFI, IMD69}, TNFRSF14 (TNF receptor superfamily member 14) [NCBI Gene 8764] {aka ATAR, CD270, HVEA, HVEM, LIGHTR, TR2}, MTOR (mechanistic target of rapamycin kinase) [NCBI Gene 2475] {aka FRAP, FRAP1, FRAP2, RAFT1, RAPT1, SKS}
- **Diseases:** hepatocellular carcinoma (MESH:D006528), pancreatic cancer (MESH:D010190), melanoma (MESH:D008545), Ewing sarcoma (MESH:D012512), non-small cell lung cancer (MESH:D002289), lymph node metastasis (MESH:D008207), distant metastasis (MESH:D009362), anxiety (MESH:D001007), gynecological malignancies (MESH:D005833), acute lymphoblastic leukemia (MESH:D054198), head and neck epithelial carcinoma (MESH:D006258), head and neck squamous cell carcinoma (MESH:D000077195), rectal cancer (MESH:D012004), ESCC (MESH:D000077277), prostate cancer (MESH:D011471), cutaneous melanoma (MESH:C562393), Cancer (MESH:D009369), nasopharyngeal cancer (MESH:D009303), inflammation (MESH:D007249), bladder cancer (MESH:D001749), hematologic (MESH:D006402), reperfusion injury (MESH:D015427), clear cell renal cell carcinoma (MESH:D002292), osteosarcoma (MESH:D012516), cervical cancer (MESH:D002583), ovarian and cervical cancers (MESH:D010051), cerebral ischemia (MESH:D002545), leukemia (MESH:D007938), lymph node (MESH:D000072717), carcinogenesis (MESH:D063646), lung cancer (MESH:D008175), MF (MESH:C567116), breast cancer (MESH:D001943), EC (MESH:D004938), primary cutaneous T-cell lymphoma (MESH:D016410), triple-negative breast cancer (MESH:D064726), infarct (MESH:D007238), colon cancer (MESH:D015179), gastrointestinal stromal tumors (MESH:D046152), gastric adenocarcinoma (MESH:D013274)
- **Chemicals:** lipid (MESH:D008055), paclitaxel (MESH:D017239), cisplatin (MESH:D002945), Vorinostat (MESH:D000077337), staurosporine (MESH:D019311), duloxetine (MESH:D000068736), CO2 (MESH:D002245), Methotrexate (MESH:D008727), disulfiram (MESH:D004221), CHIR-99021 (MESH:C473711), TRIzol (MESH:C411644), capivasertib (MESH:C575618), ATP (MESH:D000255), PF-4708671 (MESH:C552719), Tubeimoside V (MESH:C514048), BH4 (MESH:C003402), Camptothecin (MESH:D002166), dihydrobiopterin (MESH:C017226), BH2 (-)
- **Species:** Homo sapiens (human, species) [taxon 9606]
- **Mutations:** (AUC) at 1, -197 A/G, G12D, A/G
- **Cell lines:** MDA-MB-231 — Homo sapiens (Human), Breast adenocarcinoma, Cancer cell line (CVCL_0062), KYSE-30 — Homo sapiens (Human), Esophageal squamous cell carcinoma, Cancer cell line (CVCL_1351), KYSE-180 — Homo sapiens (Human), Esophageal squamous cell carcinoma, Cancer cell line (CVCL_1349), NE-2 — Mus musculus (Mouse), Mouse teratocarcinoma, Cancer cell line (CVCL_3554)

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC12311635/full.md

## Figures

12 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12311635/full.md

## References

72 references — full list in the complete paper: https://tomesphere.com/paper/PMC12311635/full.md

---
Source: https://tomesphere.com/paper/PMC12311635