Dataset Information

Single sample scoring of hepatocellular carcinoma: A study based on data mining.

ABSTRACT: Hepatocellular carcinoma (HCC) is a high mortality malignancy and the second leading cause of cancer-related deaths. Because the immune system plays a dual role by assisting the host barrier and tumor progression, there are complex interactions with considerable prognostic significance. Herein, we performed single-sample gene set enrichment (ssGSEA) to explore the tumor microenvironment (TME) and quantify the tumor-infiltrating immune cell (TIIC) subgroups of immune responses based on the HCC cohort of The Cancer Genome Atlas (TCGA) database. We evaluate molecular subpopulations, survival, function, and expression differential associations, as well as reveal potential targets, and biomarkers for immunotherapy. We combined the TME score and the 29 immune cell types in the low, medium, and high immunity groups. The stromal score, immune score, and ESTIMATE score were positively correlated with immune activity but negatively correlated with the tumor purity. There were 23 human leukocyte antigen (HLA)-related genes that were significantly different. However, KIAA1429 was not significant among the different immunity groups. Besides, programmed death-ligand 1 (PD-L1) and cytotoxic T-lymphocyte-associated antigen 4 (CTLA-4) expression increased with the increase of immune activity. This may provide valuable information for HCC immunotherapy. We also found that there was no significant difference in naïve B cells, macrophages M1, activated mast cells, resting natural killer (NK) cells, and T cells gamma delta among the different immunity groups. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that the differential proteins were mainly enriched in alpha-linolenic acid (ALA) metabolism, cytokine-cytokine receptor interaction, glycosaminoglycan biosynthesis-heparan sulfate/heparin, glycosphingolipid biosynthesis-ganglio series and proteasome. Our findings provide a deeper understanding of the immune scene, uncovering remarkable immune infiltration patterns of various subtypes of HCC using ssGSEA. This study advances the understanding of immune response and provides a basis for research to enhance immunotherapy.

SUBMITTER: Zhu D

