scholarly journals Comprehensive analysis of immune-related prognostic genes in the tumour microenvironment of hepatocellular carcinoma

BMC Cancer ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Weike Gao ◽  
Luan Li ◽  
Xinyin Han ◽  
Siyao Liu ◽  
Chengzhen Li ◽  
...  

Abstract Background The mortality rate of hepatocellular carcinoma (HCC) remains high worldwide despite surgery and chemotherapy. Immunotherapy is a promising treatment for the rapidly expanding HCC spectrum. Therefore, it is necessary to further explore the immune-related characteristics of the tumour microenvironment (TME), which plays a vital role in tumour initiation and progression. Methods In this research, 866 immune-related differentially expressed genes (DEGs) were identified by integrating the DEGs of samples from The Cancer Genome Atlas (TCGA)-HCC dataset and the immune-related genes from databases (InnateDB; ImmPort). Afterwards, 144 candidate prognostic genes were defined through weighted gene co-expression network analysis (WGCNA). Results Seven immune-related prognostic DEGs were identified using the L1-penalized least absolute shrinkage and selection operator (LASSO) Cox proportional hazards (PH) model, and the ImmuneRiskScore model was constructed on this basis. The prognostic index of the ImmuneRiskScore model was then validated in the relevant dataset. Patients were divided into high- and low-risk groups according to the ImmuneRiskScore. Differences in the immune cell infiltration of patients with different ImmuneRiskScore values were clarified, and the correlation of immune cell infiltration with immunotherapy biomarkers was further explored. Conclusion The ImmuneRiskScore of HCC could be a prognostic marker and can reflect the immune characteristics of the TME. Furthermore, it provides a potential biomarker for predicting the response to immunotherapy in HCC patients.

2022 ◽  
Author(s):  
Yang Bu ◽  
Kejun Liu ◽  
Yiming Niu ◽  
Ji Hao ◽  
Lei Cui ◽  
...  

Abstract Background: Glucose-6-phosphate dehydrogenase (G6PD) plays an important role in the metabolic and immunological aspects of tumors. In hepatocellular carcinoma (HCC), the alteration of tumor microenvironment influences recurrence and metastasis. We extracted G6PD-related data from public databases of HCC tissues and used a bioinformatics approach to explore the correlation between G6PD expression and clinicopathological features and prognosis of immune cell infiltration in HCC.Methods: We extract G6PD expression information from TCGA and GEO databases in liver cancer tissues and normal tissues, validated by immunohistochemistry, and the correlation between G6PD expression and clinical features is analyzed, and the clinical significance of G6PD in liver cancer is assessed by Kaplan-Meier, Cox regression and prognostic line graph models. Functional enrichment analysis is performed by protein-protein interaction (PPI) network, GO/KEGG, GSEA and G6PD-associated differentially expressed genes (DEGs). TIMER and ssGSEA packages are used to assess the correlation between expression and the level of immune cell infiltration.Results: Our results show that G6PD expression is significantly upregulated in hepatocellular carcinoma tissues (P < 0.001). G6PD expression is associated with histological grade, pathological stage, T-stage, vascular infiltration and AFP level (P < 0.05); HCC patients in the low G6PD expression group had longer overall survival and better prognosis compared with the high G6PD expression group (P < 0.05). The level of G6PD expression also affects the levels of macrophages, unactivated dendritic cells, B cells, and follicular helper T cells in the tumor microenvironment.Conclusion: High expression of G6PD is a potential biomarker for poor prognosis of hepatocellular carcinoma, and G6PD may be a target for immunotherapy of HCC.


2021 ◽  
Vol 12 ◽  
Author(s):  
Zhuomao Mo ◽  
Daiyuan Liu ◽  
Dade Rong ◽  
Shijun Zhang

