AUTHOR=Huo Junyu , Wu Liqun , Zang Yunjin TITLE=A Prognostic Model of 15 Immune-Related Gene Pairs Associated With Tumor Mutation Burden for Hepatocellular Carcinoma JOURNAL=Frontiers in Molecular Biosciences VOLUME=7 YEAR=2020 URL=https://www.frontiersin.org/journals/molecular-biosciences/articles/10.3389/fmolb.2020.581354 DOI=10.3389/fmolb.2020.581354 ISSN=2296-889X ABSTRACT=Introduction

Tumor mutation burden (TMB) is an emerging biomarker for immunotherapy of hepatocellular carcinoma (HCC), but its value for clinical application has not been fully revealed.

Materials and Methods

We used the Wilcox test to identify the differentially expressed immune-related genes (DEIRGs) in groups with high and low TMB as well as screened the immune gene pairs related to prognosis using univariate Cox regression analysis. A LASSO Cox regression prognostic model was developed by combining The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) with the International Cancer Genome Consortium (ICGC) LIRI-JP cohort, and internal (TCGA, ICGC) and external (GSE14520) validation analyses were conducted on the predictive value of the model. We also explored the relationship between the prognostic model and tumor microenvironment via the ESTIMATE algorithm and performed clinical correlation analysis by the chi-square test, revealing its underlying molecular mechanism with the help of Gene Set Enrichment Analysis (GSEA).

Results

The prognostic model consisting of 15 immune gene pairs showed high predictive value for short- and long-term survival of HCC in three independent cohorts. Based on univariate multivariate Cox regression analysis, the prognostic model could be used to independently predict the prognosis in each independent cohort. The immune score, stromal score, and estimated score values were lower in the high-risk group than in the low-risk group. As shown by the chi-square test, the prognostic model exhibited an obvious correlation with the tumor stage [American Joint Committee on Cancer tumor–node–metastasis (AJCC-TNM) (p < 0.001), Barcelona Clinic Liver Cancer (BCLC) (p = 0.003)], histopathological grade (p = 0.033), vascular invasion (p = 0.009), maximum tumor diameter (p = 0.013), and background of liver cirrhosis (p < 0.001). GSEA revealed that the high-risk group had an enrichment of many oncology features, including the cell cycle, mismatch repair, DNA replication, RNA degradation, etc.

Conclusion

Our research developed and validated a reliable prognostic model associated with TMB for HCC, which may help to further enrich the therapeutic targets of HCC.