A Pyroptosis-Based Prognostic Model for Immune Microenvironment Estimation of Hepatocellular Carcinoma

Background. Hepatocellular carcinoma (HCC), an aggressive malignant tumor, has a high incidence and unfavorable prognosis. Recently, the synergistic effect of pyroptosis in antitumor therapy and regulation of tumor immune microenvironment has made it possible to become a novel therapeutic method, but its potential mechanism still needs further exploration. Methods. Differentially expressed genes with prognostic value in Liver Hepatocellular Carcinoma Project of The Cancer Genome Atlas (TCGA-LIHC) cohort were screened and incorporated into the risk signature by Cox proportional hazards regression model and least absolute shrinkage and selection operator. Kaplan-Meier (KM) curves and receiver operating characteristic (ROC) curves were applied to conduct survival comparisons and estimate prediction ability. The dataset of Liver Cancer-RIKEN, Japan Project from International Cancer Genome Consortium (ICGC-LIRI-JP) cohort was used to verify the reliability of the signature. Correlation analysis between clinicopathological characteristics, immune infiltration, drug sensitivities, and risk scores was conducted. Functional annotation analyses were performed for the genes differentially expressed between high-risk and low-risk groups. Results. A risk signature consisting of 6 pyroptosis-related genes in HCC was developed and validated. KM curves and ROC curves revealed its considerable predictive accuracy. Higher risk scores meant more advanced grade, higher alpha-fetoprotein level, and stronger invasive ability. Overexpressed genes in high-risk population were more enriched in the immune-associated pathways, and these patients might be more sensitive to immune checkpoint inhibitors instead of Sorafenib. Intriguingly, 6 identified genes were promising to be prognostic biomarkers and therapeutic targets of HCC. Conclusions. The signature may have crucial clinical significance in predicting survival prognosis, immune infiltration, and drug efficacy based on pyroptosis-related genes.

Download Full-text

Prognostic Modeling of Patients with Metastatic Melanoma Based on Tumor Immune Microenvironment Characteristics

10.21203/rs.3.rs-730205/v1 ◽

2021 ◽

Author(s):

Jing Liu ◽

Ting Ye ◽

Xue fang Zhang ◽

Yong jian Dong ◽

Wen feng Zhang ◽

...

Keyword(s):

High Risk ◽

Metastatic Melanoma ◽

Prognostic Model ◽

Cox Regression ◽

Risk Group ◽

Low Risk ◽

Risk Scores ◽

Immune Microenvironment ◽

Immune Infiltration ◽

Tumor Immune Microenvironment

Abstract Most of the malignant melanomas are already in the middle and advanced stages when they are diagnosed, which is often accompanied by the metastasis and spread of other organs.Besides, the prognosis of patients is bleak. The characteristics of the local immune microenvironment in metastatic melanoma have important implications for both tumor progression and tumor treatment. In this study, data on patients with metastatic melanoma from the TCGA and GEO datasets were selected for immune, stromal, and estimate scores, and overlapping differentially expressed genes (DEGs) were screened. A nine-IRGs prognostic model (ALOX5AP, ARHGAP15, CCL8, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22) was established by univariate COX regression, LASSO and multivariate COX regression. Receiver operating characteristic (ROC) curves were used to test the predictive accuracy of the model. Immune infiltration was analyzed by using CIBERSORT, Xcell and ssGSEA in high-risk and low-risk groups. The immune infiltration of the high-risk group was significantly lower than that of the low-risk group. Immune checkpoint analysis revealed that the expression of PDCD1, CTLA4, TIGIT, CD274, HAVR2 and LAG3 were significantly different in groups with different levels of risk scores. WGCNA analysis found that the yellow-green module contained seven genes (ALOX5AP, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22) from the nine-IRG prognostic model, of which the yellow-green module had the highest correlation with risk scores. The results of GO and KEGG suggested that the genes in the yellow-green module were mainly enriched in immune-related biological processes. Finally, we analyzed the prognostic ability and expression characteristics of ALOX5AP, ARHGAP15, CCL8, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22 in metastatic melanoma. Overall, a prognostic model for metastatic melanoma based on the characteristics of the tumor immune microenvironment was established, which was helpful for further studies.It could function well in helping people to understand the characteristics of the immune microenvironment in metastatic melanoma and to find possible therapeutic targets.

