scholarly journals Identification and analysis of immune-related subtypes of hepatocellular carcinoma

2020 ◽  
pp. 153537022097013
Author(s):  
Qimeng Wang ◽  
Jin Huang ◽  
Huihua Zhang ◽  
Huan Liu ◽  
Min Yu

Hepatocellular carcinoma is a malignance that remains difficult to cure. Immunotherapy has shown its potential application in a variety of refractory malignancies. Due to the complexity of immune microenvironment of hepatocellular carcinoma, the efficacy of immunotherapy for hepatocellular carcinoma is not as effective as expected. Expression data of hepatocellular carcinoma from the TCGA and ICGC databases were used for classification and verification of hepatocellular carcinoma subtypes. The immune-related functions and pathways were identified via gene set enrichment analysis, while the sections denoting the subsets of the immune cells were estimated using the CIBERSORT algorithm. Immunity low (Immunity_L), immunity medium (Immunity_M), and immunity high (Immunity_H) were specified as the three immune-related subtypes of hepatocellular carcinoma. The quantity of stromal and immune cells was the most substantial in Immunity_H, compared to the other subtypes. Interestingly, the proportion of M0 macrophages decreased from Immunity_L to Immunity_H, while the proportion of CD8 T cells increased. Furthermore, the HLA genes expression levels, as well as those of six immune checkpoint genes were substantially lower in Immunity_L than in Immunity_H. Functional analysis was performed for 1512 differentially expressed genes between Immunity_L and Immunity_H. Finally, the PPI network was constructed with 118 nodes. The highest connectivity degree nodes were B2M, HLA-DRA, and HLA-DRB1. The above results were verified in ICGC-JP and ICGC-FR databases with a consistent trend. In this study, we divided hepatocellular carcinoma into three subtypes and explored the immune-related characteristics of these subtypes. These results may provide new insights for immunotherapy of hepatocellular carcinoma.

2022 ◽  
Vol 11 ◽  
Author(s):  
Fahui Liu ◽  
Jiadong Liang ◽  
Puze Long ◽  
Lilan Zhu ◽  
Wanyun Hou ◽  
...  

Hepatocellular carcinoma (HCC) is one of the common malignant tumors. The prognosis and five-year survival rate of HCC are not promising due to tumor recurrence and metastasis. Exploring markers that contribute to the early diagnosis of HCC, markers for prognostic evaluation of HCC patients, and effective targets for treating HCC patients are in the spotlight of HCC therapy. Zinc Finger CCHC-Type Containing 17 (ZCCHC17) encodes the RNA binding protein ZCCHC17, but its role in HCC is still unclear. Here, 90 paraffin-embedded specimens combined with bioinformatics were used to comprehensively clarify the value of ZCCHC17 in the diagnosis and prognosis of HCC and its potential functions. Paraffin-embedded specimens were used to assess ZCCHC17 protein expression and its correlation with prognosis in 90 HCC patients. the public data sets of HCC patients from TCGA, ICG, and GEO databases were also used for further analysis. It was found that protein and mRNA levels of ZCCHC17 in HCC tissues were significantly higher than those in normal tissues. The abnormally high expression may be related to the abnormal DNA methylation of ZCCHC17 in tumor tissues. The high expression of ZCCHC17 is related to AFP, histologic grade, tumor status, vascular invasion, and pathological stage. Multi-data set analysis showed that patients with high ZCCHC17 expression had a worse prognosis, and multivariate cox regression analysis showed an independent prognostic significance of ZCCHC17. The results of functional analysis, including Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Set Enrichment Analysis (GSEA), indicate that ZCCHC17 is mainly involved in immune regulation. Subsequently, further single-sample gene set enrichment analysis (ssGSEA) showed that the expression of ZCCHC17 was related to the infiltration of immune cells. Importantly, we also analyzed the relationship between ZCCHC17 and immune checkpoint genes, tumor mutation burden (TMB), microsatellite instability (MSI) and TP53 status in HCC patients and evaluated the role of ZCCHC17 in cancer immunotherapy. In summary, ZCCHC17 is a novel marker for the diagnosis and prognostic evaluation of HCC. Concurrently, it regulates immune cells in the tumor microenvironment (TME) of HCC patients, which has a specific reference value for the immunotherapy of HCC.


