Zhu Liping, Wang Zhiqiang, Sun Yilan, Giamas Georgios, Stebbing Justin, Yu Zhentao, Peng Ling
Department of Medical Oncology, Shouguang Hospital of Traditional Chinese Medicine, Shouguang, China.
Department of Urology, Shouguang Hospital of Traditional Chinese Medicine, Shouguang, China.
Front Oncol. 2021 Dec 22;11:778637. doi: 10.3389/fonc.2021.778637. eCollection 2021.
Alternative splicing (AS) is a gene regulatory mechanism that drives protein diversity. Dysregulation of AS is thought to play an essential role in cancer initiation and development. This study aimed to construct a prognostic signature based on AS and explore the role in the tumor immune microenvironment (TIME) in lung adenocarcinoma.
We analyzed transcriptome profiling and clinical lung adenocarcinoma data from The Cancer Genome Atlas (TCGA) database and lists of AS-related and immune-related signatures from the SpliceSeq. Prognosis-related AS events were analyzed by univariate Cox regression analysis. Gene set enrichment analyses (GSEA) were performed for functional annotation. Prognostic signatures were identified and validated using univariate and multivariate Cox regression, LASSO regression, Kaplan-Meier survival analyses, and proportional hazards model. The context of TIME in lung adenocarcinoma was also analyzed. Gene and protein expression data of Cyclin-Dependent Kinase Inhibitor 2A (CDKN2A) were obtained from ONCOMINE and Human Protein Atlas. Splicing factor (SF) regulatory networks were visualized.
A total of 19,054 survival-related AS events in lung adenocarcinoma were screened in 1,323 genes. Exon skip (ES) and mutually exclusive exons (ME) exhibited the most and fewest AS events, respectively. Based on AS subtypes, eight AS prognostic signatures were constructed. Patients with high-risk scores were associated with poor overall survival. A nomogram with good validity in prognostic prediction was generated. AUCs of risk scores at 1, 2, and 3 years were 0.775, 0.736, and 0.759, respectively. Furthermore, the prognostic signatures were significantly correlated with TIME diversity and immune checkpoint inhibitor (ICI)-related genes. Low-risk patients had a higher StromalScore, ImmuneScore, and ESTIMATEScore. AS-based risk score signature was positively associated with CD8+ T cells. CDKN2A was also found to be a prognostic factor in lung adenocarcinoma. Finally, potential functions of SFs were determined by regulatory networks.
Taken together, our findings show a clear association between AS and immune cell infiltration events and patient outcome, which could provide a basis for the identification of novel markers and therapeutic targets for lung adenocarcinoma. SF networks provide information of regulatory mechanisms.
可变剪接(AS)是一种驱动蛋白质多样性的基因调控机制。AS失调被认为在癌症的发生和发展中起重要作用。本研究旨在构建基于AS的预后特征,并探讨其在肺腺癌肿瘤免疫微环境(TIME)中的作用。
我们分析了来自癌症基因组图谱(TCGA)数据库的转录组图谱和临床肺腺癌数据,以及来自SpliceSeq的AS相关和免疫相关特征列表。通过单变量Cox回归分析来分析与预后相关的AS事件。进行基因集富集分析(GSEA)以进行功能注释。使用单变量和多变量Cox回归、LASSO回归、Kaplan-Meier生存分析和比例风险模型来识别和验证预后特征。还分析了肺腺癌中TIME的情况。从ONCOMINE和人类蛋白质图谱获得细胞周期蛋白依赖性激酶抑制剂2A(CDKN2A)的基因和蛋白质表达数据。对剪接因子(SF)调控网络进行可视化。
在1323个基因中筛选出了肺腺癌中总共19054个与生存相关的AS事件。外显子跳跃(ES)和互斥外显子(ME)分别表现出最多和最少的AS事件。基于AS亚型,构建了8个AS预后特征。高风险评分的患者总生存期较差。生成了在预后预测中具有良好有效性的列线图。1年、2年和3年风险评分的AUC分别为0.775、0.736和0.759。此外,预后特征与TIME多样性和免疫检查点抑制剂(ICI)相关基因显著相关。低风险患者具有较高的基质评分、免疫评分和估计评分。基于AS的风险评分特征与CD8 + T细胞呈正相关。还发现CDKN2A是肺腺癌的一个预后因素。最后,通过调控网络确定了SF的潜在功能。
综上所述,我们的研究结果表明AS与免疫细胞浸润事件和患者预后之间存在明显关联,这可为肺腺癌新型标志物和治疗靶点的鉴定提供依据。SF网络提供了调控机制的信息。