加载中…
个人资料
  • 博客等级:
  • 博客积分:
  • 博客访问:
  • 关注人气:
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
正文 字体大小:

有效群体大小Ne

(2016-01-13 10:45:20)
分类: 分子水平的自然选择

群体遗传学中,美国遗传学家Sewall Wright在两篇标志性的论文(Wright 1931, 1938)中引入了有效群体大小(effective population size,又做有效种群大小)这一概念。他定义其为“在一个理想群体中,在随机遗传漂变影响下,能够产生相同的等位基因分布或者等量的同系繁殖的个体数量。”("the number of breeding individuals in an idealized population that would show the same amount of dispersion of allele frequencies under random genetic drift or the same amount of inbreeding as the population under consideration".)有效群体大小是很多群体遗传学模型中的基本参量,通常小于绝对的种群大小(N)。参见小种群大小。

定义

有效群体大小可能有两种定义方式,即方差有效群体大小(variance effective size)和近交有效群体大小(inbreeding effective size)。两种定义由F-统计衍生而来并紧密相连。

方差有效群体大小

在Wright-Fisher理想种群模型中,在给定上一代等位基因频率http://upload.wikimedia.org/math/e/2/e/e2ecb3ae64e76e07814ced68d05588a3.png的条件方差为:

http://upload.wikimedia.org/math/7/6/c/76c5389fb0422ac919f936ec317d6b10.png

http://upload.wikimedia.org/math/8/d/9/8d9c307cb7f3c4a32822a51922d1ceaa.png求解,可得到

http://upload.wikimedia.org/math/b/d/4/bd44f32c302040e15a56d40bfa55323f.png

近交有效群体大小

另一种定义,近交有效群体大小通过近交系数相邻两代之间的变化来计算,此时Ne定义为和理想群体的具有相同近交变化的种群大小。此处文字参照Kempthorne (1957)。

对于理想群体,近交系数依以下递归方程式计算

http://upload.wikimedia.org/math/8/2/f/82faf4c542db5a5a51bf8fced9110004.png 

利用回归交配指数(1 − F)代替近交指数得到近似递归方程式

http://upload.wikimedia.org/math/e/b/d/ebd0d4cf3ab31eac5b66d9e69c0439f7.png 

每一世代的差值是

http://upload.wikimedia.org/math/c/7/7/c7774f7fc24601260100e1c1e03cfec4.png 

通过解

http://upload.wikimedia.org/math/1/7/a/17a078928625a9f5c3fce3d860a1e9e0.png

可得出近交有效群体大小为

http://upload.wikimedia.org/math/0/5/1/0513408ed1c07ba82b43fcccda9d8484.png

,尽管研究者很少直接使用这个方程式。

举例

种群大小的变化

种群大小随时间变化。假设有t个互相不重叠的世代,那么有效种群大小由种群的调和平均数得出:

http://upload.wikimedia.org/math/1/c/b/1cba9c5e7efd0d2a582234e56052bb08.png

例如,若6个世代(t = 6)的种群大小N = 10, 100, 50, 80, 20, 500,则有效种群大小为:

http://upload.wikimedia.org/math/2/8/0/280889d5eb684857a4186429442fed87.png http://upload.wikimedia.org/math/4/6/b/46b13ded38adc5b22e46d21a2816b31d.png
http://upload.wikimedia.org/math/5/8/2/58243bbefac31896980fea5de83901a6.png
http://upload.wikimedia.org/math/4/a/a/4aa45b49b8130daeedda628c5fe56e53.png
http://upload.wikimedia.org/math/a/a/4/aa40858b3b39e82656164d7719f9b312.png http://upload.wikimedia.org/math/1/5/c/15c409acb0b556c5ae52fd034d9b8e27.png

注意调和平均数的值小于算术平均数的大小。后者在本例中为126.7。

尤其需要注意的是种群瓶颈效应,即有效群体大小受最低种群数量的严重影响。

另外一种计算家族有效种群大小Ne的方法为Ne = (4N)/(Vk 2) 其中Vk是种群大小的方差。家族大小变化增大是不利的,因为这将导致家族的近亲交配。

雌雄异体

对于雌雄异体的种群,即没有自体受精的情况,则

http://upload.wikimedia.org/math/d/3/5/d35751f5330cd7050da552c254f9be2d.png

或者更普遍地,

http://upload.wikimedia.org/math/c/a/5/ca5fb2e6bb935beb8ccaf23a31bbe536.png

