《群体遗传数学模型PPT课件.ppt》由会员分享，可在线阅读，更多相关《群体遗传数学模型PPT课件.ppt（35页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、关于群体遗传的数学模型第一张，PPT共三十五页，创作于2022年6月n n一.问题与背景n n 1.问题n n 生物的细微变异在进化中很重要，可以为物种进化提供原始材料。n n但变异的个体在群体中数量很少。n n它会不会在群体的随机交配的过程中逐渐减弱直至消失？n n 人类的遗传病（如色盲）在群体中会不会由于一代一代地遗传而使患者越来越多？n n 第二张，PPT共三十五页，创作于2022年6月n n 农作物的新品种主要是用杂交的方法培育的。n n即选择某种作物具有互补性状的两个品种杂交，在它们的后代中选择培育出合乎要求的新品种。n n由于杂交后代性状分离的不稳定性，使得它不可能杂交后直接用于农
2、业生产。n n多少代以后就可以达到生产所需要的稳定性呢？第三张，PPT共三十五页，创作于2022年6月n n 2.背景：孟德尔遗传学的要点n n 1.生物体的性状是由一对（多对）基因控制的。n n染色体是基因的载体。基因在染色体上的位置称为基因座，处于同一基因座的不同的基因为等位基因。n n等位基因对性状的控制有显性和隐性之分。第四张，PPT共三十五页，创作于2022年6月n n 2.处于同一基因座上的两个基因的组合形式称为基因型。n n生物体表现出来的特征为表现型。n n一般表现型反映显性基因所控制的性状。n n 3.遗传的过程是由亲代每基因座上的等位基因分离并随机进入到一个配子中去n n在
3、交配过程中雌雄配子将随机结合形成新的合子而传入下一代。第五张，PPT共三十五页，创作于2022年6月n n二.假设n n 1.雌性和雄性在遗传上是对等的。n n 2.没有迁移、突变和选择。n n 3.群体足够大。n n 4.随机交配，控制有关性状的雌雄配子完全随机结合。第六张，PPT共三十五页，创作于2022年6月 n n三.模型n n 1.群体的遗传结构：一个基因座，两个等位基因 A，a。三个基因型 AA，Aa，aa。n n 基因型频率：f1，f2，f3，(f1+f2+f3=1).n n 基因频率:p=f1+f2/2,q=f3+f2/2.n n 群体的遗传结构:F=(f1，f2，f3)T.n
4、 n 基本群体:n n F1=(1,0,0)T,F2=(0,1,0)T,F3=(0,0,1)T.n n F=f1F1+f2F2+f3F3.第七张，PPT共三十五页，创作于2022年6月 n n 2.2.遗传过程的基本平衡关系遗传过程的基本平衡关系n n 由孟德尔遗传学的理论可知由孟德尔遗传学的理论可知,有有n n AAAAAA=AA,AAAA=AA,AAaa=Aa,aaaa=Aa,aaaa=aa,aa=aa,n n AA AAAa=(1/2)AA+(1/2)Aa,Aa=(1/2)AA+(1/2)Aa,n n aa aaAa=(1/2)Aa+(1/2)aa,Aa=(1/2)Aa+(1/2)aa,
5、n n Aa AaAa=Aa=(1/4)AA+(1/2)Aa+(1/4)aa(1/4)AA+(1/2)Aa+(1/4)aan n 即即,有如下的平衡关系有如下的平衡关系n n F F1 1F F1 1=F F1 1,F F1 1F F3 3 =F F2 2,F F3 3F F3 3=F F3 3,n n F F1 1 F F2 2=(1/2)=(1/2)F F1 1+(1/2)+(1/2)F F2 2,n n F F3 3 F F2 2=(1/2)=(1/2)F F2 2+(1/2)+(1/2)F F3 3,n n F F2 2 F F2 2=(1/4)(1/4)F F1 1+(1/2)+(1
6、/2)F F2 2+(1/4)+(1/4)F F3 3.第八张，PPT共三十五页，创作于2022年6月n n 3.群体遗传的模型:n n F(n)F(n)=(f=(f1 1(n),f(n),f2 2(n),f(n),f3 3(n)(n)T T:第 n 代群体的结构.n n 10.基本遗传过程:n n F1与F(n)交配的后代 F(n+1)F(n+1)=F F1 1F(n)F(n)=F F1 1ff1 1(n)(n)F F1 1+f+f2 2(n)(n)F F2 2+f+f3 3(n)(n)F F3 3 n n =f f1 1(n)(n)F F1 1 F F1 1+f+f2 2(n)(n)F F
7、1 1 F F2 2+f+f3 3(n)(n)F F1 1 F F3 3 n n =F F1 1 F F1 1 F F1 1 F F2 2 F F1 1 F F3 3 F(n)F(n)n n =F F1 1 (F F1 1+F F2 2)/2)/2 F F2 2 F(n)F(n)=M=M1 1 F(n)F(n)第九张，PPT共三十五页，创作于2022年6月n n F F2 2与与F(n)F(n)交配的后代交配的后代 n n F(n+1)F(n+1)=F F2 2F(n)F(n)=M=M2 2 F(n)F(n)n nF F3 3与与F(n)F(n)交配的后代交配的后代 n n F(n+1)F(n
8、+1)=F F3 3F(n)F(n)=M=M3 3 F(n)F(n)第十张，PPT共三十五页，创作于2022年6月n n 20.随机交配群体的遗传模型:一个结构为F的群体与第n世代的群体 G(n)G(n)随机交配的后代.n n G(n+1)G(n+1)=F FG(nG(n)=(f=(f1 1F F1 1+f+f2 2F F2 2+f+f3 3F F3 3)G(n)G(n)n n =f f1 1 F F1 1 G(n)G(n)+f+f2 2 F F2 2 G(n)G(n)+f+f3 3 F F3 3 G(n)G(n)n n =f=f1 1MM1 1 G(n)G(n)+f+f2 2MM2 2 G(
9、n)G(n)+f+f3 3MM3 3 G(n)G(n)n n =(f=(f1 1MM1 1+f+f2 2MM2 2+f+f3 3MM3 3)G(n)G(n)=MMr r G(n)G(n)第十一张，PPT共三十五页，创作于2022年6月n n 30.自交群体的遗传模型.自交群体指每个个体自交所产生的后代群体.n n F(n+1)F(n+1)=F(n)F(n)F(n)F(n)n n =f=f1 1(n)(n)F F F F1 1+f f2 2(n)(n)F F2 2 F F2 2+f f 3 3(n)(n)F F3 3 F F3 3nn =F F1 1 F F1 1/4+/4+F F2 2/2+/
10、2+F F3 3/4/4 F F3 3 F(n)F(n)n n =M=MS S F(n)F(n)第十二张，PPT共三十五页，创作于2022年6月n n四.模型的分析n n 1.HardyWeinberg 平衡(HW平衡)n n 在前述假设下,一个随机交配的有性繁殖的群体中，等位基因频率保持不变，基因型频率至多经过一个世代也将保持不变。n n 一个基因座，两个等位基因。n n 群体 F=(f1,f2,f3)T;n n f1,f2,f3 基因型 AA,Aa,aa 的频率;n n p=f1+f2/2,q=f3+f2/2 基因 A,a 的频率.第十三张，PPT共三十五页，创作于2022年6月 n n群
11、体自身的随机交配,有模型:n n F(n+1)F(n+1)=M=Mr r(n)(n)F(n)F(n)n n其中其中 n n p(n)=fp(n)=f1 1(n)+f(n)+f2 2(n)/2,(n)/2,n n q(n)=f q(n)=f3 3(n)+f(n)+f2 2(n)/2.(n)/2.第十四张，PPT共三十五页，创作于2022年6月n n 1.证明 p(n+1)=p(n),q(n+1)=q(n)n n由模型可以得到n nf f1 1(n+1)=p(n)f(n+1)=p(n)f1 1(n)+p(n)f(n)+p(n)f2 2(n)/2(n)/2n n =p(n)f =p(n)f1 1(n
12、)+f(n)+f2 2(n)/2=(n)/2=p(n)p(n)2 2 n nf f3 3(n+1)=q(n)f(n+1)=q(n)f3 3(n)+q(n)f(n)+q(n)f2 2(n)/2(n)/2n n =q(n)f =q(n)f3 3(n)+f(n)+f2 2(n)/2(n)/2=q(n)q(n)2 2 n nf f2 2(n+1)=(n+1)=q(n)q(n)f f1 1(n)+p(n)+q(n)f(n)+p(n)+q(n)f2 2(n)/2+p(n)f(n)/2+p(n)f3 3(n)(n)n n =q(n)fq(n)f1 1(n)+f(n)+f2 2(n)/2+p(n)f(n)/2
13、+p(n)f3 3(n)+f(n)+f2 2(n)/2(n)/2n n =2p(n)q(n)=2p(n)q(n)n n p(n+1)=fp(n+1)=f1 1(n+1)+f(n+1)+f2 2(n+1)/2=p(n),(n+1)/2=p(n),n n q(n+1)=f q(n+1)=f3 3(n+1)+f(n+1)+f2 2(n+1)/2=q(n).(n+1)/2=q(n).第十五张，PPT共三十五页，创作于2022年6月n n 2.证明基因型频率不变(至多经过一个世代)n n 矩阵 Mr 的列和为 1,因此它有全 1 的左特征向量 1=(1,1,1)=(1,1,1)T T,且有 1T Mr=
14、1T.n n 由此可知n n矩阵有特征值=1,并且只有这一个正特征根n n它所对应的右特征向量为F*0,n n (PerronFrobenious 定理)n n由 MrF*=F*不难得到 F*=(p2 2pq q2)T就是模型的平衡态。n n对于任意的初始状态F(0)F(0)=(f=(f1 1(0),f(0),f2 2(0),f(0),f3 3(0)(0)T T,容易算出 F(1)=Mr F(0)=F*.n n 因此,至多经过一个世代,群体的基因型频率将处于平衡状态不变。第十六张，PPT共三十五页，创作于2022年6月n n例 1.MN血型：n n 血红细胞膜上的特异抗原类型，分别对抗体M和抗
15、体N成阳性反应，这一性状分别由一对等位基因 LM 和 LN 控制，呈并显性。n n调查1279人三种基因型频数的比例为n n LMLM:LMLN:LNLN=363:634:282n n由此可以算出群体的遗传结构为：n n (f1 f2 f3)=(0.284，0.496，0.220)n n群体中基因 LM 和 LN 的频率分别为n n p=0.532,q=0.468n n p2=0.283,2pq=0.498,q2=0.219n n与群体的遗传结构基本相同。n n群体为HW群体。第十七张，PPT共三十五页，创作于2022年6月n n例 2.地中海贫血病，n n由一对等位基因Tm、Tn控制。n n
16、调查意大利群体10000人三种基因型的频数为 n n TmTm:TmTn:TnTn=4:400:9596n n基因 Tm 和Tn 的频率分别为n n p=0.0204，q=0.9796n n p p2 2=0.00042=0.00042，2pq=0.03997,q2pq=0.03997,q2 2=0.95961=0.95961。n n与群体的遗传结构基本相同。n n群体为HW群体第十八张，PPT共三十五页，创作于2022年6月n n问题：P175 8 民民 族族 O(%)O(%)A(%A(%）B(%)B(%)AB(%)AB(%)汉汉 30.86 30.86 31.3131.31 28.0628
17、.06 9.779.77 维吾尔维吾尔 27.5027.50 29.2229.22 31.9231.92 11.3611.36 壮壮 47.2847.28 21.2521.25 27.5727.57 3.903.90 回回 35.9435.94 27.2327.23 28.3428.34 8.498.49 哈萨克哈萨克 37.9737.97 22.8222.82 29.8329.83 9.389.38 锡伯锡伯 24.4224.42 25.0025.00 40.1240.12 10.4610.46 乌兹别克乌兹别克 25.5825.58 25.5825.58 38.7638.76 10.081
18、0.08 柯尔克孜柯尔克孜 34.6834.68 18.5418.54 39.5239.52 7.267.26 塔克尔塔克尔 21.6221.62 40.5440.54 35.1435.14 2.702.70 白白 31.4031.40 34.0034.00 23.4023.40 11.2011.20 傣傣 40.4440.44 22.0822.08 29.5929.59 7.897.89 景颇景颇 37.8137.81 34.8334.83 20.3920.39 6.976.97第十九张，PPT共三十五页，创作于2022年6月n n 2.杂交育种的稳定性n n 作物的品种：同型合子的基因型n
19、 n 育种的过程：n n双亲杂型合子 分离后代 优良组合 分离后代 优良组合第二十张，PPT共三十五页，创作于2022年6月n n问题：杂交后要自交多少世代就可以形成品种？n n F(n)F(n)=M=MS SF(n-1)F(n-1)=M=MS Sn-n-1 1F(1)F(1)=M=MS Sn nF F2 2.n n =1 =1 u u1 1=(1,0,-1)=(1,0,-1)T Tn n =1 =1 u u2 2=(0,0,1)=(0,0,1)T Tn n =0.5 =0.5 u u3 3=(1,-2,1)=(1,-2,1)T T.n n令n n则有分解式MS=U-1U,MSn=U-1nU.
20、第二十一张，PPT共三十五页，创作于2022年6月n n由此不难算出n n可得 n nF(n)F(n)=(0.5-0.5=(0.5-0.5n+1n+1,0.5,0.5n n,0.5-0.5,0.5-0.5n+1n+1)T T.n nF(6)F(6)=(0.4921875,0.015625,0.4921875)=(0.4921875,0.015625,0.4921875)T T.n n一般杂交后至少要经过六代的选择、自交才有可能作为一个稳定的品系参加品比试验。第二十二张，PPT共三十五页，创作于2022年6月n n五.状态转移的马尔可夫链模型n n 1.遗传模型的特征与马尔可夫链n n 10.研
21、究的对象具有几个可能的不同状态。n n 20.各状态之间可能会发生随机转移的现象。n n 30.状态转移时将来的状态只与现在的状态有关，而独立与过去的状态。n n 我们称描述这一类状态随机转移问题的模型为马尔可夫链模型。第二十三张，PPT共三十五页，创作于2022年6月n n2.状态转移的流程图：n n 用图示来描述状态及其相互转移的关系。n n F1 F2 1/21/2 nn n n 1/21/2AAAaaa F2 F2 1/4 1/4 1/2AAAaaa第二十四张，PPT共三十五页，创作于2022年6月F1 F 1/2 1 1 1/2AAAaaaF2 F 1/4 1/2 1/2 1/2 1
22、/2 1/2 AaaaAA F3 F 1/2 1 1/2 1Aa aa AA 第二十五张，PPT共三十五页，创作于2022年6月 n n随机交配群体n n p p1212=q p=q p2323=q/2=q/2n n p p2121=p/2 p=p/2 p3232=p=pn n n n p p1111=p p=p p2222=1/2 p=1/2 p3333=q=qF2:Aa F1:AAF3:aa第二十六张，PPT共三十五页，创作于2022年6月n n其中，fi 表示参与转移的状态的频率，pij 表示从状态 j 向状态 i 转移的概率。n n如果用 F(n)=(f1(n),f2(n),f3(n)
23、T 表示第 n 代群体的遗传结构,n n它实际上描述了一个个体在群体中处于各种状态的概率,n n我们称之为状态向量.第二十七张，PPT共三十五页，创作于2022年6月n n矩阵 M 则描述了状态随机转移的状况,n n称之为转移矩阵.n nM 非负且列和等于1。n n流程图所描述的状态转移的过程就可以用模型 F(n+1)=M F(n)来描述.n n由此得到的状态向量序列F(n)称之为马尔可夫链.模型就称为马氏链模型.第二十八张，PPT共三十五页，创作于2022年6月n n 4.4.模型的分析模型的分析n n 1 10 0.正则链正则链.对于转移矩阵对于转移矩阵 M,M,若存在正整数若存在正整数
24、k,k,使得使得 MMk k 0(p 0(pij ij 0),0),则称则称 M M 是正则的是正则的,M,M所描述的链为所描述的链为正则正则链链.我们有我们有:n n (1)M (1)M 有唯一的不动点向量有唯一的不动点向量 =(=(1 1,n n)T T,使得使得 M M =.n n (2),(2),其中其中 1 1=(1,1)=(1,1)T T.n n n n (3).(3).n n n n正则链的状态向量一定趋向于一个稳定的分布,并且这个分布与状态向量的初始取值无关.第二十九张，PPT共三十五页，创作于2022年6月n n 20.吸收链.若状态转移矩阵有形式n n其中 Ir 为 r 阶
25、单位阵,O为零阵,n n R 的行和均为正.n n则称由 M 给出的马氏链为吸收链.n n矩阵 Ir所对应的状态 S1 称为吸收态,n n其余的状态 S2 为非吸收态.n n 第三十张，PPT共三十五页，创作于2022年6月n n（1）例第三十一张，PPT共三十五页，创作于2022年6月n n(2)分析n n 10lim S(n)T=(S1T,0T)n n 20 称 F=(I Q)-1 为吸收链的基本矩阵.n nF=(fij),其中 fij 为从非吸收状态 si 在被吸收前转移到非吸收状态 sj 的平均转移次数.n n 30 令 B=F RT=(bij),则给出了从非吸收状态 fi出发,被吸收
26、状态 fj 吸收的概率.第三十二张，PPT共三十五页，创作于2022年6月n nF=(I Q)-1=(1-)-1=2n n从非吸收状态S3出发,平均转移2 次,就会被吸收状态S1,S2 吸收.n nB=F RT=2(1/4 1/4)n n =(1/2 1/2)n n从非吸收状态S3出发,被吸收状态S1,S2 吸收的概率分别为 1/2.第三十三张，PPT共三十五页，创作于2022年6月n n问题:n n 1.乒乓球比赛，每胜一球得一分，先得三分者胜一局。建模描述一局乒乓球的比赛过程。n n 2.在上面的乒乓球比赛中，如果双方打成 2:2 平局，则在后面的比赛中需要一方连续胜两球才能赢得此局。是给出模型描述一局乒乓球比赛的过程。第三十四张，PPT共三十五页，创作于2022年6月感感谢谢大大家家观观看看第三十五张，PPT共三十五页，创作于2022年6月