Identification of key snoRNAs serves as biomarkers for hepatocellular carcinoma by bioinformatics methods.

Background: Hepatocellular carcinoma (HCC) is a common malignancy with high mortality and poor prognosis due to a lack of predictive markers. However, research on small nuclear RNAs (snoRNAs) in HCC were very little. This study aimed to identify a potential diagnostic and prognostic snoRNA signature for HCC.

Methods: HCC datasets from the cancer genome atlas (TCGA) and international cancer genome consortium (ICGC) cohorts were used. Differentially expressed snoRNA (DEs) were identified using the limma package. Based on the DEs, diagnostic and prognostic models were established by the least absolute shrinkage and selection operator (LASSO) regression and COX analysis, and Kaplan-Meier (K-M) survival analysis and receiver operating characteristic (ROC) curve analysis were conducted to evaluate the efficiency of signatures. Gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) were used to analyze the risk score and further explore the potential correlation between the risk groups and tumor immune status in TCGA. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to determine the functions of key snoRNAs.

Results: We constructed a 6-snoRNAs signature which could classify patients into high- or low-risk groups and found that patients in the high-risk group had a worse prognosis than those in the low-risk group and were significantly involved in p53 processes. Tumor immune status analysis revealed that CTLA4 and PDCD1 (PD1) were highly expressed in the high-risk group, which responded to PD1 inhibitor therapy. Additionally, a 25-snoRNAs diagnostic signature was constructed with an area under the curve (AUC) of 0.933 for distinguishing HCCs from normal controls. Finally, 3 key snoRNAs (SNORA11, SNORD124, and SNORD46) were identified with both diagnostic and prognostic efficacy, some of which were closely related to the spliceosome and Notch signaling pathways.

Conclusions: Our study identified 6 snoRNAs that may serve as novel prognostic models and 3 key snoRNAs with both diagnostic and prognostic efficacy for HCC.