Ka/Ks表示的是非同义替换(Ka)和同义替换(Ks)之间的比例。这个比例可以判断是否有选择压力作用于这个蛋白质编码基因。
同义/非同义突变概念
同义突变表示氨基酸没有改变,氨基酸没有改变就不会影响蛋白质的结构的与功能,这样生物活性没有受到损失,也就不会被自然选择掉。因此这些突变碱基就会保存下来。其保留下来的碱基等于其突变率。
非同义突变:表示氨基酸发生了改变,这样就影响了蛋白质的结构与功能,这样生物的生存能下降,有的则死亡,被自然选择掉。这样有些突变没有被保留下来,这样保留下来的突变碱基低于实际的突变率。
相关参数的意义
一般认为,同义突变不受自然选择,而非同义突变则受到自然选择作用。在进化分析中,了解同义突变和非同义突变发生的速率是很有意义的。常用的参数有以下几种:同义突变频率 (Ks)、非同义突变频率(Ka)、非同义突变率与同义突变率的比值(Ka/Ks)。如果Ka/Ks>1,则认为有正选择效应。如果Ka /Ks=1,则认为存在中性选择。如果Ka/Ks<1,则认为有纯化选择作用。
计算方法:
Ks = 同义突变SNP数/同义位点数, 即同义突变率
Ka = 非同义突变SNP数/非同义位点数, 即非同义突变率
同义突变SNP数= Σ同义SNP
非同义突变SNP数= Σ非同义SNP
同义位点数= Σ同义位点
非同义位点数= Σ非同义位点
Ka>>Ks或者Ka/Ks >> 1,基因受正选择(positive selection)
Ka=Ks或者Ka/Ks =1,基因中性进化(neutral evolution)
Ka<<Ks或者Ka/Ks << 1,基因受纯化选择(purify selection)
选择压力ka/ks计算软件列表:
PAML-yn00
Kaks_calculator
K-estimator
以上就是小编今天跟大家分享的内容,你对进化选择压力有什么独特见解?欢迎来评论区留言讨论!
1. 文章越来越难发?是你没发现新思路,基因家族分析发2-4分文章简单快速,学习链接:基因家族分析实操课程、基因家族文献思路解读
2. 转录组数据理解不深入?图表看不懂?点击链接学习深入解读数据结果文件,学习链接:转录组(有参)结果解读;转录组(无参)结果解读
3. 转录组数据深入挖掘技能-WGCNA,提升你的文章档次,学习链接:WGCNA-加权基因共表达网络分析
4. 转录组数据怎么挖掘?学习链接:转录组标准分析后的数据挖掘、转录组文献解读
5. 微生物16S/ITS/18S分析原理及结果解读、OTU网络图绘制、cytoscape与网络图绘制课程
6. 生物信息入门到精通必修基础课:linux系统使用、perl入门到精通、perl语言高级、R语言画图
7. 医学相关数据挖掘课程,不用做实验也能发文章:TCGA-差异基因分析、GEO芯片数据挖掘、 GEO芯片数据不同平台标准化 、GSEA富集分析课程、TCGA临床数据生存分析、TCGA-转录因子分析、TCGA-ceRNA调控网络分析
8.其他,二代测序转录组数据自主分析、NCBI数据上传、二代测序数据解读、
如果觉得我的文章对您有用,请随意打赏。你的支持将鼓励我继续创作!