Download Full-text

Identification of Aging-Related Genes Associated With Clinical and Prognostic Features of Hepatocellular Carcinoma

Frontiers in Genetics ◽

10.3389/fgene.2021.661988 ◽

2021 ◽

Vol 12 ◽

Author(s):

Xingte Chen ◽

Lei Wang ◽

Liang Hong ◽

Zhixiong Su ◽

Xiaohong Zhong ◽

...

Keyword(s):

Hepatocellular Carcinoma ◽

Regression Analysis ◽

High Risk ◽

Risk Score ◽

Cancer Genome ◽

Gene Set Enrichment Analysis ◽

Differentially Expressed ◽

Survival Prediction ◽

Lasso Regression ◽

Prognostic Features

Background: Aging is a well-studied concept, but no studies have comprehensively analyzed the association between aging-related genes (AGs) and hepatocellular carcinoma (HCC) prognosis.Methods: Gene candidates were selected from differentially expressed genes and prognostic genes in The Cancer Genome Atlas (TCGA) database. A gene risk score for overall survival prediction was established using the least absolute shrinkage and selection operator (LASSO) regression analysis, and this was validated using data from the International Cancer Genome Consortium (ICGC) database. Functional analysis was conducted using gene ontology enrichment, Kyoto Encyclopedia of Genes and Genomes analysis, gene set enrichment analysis, and immune microenvironment and tumor stemness analyses.Results: Initially, 72 AGs from the TCGA database were screened as differentially expressed between normal and tumor tissues and as genes associated with HCC prognosis. Then, seven AGs (POLA1, CDK1, SOCS2, HDAC1, MAPT, RAE1, and EEF1E1) were identified using the LASSO regression analysis. The seven AGs were used to develop a risk score in the training set, and the risk was validated to have a significant prognostic value in the ICGC set (p < 0.05). Patients with high risk scores had lower tumor differentiation, higher stage, and worse prognosis (all p < 0.05). Multivariate Cox regression analyses also confirmed that the risk score was an independent prognostic factor for HCC in both the TCGA and ICGC sets (all p < 0.05). Further analysis showed that a high risk score was correlated with the downregulation of metabolism and tumor immunity.Conclusion: The risk score predicts HCC prognosis and could thus be used as a biomarker not only for predicting HCC prognosis but also for deciding on treatment.

Download Full-text

An m6A-Related Prognostic Biomarker Associated With the Hepatocellular Carcinoma Immune Microenvironment

Frontiers in Pharmacology ◽

10.3389/fphar.2021.707930 ◽

2021 ◽

Vol 12 ◽

Author(s):

Yingxi Du ◽

Yarui Ma ◽

Qing Zhu ◽

Tongzheng Liu ◽

Yuchen Jiao ◽

...

Keyword(s):

Hepatocellular Carcinoma ◽

Immune Checkpoint ◽

Immune Cell ◽

Cancer Genome ◽

The Cancer Genome Atlas ◽

Risk Scores ◽

Survival Prediction ◽

Consensus Clustering ◽

Immune Microenvironment ◽

Cluster 2

