单细胞分析案例拆解

用 GWAS、eQTL 和单细胞分析,解决“位点附近多个基因,先验证谁”的问题

从位点附近的一串基因,收敛到知道为什么优先、在哪类细胞验证的候选路线。

作者:发布于 2026-10-13更新于 2026-10-10约 8 分钟阅读GWAS × QTL × 单细胞
从GWAS位点、QTL共定位到单细胞类型和候选靶基因优先级的科学示意图
案例说明:本文是一个脱敏的方法型案例,保留了真实项目中常见的决策矛盾与完整分析路径,但不披露具体疾病、客户单位、位点、基因名和未公开结果。图中候选与证据强弱均为结构化示意,不构成临床结论。

一个团队拿到 GWAS 结果后,通常很快就会遇到比“有没有显著位点”更难的问题:同一位点附近有多个基因,文献各自都能找到一些依据,究竟先验证谁?这个候选应放在哪类细胞里研究?下一步实验应该测什么?

这个案例的原始材料并不少:GWAS 汇总统计、相关组织的 eQTL/pQTL 资源、单细胞表达矩阵,以及若干已经整理过的候选基因。真正的困难是,现有结果停留在“位点附近有什么”,还没有变成“哪条证据支持哪个候选、它可能在哪种细胞状态中发挥作用、先做哪一步验证”。

问题并不是候选太少,而是每个候选都像“有道理”

最靠近 lead SNP 的基因,未必是遗传效应真正影响的基因。一个非编码变异可能通过远端调控改变另一个基因的表达;一个位点也可能包含多个相互独立的遗传信号。如果直接按物理距离、最小 P 值或文献热度排序,得到的通常是一张很长的名单,而不是一个可以执行的研究计划。

我们先把项目问题收窄为一句话:哪些候选同时获得遗传定位、分子调控和疾病相关细胞场景的支持,并且能被设计成下一步可检验的实验?

GWAS、QTL与单细胞联合分析从数据统一、精细定位、共定位、细胞定位到实验验证的流程图
图 1|分析不是把三类数据简单求交集,而是让每一层证据回答不同问题:遗传信号在哪里、可能影响哪个分子、又可能在哪类细胞中发生。

第一步:先把遗传信号本身做对

我们首先统一基因组版本、染色体位置、参考/替代等位基因和效应方向,并检查 GWAS、QTL 与 LD 参考之间是否一致。A/T、C/G 这类回文位点若缺少可靠频率信息,需要谨慎处理;LD 参考还应尽量匹配研究人群祖源。否则,后面的共定位图即使“看起来很像”,也可能只是等位基因翻转或 LD 结构不匹配造成的假象。

随后在每个重点区域识别条件独立信号,并使用能够处理多信号的精细定位框架。以 SuSiE-RSS 为例,它可以基于汇总统计和 LD 估计每个变异的后验纳入概率(PIP),并形成可信集。这里的目标不是武断地宣布某一个 SNP 就是因果变异,而是把不确定性从“整个位点”收敛到若干有概率解释的变异集合。

关键检查:可信集大小、纯度和稳定性要和 LD 参考、样本量、区域边界一起解释。lead SNP 的 P 值最小,不等于它在多变量条件下拥有最高的因果后验概率。

第二步:用共定位判断“同一区域”是不是“同一信号”

接下来,我们把疾病 GWAS 与相关组织的 eQTL、必要时的 pQTL 放在同一位点比较。这里不能只看“GWAS 显著、eQTL 也显著”,因为两个性状可能在同一片区域分别由不同变异驱动。共定位分析关注的是:两组关联模式是否支持共享遗传信号。

对于存在多个信号的区域,我们优先采用支持多因果变异假设的流程,例如先用 SuSiE 分解信号,再通过 coloc.susie 比较对应成分。解释时同时看“共享信号”的支持与“两个不同信号”的支持,并对先验设定、组织选择和区域边界做敏感性检查。

如果某个基因只在一个方便获取、但与疾病生物学距离较远的组织中出现 eQTL 支持,我们不会把它直接升为最高优先级。相反,会把“组织是否相关”和“是否存在更贴近疾病场景的 QTL 资源”写进证据缺口。

第三步:把候选放回具体细胞与状态

单细胞数据的作用,不只是给候选基因画一张 FeaturePlot。我们会先检查细胞注释依据、患者组成、疾病状态和批次,再回答三个层次的问题:

  1. 表达定位:候选主要在哪些细胞类型或亚群中表达?
  2. 状态关联:在疾病相关状态、刺激条件或进展阶段中,表达或基因程序是否发生一致变化?
  3. 患者层级复核:这种变化能否在多个供体中重复,而不是由某一个样本或细胞数量差异驱动?

