Comprehensive Analysis of the Relationship Between Competitive Endogenous RNA (cerna) Networks and Tumor Infiltrating-Cells in Hepatocellular Carcinoma.

Jun Zhu,Liang Wang,Yifan Zhou,Jun Hao,Shuai Wang,Lei Liu,Jipeng Li

DOI: https://doi.org/10.21037/jgo-20-555

2020-01-01

Journal of Gastrointestinal Oncology

Abstract:BACKGROUND:The innovation of immune checkpoint blockade (ICB) represents a promising shift in the treatment of advanced hepatocellular carcinoma (HCC). However, response to ICB has varied largely due to the high tumor heterogeneity and complex tumor microenvironment (TME). The competitive endogenous RNA (ceRNA) network also plays an important role in tumor occurrence and progression, but its relation with tumor-infiltrating immune cells (TICs) remains largely unexplored in HCC. The overriding objective of our study was thus to construct a prognosis-related risk model and to further evaluate the relationship between ceRNA networks and TICs.METHODS:Differentially expressed gene (DEG) analysis was performed to identify the differentially expressed RNAs. Lasso and multivariable Cox regression analyses were used to construct risk models, which were assessed by the area under the receiver operating characteristic curve (AUC of ROC) and Kaplan-Meier (K-M) curves. Then, a single-sample gene set enrichment analysis (ssGSEA) algorithm was adopted to dissect the TICs in HCC samples. Nomograms were constructed and calibration curves were used to verify the discrimination and accuracy of the nomograms. Finally, integration analysis was performed to validate the correlation of ceRNA and TICs.RESULTS:In the study, 7 differentially expressed RNAs [5 messenger RNA s (mRNAs) and 2 micro RNAs (miRNAs)] were incorporated to construct a ceRNA risk model. The AUC of the 1-, 3-, and 5-year overall survival (OS) were 0.784, 0.685, and 0.691 respectively. Likewise, 7 types TICs were in the TICs signature model and the AUC of the 1-, 3-, and 5-year OS were 0.706, 0.731, and 0.721 respectively. The integration analysis showed that 7 pairs of mRNA-TICs and 1 pair of miRNA-TICs had a close relation (all correlation coefficients >0.2, P<0.001).CONCLUSIONS:Through constructing two risk models based on ceRNA network and TICs, we identified the hub RNAs and key TICs in the progression and prognosis of HCC, and further explored the relationship between ceRNA and TME. Importantly, targeting these hub RNAs may facilitate the remodeling of the TME and be a potential therapeutic alternative to enhancing the response to ICB, thus improving the prognosis of HCC patients.