Background: N6-methyladenosine (m6A) is related to the progression of multiple cancers. However, the underlying influences of m6A-associated genes on the tumor immune microenvironment in hepatocellular carcinoma (HCC) remain poorly understood. Therefore, we sought to construct a survival prediction model using m6A-associated genes to clarify the molecular and immune characteristics of HCC.Methods: HCC case data were downloaded from The Cancer Genome Atlas (TCGA). Then, by applying consensus clustering, we identified two distinct HCC clusters. Next, four m6A-related genes were identified to construct a prognostic model, which we validated with Gene Expression Omnibus (GEO) and International Cancer Genome Consortium (ICGC) datasets. Additionally, the molecular and immune characteristics in different subgroups were analyzed.Results: m6A RNA methylation regulators were differentially expressed between HCC and normal samples and linked with immune checkpoint expression. Using consensus clustering, we divided HCC samples into two subtypes with distinct clinical features. Cluster 2 was associated with unfavorable prognosis, higher immune checkpoint expression and immune cell infiltration levels. In addition, the immune and carcinogenic signaling pathways were enriched in cluster 2. Furthermore, we constructed a risk model using four m6A-associated genes. Patients with different risk scores had distinct survival times, expression levels of immunotherapy biomarkers, TP53 mutation rates, and sensitivities to chemotherapy and targeted therapy. Similarly, the model exhibited an identical impact on overall survival in the validation cohorts.Conclusion: The constructed m6A-based signature may be promising as a biomarker for prognostics and to distinguish immune characteristics in HCC.

Download Full-text

Prognostic modeling of patients with metastatic melanoma based on tumor immune microenvironment characteristics

Mathematical Biosciences and Engineering ◽

10.3934/mbe.2022067 ◽

2021 ◽

Vol 19 (2) ◽

pp. 1448-1470

Author(s):

Jing Liu ◽

◽

Xuefang Zhang ◽

Ting Ye ◽

Yongjian Dong ◽

...

Keyword(s):

High Risk ◽

Metastatic Melanoma ◽

Prognostic Model ◽

Cox Regression ◽

Risk Group ◽

Low Risk ◽

Risk Scores ◽

Immune Microenvironment ◽

Immune Infiltration ◽

Tumor Immune Microenvironment

<abstract> <p>Most of the malignant melanomas are already in the middle and advanced stages when they are diagnosed, which is often accompanied by the metastasis and spread of other organs. Besides, the prognosis of patients is bleak. The characteristics of the local immune microenvironment in metastatic melanoma have important implications for both tumor progression and tumor treatment. In this study, data on patients with metastatic melanoma from the TCGA and GEO datasets were selected for immune, stromal, and estimate scores, and overlapping differentially expressed genes were screened. A nine-IRGs prognostic model (ALOX5AP, ARHGAP15, CCL8, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22) was established by univariate COX regression, LASSO and multivariate COX regression. Receiver operating characteristic curves were used to test the predictive accuracy of the model. Immune infiltration was analyzed by using CIBERSORT and Xcell in high-risk and low-risk groups. The immune infiltration of the high-risk group was significantly lower than that of the low-risk group. Immune checkpoint analysis revealed that the expression of PDCD1, CTLA4, TIGIT, CD274, HAVR2 and LAG3 demonstrated the visible difference in groups with different levels of risk scores. WGCNA analysis found that the yellow-green module contained seven genes from the nine-IRG prognostic model, and the yellow-green module had the highest correlation with risk scores. The results of GO and KEGG suggested that the genes in the yellow-green module were mainly enriched in immune-related biological processes. Finally, the expression characteristics of ALOX5AP, ARHGAP15, CCL8, FCER1G, GBP4, HCK, MMP9, RARRES2 and TRIM22 were analyzed between metastatic melanoma and normal samples. Overall, a prognostic model for metastatic melanoma based on the tumor immune microenvironment characteristics was established, which left plenty of space for further studies. It could function well in helping people to understand characteristics of the immune microenvironment in metastatic melanoma.</p> </abstract>

Download Full-text

Differentially expressed liver exosome-related genes as a candidate prognostic biomarker for hepatocellular carcinoma

10.21203/rs.3.rs-105334/v1 ◽

2020 ◽

Author(s):

Bangyou Zuo ◽

Haitao Zhao ◽

Jin Bian ◽

Junyun Long ◽

Xu Yang ◽

...

Keyword(s):

Hepatocellular Carcinoma ◽

High Risk ◽

Prognostic Model ◽

