基因调控网络推断与共表达分析11 / 16
先验知识整合:TF–靶基因数据库、ChIP-seq 与 ATAC
为什么需要先验
表达数据本身只含关联信息;在真核生物中,DREAM5 显示纯表达数据驱动的推断在酵母上表现很有限。而关于“谁可能调控谁”,生物学已经积累了大量独立证据。把这些证据作为先验,可以缩小搜索空间、为方向提供依据,并减少间接边。
先验的来源
| 来源 | 证据类型 | 代表资源 | 主要局限 |
|---|---|---|---|
| 文献整理 | 小规模实验验证的调控关系 | RegulonDB(大肠杆菌)、TRRUST v2(人和小鼠,Han 等, 2018) | 偏向研究充分的转录因子;缺乏细胞类型信息 |
| 整合资源 | 多来源证据合并并分级 | DoRothEA(Garcia-Alonso 等, 2019,按证据强度分 A–E 级)、CollecTRI(Müller-Dott 等, 2023) | 不同证据的可比性;版本更新 |
| ChIP-seq | 转录因子在全基因组的结合位点 | ENCODE、ChIP-Atlas | 结合不等于调控;依赖抗体与细胞系 |
| 基序 | 调控区中的转录因子结合基序 | cisTarget 数据库、JASPAR 等基序库 | 同家族转录因子基序相似;假阳性多 |
| 染色质可及性 | ATAC-seq / scATAC-seq 显示的开放区域,与基序结合可推断候选结合 | SCENIC+、CellOracle 的基础 GRN | 开放不等于结合;增强子–基因配对不确定 |
整合策略
- 限定候选集:只允许转录因子作为调控因子(GENIE3、GRNBoost2 的做法);或只保留有基序/ChIP 支持的边(SCENIC 的剪枝步骤)。简单可靠,但先验中缺失的真边永远找不到。
- 作为软约束:在回归中对先验边减小惩罚(Inferelator 一类方法;Greenfield 等, 2013, Bioinformatics 研究了如何在先验有错时稳健地加入先验),或在贝叶斯网络中给结构加先验概率。数据仍可推翻先验。
- 估计转录因子活性:转录因子的 mRNA 水平常不反映其活性。把先验靶基因集合当作“读数”,用靶基因的表达估计转录因子活性(DoRothEA、CollecTRI 的主要用途),再用活性而非 mRNA 做下游分析或网络推断。
- 事后注释:先用数据推断网络,再检查高分边在先验中的富集,作为合理性检验。
代码示例的“回归与先验整合”单元模拟了一个“一半真边、一半假边”的先验:样本少()时,先验把 LASSO 路径的 AUPR 从 0.55 提高到 0.64;样本多()时,仅数据 AUPR 为 0.80,加入先验反而降到 0.70;两种情形下先验中错误边的得分都被抬高。先验的价值取决于数据强弱与先验质量。
循环论证:如果先验来自某数据库,评估时就不能再用同一数据库(或与之高度重叠的资源)作金标准,否则性能被系统性高估。BEELINE 使用的真实数据金标准正包括 ChIP-seq 与 DoRothEA、TRRUST 等非细胞类型特异的调控网络,比较方法时要留意每种方法是否已经“见过”这些边。此外,数据库偏向研究充分的基因,依赖先验的方法会继承这种偏倚(见数据库研究偏倚)。