Abstract
BACKGROUND: The important regulatory role of competitive endogenous RNAs (ceRNAs) in hepatocellular carcinoma (HCC) has been confirmed. Tumor infiltrating lymphocytes (TILs) are of great significance to tumor outcome and prognosis. This study will systematically analyze the key factors affecting the prognosis of HCC from the perspective of ceRNA and TILs. METHODS: The Cancer Genome Atlas (TCGA) database was used for transcriptome data acquisition of HCC. Through the analysis of the Weighted Gene Co-expression Network Analysis (WCGNA), the two modules for co-expression of the disease were determined, and a ceRNA network was constructed. We used Cox regression and LASSO regression analysis to screen prognostic factors and constructed a risk score model. The Gene Expression Omnibus (GEO) was used to validate the model. The Kyoto Encyclopedia of Genes and Genomes (KEGG) was used for mRNAs functional analysis. The cell composition of TILs was analyzed by the CIBERSORT algorithm, and Pearson correlation analysis was utilized to explore the correlation between TILs and prognostic factors. RESULTS: We constructed a ceRNA regulatory network composed of 67 nodes through WGCNA, including 44 DElncRNAs, 19 DEGs, and 4 DEmiRNAs. And based on the expression of 4 DEGs in this network (RRM2, LDLR, TXNIP, and KIF23), a prognostic model of HCC with good specificity and sensitivity was developed. CIBERSORT analyzed the composition of TILs in HCC tumor tissues. Correlation analysis showed that RRM2 is significantly correlated with T cells CD4 memory activated, T cells CD4 memory resting, T cells CD8, and T cells follicular helper, and TXNIP is negatively correlated with B cells memory. CONCLUSIONS: In this study, a ceRNA with prognostic value in HCC was created, and a prognostic risk model for HCC was constructed based on it. This risk score model is closely related to TILs and is expected to become a potential therapeutic target and a new predictive indicator.