Cox Regression ◽

Risk Group ◽

Cell Communication ◽

High Risk Group ◽

Cancer Genome ◽

The Cancer Genome Atlas ◽

Differentially Expressed

Abstract Background The function of exosome includes cell-to-cell communication, neovascularization, and metastasis of cancer cell and drug resistance, which plays an important part in the occurrence and progression of hepatocellular carcinoma (HCC). Because the mechanism in this area is less studied, our goal is to identify exosome-related genes in HCC, establish a reliable prognostic model for liver cancer patients, and explore its underlying mechanisms. Methods The exoRbase database and The Cancer Genome Atlas (TCGA) database were used to analyze differentially expressed genes (DEGs). Cox regression and LASSO analysis were applied to determine DEGs closely related to overall survival (OS). Then the exosome-related prognostic model was constructed in TCGA and validated in the database of International Cancer Genome Consortium (ICGC). Nomogram graph was performed to predict the survival. CIBERSORT was used to estimate the score of different type of immune cells. DEGs related to immunotherapy are used to predict the effect of immunotherapy. Results 48 exosome-related DEGs were obtained and five genes (XPO1, IFI30, FBXO16, CALM1, MORC3) among them were selected to construct predictive model. Then we divided the HCC patients into low-risk and high-risk groups by the best cut-off value according to the X-tile software. The high-risk related to exosome were significantly associated with a poor prognosis. Moreover, the features related to exosome could positively regulate immune response. At the same time, the proportion of T cell regulatory factors (Tregs) and macrophages M2 is higher in the high-risk group, and high-risk group exhibited higher expression of immune checkpoint molecular including PD-L1, PD-L2, TIGIT, and IDO1. Conclusions Overall, our research showed that markers related to exosomes were potential biomarkers for the prognosis of HCC, providing an immunological perspective for the development of precision treatment.

Download Full-text

Exploration of Tumor Mutation Burden Combined with Immune Infiltrates in the Prognosis of Lung Adenocarcinoma

10.21203/rs.3.rs-219997/v1 ◽

2021 ◽

Author(s):

Zhenyu Zhao ◽

Boxue He ◽

Qidong Cai ◽

Pengfei Zhang ◽

Xiong Peng ◽

...

Keyword(s):

High Risk ◽

Lung Adenocarcinoma ◽

Clinical Application ◽

Risk Group ◽

High Risk Group ◽

Immune Microenvironment ◽

Tumor Mutation Burden ◽

Immune Signature ◽

Immune Infiltration ◽

Mutation Burden

Abstract Background: Lung adenocarcinoma (LUAD) accounts for a majority of cancer-related deaths worldwide annually. A recent study shows that immunotherapy is an effective method of LUAD treatment, and tumor mutation burden (TMB) was associated with the immune microenvironment and affected the immunotherapy. Exploration of the gene signature associated with tumor mutation burden and immune infiltrates in predicting prognosis in lung adenocarcinoma in this study, we explored the correlation of TMB with immune infiltration and prognosis in LUAD.Materials and Methods: In this study, we firstly got mutation data and LUAD RNA-Seq data of the LUAD from The Cancer Genome Atlas (TCGA), and according to the TMB we divided the patients into high/low-TMB levels groups. The gene ontology (GO) pathway enrichment analysis and KOBAS-Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analysis were utilized to explore the molecular function of the differentially expressed genes (DEGs) between the two groups. The function enrichment analyses of DEGs were related to the immune pathways. Then, the ESTIMATE algorithm, CIBERSORT, and ssGSEA analysis were utilized to identify the relationship between TMB subgroups and immune infiltration. According to the results, Venn analysis was utilized to select the immune-related genes in DEGs. Univariate and Lasso Cox proportional hazards regression analyses were performed to construct the signature which positively associated with the immune infiltration and affected the survival. Finally, we verified the correlation between the signature and immune infiltration. Result: The exploration of the immune infiltration suggested that high-TMB subgroups positively associated with the high level of immune infiltration in LUAD patients. According to the TMB-related immune signature, the patients were divided into High/Low-risk groups, and the high-risk group was positively associated with poor prognostic. The results of the PCA analysis confirmed the validity of the signature. We also verified the effectiveness of the signature in GSE30219 and GSE72094 datasets. The ROC curves and C-index suggested the good clinical application of the TMB-related immune signature in LUAD prognosis. Another result suggested that the patients of the high-risk group were positively associated with higher TMB levels, PD-L1expression, and immune infiltration levels.Conclusion: In conclusion, our signature provides potential biomarkers for studying aspects of the TMB in LUAD such as TMB affected immune microenvironment and prognosis. This signature may provide some biomarkers which could improve the biomarkers of PD-L1 immunotherapy response and were inverted for the clinical application of the TMB in LUAD. LUAD male patients with higher TMB-levels and risk scores may benefit from immunotherapy. The high-risk patients along with higher PD-L1 expression of the signature may suitable for immunotherapy and improve their survival by detecting the TMB of LUAD.

