2019-07-06

Ka/Ks

在遗传学中,Ka/Ks或者dN/dS表示的是异意替换(Ka)和同意替换(Ks)之间的比例。这个比例可以判断是否有选择压力作用于这个蛋白质编码基因。
  不导致氨基酸改变的核苷酸变异我们称为同义突变,反之则称为非同义突变。一般认为,同义突变不受自然选择,而非同义突变则受到自然选择作用。在进化分析中,了解同义突变和非同义突变发生的速率是很有意义的。常用的参数有以下几种:同义突变频率(Ks)、非同义突变频率(Ka)、非同义突变率与同义突变率的比值(Ka/Ks)。如果Ka/Ks>1,则认为有正选择效应。如果Ka/Ks=1,则认为存在中性选择。如果Ka/Ks<1,则认为有纯化选择作用。

Ks = 同义突变SNP数/同义位点数
即同义突变率
Ka = 非同义突变SNP数/非同义位点数
即非同义突变率
同义突变SNP数= Σ同义SNP
非同义突变SNP数= Σ非同义SNP
同义位点数= Σ同义位点
非同义位点数= Σ非同义位点

uKa>>Ks或者Ka/Ks >> 1,基因受正选择(positive selection)
uKa=Ks或者Ka/Ks =1,基因中性进化(neutral evolution)
uKa<<Ks或者Ka/Ks << 1,基因受纯化选择(purify selection)

检测序列的功能性(funcional or pseudo)
筛选正在快速进化的基因(rapid evolution)
Ks可以反映事件发生的时间(age)

分子进化领域常用软件
系统进化树构建软件列表:
Phylip
Clustalw
PAML-Codml
其他
选择压力ka/ks计算软件列表:
PAML-yn00
Kaks_calculator
K-estimator
其他
snp搜索软件列表:
polyphred
SNPdetector
BGI-Variation analysis

非同义替换率(氨基酸改变,dn)与同义替换率的(氨基酸不改变,ds)的比值(dn/ds)也经常被用于分化分析。dn/ds的比值为1表示所研究的基因在中性选择(neutral selection)下进化,小于0. 25意味着纯化选择(purifying selection)下进化,当比值大于1时则被认为进行正向选择(positive selection)下的进化(Hurst et al, 2002; Swanson et al 2003)。

对于研究蛋白编码序列突变的一种简单而有效的分类方法是将替换分成同义
替换(Synonymous substitution) 和错义替换(Non-synonymous substitution)。同义替换是指那些可以引
起所编码的氨基酸发生不改变的替换,一般认为这样的替换不会受到选择的压力或者受到的选择作用
比较小;非同义替换是指那些可以改变所编码的氨基酸的替换,这样的替换有时候会导致新的功能。
依据密码子的简并性(degeneracy)可以将核苷酸位点分成两类:同义替换位点(Synonymous site)和错
义替换位点(Non-synonymous site)。同义替换率和错义替换率定义为:每代或者每年在每个可能的同
义(错义)位点上实际发生的同义(错义)替换数目(rS, rN)。然而由于对序列分化的时间不能确
定,因此同义替换率(Synonymous substitution, Ks 或 dS)和错义替换率(Non-synonymous substitutionrate, Ka 或dN)可以定义为:在两序列分化至今的 t 年里每个可能的同义(错义)位点上实际发生
的同义(错义)替换数目。因此有:Ka = 2rNt,Ks = 2rSt 。通过比较错义替换率和同义替换率的相对
比值可以确定这个基因在进化中受到的选择压力。Ka, Ks 的之间的比值是已经为人们所接受和广泛
应用的表现进化动力的指标。

Similar to dn / ds ratios,
the rate of accumulation of non-synonymous polymorphism
(pN) scaled by the rate of synonymous polymorphism
(pS) provides a glimpse on the selective forces
driving the evolution of a protein-coding sequence.
Thus, genes with a high pN / pS (i.e. >1) ratio are likely
to be evolving under the influence of positive selection

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 203,271评论 5 476
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,275评论 2 380
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,151评论 0 336
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,550评论 1 273
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,553评论 5 365
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,559评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,924评论 3 395
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,580评论 0 257
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,826评论 1 297
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,578评论 2 320
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,661评论 1 329
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,363评论 4 318
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,940评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,926评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,156评论 1 259
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 42,872评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,391评论 2 342

推荐阅读更多精彩内容