选择的群体遗传学原理.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,西北农林科技大学动科学院,*,Click to edit Master title style,Click to edit Master text styles,Second level,Third level,Fourth level,Fifth level,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第二章选择的群体遗传学原理,二、平衡群体各基因型频率变化,三、外来因素的干扰,从群体遗传学的角度，,Halden J B,把导致种群发生遗传演变的各种因素归纳为系统过程和分散过程两类。,第二节 系统过程,一、概念,系统过程（,Systematic processes,）,:,导致基因频率向特定方向变化的过程叫系统过程，导致这一演变过程的力量，称为系统压力,（,Systematic pressure,）。,导致系统过程的原因有迁移、突变和选择。系统过程是能够预见基因频率改变方向和数量的遗传改变，在任何规模的群体中都发挥作用。,在两性个体数相等，而且两性配子完全随机结合的标准条件下，如果群体规模为N，常染色体上某个特定基因频率为q，那么，作为漂变导致的一代间变化之数量标志的基因频率方差为：,用概率论的语言来陈述，这处在一代之间特定基因的个数由i转变为j，基因频率相应由q转变为q的情况，可称为“从状态i转变为状态j”，前述关于Pr（j）的公式，也可称为由状态i到状态j的“转换概率”，并代以Pij。,又，根据标准积分 ，可知等式左边为,突变（mutation）包括基因突变和染色体突变。,pi第i亚群中的基因A的频率；,这些分析说明，对于遗传资源保护和利用的实践来说，比之于期待或诱发新突变，倒不如收集和评价种群在长期演化过程中积累起来的既有突变可能更有实际意义。,适应度(adaptive value)是比较群体中各种基因型（以个体平均留种子女数为标准）生存适应力的相对指标。,如果突变率u=10-5，在没有其它进化压力的前提下，一个频率为0.,因而，下一代基因频率偏离亲代一定范围（的概率为：,因此，在现存的产业体制下随机保存有利的非频发突变是不可能的。,所以，新基因a在其出现后第一代消失的概率为,=u(1-u)u+(1-u)2u+(1-u)3u+(1-u)nq0,除非流失或固定，下一代将以第一代漂变的结果为起点继续漂变，与第一代的区别仅仅是公式（1）中的q值。,如果它在群中留种数为K，那么，突变基因在下一代消失的概率为 。,分群与漂变机制的比较,换言之，遗传漂变乃是两代间基因频率由于机误而发生的变迁，或者说，基因频率由于群体规模的限制而与理论上应有的频率不相一致的现象。,证明，其分布近似于普哇松分布（poisson distribution）。,这一过程是由于基因频率在世代交替过程中随机抽样误差，即遗传漂变（genetic drift）所造成。,在漂变的作用下，下一代该基因在群体中的个数就有2N+1种可能：0，1，2j；,二、迁移与分群,（一）迁移所导致的基因频率变化,所谓迁移实际上是从一个种群的角度来观察混群现象。（迁入和迁出）。,1.,一般规律,设：,q,0,某基因在原有群体的频率，,q,m,同一基因在迁入群体中的频率，,m,迁移率，即迁入个体占混合群的比例，,q,1,混群后下一代的基因频率。那么，显然有,q,1,=(1-m)q,0,+mq,m,=q,0,+m(q,m,-q,0,),基因频率的变化量是 ,q=q,1,-q,0,=mq,m,-q,0,所以，迁移导致的基因频率变化决定于原群体与迁入群基因频率之差和迁移率之乘积。,2.,两个特例,（,1,）迁移率稳定的连续迁移,已知：,q,1,=(1-m)q,0,+mq,m,则：,q,n,=(1-m)q,n-1,+mq,m,如果,m,值稳定，在连续迁移的情况下，群体原有基因频率的权，每代都要以（,1-m,）为比例下降，因而，,n,代后的基因频率是：,q,n,=(1-m),n,q,0,+1-(1-m),n,q,m,n,代间，基因频率的总改变量为：,如果,m,值较高，原有种群的基因会迅速被取代。这种情况在实践上是导致种群特性迅速消失的常见情形之一。,一个典型事例是级进杂交，其迁移率,m=0.5,。级进杂交,3,代后，群体的基因频率为,q,3,=0.5,3,q,0,+(1-0.5,3,)q,m,=0.125q,0,+0.875q,m,3,代间基因频率总改变量为,虽然杂交决策者关注和强调的是原种群某些不利性状的逐代改进，但对任何位点而言，无论其等位基因是否有利，消失的相对速率都是一样的。,（二）分群所导致的基因频率变化,分群，指孟德尔群体的再分划。一个孟德尔群体因某些自然或社会原因划分为若干个交配受到限制的亚群，这在家畜中是常见现象。,1.,分群与漂变机制的比较,两者之间存在类似性，和漂变一样，亚群对大种群基因频率的偏差也没有确定方向；各亚群和漂变可能产生的各种衍生群一样，基因频率以大种群固有的频率为均值。,两者的区别是：亚群体基因频率的差别程度由不同的外部原因决定，不存在一致性的定量规律；亚群也不是大种群的随机样本，而是对大种群典型特征有既定偏差的一个部分。,2.,瓦隆（,Whalund,）定律,一个大种群划分为,K,个同等大小的随机交配的亚群，亚群内基因频率和基因型频率分布服从,Hardy-Weinberg,定律。,设：,q,i,第,i,亚群中基因,a,的频率（）,;,p,i,第,i,亚群中的基因,A,的频率；,大种群中,a,的频率,(),则：,K,个亚群基因频率之均数为,K,个亚群基因频率的方差为,而左边是分群后,aa,纯合子的亚群平均数，即平均隐性纯合子频率，代以,R,；,同样，代以,D,于是分群后各类基因型的频率值分别为：,这个比例关系称为瓦隆公式。显然，如果不分群，整个大种群可以随机分配，上述三种基因型的频率（相应以,D,、,H,、,R,代表）应吻合于平衡公式，两种情况造成的差值为,那么，如果亚群的规模不等，又会如何？,再设：,W,i,为各亚群规模的权，即各亚群分别在大种群中占的比例；,那么，平均基因频率应当为：,基因频率的方差则为,全群总计，,aa,纯合子频率为,其值仍然比不分群高,2,q,。说明瓦隆公式普遍适用。这个公式已从许多模拟试验中获得支持。据此，小结如下：,当种群被划分为若干个随机交配的亚群时，就种群全体而言，各种纯合子的频率将以各亚群间基因频率的方差为比例而增加，杂合子频率相应下降。这就是瓦隆定律。,所以，群体再划分虽然对基因频率和各位点上等位基因的种类都没有影响，但存在导致纯合化、减少群体杂合性的效应。,亚群,p,i,q,i,p,i,2,2p,i,q,i,q,i,2,1,0.9,0.1,0.81,0.18,0.01,2,0.8,0.2,0.64,0.32,0.04,3,0.7,0.3,0.49,0.42,0.09,4,0.5,0.5,0.25,0.5,0.25,5,0.1,0.9,0.01,0.18,0.81,种群平均,0.6,0.4,0.44,0.32,0.24,不分群,0.6,0.4,0.360.48,0.16,分群差异,0,0,+0.08,-0.16,+0.08,瓦隆定律是遗传资源多样性保护的重要依据，例如,分群的间接效应：不利于保种。,3.,亚群间的迁移,设原大群体的基因频率为 ，每个亚群每个世代与大群体的某随机样本交换比例为,m,，亚群一次交换后的基因频率由,q,变为,q,，那么,亚群间迁移导致基因频率变化，与亚群和大群体的频率离差成正比。,若每代连续迁移，导致群体一致；如无反向作用，则导致群体一致化。变化情况：,起始群体,下一代,方差,如果只有亚群间迁移，没有其他压力作用，,总群体的基因频率保持不变，总群体内部异质性标志的亚群间方差将趋于缩小，,总群体的纯合子频率下降，杂合子频率上升。,亚群间迁移的效应与分群相反,，,不利于保持纯系，有利于保种。,突变（,mutation,）包括基因突变和染色体突变。但后者导致群体遗传演化的情况一般可归纳于前者之中。纯粹由突变引起的群体演化很缓慢，但突变是一切遗传变异的根源。,（一）非频发突变,偶然在某一代发生的突变称为非频发突变。,假设某个孟德尔群体某位点原来完全由,AA,基因组成，在某个世代突变产生了一个新基因,a,。于是群体中就有了一个,Aa,杂合子，如果它有机会生殖，由于交配对象只有,AA,一种基因型，所以其后代有,AA,和,Aa,两种可能的基因型，两者的概率相等。如果它在群中留种数为,K,，那么，突变基因在下一代消失的概率为 。,而,K,，可能是,0,、,1,、,2,、,3,。据,Fisher.R.A.,证明，其分布近似于普哇松分布（,poisson distribution,）。,二、突变,如果假定群体规模不变，各家系平均留种数 ，那么，根据普哇松分布概率，留种数,K,为各情况的家系之概率为,P,k,=e,-uk,/K!=e,-2,2k/K!,因此，留种数,K,、留种数为,K,的家系之概率,P,k,，在这家系中新基因消失的概率,L,k,对应如下：,K,0,1,2,3,K,P,k,e,-2,2e,-2,2,2,/2,！,e,-2,2,3,/3,！,e,-2,2,k,/K,！,e,-2,L,k,1,1/2,1/2),2,(1/2),3,(1/2),K,所以，新基因,a,在其出现后第一代消失的概率为,在突变出现后第一代得到保留的概率为,1-0.367879=0.632121,。然后，由于同样的原因，在第一代得到保留的新基因，在第二代又有随机消失的危险，其概率是,e,-0.632121,=0.531464,。,也就是说，在群体规模不扩大的情况下，非频发突变基因中一多半将在其产生的第二代随机消失。即使第二代幸存，以后每经一个世代都有一个消失危险的关口，实可谓,“,躲得过初一，躲不过十五,”,。第,n,代后消失的概率为,I,n+1,=e-,(1-In),。,Fisher,从理论上证明过，,88.7%,的非频发突变必将在其产生后的,15,代之内随机消失。这上点后来也得到生物学验证。因此，在现存的产业体制下随机保存有利的非频发突变是不可能的。,（二）单向频发突变,单向频发突变是指某些可能重复发生并在相同内外环境中有稳定频率，而且不存在回原（,back mutation,）的突变。,设：突变方向是,A a,，突变率为,；,n,代表世代；基因,A,和,a,在第,n,代的频率为,p,n,和,q,n,；,那么，由于突变，下一代基因,a,的频率是,q,n+1,=q,n,+up,n,=q,n,+u(1-q,n,),，,q,n+1,=u(1-u)q,n,基于相同的理由,q,n,=u+(1-u)q,n-1,=u+(1-u)u+(1-u)q,n-2,=u+(1-u)+(1-u),2,q,n-2,=,=u(1-u)u+(1-u),2,u+(1-u),3,u+,(1-u),n,q,0,等式右边除末项外，是以（,1-u,）为公比的几何级数，其首项是,u,，项数是,n,，其和为,因此，,q,n,的一般式为,q,n,=1-,（,1-u,）,n,(1-q,0,),将,P,n,=1-q,n,与,P,0,=1-q,0,的关系代入，又可得,如果据以追溯自驯化以来的系统史，可近似地认为 ；,这样，因为,又因为突变率,u,是一个,10,-6,-10,-5,左右的极小值，当,n,很大时，近似的有,u=1/n,所以,在这种条件下，,q,n,和,n,分别是,q,n,=1-(1-q,n,)e,-nu,由于,u,值极小，种群基因频率在一定世代期间的变化将微乎其微；而种群完成一定程度的演变所必需的世代会非常之多。,例,.,如果突变率,u=10,-5,，在没有其它进化压力的前提下，一个频率为,0.005,的基因，,10,代后的频率是,0.005009,；在此条件下，基因频率从,0.1,增加到,0.15,，需要,5715.8,个世代。假设研究对象为牛，世代间隔平均为,5,年，共经历,5715.85=28579,年。,由此可见，由突变积累多样化的过程非常缓慢。如果仅从突变的角度来看，家畜群体现存的遗传多样性，是千百年演化的结果一旦丧失，相对于现代社会的保种需要而言，恢复是不可能的。另外，上述事实也说明，为什么人类选择所未涉及，又与适应关系不大的那些位点，以及品种演化过程中，其基因频率变化极其缓慢以至可用以标志种群起源的原因。,（三）可逆的频发突变,指存在回原的频发性突变，这种突变导致的种群演变就更慢。,1,、突变平衡值,设：某位点基因,A,和,a,的频率分别是,p,和,q,；,u,为,A a,的突变率；,v,为,a A,的回原率；,在突变体和回原导致的基因频率变化相等，亦即平衡时，,A,和,a,的频率分别为 。,则：在达到平衡时，,以上分析说明：,（,1,）在一个群体，突变和回原的平衡值纯粹是由突变率和回原率的相对大小决定，而与起始的基因频率无关；,（,2,）在任何世代，如果 ，它就要逐代减小，直到,为止，如果 ，它就要逐代升高，直到 为止；,（,3,）如果由于其它进化压力，造成了,p,对 的偏差，一旦其作用消除，群体的基因频率仍然要向平衡值回复。,2,、基因频率走向平衡的速度,因为基因,a,在一代间的频率改变量为,q=up-vq=u(1-q)-vq=u-uq-vq,=(u+v),因而有,q=-(u+v),这就是说，基因频率向平衡值回复的速率与其实际值对平衡值的偏差成正比。基因频率距率平衡值越远，一代之间向平衡值靠近的幅度越大。此外，该速率还与突变率与回原率的算术和成正比。,据比，可以分析基因频率在向平衡值回复的过程中。造成一个特定变化量所需的世代。,由于,u,和,v,数值极小，,q,更小，所以，可以把,q=-,（,u+v,）作为一个微分方程，并设,t,代表以世代为单位的时间，即,两边在,n,代期间求积,又，根据标准积分 ，可知等式左边为,例如：时，由,q,0,=0.10,提高到,q,n,=0.15,以及由,q,n,=0.50,提高,q,n,=0.55,，所需要的代数分别是：,基因频率向平衡值回复是决定种群稳定性的机制之一。然而，既便是有这种变化，也极其缓慢。,根据,q=up-vq,也可以分析一定世代之后的基因频率。,如果以,q,n,代表某个世代的基因频率，下一代的基因频率则为：,q,n+1,=q,n,+u(1-q,n,)-vq,n,=u+(1-u-v)q,n,将此作为相邻两代基因频率的一般关系，反复代入,q,n,、,q,n-1,、,q,n-2,q,0,，可得,q,n,=u+u(1-u-v)+u(1-u-v),2,+u(1-u-v),n,q,0,右边除最后一项外，是以,(1-u-v),为公比的几何级数，首项为,u,，项数是,n,。因而，其和为,又，当所历经世代极多，以至,时，,例如：当,u=210,-5,，,v=10,-5,，,q,0,=0.3,，,n=15,代时，,q,n,=0.300165,。也就是说,15,代期间突变导致的基因频率变化仅为,1.6510,-4,。,这些分析说明，对于遗传资源保护和利用的实践来说，比之于期待或诱发新突变，,倒不,如收集和评价种群在长期演化过程中积累起来的既有突变可能更有实际意义。因此需要注意：,（,1,）没有经过选择的位点，其基因频率变化非常微小，常常被忽视；,（,2,）大数法则非常重要，不能只靠某些层次或少数位点进行群体分类，目前流行的用部分位点分类，只能是一种近似估计。,（,3,）抽样要有代表性，往往抽样造成的差异大于实际差异。,突变和迁移虽然有方向性，但它们的方向对于适应或育种目的而言，是随机化的。在种群遗传演变过程中，选择是唯一的导向因素。,（一）适应度和选择系数,适应度,(adaptive value),是比较群体中各种基因型（以个体平均留种子女数为标准）生存适应力的相对指标。以最佳基因型的个体平均留种子女数为,1,，特定基因型 的平均留种子女数所相当的比例，就是该基因型 的适应度，通常以,w,表示。选择系数，通常以,S,表示，,S=1-W,。例 如,:,三、选择,可见，适应度是特定基因型的留种率和群体最佳基因型之比值。,（二）关于单个位点选择进度的归纳,假设：群体在某一位点存在,A,1,A,2,两个等位基因，其频率分别为,p,和,q,；基因型,A,1,A,1,A,1,A,2,A,2,A,2,的适应度分别为,w,1,w,2,w,3,选择系数相应为,s,1,s,2,s,3,；则群体平均适应度为,w=p,2,w,1,+2pqw,2,+q,2,w,3,=1-s,1,p,2,-2s,2,pq+-s,3,q,2,。,基因型,AA,Aa,aa,亲本个体数,25,50,25,留种子女数,70,70,14,个体平均留种子女数,2.8,1.4,0.56,适应度,1,1.4/2.8=0.5,0.56/2.8=0.2,根据适应度和选择系数的定义，三种基因型的适应度中必然有一个为,1,，相应的选择系数为,0,，为了分析一般规律在此不代入这两个数值。相邻两代基因频率的变化量为,这可以作为单个位点选择时，相邻两代基因频率改变量的一般公式。,S,1,、,S,2,、,S,3,的具体情况决定着不同的选择模式和不同选择进度。,（三）关于多位点选择进度的归纳,单个位点选择进度的归纳尚未包含有关选择影响种群遗传演变的全部方面。虽然种群的遗传演变归根结底是各有关位点基因频率的变化，但作为选择单位的不是基因、染色体或配子，而是个体。以下几个方面是单个位点以外的有关选择的重要问题。,1.,在特定环境背景中，在某个或某几个位点存在有利基因型的个体，在一次选择中得以保留，它在其他位点不利基因也随之保留；反之，具有某种突出不利特征个体的淘汰，其携带的有利基因也随之淘汰。位点间这种效应的掩盖与同一位点上隐性基因的潜存机制相似。,2.,如果存在连锁，上述现象就不仅发生在一次选择过程，而成为经常性的事件。这是人工选择中常见难题。值得注意的是，有利基因的连锁材料，是更珍贵的遗传资源。,3.,位点间的互作效应，虽然不一定稳定，但作为遗传共适应的一个方面，在一定条件下，和基因资源一样是种质资源的组成部分。在人工选择中将有利互作的各位点的基因在群体中固定，是遗传资源开拓的内容之一。,一、概念,分散过程（,Systematic processes,）,:,导致基因频率发生无方向变化的过程叫分散过程。,这一过程是由于基因频率在世代交替过程中随机抽样误差，即遗传漂变（,genetic drift,）所造成。从保种的实践角度来看，分散过程是指能预见基因频率的变化数量，但不能确知其变化方向的过程。主要环节有：受精、精卵结合、着床、发育、留种，均是自由随机的，可以看作随机抽样。,二、遗传漂变的性质与度量,第三节 分散过程,1.,遗传漂变的概念,二倍体生物在世代交替过程中，亲代两性配子的随机结合以及合子的随机成活，使遗传过程经历了两次抽样。就此，可以把一个家畜群体视为来自代配子群的一个随机样本。在这一过程中，无论是否存在外来干扰，如果群体规模较小，下一代的实际基因频率都可能由于抽样误差而偏离理论上应有的频率。这种现象称为遗传漂变。换言之，遗传漂变乃是两代间基因频率由于机误而发生的变迁，或者说，基因频率由于群体规模的限制而与理论上应有的频率不相一致的现象。,2.,遗传漂变的度量,1931,年,Wright,，,S.,根据二项分布概率论证过基因频率的随机抽样方差。因为在一个位点，可以把所有的等位基因按是否某个特定基因分为两种情况，即将特定基因频率和其它所有等位基因的累加频率分别作为二项分布的基因概率来进行分析。,遗传漂变没有确定的方向，但可以用基因频率的方差来预见其一代之间的变化程度。在两性个体数相等，而且两性配子完全随机结合的标准条件下，如果群体规模为,N,，常染色体上某个特定基因频率为,q,，那么，作为漂变导致的一代间变化之数量标志的基因频率方差为：,如果其它前提不变，而该基因 性染色体上，其方差则为,不言而喻，规模越小，基因频率的方差越大，一代间由漂变引起的基因频率改变量也越大；在基因频率可能的变化范围,0,1,之间，其值为,0.5,时，其方差最大，其值越趋近两个极端值（,0,或,1,），方差也越小。,三、出现定值改变的概率,1,、下一代基因频率转变为定值的概率,如果在一个规模为,N,的小群体中特定基因的频率为,q,，那么这个基因在群体中的个数为,i=2Nq,。在漂变的作用下，下一代该基因在群体中的个数就有,2N+1,种可能：,0,，,1,，,2,j,；,2N,；为,0,时，意味着流失，为,2N,时，意味着固定。每种可能情况的概率值，可以用二项式,（,1-q,）,+q,2N,的展开式来求出。特定基因在下一代的群体中的个数为,j,的概率，正是展开式的第,j,项：（,1,）,例如：当,N=10,，,q=0.5,时，由于漂变，下一代频率转变为,0,，,0.2,，,0.4,，,1,，以有保持,0.5,不变的概率如下所示：,频率转变,j,Pr,0.50,0,9.53710,-7,0.50.2,4,4.62110,-7,0.50.4,8,0.120,0.50.5,10,0.176,0.51.0,20,9.53710,-7,用概率论的语言来陈述，这处在一代之间特定基因的个数由,i,转变为,j,，基因频率相应由,q,转变为,q,的情况，可称为,“,从状态,i,转变为状态,j,”,，前述关于,Pr,（,j,）的公式，也可称为由状态,i,到状态,j,的,“,转换概率,”,，并代以,Pij,。,2,、下一代基因频率转变在特定范围内的概率,因为群体中一代配子群可能形成,2N+1,种样本，样本间的基因频率分布近似于以亲代原有频率为中值的正态分布，因而，一代之间基因频率转变在特定数值范围内的概率，可由正态分布概率公式来估计。,如以,q,代表亲代基因频率，,q,1,为子代基因频率，则基因频率的标准偏差为：,例如：在,N=10,，,q=0.5,的群体，由于漂变，下一代基因频率偏离,0.1,和,0.3,以内（即分别在,0.4-0.6,和,0.2-0.8,范围内）的概率分别是,0.6281,和,0.9926,；下一代基因频率落在,0.2-0.4,、,0.2-0.5,范围内的概率分别为,0.1823,和,0.4963,。,因而，下一代基因频率偏离亲代一定范围（的概率为：,下一代基因频率落在高或低于原频率（单侧）一定范围内的概率则为：,四、漂变的总体影响,个体数为,N,的群体，常染色全基因一代之间的漂变有,2N+1,种可能的结果；除非流失或固定，下一代将以第一代漂变的结果为起点继续漂变，与第一代的区别仅仅是公式（,1,）中的,q,值。也就是说，第二代的漂变，本质上是以第一代漂变造成的基因频率为均值的另一次随机抽样，而与以前的基因频率无关。以后各代的情况也类似于此。这种现象在统计学上称为,“,马尔柯夫,(Markov),过程,”,或,“,无后效过程,”,(,所谓,“,马尔柯夫过程,”,是俄国人,JI.JI.Markov,首先论证的一种随机过程。其特点是,“,当现在情况已经确定时,以后的一切统计特性就跟过去的情况无关。,”,),。连续漂变的总结果，可据以预计。,设：,N,群体规模；,Pi,j,由于漂变，群体中特定基因的个数从,i,转变为,j,的概率，亦即由状态,i,到,j,的转换概率，,i=0,，,1,，,2,，,2N,，,j=0,，,1,，,2,，,2N,。,t,世代；,f,t(q),第,t,代基因频率为,q,的概率；,当然，,q=i/2N,，或,q=j/2N,。,那么，当,N,为定值时，任何确定,i,和,j,之间的转换概率就是定值。例如，当,N=2,时，,因此，在,N,已确定时，可以用各种状态（以群体包含有特定基因的个数,i,标志）之间的转换概率值为元素建立一个矩阵,P,：,其次，如前所述，,N,为定值时，作为各世代漂变起点的基因频率值,i/2N,有,2N+1,种可能的情况。例如当,N=2,时，,5,种可能的情况是,q=0,0.25,0.50,0.75,1,。在各世代，各种情况的可能性亦即出现概率,ft(q),并不相等。以此,2N+1,个概率值可构成一个列向量：,显然，作为相邻两代特定基因频率间转换可能性参数的矩阵,P,与作为遗传漂变起点的列向量,f,t,之乘积，就是下一代基因频率各种可能数值的概率,f,t+1,，即,f,t+1,=Pf,t,