Download Full-text

Identification of Biomarkers Related to Immune Cell Infiltration in Hepatocellular Carcinoma Using Gene Co-Expression Network

Pathology & Oncology Research ◽

10.3389/pore.2021.601693 ◽

2021 ◽

Vol 27 ◽

Author(s):

Wanbang Zhou ◽

Yiyang Chen ◽

Ruixing Luo ◽

Zifan Li ◽

Guanwei Jiang ◽

...

Keyword(s):

Hepatocellular Carcinoma ◽

Network Analysis ◽

Tumor Progression ◽

Protein Interactions ◽

Immune Cells ◽

Immune Cell ◽

Immune Microenvironment ◽

Hub Genes ◽

Expression Levels ◽

Immune Infiltration

Hepatocellular carcinoma (HCC) is a common cancer with poor prognosis. Due to the lack of effective biomarkers and its complex immune microenvironment, the effects of current HCC therapies are not ideal. In this study, we used the GSE57957 microarray data from Gene Expression Omnibus database to construct a co-expression network. The weighted gene co-expression network analysis and CIBERSORT algorithm, which quantifies cellular composition of immune cells, were used to identify modules related to immune cells. Four hub genes (EFTUD2, GAPDH, NOP56, PA2G4) were identified by co-expression network and protein-protein interactions network analysis. We examined these genes in TCGA database, and found that the four hub genes were highly expressed in tumor tissues in multiple HCC groups, and the expression levels were significantly correlated with patient survival time, pathological stage and tumor progression. On the other hand, methylation analysis showed that the up-regulation of EFTUD2, GAPDH, NOP56 might be due to the hypomethylation status of their promoters. Next, we investigated the correlations between the expression levels of four hub genes and tumor immune infiltration using Tumor Immune Estimation Resource (TIMER). Gene set variation analysis suggested that the four hub genes were associated with numerous pathways that affect tumor progression or immune microenvironment. Overall, our results showed that the four hub genes were closely related to tumor prognosis, and may serve as targets for treatment and diagnosis of HCC. In addition, the associations between these genes and immune infiltration enhanced our understanding of tumor immune environment and provided new directions for the development of drugs and the monitoring of tumor immune status.

Download Full-text

A Starvation-Based 9-mRNA Signature Correlates With Prognosis in Patients With Hepatocellular Carcinoma

Frontiers in Oncology ◽

10.3389/fonc.2021.716757 ◽

2021 ◽

Vol 11 ◽

Author(s):

Dengliang Lei ◽

Yue Chen ◽

Yang Zhou ◽

Gangli Hu ◽

Fang Luo

Keyword(s):

Hepatocellular Carcinoma ◽

High Risk ◽

Liver Cancer ◽

Cancer Patients ◽

Enrichment Analysis ◽

Clinical Information ◽

Cancer Genome ◽

Gene Set Enrichment Analysis ◽

Biological Processes ◽

Gene Set Enrichment

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.

Download Full-text