摘要

为探讨RNAm6A甲基化调节因子在肺腺癌中的作用,从TCGA数据库下载肺腺癌患者的RNA表达数据和临床数据。通过limma软件包分析12种m6A调节剂的表达情况。使用Pheatmap、vioplot和corrplot软件包生成热图、小提琴图和表达相关图。采用Kaplan-Meier方法分别计算肺腺癌中12种RNAm6A调节因子的生存曲线。使用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在肺腺癌中的表达与临床分期(IVvsI,OR=3.692308)和组织浸润(T2vsT1,OR=1.776 471;T4 vs T1, OR=6.303 03)显著相关(所有p<0.05)。结论认为HNRNPC可能作为肺腺癌的独立的预后因子。