Dataset Information

Bioinformatics Analysis Using ATAC-seq and RNA-seq for the Identification of 15 Gene Signatures Associated With the Prediction of Prognosis in Hepatocellular Carcinoma

ABSTRACT:

Background

Gene expression (RNA-seq) and overall survival (OS) in TCGA were combined using chromosome accessibility (ATAC-seq) to search for key molecules affecting liver cancer prognosis.

Methods

We used the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) to analyse chromatin accessibility in the promoter regions of whole genes in liver hepatocellular carcinoma (LIHC) and then screened differentially expressed genes (DEGs) at the mRNA level by transcriptome sequencing technology (RNA-seq). We obtained genes significantly associated with overall survival (OS) by a one-way Cox analysis. The three were screened by taking intersection and further using a Kaplan–Meier (KM) for validation. A prognostic model was constructed using the obtained genes by LASSO regression analysis.The expression of these genes in hepatocellular carcinomas was then analysed. The protein expression of these genes was verified using the Human Protein Atlas(HPA) online datasets and immunohistochemistry.

Results

ATAC-seq, RNA-seq and survival analysis, combined with a LASSO prediction model, identified signatures of 15 genes (PRDX6, GCLM, HTATIP2, SEMA3F, UCK2, NOL10, KIF18A, RAP2A, BOD1, GDI2, ZIC2, GTF3C6 SLC1A5, ERI3 and SAC3D1), all of which were highly expressed in hepatocellular carcinoma. The LASSO prognostic model showed that this risk score had high predictive accuracy for the survival prognosis at 1, 3 and 5 years. A KM curve analysis showed that high expression of all 15 gene signatures was significantly associated with a poor prognosis in LIHC patients. HPA analysis of protein expression showed that PRDX6, GCLM, HTATIP2, NOL10, KIF18A, RAP2A and GDI2 were highly expressed in the hepatocellular carcinoma tissues compared with normal control tissues.

Conclusions

PRDX6, GCLM, HTATIP2, SEMA3F, UCK2, NOL10, KIF18A, RAP2A, BOD1, GDI2, ZIC2, GTF3C6, SLC1A5, ERI3 and SAC3D1 may affect the prognosis of LIHC.

SUBMITTER: Yang H

PROVIDER: S-EPMC8573251 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundAs an iron-dependent type of programmed cell death, ferroptosis plays an important role in the pathogenesis and progression of hepatocellular carcinoma (HCC). Long noncoding RNAs (lncRNAs) have been linked to the prognosis of patients with HCC in a number of studies. Nevertheless, the predictive value of lncRNAs (FRLs) associated with ferroptosis in HCC has not been fully elucidated.MethodsDownload RNA sequencing data and clinical profiles of HCC patients from The Cancer Genome Atlas (TCGA) database. The FRLs associated with prognosis were determined by Pearson's correlation analysis. After that, prognostic signature for FRLs was established using Cox and LASSO regression analyses. Meanwhile, survival analysis, correlation analysis of clinicopathological features, Cox regression, receiver operating characteristic (ROC) curve, and nomogram were used to analyze the FRL signature's predictive capacity. The relationship between signature risk score, immune cell infiltration, and chemotherapy drug sensitivity is further studied.ResultsIn total, 93 FRLs were found to be of prognostic value in patients with HCC. A five-FRL signature comprising AC015908.3, LINC01138, AC009283.1, Z83851.1, and LUCAT1 was created in order to enhance the prognosis prediction with HCC patients. The signature demonstrated a good predictive potency, according to the Kaplan-Meier and ROC curves. The five-FRL signature was found to be a risk factor independent of various clinical factors using Cox regression and stratified survival analysis. The high-risk group was shown to be enriched in tumorigenesis and immune-related pathways according to GSEA analysis. Additionally, immune cell infiltration, immune checkpoint molecules, and half-inhibitory concentrations differed considerably between risk groups, implying that this signature could be used to evaluate the clinical efficacy of chemotherapy and immunotherapy.ConclusionThe five-FRL risk signature is helpful for assessing the prognosis of HCC patients and improving therapy options, so it can be further applied clinically.

Project description:BackgroundEvidence from prevailing studies show that hepatocellular carcinoma (HCC) is among the top cancers with high mortality globally. Gene regulation at post-transcriptional level orchestrated by RNA-binding proteins (RBPs) is an important mechanism that modifies various biological behaviors of HCC. Currently, it is not fully understood how RBPs affects the prognosis of HCC. In this study, we aimed to construct and validate an RBP-related model to predict the prognosis of HCC patients.MethodsDifferently expressed RBPs were identified in HCC patients based on the GSE54236 dataset from the Gene Expression Omnibus (GEO) database. Integrative bioinformatics analyses were performed to select hub genes. Gene expression patterns were validated in The Cancer Genome Atlas (TCGA) database, after which univariate and multivariate Cox regression analyses, as well as Kaplan-Meier analysis were performed to develop a prognostic model. Then, the performance of the prognostic model was assessed using receiver operating characteristic (ROC) curves and clinicopathological correlation analysis. Moreover, data from the International Cancer Genome Consortium (ICGC) database were used for external validation. Finally, a nomogram combining clinicopathological parameters and prognostic model was established for the individual prediction of survival probability.ResultsThe prognostic risk model was finally constructed based on two RBPs (BOP1 and EZH2), facilitating risk-stratification of HCC patients. Survival was markedly higher in the low-risk group relative to the high-risk group. Moreover, higher risk score was associated with advanced pathological grade and late clinical stage. Besides, the risk score was found to be an independent prognosis factor based on multivariate analysis. Nomogram including the risk score and clinical stage proved to perform better in predicting patient prognosis.ConclusionsThe RBP-related prognostic model established in this study may function as a prognostic indicator for HCC, which could provide evidence for clinical decision making.

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.