涉及组间比较时,我们把供体视为统计重复,按细胞类型进行 pseudobulk 汇总或采用适合重复测量的模型,避免把同一个患者的上千个细胞误当作上千个独立样本。若只有表达定位、没有细胞特异 QTL 或调控证据,我们会明确写成“细胞场景支持”,而不是“细胞因果证据”。

候选基因在遗传定位、QTL共定位、细胞场景、跨数据复核和实验可验证性五条证据轴上的分层矩阵
图 2|候选排序不是一个神秘的综合分数。把证据拆开后,团队能看到候选为什么优先、还缺什么、下一步怎样补证据。

第四步:把“排名”改造成“验证路线”

这个案例最终没有交付一张从 1 排到 N 的基因表,而是形成三层候选路线:

  • 优先候选:遗传精细定位、相关组织 QTL 共定位和疾病相关细胞状态三层证据方向一致,可直接进入针对性验证设计;
  • 备选候选:其中两层支持较强,但存在组织不匹配、跨队列不稳定或多信号难以区分的问题,先补一项关键证据;
  • 暂缓候选:主要依据来自物理距离、文献热度或单一数据集表达差异,目前不足以占用首轮实验资源。

对每个优先候选,我们继续写清三件事:建议使用哪类细胞模型或组织场景,干预后观察什么分子或表型读出,以及什么结果会支持或反驳当前假设。这样,生信分析的终点不再是“图画完了”,而是下一步实验可以据此排顺序。

交付物为什么不只是几张图

交付模块 回答的问题 关键技术细节 可用于什么决策
位点与可信集报告 遗传信号究竟集中在哪里? 等位基因协调、LD 匹配、独立信号、PIP 与可信集 避免只追 lead SNP 或最近基因
QTL 共定位证据表 疾病与表达/蛋白是否可能共享信号? 多信号建模、H3/H4 比较、先验和组织敏感性 区分共享调控与区域重叠
细胞场景图谱 候选可能在哪类细胞和状态中起作用? 注释复核、供体层级比较、pseudobulk、跨队列复核 选择模型与实验材料
候选分层与验证清单 先验证谁、为什么、测什么? 支持、冲突、缺口与可验证性分项呈现 安排首轮实验和补数据顺序
表 1|数据和疾病不同,具体方法会调整,但“每个结果必须服务于一个决策”这一交付原则不变。

哪些做法看起来快,实际容易把项目带偏

  • 把 lead SNP 最近的基因直接当成靶基因;
  • 只要 GWAS 和 eQTL 的显著峰重叠,就宣布共定位;
  • 在存在多个独立信号时,仍套用单一因果变异假设;
  • 只凭某个基因在一种细胞里高表达,就把它解释成疾病驱动细胞;
  • 把每个细胞当独立样本做组间检验,忽略供体层级重复;
  • 用一个综合分数遮住证据冲突,让排名无法复核。
结论边界:精细定位给出的是概率化候选集合;共定位支持共享遗传信号,不直接证明某基因介导疾病;单细胞表达定位不等于细胞特异因果。遗传效应方向也不能直接替代药物干预方向。真正的因果与功能结论仍需要适当的实验或独立研究设计验证。

这类联合分析,什么情况下最值得做

如果你已经有 GWAS 位点、候选基因或 MR 结果,但卡在“候选太多”“不知道放在哪类细胞验证”“不同证据互相冲突”,这类分析比继续堆更多相关性图更有价值。它也适合准备功能实验、申请课题、整理文章主线,或者需要判断现有公开数据能否支撑下一步研究的问题。

把一个位点发来,我们先判断证据能不能接起来

请准备:研究疾病或性状、候选位点/基因、GWAS 数据来源、已有 QTL 或单细胞数据、目标组织,以及你最想回答的实验问题。我们会先评估数据版本、祖源、样本条件和可用公开资源,再说明能连接到哪一层证据、哪些结论目前不能做。

在线沟通遗传与单细胞联合分析需求 ↗

方法与资源参考

  1. SuSiE-RSS:基于汇总统计和 LD 的遗传精细定位
  2. coloc + SuSiE:允许位点存在多个因果变异的共定位框架
  3. GTEx Portal:组织表达与 QTL 公共资源
  4. NHGRI-EBI GWAS Catalog:GWAS 汇总统计说明与获取
  5. 单细胞研究中的 pseudoreplication 与供体层级分析

把知识用于当前任务

把方法直接跑起来

使用单细胞在线分析入口,减少环境配置和重复操作。

打开单细胞工具仍有疑问,联系客服 ↗