2021 ◽  
Vol 44 (3) ◽  
pp. E32-44
Author(s):  
Jia Shen ◽  
Ming Shu ◽  
Shujie Xie ◽  
Jia Yan ◽  
Kaile Pan ◽  
...  

Purpose: This study aimed to screen hepatitis B virus (HBV)-associated hepatocellular carcinoma (HCC)-related feature ribonucleic acids (RNAs) and to establish a prognostic model. Methods: The transcriptome expression data of HBV-associated HCC were downloaded from The Cancer Genome Atlas (TCGA) database and Gene Expression Omnibus database. Differential RNAs between HBV-associated HCC and normal controls were identified by a meta-analysis of TCGA, GSE55092 and GSE121248. Weighted gene co-expression network analysis was performed to identify key RNAs and modules. A prognostic score model was established using TCGA as a training set by Cox regression analysis and was validated in E-TABM-36 dataset. Additionally, independent prognostic clinical factors were screened, and the function of lncRNAs waspredicted through Gene Set Enrichment Analysis. Results: A total of 710 consistent differential RNAs between HBV-associated HCC and normal controls were obtained, including five lncRNAs and 705 mRNAs. An optimized combination of six differential RNAs (DSCR4, DBH, ECM1, GDAP1, MATR3 and RFC4) was selected and a prognostic score model was constructed. Kaplan-Meier analysis demonstrated that the prognosis of the high-risk and low-risk groups separated by this model was significantly different in the training set and the validation set. Gene Set Enrichment Analysis showed that the co-expression genes of DSCR4 were significantly correlated with neuroactive ligand receptor interactionpathway. Conclusion: A prognostic model based on DSCR4, DBH, ECM1, GDAP1, MATR3 and RFC4 was developed that can accurately predict the prognosis of patients with HBV-associated HCC. These genes, as well as histologic grade, may serve as independent prognostic factors in HBV-associated HCC.


2021 ◽  
Vol 14 (1) ◽  
Author(s):  
Yang Shen ◽  
Li-rong Xu ◽  
Xiao Tang ◽  
Chang-po Lin ◽  
Dong Yan ◽  
...  

Abstract Background Atherosclerosis is a chronic inflammatory disease that affects multiple arteries. Numerous studies have shown the inherent immune diversity in atheromatous plaques and suggest that the dysfunction of different immune cells plays an important role in atherosclerosis. However, few comprehensive bioinformatics analyses have investigated the potential coordinators that might orchestrate different immune cells to exacerbate atherosclerosis. Methods Immune infiltration of 69 atheromatous plaques from different arterial beds in GSE100927 were explored by single-sample-gene-set enrichment analysis (presented as ssGSEA scores), ESTIMATE algorithm (presented as immune scores) and CIBERSORT algorithm (presented as relative fractions of 22 types of immune cells) to divide these plaques into ImmuneScoreL cluster (of low immune infiltration) and ImmuneScoreH cluster (of high immune infiltration). Subsequently, comprehensive bioinformatics analyses including differentially-expressed-genes (DEGs) analysis, protein–protein interaction networks analysis, hub genes analysis, Gene-Ontology-terms and KEGG pathway enrichment analysis, gene set enrichment analysis, analysis of expression profiles of immune-related genes, correlation analysis between DEGs and hub genes and immune cells were conducted. GSE28829 was analysed to cross-validate the results in GSE100927. Results Immune-related pathways, including interferon-related pathways and PD-1 signalling, were highly enriched in the ImmuneScoreH cluster. HLA-related (except for HLA-DRB6) and immune checkpoint genes (IDO1, PDCD-1, CD274(PD-L1), CD47), RORC, IFNGR1, STAT1 and JAK2 were upregulated in the ImmuneScoreH cluster, whereas FTO, CRY1, RORB, and PER1 were downregulated. Atheromatous plaques in the ImmuneScoreH cluster had higher proportions of M0 macrophages and gamma delta T cells but lower proportions of plasma cells and monocytes (p < 0.05). CAPG, CECR1, IL18, IGSF6, FBP1, HLA-DPA1 and MMP7 were commonly related to these immune cells. In addition, the advanced-stage carotid plaques in GSE28829 exhibited higher immune infiltration than early-stage carotid plaques. Conclusions Atheromatous plaques with higher immune scores were likely at a more clinically advanced stage. The progression of atherosclerosis might be related to CAPG, IGSF6, IL18, CECR1, FBP1, MMP7, FTO, CRY1, RORB, RORC, PER1, HLA-DPA1 and immune-related pathways (IFN-γ pathway and PD-1 signalling pathway). These genes and pathways might play important roles in regulating immune cells such as M0 macrophages, gamma delta T cells, plasma cells and monocytes and might serve as potential therapeutic targets for atherosclerosis.


