gsea富集分析结果怎么看_单基因联合自噬分析

最新推荐文章于 2023-08-09 11:46:59 发布

weixin_39922642

最新推荐文章于 2023-08-09 11:46:59 发布

阅读量1.1k

点赞数

文章标签： gsea富集分析结果怎么看

Identification of molecular correlations of RBM8A with autophagy in Alzheimer's disease

二. 文章思路

三. 结果解读

作者探索RBM8A在AD中的作用使用的是GSE33000数据集，样本为310AD患者 VS 157 norm，用limma包进行差异分析

A：箱线图展示RBM8A在AD患者和正常样本中的表达量。作者在AD vs norm的差异分析结果中发现RBM8A在AD中表达量显著降低(p=1.620e-19，logFC=-0.078)。之后根据AD患者RBPM8A表达量的中位数分为高低表达组。
B：AD vs norm的差异分析火山图
C：PBM8A高低表达组间的差异分析火山图。结合B中的结果取交集，得到共同的9186个差异表达基因(DEGs)
D：热图展示在AD中表达量上调/下调的DEGs中差异最显著的前25个DEGs在各样本中的表达量

图1. 分析AD中的差异表达基因

A：根据结果一中分析出的9186个与RBM8A相关的AD相关的基因，进行WGCNA(基因加权共表达网络分析)，共得到10个模块
B：10个共表达模块与临床表型进行关联性分析。可以看到棕色模块与AD的发生呈显著正相关(r=0.69，p=3e-67)；青绿色模块与AD的发生呈显著的负相关(r=-0.69，p=3e-66)(两者都正相关/负相关模块中|r|最大的)
C：分析上述两个模块中关键基因子集，并用相关图展示他们表达量之间的相关系数。一共识别出包含RBM8A在内的15个关键基因(GS>0.7，GS即基因与AD间的相关性；MM>0.9，MM即基因的模块成员度，即各基因表达量与相应模块特征基因的相关性)。图中可以看出，有12个基因与RBM8A的表达量呈显著正相关

上述分析按照WGCNA包的流程进行

图2. 对DEGs进行WGCNA的结果

、

D-E：针对一些模块中的基因分别进行GO：BP分析(B)以及KEGG分析(E)，用的是clusterProfiler这个包。图中右侧条形图表示各个条目富集的基因个数以及表达上调还是下调的比例。作者在结果中侧重于对棕色模块以及青绿色模块中基因功能富集分析结果的解读

图3.对各模块中的基因进行功能富集分析

通过GSEA分析对结果2中的功能富集分析结果进行验证，基因集分别来自MsigDB数据库中的GO：BP以及KEGG基因集

图4.GSEA分析结果

E：作者发现RBM8A的低表达会影响许多自噬通路相关的基因((FIP200，Beclin 1，NRBF2，VPS15，ATG12)，并给出了相关机制图。蓝色表示在AD中表达下调的基因，红色表示在AD中表达上调的基因

图5.RBMBA影响细胞自噬机制图

作者提取各样本15个核心基因的基因表达谱，以7：3的比例将样本分为训练集和测试集，用glmmet包进行LASSO回归。另外以GSE5281(87AD VS 74norm)和GSE483350(80AD VS 173norm)数据集为验证集

A：LASSO回归模型构建过程中的CV统计图，作者在构建模型时λ值选择是lambda.min(左侧虚线)，大小为0.0401，此时对应的预测基因有8个。作者构建的8基因AD预测模型公式如下：index=RBM8A(-2.38668779488564)+RHBDF2(2.00115481990953)+TNFRSF10B(0.817520478917702)+ACP1(-3.41028393841058)+ANKRD39(-0.279104767589027)+CA10(-0.988480656766608)+CBLN4(-0.70290303609009) + PPEF1(-1.67911758870231)
B：ROC曲线展示上述预测模型在训练集中的预测效果，AUC=0.948(越接近1表示模型预测效果越好)
C：ROC曲线预测上述模型在测试集中的预测效果，AUC=0.947(左)；ROC曲线预测上述模型在验证集中的预测效果，AUC=0.948
D：RBM8A在GSE5281数据集的AD患者中低表达
E：RBM8A在GSE48350数据集的AD患者的一些脑区中低表达(横坐标表示不同大脑组织)

图6.AD预测模型的构建以及验证

关注