Construction of oxidative phosphorylation-related prognostic risk score model in uveal melanoma

Uveal melanoma (UVM) is a malignant intraocular tumor in adults. Targeting genes related to oxidative phosphorylation (OXPHOS) may play a role in anti-tumor therapy. However, the clinical significance of oxidative phosphorylation in UVM is unclear.

The 134 OXPHOS-related genes were obtained from the KEGG pathway, the TCGA UVM dataset contained 80 samples, served as the training set, while GSE22138 and GSE39717 was used as the validation set. LASSO regression was carried out to identify OXPHOS-related prognostic genes. The coefficients obtained from Cox multivariate regression analysis were used to calculate a risk score, which facilitated the construction of a prognostic model. Kaplan-Meier survival analysis, logrank test and ROC curve using the time "timeROC" package were conducted. The immune cell frequency in low- and high-risk group was analyzed through Cibersort tool. The specific genomic alterations were analyzed by "maftools" R package. The differential expressed genes between low- or high-risk group were analyzed and performed Gene Ontology (GO) and GSEA. Finally, we verified the function of CYC1 in UVM by gene silencing in vitro.

A total of 9 OXPHOS-related prognostic genes were identified, including NDUFB1, NDUFB8, ATP12A, NDUFA3, CYC1, COX6B1, ATP6V1G2, ATP4B and NDUFB4. The UVM prognostic risk model was constructed based on the 9 OXPHOS-related prognostic genes. The prognosis of patients in the high-risk group was poorer than low-risk group. Besides, the ROC curve demonstrated that the area under the curve of the model for predicting the 1 to 5-year survival rate of UVM patients were all more than 0.88. External validation in GSE22138 and GSE39717 dataset revealed that these 9 genes could also be utilized to evaluate and predict the overall survival of patients with UVM. The risk score levels related to immune cell frequency and specific genomic alterations. The DEGs between the low- and high- risk group were enriched in tumor OXPHOS and immune related pathway. In vitro experiments, CYC1 silencing significantly inhibited UVM cell proliferation and invasion, induced cell apoptosis.

In sum, a prognostic risk score model based on oxidative phosphorylation-related genes in UVM was developed to enhance understanding of the disease. This prognostic risk score model may help to find potential therapeutic targets for UVM patients. CYC1 acts as an oncogene role in UVM.