2021 ◽  
Vol 11 ◽  
Author(s):  
Dengliang Lei ◽  
Yue Chen ◽  
Yang Zhou ◽  
Gangli Hu ◽  
Fang Luo

BackgroundHepatocellular carcinoma (HCC) is one of the world’s most prevalent and lethal cancers. Notably, the microenvironment of tumor starvation is closely related to cancer malignancy. Our study constructed a signature of starvation-related genes to predict the prognosis of liver cancer patients.MethodsThe mRNA expression matrix and corresponding clinical information of HCC patients were obtained from the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). Gene set enrichment analysis (GSEA) was used to distinguish different genes in the hunger metabolism gene in liver cancer and adjacent tissues. Gene Set Enrichment Analysis (GSEA) was used to identify biological differences between high- and low-risk samples. Univariate and multivariate analyses were used to construct prognostic models for hunger-related genes. Kaplan-Meier (KM) and receiver-operating characteristic (ROC) were used to assess the model accuracy. The model and relevant clinical information were used to construct a nomogram, protein expression was detected by western blot (WB), and transwell assay was used to evaluate the invasive and metastatic ability of cells.ResultsFirst, we used univariate analysis to identify 35 prognostic genes, which were further demonstrated to be associated with starvation metabolism through Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO). We then used multivariate analysis to build a model with nine genes. Finally, we divided the sample into low- and high-risk groups according to the median of the risk score. KM can be used to conclude that the prognosis of high- and low-risk samples is significantly different, and the prognosis of high-risk samples is worse. The prognostic accuracy of the 9-mRNA signature was also tested in the validation data set. GSEA was used to identify typical pathways and biological processes related to 9-mRNA, cell cycle, hypoxia, p53 pathway, and PI3K/AKT/mTOR pathway, as well as biological processes related to the model. As evidenced by WB, EIF2S1 expression was increased after starvation. Overall, EIF2S1 plays an important role in the invasion and metastasis of liver cancer.ConclusionsThe 9-mRNA model can serve as an accurate signature to predict the prognosis of liver cancer patients. However, its mechanism of action warrants further investigation.


2014 ◽  
Vol 13s1 ◽  
pp. CIN.S13882 ◽  
Author(s):  
Binghuang Cai ◽  
Xia Jiang

Analyzing biological system abnormalities in cancer patients based on measures of biological entities, such as gene expression levels, is an important and challenging problem. This paper applies existing methods, Gene Set Enrichment Analysis and Signaling Pathway Impact Analysis, to pathway abnormality analysis in lung cancer using microarray gene expression data. Gene expression data from studies of Lung Squamous Cell Carcinoma (LUSC) in The Cancer Genome Atlas project, and pathway gene set data from the Kyoto Encyclopedia of Genes and Genomes were used to analyze the relationship between pathways and phenotypes. Results, in the form of pathway rankings, indicate that some pathways may behave abnormally in LUSC. For example, both the cell cycle and viral carcinogenesis pathways ranked very high in LUSC. Furthermore, some pathways that are known to be associated with cancer, such as the p53 and the PI3K-Akt signal transduction pathways, were found to rank high in LUSC. Other pathways, such as bladder cancer and thyroid cancer pathways, were also ranked high in LUSC.


