Immunogenic profiling of metastatic uveal melanoma discerns a potential signature related to prognosis

Abstract

Abstract Background Uveal melanoma (UM) is an aggressive intraocular malignant tumor. The present study aimed to identify the key genes associated with UM metastasis and established a gene signature to analyze the relationship between the signature and prognosis and immune cell infiltration. Later, a predictive model combined with clinical variables was developed and validated. Methods Two UM gene expression profile chip datasets were downloaded from TCGA and GEO databases. Immune-related genes (IRGs) were obtained from IMPORT database. First, these mRNAs were intersected with IRGs, and weighted gene co-expression network analysis (WGCNA) was used to identify the co-expression of genes primarily associated with metastasis of UM. Univariate Cox regression analysis screened the genes related to prognosis. LASSO-Cox established a risk score to distinguish high-risk group and low-risk group. Then the GSEA enrichment pathway and immune cell infiltration of the two groups were compared. And combined with clinical variables, a predictive model was constructed. The time-dependent receiver operating characteristic (ROC) curve, calibration curve, and decision curve analysis (DCA) curve were used to verify the stability and accuracy of the final predictive model, and a nomogram was then drawn. Results The MEblack, MEpurple, and MEblue modules were significantly associated with the metastasis of UM patients (P value < 0.001, = 0.001, = 0.022, respectively). Four genes (UBXN2B, OTUD3, KAT8, LAMTOR2) were obtained by Pearson correlation analysis, weighted gene correlation network analysis (WGCNA), univariate Cox, and LASSO-Cox. And a novel prognostic risk score was established. Immune-related prognostic signature can well classify UM patients into high-risk and low-risk groups. Kaplan–Meier curve showed that the OS of high-risk patients was worse than that of low-risk patients. In addition, the risk score played an important role in evaluating the signaling pathway and immune cell infiltration of UM patients in high-risk and low-risk groups. Both the training set and validation set of the model showed good predictive accuracy in the degree of differentiation and calibration (e.g., 1-year overall survival: AUC = 0.930 (0.857–1.003)). Finally, a nomogram was established to serve in clinical practice. Significance UM key gene signature and prognosis predictive model might provide insights for further investigation of the pathogenesis and development of UM at the molecular level, and provide theoretical basis for determining new prognostic markers of UM and immunotherapy.