Dataset Information

Identification of latent biomarkers in connection with progression and prognosis in oral cancer by comprehensive bioinformatics analysis.

ABSTRACT:

Background

Oral cancer (OC) is a common and dangerous malignant tumor with a low survival rate. However, the micro level mechanism has not been explained in detail.

Methods

Gene and miRNA expression micro array data were extracted from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) and miRNAs (DE miRNAs) were identified by R software. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of genes and genomes (KEGG) pathway analysis were used to assess the potential molecular mechanisms of DEGs. Cytoscape software was utilized to construct protein-protein interaction (PPI) network and miRNA-gene network. Central genes were screened out with the participation of gene degree, molecular complex detection (MCODE) plugin, and miRNA-gene network. Then, the identified genes were checked by The Cancer Genome Atlas (TCGA) gene expression profile, Kaplan-Meier data, Oncomine, and the Human Protein Atlas database. Receiver operating characteristic (ROC) curve was drawn to predict the diagnostic efficiency of crucial gene level in normal and tumor tissues. Univariate and multivariate Cox regression were used to analyze the effect of dominant genes and clinical characteristics on the overall survival rate of OC patients.

Results

Gene expression data of gene expression profiling chip(GSE9844, GSE30784, and GSE74530) were obtained from GEO database, including 199 tumor and 63 non-tumor samples. We identified 298 gene mutations, including 200 upregulated and 98 downregulated genes. GO functional annotation analysis showed that DEGs were enriched in extracellular structure and extracellular matrix containing collagen. In addition, KEGG pathway enrichment analysis demonstrated that the DEGs were significantly enriched in IL-17 signaling pathway and PI3K-Akt signaling pathway. Then, we detected three most relevant modules in PPI network. Central genes (CXCL8, DDX60, EIF2AK2, GBP1, IFI44, IFI44L, IFIT1, IL6, MMP9,CXCL1, CCL20, RSAD2, and RTP4) were screened out with the participation of MCODE plugin, gene degree, and miRNA-gene network. TCGA gene expression profile and Kaplan-Meier analysis showed that high expression of CXCL8, DDX60, IL6, and RTP4 was associated with poor prognosis in OC patients, while patients with high expression of IFI44L and RSAD2 had a better prognosis. The elevated expression of CXCL8, DDX60, IFI44L, RSAD2, and RTP44 in OC was verified by using Oncomine database. ROC curve showed that the mRNA levels of these five genes had a helpful diagnostic effect on tumor tissue. The Human Protein Atlas database showed that the protein expressions of DDX60, IFI44L, RSAD2, and RTP44 in tumor tissues were higher than those in normal tissues. Finally, univariate and multivariate Cox regression showed that DDX60, IFI44L, RSAD2, and RTP44 were independent prognostic indicators of OC.

Conclusion

This study revealed the potential biomarkers and relevant pathways of OC from publicly available GEO database, and provided a theoretical basis for elucidating the diagnosis, treatment, and prognosis of OC.

SUBMITTER: Reyimu A

PROVIDER: S-EPMC8361649 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Epithelial ovarian cancer (EOC) is one of the malignancies in women, which has the highest mortality. However, the microlevel mechanism has not been discussed in detail. The expression profiles GSE27651, GSE38666, GSE40595, and GSE66957 including 188 tumor and 52 nontumor samples were downloaded from the Gene Expression Omnibus database. The differentially expressed genes (DEGs) were filtered using R software, and we performed functional analysis using the clusterProfiler. Cytoscape software, the molecular complex detection plugin and database STRING analyzed DEGs to construct protein-protein interaction network. We identified 116 DEGs including 81 upregulated and 35 downregulated DEGs. Functional analysis revealed that they were significantly enriched in the extracellular region and biosynthesis of amino acids. We next identified four bioactive compounds (vorinostat, LY-294002,trichostatin A, and tanespimycin) based on ConnectivityMap. Then 114 nodes were obtained in protein-protein interaction. The three most relevant modules were detected. In addition, according to degree ≥ 10, 14 core genes including FOXM1, CXCR4, KPNA2, NANOG, UBE2C, KIF11, ZWINT, CDCA5, DLGAP5, KIF15, MCM2, MELK, SPP1, and TRIP13 were identified. Kaplan-Meier analysis, Oncomine, and Gene Expression Profiling Interactive Analysis showed that overexpression of FOXM1, SPP1, UBE2C, KIF11, ZWINT, CDCA5, UBE2C, and KIF15 was related to bad prognosis of EOC patients. CDCA5, FOXM1, KIF15, MCM2, and ZWINT were associated with stage. Receiver operating characteristic (ROC) curve showed that messenger RNA levels of these five genes exhibited better diagnostic efficiency for normal and tumor tissues. The Human Protein Atlas database was performed. The protein levels of these five genes were significantly higher in tumor tissues compared with normal tissues. Functional enrichment analysis suggested that all the hub genes played crucial roles in citrate cycle tricarboxylic acid cycle. Furthermore, the univariate and multivariate Cox proportional hazards regression showed that ZWINT was independent prognostic indictor among EOC patients. The genes and pathways discovered in the above studies may open a new direction for EOC treatment.