Background: Generally, hepatocellular carcinoma (HCC) exists in an immunosuppressive microenvironment that promotes tumor evasion. Hypoxia can impact intercellular crosstalk in the tumor microenvironment. This study aimed to explore and elucidate the underlying relationship between hypoxia and immunotherapy in patients with HCC.Methods: HCC genomic and clinicopathological datasets were obtained from The Cancer Genome Atlas (TCGA-LIHC), Gene Expression Omnibus databases (GSE14520) and International Cancer Genome Consortium (ICGC-LIRI). The TCGA-LIHC cases were divided into clusters based on single sample gene set enrichment analysis and hierarchical clustering. After identifying patients with immunosuppressive microenvironment with different hypoxic conditions, correlations between immunological characteristics and hypoxia clusters were investigated. Subsequently, a hypoxia-associated score was established by differential expression, univariable Cox regression, and lasso regression analyses. The score was verified by survival and receiver operating characteristic curve analyses. The GSE14520 cohort was used to validate the findings of immune cell infiltration and immune checkpoints expression, while the ICGC-LIRI cohort was employed to verify the hypoxia-associated score.Results: We identified hypoxic patients with immunosuppressive HCC. This cluster exhibited higher immune cell infiltration and immune checkpoint expression in the TCGA cohort, while similar significant differences were observed in the GEO cohort. The hypoxia-associated score was composed of five genes (ephrin A3, dihydropyrimidinase like 4, solute carrier family 2 member 5, stanniocalcin 2, and lysyl oxidase). In both two cohorts, survival analysis revealed significant differences between the high-risk and low-risk groups. In addition, compared to other clinical parameters, the established score had the highest predictive performance at both 3 and 5 years in two cohorts.Conclusion: This study provides further evidence of the link between hypoxic signals in patients and immunosuppression in HCC. Defining hypoxia-associated HCC subtypes may help reveal potential regulatory mechanisms between hypoxia and the immunosuppressive microenvironment, and our hypoxia-associated score could exhibit potential implications for future predictive models.


2022 ◽  
Vol 2022 ◽  
pp. 1-24
Author(s):  
Bin-Bin Da ◽  
Shuai Luo ◽  
Ming Huang ◽  
Fei Song ◽  
Rong Ding ◽  
...  

It has been demonstrated that the inflammatory response influences cancer development and can be used as a prognostic biomarker in various tumors. However, the relevance of genes associated with inflammatory responses in hepatocellular carcinoma (HCC) remains unknown. The Cancer Genome Atlas (TCGA) database was analyzed using weighted gene coexpression network analysis (WGCNA) and differential analysis to discover essential inflammatory response-related genes (IFRGs). Cox regression studies, both univariate and multivariate, were employed to develop a prognostic IFRGs signature. Additionally, Gene Set Enrichment Analysis (GSEA) was used to deduce the biological function of the IFRGs signature. Finally, we estimated immune cell infiltration using a single sample GSEA (ssGSEA) and x-cell. Our results revealed that, among the major HCC IFRGs, two (DNASE1L3 and KLKB1) were employed to create a predictive IFRG signature. The IFRG signature could correctly predict overall survival (O.S) as per Kaplan-Meier time-dependent roc curves analysis. It was also linked to pathological tumor stage and T stage and might be used as a prognostic predictor in HCC. GSEA analysis concluded that the IFRG signature might influence the immune response in HCC. Immunological cell infiltration and immune checkpoint molecule expression differed in the high-risk and low-risk groups. As a result of our findings, DNASILE may play a role in the tumor microenvironment. However, more research is necessary to confirm the role of DNASE1L3 and KLKB1.


2021 ◽  
Vol 11 ◽  
Author(s):  
Xiaoxi Shi ◽  
Yuanlin Liu ◽  
Shuai Cheng ◽  
Haidi Hu ◽  
Jian Zhang ◽  
...  

