用途
GWAS SNP Interpretation Skill 用于解读 GWAS（全基因组关联分析）研究中的遗传变异（SNP），通过汇总来自多个来源的证据，提供全面的临床与生物学背景。典型使用场景包括：解读 rs7903146（TCF7L2 糖尿病变异）、查询 rs429358（APOE 阿尔茨海默病变异）关联的疾病、了解 rs1801133（MTHFR 变异）的临床意义，以及判断 rs12913832（眼睛颜色变异）是否位于任何精细定位位点中。
主要能力
该技能通过以下方式提供 SNP 的全面解读：
- SNP 注释：获取变异的基本信息，包括基因组坐标、等位基因、功能后果和映射基因。
- 关联发现：查找所有具有统计显著性的 GWAS 性状/疾病关联。
- 精细定位证据：识别该变异所属的可信集（精细定位的因果位点）。
- 基因映射：使用 Locus-to-Gene（L2G）预测识别可能的因果基因。
- 临床摘要：将证据汇总为可操作的临床意义。
工作流程为：用户输入（如 rs7903146）→ SNP 查找（gwas_get_snp_by_id，获取位置、后果、MAF）→ 关联搜索（gwas_get_associations_for_snp，查找所有性状/疾病关联）→ 精细定位（可选，OpenTargets_get_variant_credible_sets，获取可信集成员关系）→ 基因预测（提取因果基因的 L2G 分数，内嵌于可信集）→ 临床摘要（汇总证据，识别关键性状与基因）→ 输出：综合解读报告。
数据来源
- GWAS Catalog（EMBL-EBI）：提供 SNP 注释（功能后果、映射基因、群体频率）与关联（P 值、效应量、研究元数据），覆盖 350,000+ 篇文献、670,000+ 条关联。
- Open Targets Genetics：提供精细定位（来自 SuSiE、FINEMAP 方法的统计可信集）、L2G 预测（基于机器学习的基因优先级排序）、共定位（因果基因的 QTL 证据），覆盖 UK Biobank、FinnGen 等大型队列。
输入参数
必填：
- rs_id（str）：dbSNP rs 标识符，格式为 "rs" + 数字（如 "rs7903146"），必须是 GWAS Catalog 中的有效 rsID。
可选：
- include_credible_sets（bool，默认 True）：是否查询精细定位数据。True 表示完整解读（较慢，约 10-30 秒）；False 表示仅快速返回关联（约 2-5 秒）。
- p_threshold（float，默认 5e-8）：全基因组显著性阈值。
- max_associations（int，默认 100）：最多检索的关联数量。
输出格式
返回 SNPInterpretationReport，包含：
- SNP 基本信息：如 rs_id、染色体、位置、参考等位基因、替代等位基因、功能后果、映射基因、MAF。
- 性状关联：每条记录包含性状、P 值、beta（效应量）、study_id、pubmed_id、效应等位基因。
- 可信集（精细定位）：包含 study_id、性状、精细定位方法（如 SuSiE-inf）、P 值、预测基因（含 L2G 分数）、区域。
- 临床意义：全基因组显著关联的性状/疾病列表、被识别的精细定位位点数量、预测的因果基因。
使用方式
技能通过调用以下工具完成：
- GWAS Catalog 工具：gwas_get_snp_by_id（获取 SNP 注释）、gwas_get_associations_for_snp（获取所有性状关联）。
- Open Targets 工具：OpenTargets_get_variant_info（获取变异详情与群体频率）、OpenTargets_get_variant_credible_sets（获取精细定位可信集与 L2G）。
平台特定的示例参见 QUICK_START.md（资料未提供其内容）。
解读指南
P 值显著性水平：
- p < 5e-8：全基因组显著（强证据）
- p < 5e-6：提示性（中等证据）
- p < 0.05：名义（弱证据）
L2G 分数解读：
-
0.5：高置信度因果基因
- 0.1-0.5：中等置信度
- < 0.1：低置信度
临床可操作性：
- 高：多个全基因组显著关联 + 位于可信集中 + 高 L2G 分数
- 中：有全基因组显著关联但精细定位有限
- 低：提示性关联或重复证据有限
适用场景
适用于需要综合多个数据库证据解读特定 SNP 的临床或生物学背景、查找性状/疾病关联、进行精细定位与因果基因预测的场景，例如解读 TCF7L2、APOE、MTHFR 等变异。
限制
- 变异 ID 转换：OpenTargets 需要 chr_pos_ref_alt 格式，可能需要查找等位基因。
- 群体特异性：关联可能因祖源而异。
- 效应量：beta 值依赖于研究（不同表型尺度不同）。
- 因果性：关联不能证明因果关系；精细定位可提高置信度。
- 时效性：数据反映已发表的 GWAS，最新研究可能未包含。
最佳实践
- 使用完整解读：做临床决策时启用 include_credible_sets=True。
- 检查多个变异：查看同一位点上的其他变异。
- 验证群体：考虑祖源特异性效应量。
- 审阅文献：查看原始研究以了解背景。
- 整合证据：结合功能数据、eQTLs、pQTLs。
技术说明
性能：快速模式（无可信集）约 2-5 秒；完整模式（含可信集）约 10-30 秒；瓶颈为 OpenTargets GraphQL API 速率限制。
错误处理：无效 rs_id 返回错误消息；无关联返回空列表并附说明；API 失败时优雅降级（返回部分结果）。
相关技能：基因功能分析、疾病本体查询、PubMed 文献检索、变异效应预测。
版本：1.0.0，最后更新 2026-02-13；要求 ToolUniverse 版本 >= 1.0.0；需要工具：gwas_get_snp_by_id、gwas_get_associations_for_snp、OpenTargets_get_variant_credible_sets。