2021 ◽  
Vol 11 ◽  
Author(s):  
Yuan Nie ◽  
Mei-chun Jiang ◽  
Cong Liu ◽  
Qi Liu ◽  
Xuan Zhu

BackgroundsTumor microenvironment (TME) plays a crucial role in the initiation and progression of Hepatocellular Carcinoma (HCC), especially immune infiltrates. However, there is still a challenge in understanding the modulation of the immune and stromal components in TME, especially TME related genes.MethodsThe proportion of tumor-infiltrating immune cells (TICs) and the immune and stromal scores in 374 HCC patients from The Cancer Genome Atlas (TCGA) database were determined using CIBERSORT and ESTIMATE computational methods. The final screened genes were confirmed by the PPI network and univariate Cox regression of the differentially expressed genes based on different immune or stromal scores. The correlation between the expression levels of the final gene interactions and the clinical characteristics was based on TCGA database and local hospital data. Gene set enrichment analysis (GSEA) and the effect of CXCL5 expression on TICs were conducted.ResultsThere were correlations between the expression of CXCL5 and survival of HCC patients and TMN classification both in TCGA database and local hospital data. The immune-related activities were enriched in the high-expression group; however, the metabolic pathways were enriched in the low-expression group. The result of CIBERSORT analyzing had indicated that CXCL5 expression were correlated with the proportion of NK cells activated, macrophages M0, Mast cells resting, Neutrophils.ConclusionsCXCL5 was a potential prognostic marker for HCC and provides clues regarding immune infiltrates, which offers extra insight for therapeutics of HCC, however, more independent cohorts and functional experiments of CXCL5 are warranted.


2020 ◽  
Author(s):  
Chen Xu ◽  
Ling-bing Meng ◽  
Yu Xiao ◽  
Yong Qiu ◽  
Ying-jue Du ◽  
...  

Abstract Background Osteoarthritis (OA) is a chronic, progressive, inflammatory, degenerative disease, which has become an osteoarthropathy that seriously affects physical health and quality of life of elderly people. However, the etiology and pathogenesis of OA remains unclear. Therefore, the study purposed to utilize bioinformatics technology to perform identification and functional enrichment analysis of differentially expressed genes in osteoarthritis. Method The main methods of this study consist of access to microarray data (GSE82107 and GSE55235), identification of differently expressed genes (DEGs) by GEO2R between OA and normal synovium samples, enrichment analysis of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) by Gene Set Enrichment Analysis (GSEA), construction and analysis of protein-protein interaction (PPI) network, significant module and hub genes. Result A total of 300 DEGs were identified, consisting of 64 up-regulated genes and 11 down-regulated genes in OA samples compared to normal synovium tissues. Gene set enrichment analysis of DEGs provided a comprehensive overview of some major pathophysiological mechanisms in OA: cellular response to hydrogen peroxide, P53 signaling pathway and so on. The study also built the PPI network, and a total of 10 key genes were identified: CYR61, PENK, GOLM1, DUSP1, ATF3, STC2, FOSB, PRSS23, TF, and TNC. Conclusion DEGs exists between OA patients and normal cartilage tissue, which may be involved in the related mechanism of OA development, especially cellular response to hydrogen peroxide and CYR61.


Cells ◽  
2021 ◽  
Vol 10 (3) ◽  
pp. 488
Author(s):  
Yoshihisa Tokumaru ◽  
Masanori Oshi ◽  
Ankit Patel ◽  
Wanqing Tian ◽  
Li Yan ◽  
...  

The colon adenoma–carcinoma sequence is a multistep genomic-altering process that occurs during colorectal cancer (CRC) carcinogenesis. Organoids are now commonly used to model both non-cancerous and cancerous tissue. This study aims to investigate how well organoids mimic tissues in the adenoma–carcinoma sequence by comparing their transcriptomes. A total of 234 tissue samples (48 adenomas and 186 CRC) and 60 organoid samples (15 adenomas and 45 CRC) were analyzed. We found that cell-proliferation-related gene sets were consistently enriched in both CRC tissues and organoids compared to adenoma tissues and organoids by gene set enrichment analysis (GSEA). None of the known pathways in the colon adenoma–carcinoma sequence were consistently enriched in CRC organoids. There was no enrichment of the tumor microenvironment-related gene sets in CRC organoids. CRC tissues enriched immune-response-related gene sets, whereas CRC organoids did not. The proportions of infiltrating immune cells were different between tissues and organoids, whereas there was no difference between cancer and adenoma organoids. The amounts of cancer stem cells and progenitor cells were not different between CRC and adenoma organoids, whereas a difference was noted between CRC and adenoma tissues. In conclusion, we demonstrated that organoids model only part of the adenoma–carcinoma sequence and should be used with caution after considering their limitations.