其中D表示雌雄异体特性,可取0(非雌雄异体)或者1(雌雄异体)。

N很大时,Ne近似等于N,因此这通常不重要,可以忽略,即:

http://upload.wikimedia.org/math/5/c/0/5c0cfa9b79aa81706023f02aaeda7ea1.png

非费雪式性别比

当一个种群性别比不符合费雪式(Fisherian)的1:1的比值,有效种群大小为:

http://upload.wikimedia.org/math/d/d/b/ddb4e291b6a32810f66c165deb3b19b0.png

其中Nm是雄性数量而Nf为雌性数量。例如,若有100个个体的群体中有80个雄性和20个雌性,则

http://upload.wikimedia.org/math/a/a/4/aa40858b3b39e82656164d7719f9b312.png http://upload.wikimedia.org/math/c/1/1/c11f97b7dd017c321e2b7b0d3d7833d8.png
http://upload.wikimedia.org/math/8/6/4/8645cb09f2800ff4b3a3018b9225c9a5.png
http://upload.wikimedia.org/math/b/e/0/be0e1766dc69c620c17edf8871535b21.png

同样,这里有效种群大小Ne小于绝对种群大小N

对下一代的不均等贡献

如果种群大小保持恒定,每个个体必须平均向下一代贡献两个配子。一个理想种群假设这遵循泊松分布(Poisson distribution),此时每个个体贡献配子数k的方差也等于贡献配子的平均数,即2:

http://upload.wikimedia.org/math/d/c/0/dc0100b5cc28afd668bf80afda1f664e.png

然而,在自然种群中,方差通常大于此:

http://upload.wikimedia.org/math/8/1/c/81c2986db311c2469b831d09fbf002bd.png

则有效种群数量为

http://upload.wikimedia.org/math/4/a/b/4abbe244a16456e7c382e689441b965a.png

注意如果k的方差小于2,则Ne大于N。生育能力造成的遗传偏差通常造成Ne值下降。

世代重叠及种群和年龄结构的关系

当生物的生存时间跨越多于一个生育季节,有效种群大小需要考虑进物种的生命表(life table)。

单倍体

假设单倍体种群具有离散的年龄结构。例如一种生物可以再若干个离散的繁殖季节中存活。此外,定义如下年龄结构特征:

http://upload.wikimedia.org/math/8/6/5/865c0c0b4ab0e063e5caa3387c1a8741.png的费雪氏繁殖值(Fisher's reproductive value),
http://upload.wikimedia.org/math/8/6/5/865c0c0b4ab0e063e5caa3387c1a8741.png的几率,以及
http://upload.wikimedia.org/math/d/3/7/d379e006b6927ba107781f0273eaf82a.png 每个繁殖季节的新生个体数。

世代间隔可计算为生于个体的平均年龄:

http://upload.wikimedia.org/math/9/c/a/9ca2834f7e7e9e78242b41cc473add17.png

那么,近交有效种群大小为

http://upload.wikimedia.org/math/6/c/e/6ce842898c49f02308f674bb46e7a393.png (Felsenstein 1971)

二倍体

类似地,对于具有离散年龄结构的二倍体种群的近交有效种群大小也可以计算。这首先由Johnson(1977)提出,但本文的符号更接近Emigh和Pollak(1979)。

假定生命表的基本参量和上述单倍体情况相同,但区分雌雄性,如用N0ƒ N0m分别表示新生的雌性和雄性个体数量(注意小写ƒ代表雌性而大写F表示近交)。则近交有效种群大小为

http://upload.wikimedia.org/math/b/2/a/b2a0c17f9d2a4368675534a94b24daeb.png

参见

  • Minimum viable population

参考文献

  • Emigh TH, Pollak E. Fixation probabilities and effective population numbers in diploid populations with overlapping generations. Theoretical Population Biology. 1979, 15: 86–107.
  • Felsenstein J. Inbreeding and variance effective numbers in populations with overlapping generations. Genetics. 1971, 68: 581–597.
  • Johnson DL. Inbreeding in populations with overlapping generations. Genetics. 1977, 87: 581–591.
  • Kempthorne O. An Introduction to Genetic Statistics. Iowa State University Press. 1957, [1969].
  • Wright S. Evolution in Mendelian populations. Genetics. 1931, 16: 97–159. reprint
  • Wright S. Size of population and breeding structure in relation to evolution. Science. 1938, 87: 430–431.

外部链接

0

阅读 收藏 喜欢 打印举报/Report
前一篇:减数分裂
后一篇:新年目标
  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有