clusterProfiler、GSEA 与 GSVA
clusterProfiler 常用于 GO、KEGG 过表达分析,输入是差异基因列表;GSEA 使用全部基因排序,能发现整体小幅协同变化的通路;GSVA 和 ssGSEA 为每个样本算通路得分,便于分组比较或关联表型。
15个分析 Skill · Codar 生信插件
网络与系统生物学把单个基因的变化放到通路、模块和调控关系里理解,典型输入是差异基因列表、表达矩阵或时间序列数据。这个插件在 Codar 里覆盖通路与功能富集、GSEA 与样本级基因集评分、WGCNA 共表达网络、蛋白互作网络、基因调控网络推断、差异网络、网络药理学、SBML 建模和网络可视化。
确定手上是基因列表、完整表达矩阵还是时间序列,这决定能用富集、共表达还是调控推断。
统一基因 ID,过滤低表达和低变异基因,检查离群样本;时序数据要确认时间点和重复。
按输入类型选择过表达分析或 GSEA,必要时对每个样本计算通路评分,用于后续比较。
做 WGCNA 模块、蛋白互作网络或调控网络推断,并比较不同条件下网络连接的变化。
计算度、中介中心性等指标并做社区划分,识别枢纽节点和功能模块,再对模块做富集。
输出网络图、模块-性状热图、富集图、节点表和报告,附可复现脚本。
clusterProfiler 常用于 GO、KEGG 过表达分析,输入是差异基因列表;GSEA 使用全部基因排序,能发现整体小幅协同变化的通路;GSVA 和 ssGSEA 为每个样本算通路得分,便于分组比较或关联表型。
通过基因间相关性构建加权共表达网络并划分模块,再把模块特征基因与性状关联。软阈值通常按无标度拓扑拟合度选择,常用 0.8 左右作为参考。作者建议样本数不少于 15 个,样本过少时模块不稳定。
STRING 汇总实验、共表达和文献挖掘等来源的蛋白互作,常用综合分数 0.4 或 0.7 作为中、高置信阈值。Cytoscape 用于网络展示和拓扑分析,MCODE、cytoHubba 等插件常用来找子网络和枢纽节点。
用于从表达数据推断基因调控关系。GENIE3 基于随机森林,ARACNe 基于互信息并剔除间接关联,SCENIC 在共表达基础上加入转录因子结合基序验证,常用于单细胞数据。推断出的调控边是统计关联,需要实验验证因果方向。
共 15 项。启用插件后,Codar 会按任务自动选用,不需要逐个记住名称。
在 Codar 左上角选择 Codar 生信。
在插件页找到「网络与系统生物学」,安装或启用。
用一句话写清数据和目标,确认分析方案后开始执行。
帮我对这份表达矩阵做 WGCNA,找出与临床性状最相关的模块,并对模块基因做通路富集。
GO/KEGG 过表达分析只看按阈值筛出来的差异基因,结果受阈值影响较大;GSEA 使用全部基因按差异程度的排序,不需要设阈值,能发现多数基因小幅同向变化的通路。差异基因很少或很多时,GSEA 往往更稳健,两者也可以结合报告。
WGCNA 作者在官方说明中建议样本数不少于 15 个,样本更多时模块更稳定。样本过少时,相关系数估计噪声大,模块划分容易随参数改变。样本有限时,可以考虑先做差异分析再用 GSEA 等方法补充,而不是强行构建共表达网络。
网络药理学基于数据库中的成分、靶点和通路关系做预测,结果受数据库收录质量影响,靶点预测本身有较高假阳性。它适合用来提出假设和缩小验证范围,不能单独作为作用机制的证据,关键靶点通常需要分子对接或细胞、动物实验进一步验证。