Dataset Information

Screening and verifying key genes with poor prognosis in colon cancer through bioinformatics analysis

ABSTRACT:

Background

Colon cancer (CC) is one of the tumors with high morbidity and mortality in the world, and has a trend of younger generation. The molecular level of CC has not been fully elaborated. The purpose of this study is to screen and identify important genes with poor prognosis and their mechanisms at different levels.

Methods

GSE74602 and GSE10972 gene expression profiles were downloaded from the Gene Expression Omnibus (GEO) database. There were 58 normal tissues and 58 CC tissues. Differentially expressed genes (DEGs) were screened out by using the GEO2R tool and Venn diagram. Then, the DAVID online database was used to perform the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Six hub genes with the highest correlation were screened out after the modular analysis of the protein-protein interaction (PPI) network by using Cytoscape’s MCODE plug-in. Finally, the overall survival of key hub genes and potential pathways were verified in GEPIA and UALCAN database.

Results

A total of 78 up-regulated DEGs were enriched in the mitotic nuclear division, cell division, cell proliferation, anaphase-promoting complex-dependent catabolic process and G2/M transition of the mitotic cell cycle. In total, 130 down-regulated DEGs were enriched in muscle contraction, bicarbonate transport, cellular response to zinc ion, negative regulation of growth, negative regulation of leukocyte apoptotic process and one-carbon metabolic process. CDK1, CCNB1, CDC20, AURKA, CCNA2 and TOP2A were the top six hub genes, mainly enriched in cell cycle pathways. Among them, CCNB1, CDK1, CDC20, CCNA2 were enriched in the G2/M phase. GEPIA and UALCAN database confirmed that CCNA2 and CCNB1 had a significant relationship with the poor prognosis of CC patients. Meanwhile, there was a positive correlation between the two.

Conclusions

Screening out genes with abnormal expression in CC help understand the initiation and progression of CC at the molecular level and explore candidate biomarkers for diagnosis, treatment and prognosis.

SUBMITTER: Dong B

PROVIDER: S-EPMC8797306 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Colorectal cancer (CRC) is a life-threatening disease with a poor prognosis. Therefore, it is crucial to identify molecular prognostic biomarkers for CRC. The present study aimed to identify potential key genes that could be used to predict the prognosis of patients with CRC. Three CRC microarray datasets (GSE20916, GSE73360 and GSE44861) were downloaded from the Gene Expression Omnibus (GEO) database, and one dataset was obtained from The Cancer Genome Atlas (TCGA) database. The three GEO datasets were analyzed to detect differentially expressed genes (DEGs) using the BRB-ArrayTools software. Functional and pathway enrichment analyses of these DEGs were performed using the Database for Annotation, Visualization and Integrated Discovery tool. A protein-protein interaction (PPI) network of DEGs was constructed, hub genes were extracted, and modules of the PPI network were analyzed. To investigate the prognostic values of the hub genes in CRC, data from the CRC datasets of TCGA were used to perform the survival analyses based on the sample splitting method and Cox regression model. Correlation among the hub genes was evaluated using Spearman's correlation analysis. In the three GEO datasets, a total of 105 common DEGs were identified, including 51 down- and 54 up-regulated genes in CRC compared with normal colorectal tissues. A PPI network consisting of 100 DEGs and 551 edges was constructed, and 44 nodes were identified as hub genes. Among these 44 genes, the four hub genes TIMP metallopeptidase inhibitor 1 (TIMP1), solute carrier family 4 member 4 (SLC4A4), aldo-keto reductase family 1 member B10 (AKR1B10) and ATP binding cassette subfamily E member 1 (ABCE1) were associated with overall survival (OS) in patients with CRC. Three significant modules were extracted from the PPI network. The hub gene TIMP1 was present in Module 1, ABCE1 was involved in Module 2 and SLC4A4 was identified in Module 3. Univariate analysis revealed that TIMP1, SLC4A4, AKR1B10 and ABCE1 were associated with the OS of patients with CRC. Multivariate analysis demonstrated that SLC4A4 may be an independent prognostic factor associated with OS. Furthermore, the results from correlation analysis revealed that there was no correlation between TIMP1, SLC4A4 and ABCE1, whereas AKR1B10 was positively correlated with SLC4A4. In conclusion, the four key genes TIMP1, SLC4A4, AKR1B10 and ABCE1 associated with the OS of patients with CRC were identified by integrated bioinformatics analysis. These key genes may be used as prognostic biomarkers to predict the survival of patients with CRC, and may therefore represent novel therapeutic targets for CRC.