BackgroundCancer stem cells (CSCs) have been proven to influence drug resistance, recurrence, and metastasis in tumors. Our study aimed to identify stemness-related prognostic biomarkers for new therapeutic strategies in adrenocortical carcinoma.MethodsRNA-seq data and clinical characteristics were downloaded from The Cancer Genome Atlas (TCGA). The stemness indexes, mDNAsi and mRNAsi, were calculated to classify all samples into low-score and high-score groups. Two algorithms, based on the R language, ESTIMATE and single-sample Gene Set Enrichment Analysis (ssGSEA) were used to assess the immune cell infiltration states of adrenocortical carcinoma patients. Weighted Gene Co-expression Network Analysis (WGCNA) was used to find genes that were related to the stemness of cancer. By bioinformatics methods, the correlations between biomarkers capable of predicting immune checkpoint inhibitors (ICIs) responses and stemness of cancer were explored.ResultsHigh-mRNAsi predicted shorter overall survival (OS) and a higher metastatic trend in adrenocortical carcinoma (ACC) patients. Compared with the low-mRNAsi group, the high-mRNAsi group had a lower ImmuneScore and StromalScroe. Twenty-two stemness-related prognostic genes were obtained by WGCNA, which focused on the function of the cell cycle and cell mitosis. Immune cell infiltration, especially CD8+T cell, increased in the low-mRNAsi group compared with the high-mRNAsi group. Lower expression of PD-L1, CTLA-4, and TIGHT was evaluated in the high-mRNAsi group.ConclusionsACC patients with high-mRNAsi have poor prognosis and less immune cell infiltration. Combined with the finding of lower expression of CTLA-4, TIGHT, and PD-L1 in the high-mRNAsi group, we came to the conclusion that stemness index is a potential biomarker to predict the effectiveness of ICIs.


2021 ◽  
Vol 12 ◽  
Author(s):  
Shiyan Yang ◽  
Yajun Cheng ◽  
Xiaolong Wang ◽  
Ping Wei ◽  
Hui Wang ◽  
...  

Background: Globally, hepatocellular carcinoma (HCC) is the sixth most frequent malignancy with a high incidence and a poor prognosis. Immune cell infiltration (ICI) underlies both the carcinogenesis and immunogenicity of tumors. However, a comprehensive classification system based on the immune features for HCC remains unknown.Methods: The HCC dataset from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) cohorts was used in this study. The ICI patterns of 571 patients were characterized using two algorithms: the patterns were determined based on the ICI using the ConsensusClusterPlus package, and principal component analysis (PCA) established the ICI scores. Differences in the immune landscape, biological function, and somatic mutations across ICI scores were evaluated and compared, followed by a predictive efficacy evaluation of ICI scores for immunotherapy by the two algorithms and validation using an external immunotherapy cohort.Results: Based on the ICI profile of the HCC patients, three ICI patterns were identified, including three subtypes having different immunological features. Individual ICI scores were determined; the high ICI score subtype was characterized by enhanced activation of immune-related signaling pathways and a significantly high tumor mutation burden (TMB); concomitantly, diminished immunocompetence and enrichment of pathways associated with cell cycle and RNA degradation were found in the low ICI score subtype. Taken together, our results contribute to a better understanding of an active tumor and plausible reasons for its poor prognosis.Conclusion: The present study reveals that ICI scores may serve as valid prognostic biomarkers for immunotherapy in HCC.


2021 ◽  
Author(s):  
Xinyu Gu ◽  
Haibo Zhou ◽  
Qingfei Chu ◽  
Qiuxian Zheng ◽  
Jing Wang ◽  
...  

