RESUMO
BACKGROUND: Melanomas are malignant tumors that can occur in different body parts or tissues such as the skin, mucous membrane, uvea, and pia mater. Long non-coding RNAs (lncRNAs) are key factors in the occurrence and development of many malignant tumors, and are involved in the prognosis of some patients. AIM: To identify autophagy-related lncRNAs in melanoma that are crucial for the diagnosis, treatment, and prognosis of melanoma patients. METHODS: We retrieved transcriptome expression profiles and clinical information of 470 melanoma patients from The Cancer Genome Atlas (TCGA) database. Then, we identified autophagy-related genes in the Human Autophagy Database. Using R, coexpression analysis of lncRNAs and autophagy-related genes was conducted to obtain autophagy-related lncRNAs and their expression levels. We also performed univariate and multivariate Cox proportional risk analyses on the obtained datasets, to systematically evaluate the prognostic value of autophagy-related lncRNAs in melanoma. Fifteen autophagy-related lncRNAs were identified and an autophagy-related prognostic signature for melanoma was established. The Kaplan-Meier and univariate and multivariate Cox regression analyses were used to calculate risk scores. Based on the risk scores, melanoma patients were randomly divided into high- and low-risk groups. Receiver operating characteristic curve analysis, dependent on time, was performed to assess the accuracy of the prognostic model. At the same time, we also downloaded the melanoma data sets GSE65904, GSE19234, and GSE78220 from the GENE EXPRESSION OMNIBUS database for model verification. Finally, we performed Gene Set Enrichment Analysis functional annotation, which showed that the low and the high-risk groups had different enriched pathways. RESULTS: The co-expression network for autophagy-related genes was constructed using R, and 936 lncRNAs related to autophagy were identified. Then, 52 autophagy-related lncRNAs were significantly associated with TCGA melanoma patients' survival by univariate Cox proportional risk analysis (P < 0.01). Further, the 52 autophagy-related lncRNAs mentioned above were analyzed by multivariate Cox analysis with R. Fifteen lncRNAs were selected: LINC01943, AC090948.3, USP30-AS1, AC068282.1, AC004687.1, AL133371.2, AC242842.1, PCED1B-AS1, HLA-DQB1-AS1, AC011374.2, LINC00324, AC018553.1, LINC00520, DBH-AS1, and ITGB2-AS1. The P values in all survival analyses using these 15 lncRNAs were < 0.05. These lncRNAs were used to build a risk model based on the risk score. Negative correlations were observed between risk scores and overall survival rate in melanoma patients over time. Additionally, the melanoma risk curve and scatter plot analyses showed that the death number increased along with the increase in the risk score. Overall, we identified and established a new prognostic risk model for melanoma using 15 autophagy-related lncRNAs. The risk model constructed with these lncRNAs can help and guide melanoma patient prognosis predictions and individualized treatments in the future. CONCLUSION: Overall, the risk model developed based on the 15 autophagy-related lncRNAs can have important prognostic value and may provide autophagy-related clinical targets for melanoma treatment.
