Integrated bioinformatics analysis of key genes involved in progress of colon cancer.
Ontology highlight
ABSTRACT: Background
Colon cancer is one of most malignant cancers around worldwide. Nearly 20% patients were diagnosed at colon cancer with metastasis. However, the lack of understanding regarding its pathogenesis brings difficulties to study it.Methods
In this study, we acquired high-sequence data from GEO dataset, and performed integrated bioinformatic analysis including differently expressed genes, gene ontology and Kyoto Encyclopedia of Genes and Genomes pathways analysis, protein-protein analysis, survival analysis to analyze the development of colon cancer.Results
By comparing the colon cancer tissues with normal colon tissues, 109 genes were dysregulated; among them, 83 genes were downregulated and 26 genes were upregulated. Two clusters were founded based on the STRING database and MCODE plugin of cytoscape software. Then, six genes with prognostic value were filtered out in UALCAN website.Conclusion
We found that SPP1, VIP, COL11A1, CA2, ADAM12, INHBA could provide great significant prognostic value for colon cancer.
SUBMITTER: Yang H
PROVIDER: S-EPMC6465657 | biostudies-literature | 2019 Apr
REPOSITORIES: biostudies-literature
Publications
Integrated bioinformatics analysis of key genes involved in progress of colon cancer.
Yang Haojie H Wu Jiong J Zhang Jingjing J Yang Zhigang Z Jin Wei W Li Ying Y Jin Lei L Yin Lu L Liu Hua H Wang Zhenyi Z
Molecular genetics & genomic medicine 20190211 4
<h4>Background</h4>Colon cancer is one of most malignant cancers around worldwide. Nearly 20% patients were diagnosed at colon cancer with metastasis. However, the lack of understanding regarding its pathogenesis brings difficulties to study it.<h4>Methods</h4>In this study, we acquired high-sequence data from GEO dataset, and performed integrated bioinformatic analysis including differently expressed genes, gene ontology and Kyoto Encyclopedia of Genes and Genomes pathways analysis, protein-pro ...[more]
Similar Datasets
Project description:Background: Valvular heart disease is obtaining growing attention in the cardiovascular field and it is believed that calcific aortic valve disease (CAVD) is the most common valvular heart disease (VHD) in the world. CAVD does not have a fully effective treatment to delay its progression and the specific molecular mechanism of aortic valve calcification remains unclear. Materials and Methods: We obtained the gene expression datasets GSE12644 and GSE51472 from the public comprehensive free database GEO. Then, a series of bioinformatics methods, such as GO and KEGG analysis, STING online tool, Cytoscape software, were used to identify differentially expressed genes in CAVD and healthy controls, construct a PPI network, and then identify key genes. In addition, immune infiltration analysis was used via CIBERSORT to observe the expression of various immune cells in CAVD. Results: A total of 144 differential expression genes were identified in the CAVD samples in comparison with the control samples, including 49 up-regulated genes and 95 down-regulated genes. GO analysis of DEGs were most observably enriched in the immune response, signal transduction, inflammatory response, proteolysis, innate immune response, and apoptotic process. The KEGG analysis revealed that the enrichment of DEGs in CAVD were remarkably observed in the chemokine signaling pathway, cytokine-cytokine receptor interaction, and PI3K-Akt signaling pathway. Chemokines CXCL13, CCL19, CCL8, CXCL8, CXCL16, MMP9, CCL18, CXCL5, VCAM1, and PPBP were identified as the hub genes of CAVD. It was macrophages that accounted for the maximal proportion among these immune cells. The expression of macrophages M0, B cells memory, and Plasma cells were higher in the CAVD valves than in healthy valves, however, the expression of B cells naïve, NK cells activated, and macrophages M2 were lower. Conclusion: We detected that chemokines CXCL13, CXCL8, CXCL16, and CXCL5, and CCL19, CCL8, and CCL18 are the most important markers of aortic valve disease. The regulatory macrophages M0, plasma cells, B cells memory, B cells naïve, NK cells activated, and macrophages M2 are probably related to the occurrence and the advancement of aortic valve stenosis. These identified chemokines and these immune cells may interact with a subtle adjustment relationship in the development of calcification in CAVD.
| S-EPMC8144713 | biostudies-literature
Project description:Wilms tumor (WT) is one of the most common types of pediatric solid tumors; however, its molecular mechanisms remain unclear. The present study aimed to identify key genes and microRNAs (miRNAs), and to predict the underlying molecular mechanisms of WT using integrated bioinformatics analysis. Original gene expression profiles were downloaded from the Gene Expression Omnibus (GEO; accession, GSE66405) and The Cancer Genome Atlas (TCGA) databases. Similarly, miRNA expression patterns were downloaded from GEO (accession, GSE57370) and TCGA. R version 3.5.0 software was used to identify differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs) using the limma and edgeR packages. Kyoto Encyclopedia of Genes and Genomes pathway and Gene Ontology analyses were performed to examine the biological functions of the DEGs. Additionally, a protein-protein interaction (PPI) network was constructed to screen hub gene modules using Cytoscape software. By predicting target genes of the DEMs and integrating them with DEGs, the present study constructed a miRNA-mRNA regulatory network to predict the possible molecular mechanism of WT. Expression of hub genes was validated using the Oncomine database. A total of 613 genes and 29 miRNAs were identified to be differentially expressed in WT. By constructing a PPI network and screening hub gene modules, 5 upregulated genes, including BUB1 mitotic checkpoint serine/threonine kinase, BUB1B mitotic checkpoint serine/threonine kinase B, cell division cycle protein 45, cyclin B2 and pituitary tumor-transforming 1. These genes were identified to be associated with the cell cycle pathway, which suggested that these genes may serve important roles in WT. In addition, a miRNA-mRNA regulatory network was constructed and comprised 16 DEMs and 19 DEGs. In conclusion, key genes, miRNAs and the mRNA-miRNA regulatory network identified in the present study may improve understanding of the underlying molecular mechanisms in the occurrence and development of WT, and may aid the identification of potential biomarkers and therapeutic targets.
| S-EPMC6755433 | biostudies-literature
Project description:Breast cancer (BC) is the leading malignancy in women worldwide, yet relatively little is known about the genes and signaling pathways involved in BC tumorigenesis and progression. The present study aimed to elucidate potential key candidate genes and pathways in BC. Five gene expression profile data sets (GSE22035, GSE3744, GSE5764, GSE21422 and GSE26910) were downloaded from the Gene Expression Omnibus (GEO) database, which included data from 113 tumorous and 38 adjacent non?tumorous tissue samples. Differentially expressed genes (DEGs) were identified using t?tests in the limma R package. These DEGs were subsequently investigated by pathway enrichment analysis and a protein?protein interaction (PPI) network was constructed. The most significant module from the PPI network was selected for pathway enrichment analysis. In total, 227 DEGs were identified, of which 82 were upregulated and 145 were downregulated. Pathway enrichment analysis results revealed that the upregulated DEGs were mainly enriched in 'cell division', the 'proteinaceous extracellular matrix (ECM)', 'ECM structural constituents' and 'ECM?receptor interaction', whereas downregulated genes were mainly enriched in 'response to drugs', 'extracellular space', 'transcriptional activator activity' and the 'peroxisome proliferator?activated receptor signaling pathway'. The PPI network contained 174 nodes and 1,257 edges. DNA topoisomerase 2?a, baculoviral inhibitor of apoptosis repeat?containing protein 5, cyclin?dependent kinase 1, G2/mitotic?specific cyclin?B1 and kinetochore protein NDC80 homolog were identified as the top 5 hub genes. Furthermore, the genes in the most significant module were predominantly involved in 'mitotic nuclear division', 'mid?body', 'protein binding' and 'cell cycle'. In conclusion, the DEGs, relative pathways and hub genes identified in the present study may aid in understanding of the molecular mechanisms underlying BC progression and provide potential molecular targets and biomarkers for BC.
| S-EPMC5983982 | biostudies-literature
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.
| S-EPMC6625394 | biostudies-literature
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.
| S-EPMC6924121 | biostudies-literature
Project description:Colon cancer is one of the most malignant cancers worldwide. Nearly 20% of all colon cancer patients are diagnosed at stage IV (metastasis). However, further study of colon cancer is difficult due to a lack of understanding of its pathogenesis. In this study, we acquired high-throughput sequence data from TCGA datasets and performed integrated bioinformatic analysis including differential gene expression analysis, gene ontology and KEGG pathways analysis, protein-protein analysis, survival analysis, and multivariate Cox proportional hazards regression analysis in order to identify a panel of key candidate genes involved in the metastasis of colon cancer. We then constructed a prognostic signature based on the expression of REG1B, TGM6, NTF4, PNMA5, and HOXC13 which could provide significant prognostic value for colon cancer.
| S-EPMC6536975 | biostudies-literature
Project description: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.
| S-EPMC8797306 | biostudies-literature
Project description:Background The purpose of the present study was to investigate the molecular mechanisms of tamoxifen resistance in breast cancer and to identify potential targets for antitamoxifen resistance. Methods Differentially expressed genes (DEGs) in tamoxifen-resistant and tamoxifen-sensitive breast cancer cells were assessed using the GSE67916 dataset acquired from the Gene Expression Omnibus database. Gene ontology (GO) and pathway enrichment analyses were applied to investigate the functions and pathways of the DEGs. Subsequently, the protein-protein interaction (PPI) network was constructed using the Search Tool for the Retrieval of Interacting Genes (STRING), and subnetworks were further analyzed by Molecular Complex Detection (MCODE). The PPI network and subnetworks were visualized using Cytoscape software. Results In total, 438 DEGs were identified, of which 300 were upregulated and 138 were downregulated. The DEGs were significantly enriched in the protein binding, cellular response to estradiol stimulus, and immune response GO terms while the most significant pathways included the mitogen-activated protein kinase (MAPK) signaling pathway in cancer. The PPI network of DEGs was constructed with 288 nodes and 629 edges, and 2 subnetworks were screened out from the entire network. Conclusions A number of significant hub DEGs were identified based on their degree of connectivity in the PPI network, , included MAPK1 (node degree 36), ESR1 (node degree 27), SMARCA4 (node degree 27), RANBP2 (node degree 25), and PRKCA (node degree 21). These critical hub genes were found to be related to tamoxifen resistance in breast cancer. The results of this study further the understanding of tamoxifen resistance at the molecular level and identify potential therapeutic targets for tamoxifen-resistant breast cancer.
| S-EPMC8798269 | biostudies-literature
Project description:Different drug combinations including irinotecan remain some of the most important therapeutic modalities in treating colorectal cancer (CRC). However, chemotherapy often leads to the acquisition of cancer drug resistance. To bridge the gap between in vitro and in vivo models, we compared the mRNA expression profiles of CRC cell lines (HT29, HTC116, and LoVo and their respective irinotecan-resistant variants) with patient samples to select new candidate genes for the validation of irinotecan resistance. Data were downloaded from the Gene Expression Omnibus (GEO) (GSE42387, GSE62080, and GSE18105) and the Human Protein Atlas databases and were subjected to an integrated bioinformatics analysis. The protein-protein interaction (PPI) network of differently expressed genes (DEGs) between FOLFIRI-resistant and -sensitive CRC patients delivered several potential irinotecan resistance markers: NDUFA2, SDHD, LSM5, DCAF4, COX10 RBM8A, TIMP1, QKI, TGOLN2, and PTGS2. The chosen DEGs were used to validate irinotecan-resistant cell line models, proving their substantial phylogenetic heterogeneity. These results indicated that in vitro models are highly limited and favor different mechanisms than in vivo, patient-derived ones. Thus, cell lines can be perfectly utilized to analyze specific mechanisms on their molecular levels but cannot mirror the complicated drug resistance network observed in patients.
| S-EPMC9312838 | biostudies-literature
Project description:Purpose:Triple-negative breast cancer refers to breast cancer that does not express estrogen receptor (ER), progesterone receptor (PR), or human epidermal growth factor receptor 2 (Her2). This study aimed to identify the key pathways and genes and find the potential initiation and progression mechanism of triple-negative breast cancer (TNBC). Methods:We downloaded the gene expression profiles of GSE76275 from Gene Expression Omnibus (GEO) datasets. This microarray Super-Series sets are composed of gene expression data from 265 samples which included 67 non-TNBC and 198 TNBC. Next, all the differentially expressed genes (DEGs) with p<0.01 and fold change ?1.5 or ?-1.5 were identified. Result:56 upregulated and 151 downregulated genes were listed and the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) enrichment analysis was performed. These significantly changed genes were mainly involved in the biological process termed prostate gland morphogenesis, inner ear morphogenesis, cell maturation, digestive tract morphogenesis, autonomic nervous system development, monovalent inorganic anion homeostasis, neural crest cell development, regulation of dendrite extension and glial cell proliferation, immune system process termed T cell differentiation, regulation of immune response, and macrophage activation. Genes are mainly involved in the KEGG pathway termed Oocyte meiosis. All DEGs underwent survival analysis using datasets from The Cancer Genome Atlas (TCGA) integrated by cBioPortal, of which amplification of SRY-related HMG-box 8 (SOX8), androgen receptor (AR), and Chromosome 9 Open Reading Frame 152 (C9orf152) were significantly negative while Nik Related Kinase (NRK) and RAS oncogene family 30 (RAB30) were positively correlated to the life expectancy (p<0.05). Conclusions:In conclusion, these pathways and genes identified could help understanding the mechanism of development of TNBC. Besides, SOX8, AR, C9orf152, NRK and RAB30, and other key genes and pathways might be promising targets for the TNBC treatment.
| S-EPMC6098886 | biostudies-literature