生物进化中序列正向选择位点筛选方法的模拟研究

目的探寻生物进化中序列正向选择位点筛选的最佳方法。方法以流感病毒H3亚型HA1序列的进化树分析结果作为产生模拟序列的进化模式,按照不同的进化参数模拟序列数为8、16、32、64、128的样本,分别采用计数模型、固定效应似然比检验模型及随机效应模型筛选模拟序列的正向选择位点。结果随着分析序列数的增多、进化强度的增强,三种方法筛选正向选择位点的灵敏度均增高;同序列数、同进化强度的筛选中,计数模型、随机效应模型的灵敏度显著低于固定效应似然比检验模型;三种模型的运算速度由快至慢分别为计数模型、固定效应似然比检验模型、随机效应模型。结论对于高进化速度、大序列样本的流感病毒正向选择位点筛选,固定效应模型比其他两类方法更加适宜。