Identifying Key Epigenetic Modification-Related Genes for Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma and Cellular Validation.

Abstract

Introduction: Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma (CESC) is a highly prevalent female malignancy. As the epigenomic characteristics of immune cells and cancer cells can serve as predictive indicators for the response to cancer immunotherapy, analysis of epigenetically modified genes (EpiGenes) could contribute to CESC treatment.

Methods: The ssGSEA algorithm was employed to compute EpiGenes scores. Core genes that exhibited significant module association and a close correlation with EpiGenes scores were identified via the WGCNA package. Univariate Cox proportional hazards regression was performed on the core genes using the survival package, followed by gene set reduction via LASSO Cox regression. Ultimately, key genes were determined through multivariate Cox regression to establish a RiskScore model. Further, the optimal risk cutoff was determined using the survminer package to stratify CESC patients into high- and low-risk subgroups. For enrichment analysis, clusterProfiler and GSEA were utilized. Immune infiltration across risk groups was evaluated via ssGSEA, the MCPcounter algorithm, and the ESTIMATE algorithm. TIDE was employed to compare immunotherapeutic responses between the risk groups, while the pRRophetic software was utilized to predict patients' chemotherapeutic drug sensitivity. The biomarkers identified were validated by performing in vitro experiments.

Results: CEP78, DOCK7, DPY19L4, and POM121 were identified by computational analyses as the key genes for CESC and further validated through in vitro experiments. Pathway enrichment analysis revealed predominant enrichment in immune-related pathways in the high-risk group, whereas the low-risk group was more enriched in energy and metabolic pathways. A significant negative correlation was observed between CD8+ T cell abundance and RiskScore, with higher ESTIMATEScores and StromalScores in high-risk patients. Notably, the high-risk group also demonstrated lower potential sensitivity to immunotherapy but more active responsiveness to a broader spectrum of chemotherapeutic agents.

Discussion: Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that module genes are significantly enriched in cell cycle regulatory pathways, and these genes, in conjunction with Human Papillomavirus (HPV) infection-induced cell cycle dysregulation, jointly participate in CESC pathogenesis, providing a mechanistic basis for understanding the disease.

Conclusion: This study provided novel theoretical evidence for immunotherapy and chemotherapy selection in the management of CESC.