Project description:The aim of the present study was to identify potential key genes associated with the progression and prognosis of colorectal cancer (CRC). Differentially expressed genes (DEGs) between CRC and normal samples were screened by integrated analysis of gene expression profile datasets, including the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was conducted to identify the biological role of DEGs. In addition, a protein?protein interaction network and survival analysis were used to identify the key genes. The profiles of GSE9348, GSE22598 and GSE113513 were downloaded from the GEO database. A total of 405 common DEGs were identified, including 236 down? and 169 upregulated. GO analysis revealed that the downregulated DEGs were mainly enriched in 'detoxification of copper ion' [biological process, (BP)], 'oxidoreductase activity, acting on CH?OH group of donors, NAD or NADP as acceptor' [molecular function, (MF)] and 'brush border' [cellular component, (CC)]. Upregulated DEGs were mainly involved in 'nuclear division' (BP), 'snoRNA binding' (MF) and 'nucleolar part' (CC). KEGG pathway analysis revealed that DEGs were mainly involved in 'mineral absorption', 'nitrogen metabolism', 'cell cycle' and 'caffeine metabolism'. A PPI network was constructed with 268 nodes and 1,027 edges. The top one module was selected, and it was revealed that module?related genes were mainly enriched in the GO terms 'sister chromatid segregation' (BP), 'chemokine activity' (MF), and 'condensed chromosome (CC)'. The KEGG pathway was mainly enriched in 'cell cycle', 'progesterone?mediated oocyte maturation', 'chemokine signaling pathway', 'IL?17 signaling pathway', 'legionellosis', and 'rheumatoid arthritis'. DNA topoisomerase II?? (TOP2A), mitotic arrest deficient 2 like 1 (MAD2L1), cyclin B1 (CCNB1), checkpoint kinase 1 (CHEK1), cell division cycle 6 (CDC6) and ubiquitin conjugating enzyme E2 C (UBE2C) were indicated as hub genes. Furthermore, survival analysis revealed that TOP2A, MAD2L1, CDC6 and CHEK1 may serve as prognostic biomarkers in CRC. The present study provided insights into the molecular mechanism of CRC that may be useful in further investigations.

Project description:Hepatocellular carcinoma (HCC) is a malignant tumor with unsatisfactory prognosis. The abnormal genes expression is significantly associated with initiation and poor prognosis of HCC. The aim of the present study was to identify molecular biomarkers related to the initiation and development of HCC via bioinformatics analysis, so as to provide a certain molecular mechanism for individualized treatment of hepatocellular carcinoma.Three datasets (GSE101685, GSE112790, and GSE121248) from the GEO database were used for the bioinformatics analysis. Differentially expressed genes (DEGs) of HCC and normal liver samples were obtained using GEO2R online tools. Gene ontology term and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway analysis were conducted via the Database for Annotation, Visualization, and Integrated Discovery online bioinformatics tool. The protein-protein interaction (PPI) network was constructed by the Search Tool for the Retrieval of Interacting Genes database and hub genes were visualized by Cytoscape. Survival analysis and RNA sequencing expression were conducted by UALCAN and Gene Expression Profiling Interactive Analysis.A total of 115 shared DEGs were identified, including 30 upregulated genes and 85 downregulated genes in HCC samples. P53 signaling pathway and cell cycle were the major enriched pathways for the upregulated DEGs whereas metabolism-related pathways were the major enriched pathways for the downregulated DEGs. The PPI network was established with 105 nodes and 249 edges and 3 significant modules were identified via molecular complex detection. Additionally, 17 candidate genes from these 3 modules were significantly correlated with HCC patient survival and 15 of 17 genes exhibited high expression level in HCC samples. Moreover, 4 hub genes (CCNB1, CDK1, RRM2, BUB1B) were identified for further reanalysis of KEGG pathway, and enriched in 2 pathways, the P53 signaling pathway and cell cycle pathway.Overexpression of CCNB1, CDK1, RRM2, and BUB1B in HCC samples was correlated with poor survival in HCC patients, which could be potential therapeutic targets for HCC.

Project description:Colorectal cancer (CRC) is known to display a high risk of metastasis and recurrence. The main objective of our investigation was to shed more light on CRC pathogenesis by screening CRC datasets for the identification of key genes and signaling pathways, possibly leading to new approaches for the diagnosis and treatment of CRC. We downloaded the colorectal cancer datasets from the Gene Expression Omnibus (GEO) database site. We used GEO2R to screen for differentially expressed genes (DEGs) of which those with a fold change >1 were considered as up-regulated and those with a fold change <-1 were considered as down-regulated on the basis of a P < .05. "Gene ontology (GO)" and "Kyoto Encyclopedia of Genes and Genomes (KEGG)" data were analyzed by the "DAVID" software. The online search tool "STRING" was used to search for interacting genes or proteins and we used Cytoscape (v3.8.0) to generate a PPI network map and to identify key genes. Finally, survival analysis and stage mapping of key genes were performed using "GEPIA" with the aim of elucidating their potential impact on CRC. Our study revealed 120 intersecting genes of which 55 were up- and 65 were downregulated, respectively. GO analysis revealed that these genes were involved in cell proliferation, exosome secretion, G2/M transition, cytosol, protein binding, and protein kinase activity. KEGG pathway analysis showed that these genes were involved in cell cycle and mineral absorption. The Cytoscape PPI map showed 17 nodes and 262 edges, and 10 hub genes were identified by top 10 degrees. Survival analysis demonstrated that the AURKA, CCNB1, and CCNA2 genes were strongly associated with the survival rate of CRC patients. In addition, CCNB1, CCNA2, CDK1, CKS2, MAD2L1, and DLGAP5 could be correlated to pathological CRC staging. In this research, we identified key genes that may explain the molecular mechanism of occurrence and progression of CRC but may also contribute to an improvement in the clinical staging and prognosis of CRC patients.