Abstract Background: 5-Methylcytosine (m5C) plays essential roles in hepatocellular carcinoma (HCC), but the association between m5C regulation and immune cell infiltration in HCC has not yet been clarified.Methods: In this study, we analysed 371 patients with HCC from The Cancer Genome Atlas (TCGA) database, and the expression of 13 m5C regulators was investigated. Additionally, gene set variation analysis (GSVA), unsupervised clustering analysis, single-sample gene set enrichment analysis (ssGSEA), correlation analysis, and immunohistochemical (IHC) staining were performed.Results: Among the 371 patients, 41 had mutations in m5C regulators, the frequency of which was 11.26%. Then, we identified three m5C modification patterns that had obvious tumour microenvironment (TME) cell infiltration characteristics. Cluster-1 had an immune rejection phenotype; Cluster-2 had an immunoinflammatory phenotype; and Cluster-3 had an immune desert phenotype. In addition, we found that DNMT1 was highly expressed in tumour tissues compared with normal tissues in a tissue microarray (TMA) and that it was positively correlated with many TME-infiltrating immune cells. High expression of the m5C regulator DNMT1 was related to a poor prognosis in patients with HCC. Furthermore, we developed three Immu-clusters that were consistent with the immune characteristics of the m5C methylation modification patterns. We also discovered differences in the levels of immune cells and expression of chemokines and cytokines among the three Immu-clusters.Conclusions: Our work revealed the association between m5C modification and immune regulators in the TME. These findings also suggest that DNMT1 has great potential as a prognostic biomarker and therapeutic target for HCC.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Qiuxian Zheng ◽  
Qin Yang ◽  
Jiaming Zhou ◽  
Xinyu Gu ◽  
Haibo Zhou ◽  
...  

Abstract Background Hepatocellular carcinoma (HCC) has a poor prognosis and has become the sixth most common malignancy worldwide due to its high incidence. Advanced approaches to therapy, including immunotherapeutic strategies, have played crucial roles in decreasing recurrence rates and improving clinical outcomes. The HCC microenvironment is important for both tumour carcinogenesis and immunogenicity, but a classification system based on immune signatures has not yet been comprehensively described. Methods HCC datasets from The Cancer Genome Atlas (TCGA), the Gene Expression Omnibus (GEO), and the International Cancer Genome Consortium (ICGC) were used in this study. Gene set enrichment analysis (GSEA) and the ConsensusClusterPlus algorithm were used for clustering assessments. We scored immune cell infiltration and used linear discriminant analysis (LDA) to improve HCC classification accuracy. Pearson's correlation analyses were performed to assess relationships between immune signature indices and immunotherapies. In addition, weighted gene co-expression network analysis (WGCNA) was applied to identify candidate modules closely associated with immune signature indices. Results Based on 152 immune signatures from HCC samples, we identified four distinct immune subtypes (IS1, IS2, IS3, and IS4). Subtypes IS1 and IS4 had more favourable prognoses than subtypes IS2 and IS3. These four subtypes also had different immune system characteristics. The IS1 subtype had the highest scores for IFNγ, cytolysis, angiogenesis, and immune cell infiltration among all subtypes. We also identified 11 potential genes, namely, TSPAN15, TSPO, METTL9, CD276, TP53I11, SPINT1, TSPO, TRABD2B, WARS2, C9ORF116, and LBH, that may represent potential immunological biomarkers for HCC. Furthermore, real-time PCR revealed that SPINT1, CD276, TSPO, TSPAN15, METTL9, and WARS2 expression was increased in HCC cells. Conclusions The present gene-based immune signature classification and indexing may provide novel perspectives for both HCC immunotherapy management and prognosis prediction.


2021 ◽  
Vol 41 (2) ◽  
Author(s):  
Junhui Chen ◽  
Jie Yang ◽  
Qingchun Xu ◽  
Zhenyu Wang ◽  
Jun Wu ◽  
...  

Abstract Liver hepatocellular carcinoma (LIHC) is one of the most frequently occurring primary malignant liver tumors and seriously harms people’s health in the world. Methylenetetrahydrofolate dehydrogenase 1-like (MTHFD1L) has been shown to be associated with colon cancer cell proliferation, colony formation and invasion. In the present study, a total of 370 LIHC and 51 normal samples data were downloaded from The Cancer Genome Atlas (TCGA) database. Bioinformatics and immunohistochemistry (IHC) analysis showed that MTHFD1L is highly expressed in liver tumors. Correlation analysis suggested the differences of vital status between high- and low-expression MTHFD1L groups of LIHC. Univariate and multivariate Cox proportional hazards regression were performed to identify the relationship between clinical characteristics and overall survival (OS). In addition, to explore whether MTHFD1L has an effect on the immune infiltration of LIHC. The correlation between MTHFD1L expression and 24 immune cells were analyzed by ImmuneCellAI database. Furthermore, we combined three databases CIBERSORT, TIMER and ImmuneCellAI to do a comprehensive validation and determined that dendritic cells (DCs) resting, macrophage M0 and macrophage M2 closely related to the expression of MTHFD1L. The results showed that MTHFD1L was a potential prognostic biomarker for LIHC, and could help to elucidate that how the immune microenvironment promotes liver cancer development.


