摘要
为探讨RNA m^(6)A甲基化调节因子在肺腺癌中的作用,从TCGA数据库下载肺腺癌患者的RNA表达数据和临床数据。通过limma软件包分析12种m^(6)A调节剂的表达情况。使用Pheatmap、vioplot和corrplot软件包生成热图、小提琴图和表达相关图。采用Kaplan⁃Meier方法分别计算肺腺癌中12种RNA m^(6)A调节因子的生存曲线。使用Cox回归和Kaplan⁃Meier方法分析TCGA肺腺癌患者的总体存活相关的临床病理学特征。最后用Kruskal(KS)检验和logistic回归分析临床病理学特征与HNRNPC表达的关系。在肺腺癌的TCGA队列中,发现HNRNPC、WTAP、YTHDF3、FTO、ZC3H13、METTL14、METTL3、YTHDF1、YTHDF2这些基因是差异表达的。Kaplan⁃Meier生存分析显示,在这些差异表达的基因中仅仅HNRNPC和YTHDF2的表达与生存显著相关。然后,通过多因素Cox回归结果表明HNRNPC的表达在肺腺癌TCGA队列中是个独立危险因素。最后,HNRNPC在肺腺癌中的表达与临床分期(IV vs I,OR=3.692308)和组织浸润(T2 vs T1,OR=1.776471;T4 vs T1,OR=6.30303)显著相关(所有p<0.05)。结论认为HNRNPC可能作为肺腺癌的独立的预后因子。
This study aimed to investigate the role of RNA m^(6)A methylation regulators in lung adenocarcinoma.RNA expression data and clinical data of lung adenocarcinoma patients were downloaded from the TCGA database.The expression of 12 m^(6) A regulators was analyzed by the limma software package.Pheatmap,violin maps,and expression correlation maps were generated using Pheatmap,vioplot,and corrplot software packages.Kaplan⁃Meier method was used to calculate the survival curves of 12 RNA m^(6)A regulators in lung adenocarcinoma.Cox regression and Kaplan⁃Meier methods were used to correlate clinicopathological features associated with overall survival in patients with TCGA lung adenocarcinoma.Finally,the relationship between clinicopathological features and HNRNPC expression was analyzed by Kruskal(KS)test and logistic regression.In the TCGA cohort of lung adenocarcinoma,it was found that HNRNPC,WTAP,YTHDF3,FTO,ZC3H13,METTL14,METTL3,YTHDF1,and YTHDF2 were differentially⁃expressed.Kaplan⁃Meier survival analysis showed that only the expression of HNRNPC and YTHDF2 was significantly associated with survival among these differentially expressed genes.Then,results of multivariate Cox regression showed that the expression of HNRNPC was an independent risk factor in the TCGA cohort of lung adenocarcinoma.Finally,the expression of HNRNPC in lung adenocarcinoma was significantly correlated with clinical stage(IV vs I,OR=3.692308)and tissue invasion(T2 vs T1,OR=1.776471;T4 vs T1,OR=6.30303)(all p values<0.05).Conclusions cna be drawn that HNRNPC may serve as an independent prognostic factor for lung adenocarcinoma.
作者
胡滨滨
张明
HU Binbin;ZHANG Ming(Cancer Institute,Affiliated Gaoyou Hospital of Soochow University,Gaoyou 225600,Jiangsu,China)
出处
《生物信息学》
2022年第2期124-135,共12页
Chinese Journal of Bioinformatics