Project description:BackgroundEndometrial cancer (EC) is one of the female malignant tumors. Endometrial cancer predominately affects post-menopausal women. Bioinformatics analysis has been widely applied to screen and analyze genes in linkage to various types of cancer progression.MethodsDownload the gene expression profile from Gene Expression Omnibus (GEO). Calculate raw expression data according to pre-processing procedures. We performed the "limma" R language package to screen DEGs between Endometrial cancer tissue samples and normal uterus tissue samples. Enrichment of the functions and pathways was analyzed by using clusterprofiler. We utilized Search Tool for the Retrieval of Interacting Genes Database (STRING) to assess protein-protein interaction (PPI) information, and then we used plug-in Molecular Complex Detection (MCODE) to screen hub modules of PPI network in Cytoscape. We also performed functional analysis on the genes in the hub module by using clusterprofiler. Next, we utilized the "WGCNA" package in R to establish co-expression network for the DEGs. The Venn diagram was performed to overlap the gene in key module and hub PPI cluster. We validated the key genes in TCGA, GEPIA, UALCAN and Immunohistochemistry staining obtained from The Human Protein Atlas database. And then we did ROC curve analysis by SPSS. Gene set enrichment analysis (GSEA) and mutation analysis were also performed for hub genes.ResultsFunctional and pathway enrichment analysis demonstrated that the upregulated differentially expressed genes (DEGs) were significantly enriched in CXCR chemokine receptor binding, chemokine activity, chemokine receptor binding, G-protein coupled receptor binding, RAGE receptor binding, cytokine activity, microtubule binding, receptor regulator activity and microtubule motor activity, and the down-regulated genes were highly enriched in collagen binding. After using STRING software to construct PPI network, 30 prominent proteins were identified and the first two significant modules were selected. In co-expression network, 5 EC-related modules were identified. Among them, the turquoise module has the highest correlation with the EC. We further analyzed the genes in the PPI and turquoise module, and selected eleven key genes related to EC after validation of TCGA database, GEPIA, UALCAN and immunohistochemistry. Six of them had mutation significance.ConclusionsIn summary, these 11 genes may become new therapy targets for EC treatment.

Project description:BackgroundThis study aimed to identify the key genes associated with prognosis in thyroid cancer (TC), and to explore potential pathways.MethodsGSE66783, GSE58545, and GSE129562 datasets were used to identify differentially expressed genes (DEGs) between tumor and normal tissues, followed by KEGG analyses on DEGs. The protein-protein interaction (PPI) network of DEGs was subsequently constructed to find the top 10 hub genes and seed genes in the whole network. Furthermore, the mRNA expressions of hub genes and prognostic values were explored. Regarding the seed gene, pathway activity score and GSEA analyses were conducted as well.Results1) A total of 183 DEGs were consistently expressed in three datasets comprising 76 up-regulated and 107 down-regulated genes. DEGs were mainly enriched in the p53 signaling pathway, complement and coagulation cascades, and hedgehog signaling pathway. 2) The top 10 hub genes, including CCND1, TIMP1, ICAM1, MET, PLAU, LDLR, PLAUR, ITGA2, ITGA3, and LGALS3, were identified. All hub genes were highly expressed in TC compared with normal samples. 3) High expression of CCND1, TIMP1, MET, and LGALS3 statistically correlated with a favorable prognosis of patients. Poor survival was observed in patients with ITGA2 and ITGA3 high expression. There was no association between ICAM1, PLAU, and PLAUR expression and survival of patients. LGALS3 and TIMP1 were further identified as independent prognostic factors in TC. 4) Among 10 hub genes, TIMP1 was determined as the seed gene, indicating its significance in the whole network. We further found that in most of the famous cancer-related pathways, TIMP1 higher expression caused a lower pathway activity, suggesting its inhibitory effect to these pathways in TC. In addition, TIMP1 positively correlated with the p53 signaling pathway, complement, and coagulation cascades involved in TC.ConclusionThe present study provided seven prognosis-associated genes in TC and revealed several significant pathways, which contributed to elucidate the pathogenesis of TC.

