Abstract Background: Hepatocellular carcinoma (HCC) is one of the most common cancers worldwide with a high recurrence rate. Although immune-checkpoint treatment provides survival benefits for patients with HCC, the mortality of HCC is still high. Therefore, it is necessary to understand the molecular mechanism, immune microenvironment and related prognostic factors of HCC. Methods: HCC patients with both RNAseq and miRNA expression data and their clinical information were downloaded from the NCI's Genomic Data Commons. Gene expression data were analyzed to identify differentially expressed genes (DEGs). The results were further used to construct the competing endogenous RNA (ceRNA) networks related to HCC. Meanwhile, the proportions of 22 immune cell types in HCC were estimated by using CIBERSORT (https://cibersort.stanford.edu). Subsequently, two nomograms based on the expression value of ceRNAs and the percentage of immune cells were constructed to predict the prognosis of HCC patients, respectively. The calibration and time series ROC curves were used to evaluate these models. Finally, the correlation between ceRNAs and immune cells were explored to reveal the underlying molecular mechanism. Results: Gene expression data and clinical information of 367 primary HCC and 50 adjacent normal samples were downloaded from the TCGA database. These data were used to identify DEGs. A total of 12 lncRNAs, 71 mRNAs and 39 miRNAs were identified to be in ceRNA networks according to starBase v2.0 database. The results showed that stages of cancer, the expression values of CBX2, MAFG, DNM3, BUB1, ARL5B and other 23 genes were important predictors. The AUCs of the constructed nomogram in predicting the 3- and 5-year OS were 0.79 and 0.81 respectively. Stages of cancer, proportions of resting memory CD4+ T cell, M0 macrophages and resting mast cells constructed another nomogram with AUC of 0.79 and 0.81 in predicting the 3- and 5-year OS. CIBERSORT results showed that the proportions of regulatory T cells, CD8+ T cells, follicular helper T cells and M0 macrophages were slightly higher while the proportions of M2 macrophages, monocytes and activated mast cells were slightly lower in HCC than in adjacent normal samples. Integration analysis of CIBERSORT results, ceRNAs and gene set enrichment analysis of immunologic signatures identified BUB1, E2F8, H2AFZ, JPT1, LPCAT1, MAFG and PSMC3IP to be positively related with regulatory T cells, follicular helper T cells, activated memory CD4+ T cells and M0 macrophages. This is consistent with current knowledge that immune activity is suppressed in certain types of cancers. However, the underlying mechanism warrant further exploration. Conclusions: In this study, two nomograms were constructed to predict the prognosis of HCC patients, which might be useful in predicting the prognosis of HCC patients. The results suggest that the immune suppress effect of HCC may be related to BUB1, E2F8, H2AFZ, JPT1, LPCAT1, MAFG and PSMC3IP genes. These genes may serve as potential therapeutic targets. Acknowledgments: This work was supported by the National Science and Technology Major Project of China (2018ZX10302205). Citation Format: Jianxiang Shi, Chi Cui, Yaru Duan, Hua Ye, Peng Wang, Jitian Li, Shuiling Jin, Liping Dai, Jianying Zhang. Construction of ce-RNA networks and comprehensive analysis with tumor infiltrating immune cells in hepatocellular carcinoma [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 35.
Read full abstract