Dataset Information

Identification of Stemness-Related Genes for Cervical Squamous Cell Carcinoma and Endocervical Adenocarcinoma by Integrated Bioinformatics Analysis.

ABSTRACT:

Background

Invasion and metastasis of cervical cancer are the main factors affecting the prognosis of patients with cervical squamous cell carcinoma (CESC). Therefore, it is of vital importance to find novel biomarkers that are associated with CESC invasion and metastasis, which will aid in the amelioration of individualized therapeutic methods for advanced patients.

Methods

The gene expression profiles of 10 metastatic and 116 non-metastatic samples were downloaded from The Cancer Genome Atlas (TCGA), where differentially expressed genes (DEGs) were defined. Weighted gene correlation network analysis (WGCNA) was employed to identify the stemness-related genes (SRGs). Univariate and multivariate regression analyses were used to identify the most significant prognostic key genes. Differential expression analysis of transcription factor (TF) and Gene Set Variation Analysis (GSVA) were utilized to explore the potential upstream regulation of TFs and downstream signaling pathways, respectively. Co-expression analysis was performed among significantly enriched TFs, key SRGs, and signaling pathways to construct a metastasis-specific regulation network in CESC. Connectivity Map (CMap) analysis was performed to identify bioactive small molecules which might be potential inhibitors for the network. Additionally, direct regulatory patterns of key genes were validated by ChIP-seq and ATAC-seq data.

Results

DEGs in yellow module acquired via WGCNA were defined as key genes which were most significantly related to mRNAsi. A multivariate Cox regression model was constructed and then utilized to explore the prognostic value of key SRGs by risk score. Area under curve (AUC) of the receiver operating characteristic (ROC) curve was 0.842. There was an obvious co expression pattern between the TF NR5A2 and the key gene VIM (R = 0.843, p < 0.001), while VIM was also significantly co-expressed with hallmark epithelial mesenchymal transition (EMT) signaling pathway (R = 0.318, p < 0.001). Naringenin was selected as the potential bioactive small molecule inhibitor for metastatic CESC based on CMap analysis.

Conclusions

VIM positively regulated by NR5A2 affected EMT signaling pathways in metastatic CESC, and naringenin was the inhibitor for the treatment of metastatic CESC via suppressing cancer stemness. This hypothetical signaling axis and potential inhibitors provide biomarkers and novel therapeutic targets for metastatic CESC.

SUBMITTER: Guo H

PROVIDER: S-EPMC8027330 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundThe purpose of this study was to investigate the prognostic signature of necroptosis-related lncRNAs (NRLs) and explore their association with immune-related functions and sensitivity of the therapeutic drug in cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC).MethodsUCSC Xena provided lncRNA sequencing and clinical data about CESC, and a necroptosis gene list was obtained from the KEGG database. NRLs were selected by structuring a co-expression network of lncRNAs and necroptosis-related genes. To further screen lncRNAs, we used the univariate Cox regression method, Lasso regression, and multivariate Cox regression. Afterward, an NRL signature was established. We used the xCell algorithm and single-sample gene set enrichment analysis (ssGSEA) to clarify the pertinence between immune infiltration and NRL expressions in CESC patients and explored the relationship between the target lncRNAs and immune-related genes. By leveraging the GDSC database, the therapy-sensitive response of the prognostic signature was forecasted and an experimental validation was performed. We performed GSEA with the aim of recognizing the potential pathway related to the individual prognostic signature.ResultsThe two prognostic NRLs (AC009095.1 and AC005332.4) showed significant diversity and constituted the NRL signature. On the grounds of our signature, risk score was an independent element which was bound up with patient outcome (HR = 4.97 CI: 1.87-13.2, P = 0.001). The CESC patients were classified by the median risk score. Immune infiltration analysis revealed significant increases in CD4 + Tcm, eosinophils, epithelial cells, fibroblasts, NKT, plasma cells, platelets, and smooth muscle in the high-risk group (P< 0.05). Target lncRNAs also showed some correlation with NRGs. The estimated IC50 values of bicalutamide, CHIR.99021, and imatinib were lower in the high-risk group. Through the subsequent experimental validation, both AC009095.1 and AC005332.4 were significantly more highly expressed in SiHa than in Hela. AC009095.1 was expressed more highly in SiHa than in HUCEC, but the expression of AC005332.4 was reversed.ConclusionsThis study elucidated that NRLs, as a novel signature, were indispensable factors which can significantly influence the prognosis of patients with CESC and could provide novel clinical evidence to serve as a potential molecular biomarker for future therapeutic targets.

