Developing a RiskScore Model based on Angiogenesis-related lncRNAs for Colon Adenocarcinoma Prognostic Prediction

Aim: We screened key angiogenesis-related lncRNAs based on colon adenocarcinoma (COAD) to construct a RiskScore model for predicting COAD prognosis and help reveal the pathogenesis of the COAD as well as optimize clinical treatment. Background: Regulatory roles of lncRNAs in tumor progression and prognosis have been confirmed, but few studies have probed into the role of angiogenesis-related lncRNAs in COAD. Objective: To identify key angiogenesis-related lncRNAs and build a RiskScore model to predict the survival probability of COAD patients and help optimize clinical treatment. Methods: Sample data were collected from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. The HALLMARK pathway score in the samples was calculated using the single sample gene set enrichment analysis (ssGSEA) method. LncRNAs associated with angiogenesis were filtered by an integrated pipeline algorithm. LncRNA-based subtypes were classified by ConsensusClusterPlus and then compared with other established subtypes. A RiskScore model was created based on univariate Cox, least absolute shrinkage and selection operator (LASSO) regression and stepwise regression analysis. The Kaplan-Meier curve was drawn by applying R package survival. The time-dependent ROC curves were drawn by the timeROC package. Finally, immunotherapy benefits and drug sensitivity were analyzed using tumor immune dysfunction and exclusion (TIDE) software and pRRophetic package. Results: Pathway analysis showed that the angiogenesis pathway was a risk factor affecting the prognosis of COAD patients. A total of 66 lncRNAs associated with angiogenesis were screened, and three molecular subtypes (S1, S2, S3) were obtained. The prognosis of S1 and S2 was better than that of S3. Compared with the existing subtypes, the S3 subtype was significantly different from the other two subtypes. Immunoassay showed that immune cell scores of the S2 subtype were lower than those of the S1 and S3 subtypes, which also had the highest TIDE scores. We recruited 8 key lncRNAs to develop a RiskScore model. The high RiskScore group with inferior survival and higher TIDE scores was predicted to benefit limitedly from immunotherapy, but it may be more sensitive to chemotherapeutics. A nomogram designed by RiskScore signature and other clinicopathological characteristics shed light on rational predictive power for COAD treatment. Conclusion: We constructed a RiskScore model based on angiogenesis-related lncRNAs, which could serve as potential prognostic predictors for COAD patients and may offer clues for the intervention of anti-angiogenic application. Our results may help evaluate the prognosis of COAD and provide better treatment strategies.