2019 ◽  
Author(s):  
Ludwig Geistlinger ◽  
Gergely Csaba ◽  
Mara Santarelli ◽  
Marcel Ramos ◽  
Lucas Schiffer ◽  
...  

AbstractBackgroundAlthough gene set enrichment analysis has become an integral part of high-throughput gene expression data analysis, the assessment of enrichment methods remains rudimentary and ad hoc. In the absence of suitable gold standards, evaluations are commonly restricted to selected data sets and biological reasoning on the relevance of resulting enriched gene sets. However, this is typically incomplete and biased towards the goals of individual investigations.ResultsWe present a general framework for standardized and structured benchmarking of enrichment methods based on defined criteria for applicability, gene set prioritization, and detection of relevant processes. This framework incorporates a curated compendium of 75 expression data sets investigating 42 different human diseases. The compendium features microarray and RNA-seq measurements, and each dataset is associated with a precompiled GO/KEGG relevance ranking for the corresponding disease under investigation. We perform a comprehensive assessment of 10 major enrichment methods on the benchmark compendium, identifying significant differences in (i) runtime and applicability to RNA-seq data, (ii) fraction of enriched gene sets depending on the type of null hypothesis tested, and (iii) recovery of the a priori defined relevance rankings. Based on these findings, we make practical recommendations on (i) how methods originally developed for microarray data can efficiently be applied to RNA-seq data, (ii) how to interpret results depending on the type of gene set test conducted, and (iii) which methods are best suited to effectively prioritize gene sets with high relevance for the phenotype investigated.ConclusionWe carried out a systematic assessment of existing enrichment methods, and identified best performing methods, but also general shortcomings in how gene set analysis is currently conducted. We provide a directly executable benchmark system for straightforward assessment of additional enrichment methods.Availabilityhttp://bioconductor.org/packages/GSEABenchmarkeR


2021 ◽  
Vol 12 ◽  
Author(s):  
Yumei Fan ◽  
Jiajie Hou ◽  
Xiaopeng Liu ◽  
Bihui Han ◽  
Yanxiu Meng ◽  
...  

Hepatocellular carcinoma (HCC) is one of the most common malignancies and ranks as the second leading cause of cancer-related mortality worldwide. Heat shock factor 2 (HSF2) is a transcription factor that plays a critical role in development, particularly corticogenesis and spermatogenesis. However, studies examining the expression and prognostic value of HSF2 and its association with tumor-infiltrating immune cells in HCC are still rare. In the present study, we found that HSF2 expression was significantly upregulated in HCC tissues compared with normal liver tissues using the TCGA, ICGC, GEO, UALCAN, HCCDB and HPA databases. High HSF2 expression was associated with shorter survival of patients with HCC. Cox regression analyses and nomogram were used to evaluate the association of HSF2 expression with the prognosis of patients with HCC. Gene Ontology (GO) analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis and gene set enrichment analysis (GSEA) revealed that HSF2 was associated with various signaling pathways, including the immune response. Notably, HSF2 expression was significantly correlated with the infiltration levels of different immune cells using the TIMER database and CIBERSORT algorithm. HSF2 expression also displayed a significant correlation with multiple immune marker sets in HCC tissues. Knockdown of HSF2 significantly inhibited the proliferation, migration, invasion and colony formation ability of HCC cells. In summary, we explored the clinical significance of HSF2 and provided a therapeutic basis for the early diagnosis, prognostic judgment, and immunotherapy of HCC.


Sign in / Sign up

Export Citation Format

Share Document