Dataset Information

Identification of key genes associated with esophageal adenocarcinoma based on bioinformatics analysis.

ABSTRACT:

Background

Esophageal adenocarcinoma (EAC) is an aggressive malignancy and accounts for the majority of cancer-related death worldwide. It is often diagnosed at an advanced stage and entails a poor prognosis for those afflicted. The mechanisms of its pathogenesis and progress remain unclear and require urgent elucidation. This study aimed to identify specific genes and potential pathways associated with the progression and prognosis of EAC using bioinformatics analyses.

Methods

EAC microarray datasets from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases were analyzed to identify differentially expressed genes (DEGs) using bioinformatics analysis. The DEGs in TCGA were then analyzed to construct a co-expression network by weighted correlation network analysis (WGCNA), and module-clinical trait relationships were analyzed to explore the genes that associated with clinicopathological parameters of EAC. Gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analyses were performed for the cancer-related genes, and a DEG-based protein-protein interaction (PPI) network was used to extract hub genes through Cytoscape plugins. The consensus survival analysis for EAC (OSeac) was performed to identify the prognosis-related genes. The immune infiltration was evaluated by tumor immune estimation resource (TIMER) algorithms, and a risk score prognostic model was established using univariate, multivariate Cox proportional hazards regression, and lasso regression analysis.

Results

Ultimately, 190 cancer-related DEGs were identified, 6 of which were found to play vital roles in the progression of EAC, including ACTA2, BGN, CALD1, COL1A1, COL4A1, and DCN. The risk score prognostic model consisted of 6 other genes that had an important impact on the prognosis of EAC, including CLDN3, EPB41L4A, ESM1, MT1X, PAQR5, and PLAU. The area under the curve of the prognostic model for predicting the survival of patients at 1, 2, and 3 years was 0.707, 0.702, and 0.726, respectively.

Conclusions

This study identified several genes with the potential to become useful targets for the diagnosis and treatment of EAC. The 6-gene-related risk score prognostic model and nomogram based on these genes may be a reliable tool for predicting the prognosis of patients with EAC.

SUBMITTER: Qi W

PROVIDER: S-EPMC8743722 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Esophageal adenocarcinoma (EAC) is one of the most common subtypes of esophageal cancer, and is associated with a low 5-year survival rate. The present study aimed to identify key genes and pathways associated with EAC using bioinformatics analysis. The gene expression profiles of GSE92396, which includes 12 EAC samples and 9 normal esophageal samples, were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) between the EAC and normal samples were identified using the limma package in R language. Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the identified DEGs were conducted using the online analysis tool, the Database for Annotation, Visualization and Integrated Discovery. A protein-protein interaction (PPI) network of the DEGs was constructed using the Search Tool for the Retrieval of Interacting Genes (STRING) database and Cytoscape software. Finally, module analysis was conducted for the PPI network using the MCODE plug-in in Cytoscape. Of the 386 DEGs identified, the 150 upregulated genes were mainly enriched in the KEGG pathways of complement and coagulation cascades, maturity onset diabetes of the young and protein digestion and absorption; and the 236 downregulated genes were mainly enriched in amoebiasis, retinol metabolism and drug metabolism-cytochrome P450. Based on information from the STRING database, a PPI network comprising of 369 nodes and 534 edges was constructed in Cytoscape. The top 10 hub nodes with the highest degrees were determined as interleukin-8, involucrin, tissue inhibitor of metalloproteinase 1, fibronectin 1, serpin family E member 1, serpin family A member 1, cystic fibrosis transmembrane conductance regulator, secreted phosphoprotein 1, collagen type I alpha 1 chain and angiotensinogen. A total of 6 modules were detected from the PPI network that satisfied the criteria of MCODE score >4 and number of nodes >4. KEGG pathways enriched for the module DEGs were mainly within arachidonic acid metabolism, complement and coagulation cascades and rheumatoid arthritis. In conclusion, identification of these key genes and pathways may improve understanding of the mechanisms underlying the development of EAC, and may be used as diagnostic and therapeutic targets in EAC.

