1424 m⁶A regulators related to molecular feature, immune microenvironment and prognosis in lung adenocarcinoma

Zhibo Tan,Min Chen,Yujie Zhao,Feng Peng,Ying Li,Mengqi Yang,Lei Zhang,Xin Li,Pengfei Yang,Zhe Zhang,Daming Li,Zhaoming Peng,Yajie Liu

DOI: https://doi.org/10.1136/jitc-2022-sitc2022.1424

2022-01-01

Abstract:<h3>Background</h3> The dynamic changes of N⁶-methyladenosine controlled by m⁶A regulators are vital for cancer biology.¹ However, the role of m⁶A regulators in lung adenocarcinoma (LUAD) remains unclear. This study aimed to gain a deeper insight into the roles of m⁶A regulators in LUAD by bioinformatics. <h3>Methods</h3> 31 m⁶A regulators were retrieved from previous literature. The expression levels of the 31 regulators were investigated by TIMER2.0. The association between the 31 regulators and overall survival (OS) were explored via K-M plotter and Cox regression. Next, mutations status and copy number alterations of the 31 regulators were probed through cBioPortal. Moreover, the interaction between the 31 regulators were revealed via principal component analysis (PCA), protein-protein interaction (PPI), and Metascape enrichment by R programming, STRING 11.5, Cytoscape, and Metascape. Then, the immune microenvironment including immune score, specific immune cells infiltration, and expressions of immune-related gene families were probed via ESTIMATE, TIMER2.0, TIMER, and R programming. Next, LUAD samples were divided into 3 clusters based on the expression patterns of the 31 regulators by k-means method. The differences of OS, gene enrichment, immune microenvironment of the 3 clusters were compared by Log-rank test, GSVA, ESTIMATE, and CIBERSORTx. Moreover, a m⁶A Score Model was constructed using stepwise regression in this study. The predicting ability of the m⁶A Score Model was verified by GEO LUAD dataset and TCGA Pan- adenocarcinoma dataset. <h3>Results</h3> 22 out of 31 regulators showed differential expression levels between LUAD and normal lung (figure 1). 28 out of 31 regulators were linked with OS (figure 2). The frequency of mutation and copy number alterations were 26.6% and 29.1% respectively (figure 3). PCA revealed the 31 m⁶A regulators displayed almost completely distinct distribution patterns between LUAD and normal lung (figure 4). PPI and Metascape enrichment indicated the 31 m⁶A regulators displayed close interrelation (figure 5). Between high and low expression groups of 20, 20 and 29 regulators, tumor purity, immune score, immune cells infiltrations exhibited statistical difference (figure 6). Moreover, 3 clusters displayed distinct OS, gene enrichment, and immune microenvironment (figure 7). Notably, a m⁶A Score Model was built up upon 15 m⁶A regulators (figure 8). The model was associated with OS for LUAD and pan-adenocarcinoma (figure 9). Conclusions m⁶A regulators are associated with molecular features, immune microenvironment and prognosis in LUAD. m⁶A Score could be used as a reference to immunotherapy and a prognostic index for LUAD. <h3>Reference</h3> Nombela P, Miguel-López B, Blanco S. The role of m6A, m5C and ? RNA modifications in cancer: Novel therapeutic opportunities. <i>Mol Cancer</i> 2021;<b>20</b>:18.