Author(s):  
Xinyu Gu ◽  
Haibo Zhou ◽  
Qingfei Chu ◽  
Qiuxian Zheng ◽  
Jing Wang ◽  
...  

Background: 5-Methylcytosine (m5C) plays essential roles in hepatocellular carcinoma (HCC), but the association between m5C regulation and immune cell infiltration in HCC has not yet been clarified.Methods: In this study, we analysed 371 patients with HCC from The Cancer Genome Atlas (TCGA) database, and the expression of 13 m5C regulators was investigated. Additionally, gene set variation analysis (GSVA), unsupervised clustering analysis, single-sample gene set enrichment analysis (ssGSEA), correlation analysis, and immunohistochemical (IHC) staining were performed.Results: Among the 371 patients, 41 had mutations in m5C regulators, the frequency of which was 11.26%. Compared with normal hepatic tissues, the expression of m5C regulators with copy number variations (CNVs) expansion was significantly higher than that in HCC tissues. Then, we identified three m5C modification patterns that had obvious tumour microenvironment (TME) cell infiltration characteristics. The prognostic analysis of the three major m5C modification subtypes showed that Cluster-2 had a clear survival advantage over the others. In addition, we found that DNMT1 was highly expressed in tumour tissues compared with normal tissues in a tissue microarray (TMA) and that it was positively correlated with many TME-infiltrating immune cells. High expression of the m5C regulator DNMT1 was related to a poor prognosis in patients with HCC. Furthermore, we developed three distinct Immu-clusters. Importantly, mRNAs related to the transcription of growth factor β (TGF-β)/EMT pathway were significantly up-regulated in Immu-cluster 2, indicating that this cluster is considered to be the immune rejection phenotype. Immu-cluster 3 showed elevated expression of mRNAs related to immune checkpoint genes.Conclusion: Our work revealed the association between m5C modification and immune regulators in the TME. These findings also suggest that DNMT1 has great potential as a prognostic biomarker and therapeutic target for HCC.


Author(s):  
Lu Yuan ◽  
Xixi Wu ◽  
Longshan Zhang ◽  
Mi Yang ◽  
Xiaoqing Wang ◽  
...  

AbstractPulmonary surfactant protein A1 (SFTPA1) is a member of the C-type lectin subfamily that plays a critical role in maintaining lung tissue homeostasis and the innate immune response. SFTPA1 disruption can cause several acute or chronic lung diseases, including lung cancer. However, little research has been performed to associate SFTPA1 with immune cell infiltration and the response to immunotherapy in lung cancer. The findings of our study describe the SFTPA1 expression profile in multiple databases and was validated in BALB/c mice, human tumor tissues, and paired normal tissues using an immunohistochemistry assay. High SFTPA1 mRNA expression was associated with a favorable prognosis through a survival analysis in lung adenocarcinoma (LUAD) samples from TCGA. Further GeneOntology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses showed that SFTPA1 was involved in the toll-like receptor signaling pathway. An immune infiltration analysis clarified that high SFTPA1 expression was associated with an increased number of M1 macrophages, CD8+ T cells, memory activated CD4+ T cells, regulatory T cells, as well as a reduced number of M2 macrophages. Our clinical data suggest that SFTPA1 may serve as a biomarker for predicting a favorable response to immunotherapy for patients with LUAD. Collectively, our study extends the expression profile and potential regulatory pathways of SFTPA1 and may provide a potential biomarker for establishing novel preventive and therapeutic strategies for lung adenocarcinoma.


Sign in / Sign up

Export Citation Format

Share Document