系统发育树入门与实践 - SpiralKernel/bioinfo-wiki GitHub Wiki
系统发育树入门与实践
更新:2026-10-13;状态:可用;验证范围:文档核对;操作示例未在本机运行。
系统发育推断从可比较的同源序列出发,树的拓扑、枝长和支持度回答不同问题。
先判断序列是否适合比较
准备来自明确物种或样本、可比区域的序列及来源表。先理解比对相似性同源性与覆盖率:相似检索得到的命中并不自动适合放在同一棵基因树。混入旁系同源、多拷贝、片段和污染会影响解释。
多序列比对把推断同源的位置放在同一列;低质量区域、缺失和错误对齐会传递到树推断。基因树描述所分析序列的历史,不能不加判断地当作物种树。
一个学习流程
确认同源序列与标识 → 使用多序列比对工具 → 人工检查异常区域并记录处理 → 在模型下估计树 → 评估分支支持 → 根据明确外群或其他依据定根 → 连同输入比对和参数解释结果。最大似然方法在给定模型下寻找更能解释观测比对的树;搜索通常含启发式步骤。
先依照IQ-TREE官方入门下载其示例比对,运行该版本的基本建树与支持度练习,保存日志、模型选择结果、树和输入校验值。本机尚未安装运行IQ-TREE,不在此提供未经验证的跨版本命令。
图上读什么
拓扑表示分支关系;枝长通常与模型下的替换量有关,不能在没有时间校准时直接读成年份。未定根树不提供祖先方向。Bootstrap或其近似方法的支持度取决于重采样和方法,不能当作某个物种关系为真的直接概率。
比较处理前后的比对与树,检查样本ID是否一致、长枝是否来自异常序列、定根依据是否明确。单一序列区域不足以解决所有历史关系,基因重组或不同基因历史也可能造成差异。