Project description:The progression of non-small cell lung cancer (NSCLC) is linked to epithelial-mesenchymal transition (EMT), a biologic process that enables tumor cells to acquire a migratory phenotype and resistance to chemo- and immunotherapies. Discovery of novel biomarkers in NSCLC progression is essential for improved prognosis and pharmacological interventions. In the current study, we performed an integrated bioinformatics analysis on gene expression datasets of TGF-β-induced EMT in NSCLC cells to identify novel gene biomarkers and elucidate their regulation in NSCLC progression. The gene expression datasets were extracted from the NCBI Gene Expression Omnibus repository, and differentially expressed genes (DEGs) between TGF-β-treated and untreated NSCLC cells were retrieved. A protein-protein interaction network was constructed and hub genes were identified. Functional and pathway enrichment analyses were conducted on module DEGs, and a correlation between the expression levels of module genes and survival of NSCLC patients was evaluated. Prediction of interactions of the biomarker genes with transcription factors and miRNAs was also carried out. We described four protein clusters in which DEGs were associated with ubiquitination (Module 1), regulation of cell death and cell adhesions (Module 2), oxidation-reduction reactions of aerobic respiration (Module 3) and mitochondrial translation (Module 4). From the module genes, we identified ten prognostic gene biomarkers in NSCLC. Low expression levels of KCTD6, KBTBD7, LMO7, SPSB2, RNF19A, FOXA2, DHTKD1, CDH1 and PDHB and high expression level of KLHL25 were associated with reduced overall survival of NSCLC patients. Most of these biomarker genes were involved in protein ubiquitination. The regulatory network of the gene biomarkers revealed their interaction with tumor suppressor miRNAs and transcription factors involved in the mechanisms of cancer progression. This ten-gene prognostic signature can be useful to improve risk prediction and therapeutic strategies in NSCLC. Our analysis also highlights the importance of deregulation of ubiquitination in EMT-associated NSCLC progression.

Project description:Background: Cervical carcinoma is one of the most common malignant gynecological tumors and is associated with high rates of morbidity and mortality. Early diagnosis and early treatment can reduce the mortality rate of cervical cancer. However, there is still no specific biomarkers for the diagnosis and detection of cervical cancer prognosis. Therefore, it is greatly urgent in searching biomarkers correlated with the diagnosis and prognosis of cervical cancer. Results: The mRNA and microRNA expression profile datasets (GSE7803, GSE9750, GSE63514, and GSE30656) were downloaded from the Gene Expression Omnibus database (GEO). The three microarray datasets were integrated to one via integrated bioinformatics. Differentially expressed genes (DEGs) and microRNAs (DEMs) were obtained by R software. The protein-protein interaction (PPI) networks of the DEGs were performed from the STRING database and further visualized by Cytoscape software. A total of 83 DEGs and 14 DEMs were screened from the microarray expression profile datasets. The miRNAs validated to be associated with cervical cancer were obtained using HMDD online website and the target genes of DEMs were identified using the miRWalk2.0 online database. ESR1, PPP1R3C, NSG1, and TMPRSS11D were the gene targets of hsa-miR-21; the targets of hsa-miR-16 were GYS2, ENDOU, and KLF4. These targets were all downregulated in cervical cancer. Finally, we verified the expression of those targets in cervical tissues from TCGA and GTEx databases and analyzed their relationship with survival of cervical cancer patients. In the end, the expression of key genes in cervical cancer tissues was verified via experiment method, we found KLF4 and ESR1 were downregulated in tumor tissues. Conclusion: This study indicates that KLF4 and ESR1 are downregulated by the upregulated miR21 and miRNA16 in cervical cancer, respectively, using bioinformatics analysis, and the lower expression of KLF4 and ESR1 is closely related to the poor prognosis. They might be of clinical significance for the diagnosis and prognosis of cervical cancer, and provide effective targets for the treatment of cervical cancer.