变异注释与功能解释学习路线 - SpiralKernel/bioinfo-wiki-migrated GitHub Wiki

变异注释与功能解释学习路线

更新:2026-10-14;状态:可用;验证范围:人工基因案例与实际工具输出核对;具体已测范围见正文。

把已经得到的变异转换为可审查的基因、转录本、后果和候选证据表。

学习顺序

  1. 先读GFF3GTF与基因转录本外显子、正负链坐标与转录方向和CDS阅读框密码子与剪接位点,画出一个带两条转录本的基因。
  2. 用bedtools区间操作与坐标验证理解重叠、最近距离和一对多关系;不要直接把最近基因当因果基因。
  3. 读变异后果标签与影响等级、转录本选择与代表注释,再按人工基因贯通区间注释与后果提取复跑人工真值。
  4. 真数据先做GFFGTF版本核查与注释准备,根据物种和资源选择SnpEff自建数据库与后果解析或VEP缓存离线注释与CSQ。
  5. 按按基因转录本与后果提取变异保留ALT对应关系;输出进入变异注释与候选证据整理和基因组浏览器检查与候选证据。
  6. 最后运行最小注释工作流与Snakemake恢复,核对无变化重跑与缺失输出恢复。

学完应能交付什么

一份表至少有参考版本、位置、REF、ALT、基因ID、转录本ID、后果、注释版本、筛选理由和原记录标识。多转录本和多ALT应能够回溯,不压成无法解释的一行标签。

同一位点对两条转录本可能分别是外显子与内含子;负链基因的cDNA替换方向也可能与VCF不同。能解释这两个反例,比记住某软件默认参数更重要。

适用边界

本轮实际跑通人工基因SnpEff/bedtools和Snakemake;VEP、参考版本转换及特殊倍性专业流程是选型指南,未实测。注释后果不是功能实验。GWAS显著位点进入候选解释时还需LD、调控和独立证据。