Liu Jingchao, Wang Jianlong, Wu Meng, Zhang Wei, Meng Lingfeng, Wang Jiawen, Lv Zhengtong, Xia Haoran, Zhang Yaoguang, Wang Jianye
Department of Urology, Beijing Hospital, National Center of Gerontology, Institute of Geriatric Medicine, Chinese Academy of Medical Sciences, Beijing, China.
Graduate School of Peking Union Medical College, Chinese Academy of Medical Sciences, Beijing, China.
Front Med (Lausanne). 2021 Dec 23;8:757432. doi: 10.3389/fmed.2021.757432. eCollection 2021.
To comprehensively analyze N6-methyladenosine modification patterns in bladder tumors and to further systematically explore the inherent relationships between these modification patterns and multiomic tumor characteristics. A total of 901 bladder tumor samples, including 405 samples from TCGA database, 188 samples from GSE13507 and 308 samples from GSE32894, were included in this systematic analysis. The N6-methyladenosine modification patterns were identified utilizing unsupervised clustering analysis. To quantify N6-methyladenosine modification patterns, the m6Ascore of individual sample was developed using principal component analysis algorithms. Relationships among immune infiltration, tumor mutation burden, various clinical characteristics, molecular subtypes, and the m6Ascore were systematically analyzed. The guiding value of m6Ascore in immunotherapy was further validated in an external trial cohort. Genomics of Drug Sensitivity in Cancer expression references were also utilized to perform drug sensitivity analysis for patients with distinct m6A modification patterns. We determined three different N6-methyladenosine modification patterns for 901 bladder tumors. The quantitative m6Ascore of individual sample derived from N6-methyladenosine modification patterns could play a significant role in predicting overall survival, immune cell infiltration, and classic oncogene mutations. A low m6Ascore combined with high tumor mutation burden indicated better survival outcomes ( < 0.001). A higher m6Ascore also indicated a higher grade, higher T and N stage, elder ages, higher death rate, and higher PD1/PDL1/CTLA4 expressions ( < 0.01). The Basal type tended to exhibit significantly higher m6Ascores than the Luminal and Neuronal subtypes. External immunotherapy cohorts demonstrated that no difference in therapeutic effects was noted between the high and low m6Ascore groups when anti-PD1 immunotherapy was exclusively administered. When anti-PD1 and anti-CTLA4 immunotherapy were simultaneously administered, the high m6Ascore group had a significantly better prognosis than the low m6Ascore group ( < 0.001). High m6A groups were potentially sensitive to various medical treatments including Bleomycin, Bortezomib, Cisplatin, Cyclopamine, Dasatinib, Docetaxe, Rapamycin, and Vinblastine in this study. This study systematically revealed the important roles of m6A methylation modification patterns in bladder tumors. Detailed quantification of m6A modification patterns could improve our understanding of the bladder tumor microenvironments and could provide guidance for future immunotherapy strategies.
为全面分析膀胱肿瘤中的N6-甲基腺嘌呤修饰模式,并进一步系统探索这些修饰模式与多组学肿瘤特征之间的内在关系。本系统分析共纳入901例膀胱肿瘤样本,包括来自TCGA数据库的405例样本、来自GSE13507的188例样本和来自GSE32894的308例样本。利用无监督聚类分析确定N6-甲基腺嘌呤修饰模式。为量化N6-甲基腺嘌呤修饰模式,使用主成分分析算法开发了个体样本的m6A评分。系统分析了免疫浸润、肿瘤突变负荷、各种临床特征、分子亚型与m6A评分之间的关系。m6A评分在免疫治疗中的指导价值在外部试验队列中得到进一步验证。还利用癌症药物敏感性基因组学表达参考对具有不同m6A修饰模式的患者进行药物敏感性分析。我们确定了901例膀胱肿瘤的三种不同N6-甲基腺嘌呤修饰模式。源自N6-甲基腺嘌呤修饰模式的个体样本的定量m6A评分在预测总生存期、免疫细胞浸润和经典癌基因突变方面可能发挥重要作用。低m6A评分与高肿瘤突变负荷相结合表明生存结果更好(<0.001)。较高的m6A评分还表明分级较高、T和N分期较高、年龄较大、死亡率较高以及PD1/PDL1/CTLA4表达较高(<0.01)。基底型往往比管腔型和神经型亚型表现出显著更高的m6A评分。外部免疫治疗队列表明,单独给予抗PD1免疫治疗时,高m6A评分组和低m6A评分组在治疗效果上没有差异。当同时给予抗PD1和抗CTLA4免疫治疗时,高m6A评分组的预后明显优于低m6A评分组(<0.001)。在本研究中,高m6A组对包括博来霉素、硼替佐米、顺铂、环杷明、达沙替尼、多西他赛、雷帕霉素和长春碱在内的各种药物治疗可能敏感。本研究系统揭示了m6A甲基化修饰模式在膀胱肿瘤中的重要作用。对m6A修饰模式的详细量化可以提高我们对膀胱肿瘤微环境的理解,并可为未来的免疫治疗策略提供指导。