Skip to content

A 级候选的完整证据链与两处硬伤

2026-07-31 为汇报 PPT 重新定位(改为汇报课题本身)而做的全量数字核对, 查出两处会直接影响结论表述的问题。本页如实记录,已同步进 PPT(结果页 + 局限页)。

数字来源:results/ 当前产物,提取脚本 scratchpad/facts.R + facts2.R; PPT 审计脚本 audit_ppt.py 44 项回查全部 PASS。

配套:baseline 文章逐条对比 · 共定位分档与 SuSiE 降级 · 反向 MR · R9 主线


一、三个 A 级候选的完整证据链

ERMAPIFNAR1APOL1
结局糖尿病黄斑病变糖尿病黄斑病变糖尿病黄斑病变
工具变异rs11210710rs914142rs136168
变异类型基因上游变异内含子变异内含子变异
工具 F169.03,163.91,248.4
OR(95% CI)0.329(0.197–0.552)0.782(0.696–0.878)1.521(1.249–1.853)
FDR4.3e-034.9e-034.9e-03
HyPrColoc 簇后验0.8990.8370.815
coloc.abf PP.H40.9670.9400.926
糖尿病对照并发症特异并发症特异并发症特异
反向 MR(FDR)0.883 未显著0.618 未显著0.735 未显著
多效性
视网膜落位AUC 0.536 无实质富集AUC 0.757 神经节细胞AUC 0.506 无实质富集
外部复制无(黄斑病变无外部队列)
可成药性可成药,无已知药可成药,方向一致,12 个已知药可成药,无已知药

前十一行三者完全对等;差异全部出在最后三行。


二、硬伤一:ERMAP 与 APOL1 的组织落位证据不成立

问题

单细胞富集表(results/singlecell/celltype_enrichment.csv)里,ERMAP 与 APOL1 都有多条 sig = TRUE 的记录。但看 AUC 就会发现,这些「显著」没有任何实质含义

三个候选在人视网膜 18 种细胞类型中的完整数值:

基因最强细胞类型表达细胞比例AUC全部 18 类的 AUC 范围
IFNAR1伞状神经节细胞74.4%0.7570.464 – 0.757
IFNAR1视网膜神经节细胞63.8%0.711
IFNAR1小神经节细胞62.8%0.713
ERMAP视锥细胞9.3%0.5360.489 – 0.536
APOL1伞状神经节细胞1.6%0.5060.498 – 0.506

(对照:APOE 在 Müller 细胞 91.7%、AUC 0.945)

为什么会这样

单细胞图谱有数万个细胞。在这种样本量下,AUC 0.506 也能得到 p = 0sig = TRUE 反映的是统计显著性,与富集强度无关。

  • APOL1 的最高表达比例只有 1.6%,18 种细胞类型的 AUC 全部落在 0.498–0.506 —— 与随机(0.5)无法区分。
  • ERMAP 最高 9.3%,AUC 最高 0.536 —— 同样谈不上特异定位。

结论

不能声称 ERMAP 与 APOL1 在视网膜有特异定位

组织落位这一层只支持 IFNAR1(神经节细胞,74% 细胞表达,AUC 0.757)。

ERMAP 与 APOL1 需要另找组织学证据:换用 RPE / 脉络膜数据集, 或直接做免疫组化与原位杂交。

若按 FDR 显著性报出去,审稿人核对 AUC 时会立刻发现问题。

这是一个通用教训:单细胞富集分析必须同时报告效应量(AUC / 表达比例), 只报 p 值或 FDR 会系统性夸大结论。


三、硬伤二:三个 A 级候选全无外部复制

问题

三个 A 级候选全部指向糖尿病黄斑病变,而黄斑病变没有可用的外部队列

现有的外部复制覆盖情况(results/screen_R9_dm/replication_external.csv):

复制队列对数方向一致
MVP 糖网(GCST90475689)2422(92%)
Salem 2019 糖肾96(67%)
MVP 神经(GCST90475676)54(80%)
糖尿病黄斑病变0

暴露源对照(换平台不换人群)也帮不上忙:decode-rep 目前只跑了 AMD / 严格定义糖网两组结局, 没有覆盖 FinnGen R9 的糖尿病黄斑病变

结论

A 级候选目前是「内部证据完整、外部证据为零」

必须在正文与局限中明确写出。可行的补救方向:

  1. 换暴露平台做交叉验证(成本最低)——用 deCODE SomaScan 的 cis-pQTL 重跑 ERMAP / IFNAR1 / APOL1 × FinnGen R9 黄斑病变。这不需要新的结局队列, 且能同时回应 Olink 亲和力平台的表位效应质疑。 ⚠ 需先确认 deCODE 是否有这三个蛋白的适配体。
  2. 寻找外部黄斑病变端点——MVP 糖尿病眼病端点、英国生物银行自建端点,或与已发表 GWAS 合作。
  3. 退一步的写法——把 A 级定位为「待复制的优先候选」,而非「已验证靶点」。

四、这两处对全局结论的影响

层面影响
A 级三候选是否成立内部证据链(MR + 共定位 + 糖尿病对照 + 反向 MR + 多效)不受影响,五道全部通过
能否称为「靶点」不能。缺外部复制,只能称「优先候选」
三者是否等价不等价。只有 IFNAR1 同时具备组织落位与方向一致的已上市药
若只推一个IFNAR1——证据链唯一完整的一个

变异类型这一行是加分项,别漏掉

三个候选的工具变异全部是非蛋白改变型(ERMAP 上游变异、IFNAR1 与 APOL1 内含子变异), 因此表位效应风险低——这对 Olink 这类亲和力平台是必须回应的质疑。 数据一直在显著表的 gene_consequence 列里,只是以前没提取出来。

个人科研与运维文档 · 内容持续修订