Comprehensive analysis of nine m7G-related lncRNAs as prognosis factors in tumor immune microenvironment of hepatocellular carcinoma and experimental validation.

Abstract

Background: Hepatocellular carcinoma (HCC) remains the most prevalent gastrointestinal malignancy worldwide, with robust drug resistance to therapy. N7-methylguanosine (m7G) mRNA modification has been significantly related to massive human diseases. Considering the effect of m7G-modified long non-coding RNAs (lncRNAs) in HCC progression is unknown, the study aims at investigating a prognostic signature to improve clinical outcomes for patients with HCC. Methods: Two independent databases (TCGA and ICGC) were used to analyze RNAseq data of HCC patients. First, co-expression analysis was applied to obtain the m7G-related lncRNAs. Moreover, consensus clustering analysis was employed to divide HCC patients into clusters. Then, using least absolute shrinkage and selection operator-Cox regression analysis, the m7G-related lncRNA prognostic signature (m7G-LPS) was first tested in the training set and then confirmed in both the testing and ICGC sets. The expression levels of the nine lncRNAs were further confirmed via real-time PCR in cell lines, principal component analysis, and receiver operating characteristic curve. The m7G-LPS could divide HCC patients into two different risk groups with the optimal risk score. Then, Kaplan–Meier curves, tumor mutation burden (TMB), therapeutic effects of chemotherapy agents, and expressions of immune checkpoints were performed to further enhance the availability of immunotherapeutic treatments for HCC patients. Results: A total of 1465 lncRNAs associated with the m7G genes were finally selected from the TCGA database, and through the univariate Cox regression, the expression levels of 22 m7G-related lncRNAs were concerning HCC patients’ overall survival (OS). Then, the whole patients were grouped into two subgroups, and the OS in Cluster 1 was longer than that of patients in Cluster 2. Furthermore, nine prognostic m7G-related lncRNAs were identified to conduct the m7G-LPS, which were further verified. A prognostic nomogram combined age, gender, HCC grade, stage, and m7G-LPS showed strong reliability and accuracy in predicting OS in HCC patients. Finally, immune checkpoint expression, TMB, and several chemotherapy agents were remarkably associated with risk scores. More importantly, the OS of the TMB-high patients was the worst among the four groups. Conclusion: The prognostic model we established was validated by abundant algorithms, which provided a new perspective on HCC tumorigenesis and thus improved individualized treatments for patients.