Project description:BackgroundCervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) are among the most prevalent gynecologic malignancies globally. The prognosis is abysmal once cervical cancer progresses to lymphatic metastasis. Anoikis, a specialized form of apoptosis induced by loss of cell adhesion to the extracellular matrix, plays a critical role. The prediction model based on anoikis-related genes (ARGs) expression and clinical data could greatly aid clinical decision-making. However, the relationship between ARGs and CESC remains unclear.MethodsARGs curated from the GeneCards and Harmonizome portals were instrumental in delineating CESC subtypes and in developing a prognostic framework for patients afflicted with this condition. We further delved into the intricacies of the immune microenvironment and pathway enrichment across the identified subtypes. Finally, our efforts culminated in the creation of an innovative nomogram that integrates ARGs. The utility of this prognostic tool was underscored by Decision Curve Analysis (DCA), which illuminate its prospective benefits in guiding clinical interventions.ResultsIn our study, We discerned a set of 17 survival-pertinent, anoikis-related differentially expressed genes (DEGs) in CESC, from which nine were meticulously selected for the construction of prognostic models. The derived prognostic risk score was subsequently validated as an autonomous prognostic determinant. Through comprehensive functional analyses, we observed distinct immune profiles and drug response patterns among divergent prognostic stratifications. Further, we integrated the risk scores with the clinicopathological characteristics of CESC to develop a robust nomogram. DCA corroborated the utility of our model, demonstrating its potential to enhance patient outcomes through tailored clinical treatment strategies.ConclusionThe predictive signature, encompassing nine pivotal genes, alongside the meticulously constructed nomogram developed in this research, furnishes clinicians with a sophisticated tool for tailoring treatment strategies to individual patients diagnosed with CESC.

Project description:BackgroundEsophageal squamous cell carcinoma (ESCC) as the main subtype of esophageal cancer (EC) is a leading cause of cancer-related death worldwide. Despite advances in early diagnosis and clinical management, the long-term survival of ESCC patients remains disappointing, due to a lack of full understanding of the molecular mechanisms.MethodsIn order to identify the differentially expressed genes (DEGs) in ESCC, the microarray datasets GSE20347 and GSE26886 from Gene Expression Omnibus (GEO) database were analyzed. The enrichment analyses of gene ontology (GO) terms, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and Gene Set Enrichment Analysis (GSEA) were performed for the DEGs. The protein-protein interaction (PPI) network of these DEGs was constructed using the Cytoscape software based on the STRING database to select as hub genes for weighted co-expression network analysis (WGCNA) with ESCC samples from TCGA database.ResultsA total of 746 DEGs were commonly shared in the two datasets including 286 upregulated genes and 460 downregulated genes in ESCC. The DEGs were enriched in biological processes such as extracellular matrix organization, proliferation and keratinocyte differentiation, and were enriched in biological pathways such as ECM-receptor interaction and cell cycle. GSEA analysis also indicated the enrichment of upregulated DEGs in cell cycle. The 40 DEGs were selected as hub genes. The MEblack module was found to be enriched in the cell cycle, Spliceosome, DNA replication and Oocyte meiosis. Among the hub genes correlated with MEblack module, GSEA analysis indicated that DEGs of TCGA samples with DLGAP5 upregulation was enriched in cell cycle. Moreover, the highly endogenous expression of DLGAP5 was confirmed in ESCC cells. DLGAP5 knockdown significantly inhibited the proliferation of ESCC cells.ConclusionsDEGs and hub genes such as DLGAP5 from independent datasets in the current study will provide clues to elucidate the molecular mechanisms involved in development and progression of ESCC.