PROVIDER: S-EPMC8168165 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BACKGROUND:Gene set scoring provides a useful approach for quantifying concordance between sample transcriptomes and selected molecular signatures. Most methods use information from all samples to score an individual sample, leading to unstable scores in small data sets and introducing biases from sample composition (e.g. varying numbers of samples for different cancer subtypes). To address these issues, we have developed a truly single sample scoring method, and associated R/Bioconductor package singscore ( https://bioconductor.org/packages/singscore ). RESULTS:We use multiple cancer data sets to compare singscore against widely-used methods, including GSVA, z-score, PLAGE, and ssGSEA. Our approach does not depend upon background samples and scores are thus stable regardless of the composition and number of samples being scored. In contrast, scores obtained by GSVA, z-score, PLAGE and ssGSEA can be unstable when less data are available (NS?<?25). The singscore method performs as well as the best performing methods in terms of power, recall, false positive rate and computational time, and provides consistently high and balanced performance across all these criteria. To enhance the impact and utility of our method, we have also included a set of functions implementing visual analysis and diagnostics to support the exploration of molecular phenotypes in single samples and across populations of data. CONCLUSIONS:The singscore method described here functions independent of sample composition in gene expression data and thus it provides stable scores, which are particularly useful for small data sets or data integration. Singscore performs well across all performance criteria, and includes a suite of powerful visualization functions to assist in the interpretation of results. This method performs as well as or better than other scoring approaches in terms of its power to distinguish samples with distinct biology and its ability to call true differential gene sets between two conditions. These scores can be used for dimensional reduction of transcriptomic data and the phenotypic landscapes obtained by scoring samples against multiple molecular signatures may provide insights for sample stratification.

Project description:The aims of the present study were to examine the potential role of microRNA‑233‑3p (miR)‑223‑3p in the tumorigenesis of hepatocellular carcinoma (HCC), and to investigate its diagnostic accuracy and potential molecular mechanisms. The expression data of miR‑223‑3p in HCC were obtained from the Gene Expression Omnibus (GEO). Data for the precursor miR‑223 were obtained from The Cancer Genome Atlas (TCGA). The diagnostic role of miR‑223‑3p was identified by the receiver operating curve (ROC), and the diagnostic value of miR‑223‑3p in HCC was calculated from qualified reports in the literature. In addition, associated data from the GEO, TCGA and qualified experiments were pooled for comprehensive meta‑analysis. Genes, which intersected between online prediction databases, natural language processing and differentially expressed genes from TCGA were regarded as potential targets of miR‑223‑3p in HCC. The Gene Ontology enrichment analysis and the Kyoto Encyclopedia of Genes and Genomes pathways of potential targets were performed using the Database for Annotation, Visualization and Integrated Discovery. The protein‑protein interactions were mapped using the Search Tool for the Retrieval of Interacting Genes. Among 15 qualified microarray data sets from GEO, seven showed that a significantly lower level of miR‑223‑3p was present in the HCC tissues, compared with that in non‑cancerous tissues (P<0.05). In addition, five GEO data sets revealed diagnostic values of miR‑223‑3p, with an area under the curve (AUC) of >0.80 (P<0.05). The diagnostic accuracy of the precursor miR‑223 in TCGA was also calculated (AUC=0.78, P<0.05). Similarly, the precursor miR‑223 showed a higher level of downregulation in HCC tissues, compared with that in healthy controls in TCGA (P<0.001). A summary ROC was also calculated as 0.89 (95% CI, 0.85‑0.91) in the meta‑analysis. A total of 72 potential targets were extracted, mainly involved in the terms 'microRNAs in cancer', 'ATP binding' and 'prostate cancer'. Five potential target genes were considered the hub genes of miR‑223‑3p in HCC, including checkpoint kinase 1, DNA methyltransferase 1, baculoviral IAP repeat containing 5, kinesin family member 23, and collagen, type I, α1. Based on TCGA, the hub genes were significantly upregulated in HCC (P<0.05). Collectively, these results showed that miR‑223‑3p may be crucial in HCC carcinogenesis showing high diagnostic accuracy, and may be mediated by several hub genes.

Project description:BackgroundChromobox protein homolog 8 (CBX8), a transcriptional repressor, participates in many biological processes in various carcinomas. Cell differentiation, aging, and cell cycle progression are examples of such processes. It is critical to investigate CBX8 expression and its relationship with clinicopathological characteristics in liver hepatocellular carcinoma (LIHC), kidney renal clear cell carcinoma (KIRC), and ovarian cancer (OV) to investigate CBX8's potential diagnostic and prognostic values.MethodsTCGA and CPTAC databases were used to compare the data between cancer and matched normal tissues on RNA and protein expression profiles and their relevant clinical information to determine the relationship between CBX8 and clinicopathological features. Kaplan-Meier analyses were used to assess CBX8 relationship's with disease-free survival (DFS), relapse-free survival (RFS), and overall survival (OS). The multivariate Cox regression analysis was used to identify independent risk factors which affect prognosis. DNA methylation and genetic changes and their impact on prognoses were evaluated by cBioPortal and MethSurv websites. Spearman's correlation was used to determine the relationship of CBX8 expression with somatic mutation. Tumor immune estimation resource (TIMER) was adopted to investigate the relationship between CBX8 and immune cell infiltration (ICI). CBX8-relevant genes and proteins are analyzed by EnhancedVolcano and STRING databases. The gene set enrichment analysis (GSEA) was performed to investigate the potential functions of CBX8.ResultsCBX8 RNA and protein overexpression were confirmed in LIHC, KIRC, and OV (p < 0.05). High CBX8 was significantly related to the clinical features and poor prognoses. The CBX8 genetic alteration rate was 3%. DNA methylation was also associated with prognoses. CBX8 closely interacted with ICI, TMB, MSI, purity, and ploidy. GO analyses revealed that CBX8-associated genes were enriched in biological processes, cell cycle regulation, and molecular functions. KEGG analyses exhibited that CBX8 was gathered in signaling pathway regulation. GSEA revealed that cell cycle, DNA replication, and Wnt signaling pathways were differentially enriched in the high CBX8 expression group.ConclusionsCBX8 could be a potential diagnostic and prognostic biomarker for LIHC, KIRC, and OV cancers.

Project description:Several studies have demonstrated that chronic hepatitis delta virus (HDV) infection is associated with a worsening of hepatitis B virus (HBV) infection and increased risk of hepatocellular carcinoma (HCC). However, there is limited data on the role of HDV in the oncogenesis of HCC. This study is aimed at assessing the potential mechanisms of HDV-associated hepatocarcinogenesis, especially to screen and identify key genes and pathways possibly involved in the pathogenesis of HCC. We selected three microarray datasets: GSE55092 contains 39 cancer specimens and 81 paracancer specimens from 11 HBV-associated HCC patients, GSE98383 contains 11 cancer specimens and 24 paracancer specimens from 5 HDV-associated HCC patients, and 371 HCC patients with the RNA-sequencing data combined with their clinical data from the Cancer Genome Atlas (TCGA). Afterwards, 948 differentially expressed genes (DEGs) closely related to HDV-associated HCC were obtained using the R package and filtering with a Venn diagram. We then performed gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis to determine the biological processes (BP), cellular component (CC), molecular function (MF), and KEGG signaling pathways most enriched for DEGs. Additionally, we performed Weighted Gene Coexpression Network Analysis (WGCNA) and protein-to-protein interaction (PPI) network construction with 948 DEGs, from which one module was identified by WGCNA and three modules were identified by the PPI network. Subsequently, we validated the expression of 52 hub genes from the PPI network with an independent set of HCC dataset stored in the Gene Expression Profiling Interactive Analysis (GEPIA) database. Finally, seven potential key genes were identified by intersecting with key modules from WGCNA, including 3 reported genes, namely, CDCA5, CENPH, and MCM7, and 4 novel genes, namely, CDC6, CDC45, CDCA8, and MCM4, which are associated with nucleoplasm, cell cycle, DNA replication, and mitotic cell cycle. The CDCA8 and stage of HCC were the independent factors associated with overall survival of HDV-associated HCC. All the related findings of these genes can help gain a better understanding of the role of HDV in the underlying mechanism of HCC carcinogenesis.

Dataset Information

Single sample scoring of hepatocellular carcinoma: A study based on data mining.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets