Integrating network analysis with differential expression to uncover therapeutic and prognostic biomarkers in esophageal squamous cell carcinoma
Authors: Khurshid Sana , Usmani Shahabuddin , Ali Raiyan , Hamid Saira , Masoodi Tariq , Sadida Hana Q. , Ahmed Ikhlak , Khan Mohd Shahnawaz , Abeer Inara , Albalawi Ibrahim Altedlawi , Bedaiwi Ruqaiah I. , Mir Rashid , Al-Shabeeb Akil Ammira S. , Bhat Ajaz A. , Macha Muzafar A. 
Published: August 2024

Introduction: Esophageal squamous cell carcinoma (ESCC) accounts for over 90% of all esophageal tumors. However, the molecular mechanism underlying ESCC development and prognosis remains unclear, and there are still no effective molecular biomarkers for diagnosing or predicting the clinical outcome of patients with ESCC. Here, we used bioinformatics analysis to identify potential biomarkers and therapeutic targets for ESCC.

Methodology: Differentially expressed genes (DEGs) between ESCC and normal esophageal tissue samples were obtained by comprehensively analyzing publicly available RNA-seq datasets from the TCGA and GTEX. Gene Ontology (GO) annotation and Reactome pathway analysis identified the biological roles of the DEGs. Moreover, the Cytoscape 3.10.1 platform and subsidiary tools such as CytoHubba were used to visualize the DEGs’ protein-protein interaction (PPI) network and identify hub genes, Furthermore our results are validated by using Single-cell RNA analysis. Results: Identification of 2524 genes exhibiting altered expression enriched in pathways including keratinization, epidermal cell differentiation, G alpha(s) signaling events, and biological process of cell proliferation and division, extracellular matrix (ECM) disassembly, and muscle function. Moreover, upregulation of hallmarks E2F targets, G2M checkpoints, and TNF signaling. CytoHubba revealed 20 hub genes that had a valuable influence on the progression of ESCC in these patients. Among these, the high expression levels of four genes, CDK1 MAD2L1, PLK1, and TOP2A, were associated with critical dependence for cell survival in ESCC cell lines, as indicated by CRISPR dependency scores, gene expression data, and cell line metadata. We also identify the molecules targeting these essential hub genes, among which GSK461364 is a promising inhibitor of PLK1, BMS265246, and Valrubicin inhibitors of CDK1 and TOP2A, respectively. Moreover, we identified that elevated expression of MMP9 is associated with worse overall survival in ESCC patients, which may serve as potential prognostic biomarker or therapeutic target for ESCC. The single-cell RNA analysis showed MMP9 is highly expressed in myeloid, fibroblast, and epithelial cells, but low in T cells, endothelial cells, and B cells. This suggests MMP9’s role in tumor progression and matrix remodeling, highlighting its potential as a prognostic marker and therapeutic target.

Discussion: Our study identified key hub genes in ESCC, assessing their potential as therapeutic targets and biomarkers through detailed expression and dependency analyses. Notably, MMP9 emerged as a significant prognostic marker with high expression correlating with poor survival, underscoring its potential for targeted therapy. These findings enhance our understanding of ESCC pathogenesis and highlight promising avenues for treatment.