Project description:The aim of the present study was to identify the differentially expressed genes (DEGs) in esophageal squamous-cellcarcinoma (ESCC) and provide potential therapeutic targets. The microarray dataset GSE20347 was downloaded from the Gene Expression Omnibus (GEO) database, and included 17 tissue samples and 13 normal adjacent tissue samples from patients with ESCC. A total of 22,277 DEGs were identified. A heat map for the DEGs was constructed with the Morpheus online tool and the top 200 genes (100 upregulated and 100 downregulated) were selected for further bioinformatics analysis, including analysis of gene ontology (GO) terms, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, protein-protein interaction networks and Spearman's correlation tests. The results of the GO analysis indicated that the upregulated DEGs were most significantly enriched in membrane-bounded vesicles in the cellular component (CC) category, but were not significantly enriched in any GO terms of the categories biological process (BP) or molecular function (MF); furthermore, the downregulated DEGs were most significantly enriched in regulation of DNA metabolic processes, nucleotide binding and chromosomes in the categories BP, MF and CC, respectively. The KEGG analysis indicated that the downregulated DEGs were enriched in the regulation of cell cycle pathways. The top 10 hub proteins in the protein-protein interaction network were cyclin-dependent kinase 4, budding uninhibited by benzimidazoles 1, cyclin B2, heat shock protein 90AA1, aurora kinase A, H2A histone family member Z, replication factor C subunit 4, and minichromosome maintenance complex component 2, -4 and -7. These proteins are mainly involved in regulating tumor progression. The genes in the four top modules were mainly implicated in regulating cell cycle pathways. Secreted Ly-6/uPAR-related protein (SLURP) was the hub gene, and SLURP and its interacting genes were most enriched in the chromosomal part in the CC category, organelle organization in the BP category and protein binding in the MF category, and were involved in pathways including DNA replication, cell cycle and P53 signaling. The expression of SLURP-1 in fifteen patients with esophageal carcinoma was detected using quantitative polymerase chain reaction analysis, and the results indicated that SLURP-1 expression was significantly decreased in the tumor samples relative to that in normal adjacent tissues. These results suggest that several hub proteins and the hub gene SLURP-1 may serve as potential therapeutic targets, and that gene dysfunction may be involved in the tumorigenesis of ESCC.

Project description:The underlying causes of esophageal cancer (EC) are unknown. To explore the molecular mechanisms that lead to EC, gene expression profiles of large cohorts of patients with EC were obtained from The Cancer Genome Atlas and the Gene Expression Omnibus (GEO) databases (GSE5364, GSE20347 and GSE23400). The present study identified 83 upregulated and 22 downregulated genes between EC and normal tissue using R statistical software and the GEO2R web tool. The Database for Annotation, Visualization and Integrated Discovery was used to identify the associated pathways, and for functional annotation of the differentially expressed genes (DEGs). Protein-protein interactions of these DEGs were analyzed based on the Search Tool for the Retrieval of Interacting Genes database, and hub genes were visualized using Cytoscape software. An online Kaplan-Meier plotter survival analysis tool was utilized to evaluate the prognostic value of hub gene expression in patients with EC. Further analysis of an additional dataset from GEO (GSE21293) revealed that these genes were associated with infiltration and metastasis in EC. In addition, the Gene Expression Profiling Interactive Analysis tool was used to evaluate expression levels of hub genes in patients with EC for different pathological stages. The Ualcan analysis tool was used to evaluate the expression levels of hub genes for different histological types. Overall, ubiquitin conjugating enzyme E2 C, cyclin dependent kinase inhibitor 3, CDC28 protein kinase regulatory subunit 2, kinesin family member 20A (KIF20A) and RAD51 associated protein 1 (RAD51AP1) were upregulated in EC tissues compared with normal tissues, and upregulation of these genes was a poor prognostic factor for patients with EC, indicating that these genes may mediate EC cell infiltration and metastasis. Among the hub genes, KIF-20A had potential value for predicting the pathological stage of EC. KIF20A and RAD51AP1 were more informative biomarkers of esophageal squamous cell carcinoma. Further studies are required to explore the value of these genes in the treatment of EC.

Project description:BackgroundHepatocellular carcinoma (HCC) is a type of primary liver tumor with poor prognosis and high mortality, and its molecular mechanism remains incompletely understood. This study aimed to use bioinformatics technology to identify differentially expressed genes (DEGs) in HCC pathogenesis, hoping to identify novel biomarkers or potential therapeutic targets for HCC research.MethodsThe bioinformatics analysis of our research mostly involved the following two datasets: Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA). First, we screened DEGs based on the R packages (limma and edgeR). Using the DAVID database, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of DEGs were carried out. Next, the protein-protein interaction (PPI) network of the DEGs was built in the STRING database. Then, hub genes were screened through the cytoHubba plug-in, followed by verification using the GEPIA and Oncomine databases. We demonstrated differences in levels of the protein in hub genes using the Human Protein Atlas (HPA) database. Finally, the hub genes prognostic values were analyzed by the GEPIA database. Additionally, using the Comparative Toxicogenomics Database (CTD), we constructed the drug-gene interaction network.ResultsWe ended up with 763 DEGs, including 247 upregulated and 516 downregulated DEGs, that were mainly enriched in the epoxygenase P450 pathway, oxidation-reduction process, and metabolism-related pathways. Through the constructed PPI network, it can be concluded that the P53 signaling pathway and the cell cycle are the most obvious in module analysis. From the PPI, we filtered out eight hub genes, and these genes were significantly upregulated in HCC samples, findings consistent with the expression validation results. Additionally, survival analysis showed that high level gene expression of CDC20, CDK1, MAD2L1, BUB1, BUB1B, CCNB1, and CCNA2 were connected with the poor overall survival of HCC patients. Toxicogenomics analysis showed that only topotecan, oxaliplatin, and azathioprine could reduce the gene expression levels of all seven hub genes.ConclusionThe present study screened out the key genes and pathways that were related to HCC pathogenesis, which could provide new insight for the future molecularly targeted therapy and prognosis evaluation of HCC.