Construction and validation of immune-related LncRNAs classifier to predict prognosis and immunotherapy response in laryngeal squamous cell carcinoma

Abstract

Abstract Background Rapid advances in transcriptomic profiles have resulted in recognizing IRLs (immune-related long noncoding RNAs), as modulators of the expression of genes related to immune cells that mediate immune inhibition as well as immune stimulatory, indicating LncRNAs play fundamental roles in immune modulation. Hence, we establish an IRL classifier to precisely predict prognosis and immunotherapeutic efficiency in laryngeal squamous cell carcinoma (LSCC). Methods LSCC RNA-seq (RNA sequencing) datasets, somatic mutation data, and corresponding clinicopathologic information were acquired from TCGA (the Cancer Genome Atlas) and Gene Expression Omnibus (GEO) databases. Spearman correlation analysis identified LncRNAs associated with immune-related genes (IRG). Based on Lasso penalized regression and random forest (RF), we constructed an IRL classifier associated with prognosis. GEO database was utilized to validate the IRL classifier. The predictive precision and clinical application of the IRL classifier were assessed and compared to clinicopathologic features. The immune cell infiltration of LSCC was calculated via CIBERSORTx tools and ssGSEA (single-sample gene set enrichment analysis). Then, we systematically correlated the IRL classifier with immunological characteristics from multiple perspectives, such as immune-related cells infiltrating, tumor microenvironment (TME) scoring, microsatellite instability (MSI), tumor mutation burden (TMB), and chemokines. Finally, the TIDE (tumor immune dysfunction and exclusion) algorithm was used to predict response to immunotherapy. Results Based on machine learning approach, three prognosis-related IRLs (BARX1-DT, KLHL7-DT, and LINC02154) were selected to build an IRL classifier. The IRL classifier could availably classify patients into the low-risk and high-risk groups based on the different endpoints, including recurrence-free survival (RFS) and overall survival (OS). In terms of predictive ability and clinical utility, the IRL classifier was superior to other clinical characteristics. Encouragingly, similar results were observed in the GEO databases. Immune infiltration analysis displayed immune cells that are significantly richer in low-risk group, CD8 T cells and activated NK cells via CIBERSORTx algorithm as well as activated CD8 T cell via ssGSEA. Additionally, compared with the high-risk group, immune score, CD8 T effector was higher in the low-risk group, yet stromal score, score of p53 signaling pathway and TGFher in the Tx algorithm, was lower in the low-risk group. Corresponding results were confirmed in GEO dataset. Finally, TIDE analysis uncovered that the IRL classifier may be effectually predict the clinical response of immunotherapy in LSCC. Conclusion Based on BARX1-DT, KLHL7-DT, and LINC02154, the IRL classifier was established, which can be used to predict the prognosis, immune infiltration status, and immunotherapy response in LSCC patients and might facilitate personalized counseling for immunotherapy.