群体遗传的数学模型

关于群体遗传的数学模型第1页，讲稿共35页，2023年5月2日，星期三一.问题与背景

1.问题生物的细微变异在进化中很重要，可以为物种进化提供原始材料。但变异的个体在群体中数量很少。它会不会在群体的随机交配的过程中逐渐减弱直至消失？人类的遗传病（如色盲）在群体中会不会由于一代一代地遗传而使患者越来越多？

第2页，讲稿共35页，2023年5月2日，星期三农作物的新品种主要是用杂交的方法培育的。即选择某种作物具有互补性状的两个品种杂交，在它们的后代中选择培育出合乎要求的新品种。由于杂交后代性状分离的不稳定性，使得它不可能杂交后直接用于农业生产。多少代以后就可以达到生产所需要的稳定性呢？第3页，讲稿共35页，2023年5月2日，星期三

2.背景：孟德尔遗传学的要点

1.生物体的性状是由一对（多对）基因控制的。染色体是基因的载体。基因在染色体上的位置称为基因座，处于同一基因座的不同的基因为等位基因。等位基因对性状的控制有显性和隐性之分。

第4页，讲稿共35页，2023年5月2日，星期三

2.处于同一基因座上的两个基因的组合形式称为基因型。生物体表现出来的特征为表现型。一般表现型反映显性基因所控制的性状。

3.遗传的过程是由亲代每基因座上的等位基因分离并随机进入到一个配子中去在交配过程中雌雄配子将随机结合形成新的合子而传入下一代。第5页，讲稿共35页，2023年5月2日，星期三二.假设

1.雌性和雄性在遗传上是对等的。

2.没有迁移、突变和选择。

3.群体足够大。

4.随机交配，控制有关性状的雌雄配子完全随机结合。第6页，讲稿共35页，2023年5月2日，星期三

三.模型

1.群体的遗传结构：一个基因座，两个等位基因A，a。三个基因型AA，Aa，aa。基因型频率：f1，f2，f3，(f1+f2+f3=1).

基因频率:p=f1+f2/2,q=f3+f2/2.

群体的遗传结构:

F=(f1，f2，f3)T.

基本群体:

F1=(1,0,0)T,

F2=(0,1,0)T,

F3=(0,0,1)T.

F=f1×F1+f2×F2+f3×F3.第7页，讲稿共35页，2023年5月2日，星期三

2.遗传过程的基本平衡关系由孟德尔遗传学的理论可知,有

AAAA=AA,AAaa=Aa,aaaa=aa,AAAa=(1/2)AA+(1/2)Aa,aaAa=(1/2)Aa+(1/2)aa,AaAa=(1/4)AA+(1/2)Aa+(1/4)aa

即,有如下的平衡关系

F1F1=F1,

F1F3

=

F2,

F3F3=

F3,F1

F2=(1/2)F1+(1/2)F2,

F3

F2=(1/2)F2+(1/2)F3,F2

F2=(1/4)F1+(1/2)F2+(1/4)F3.第8页，讲稿共35页，2023年5月2日，星期三

3.群体遗传的模型:

F(n)=(f1(n),f2(n),f3(n))T:第n代群体的结构.10.基本遗传过程:

F1与F(n)交配的后代F(n+1)=F1F(n)=F1[f1(n)F1+f2(n)F2+f3(n)F3]=f1(n)[

F1

F1]+f2(n)[F1

F2]+f3(n)[

F1

F3]

=

[F1

F1F1

F2F1

F3]

F(n)

=[

F1

(F1+F2)/2

F2]

F(n)=M1

F(n)第9页，讲稿共35页，2023年5月2日，星期三

F2与F(n)交配的后代

F(n+1)=F2F(n)=M2

F(n)F3与F(n)交配的后代

F(n+1)=F3F(n)=M3

F(n)第10页，讲稿共35页，2023年5月2日，星期三

20.随机交配群体的遗传模型:一个结构为F的群体与第n世代的群体

G(n)

随机交配的后代.G(n+1)=FG(n)=(f1F1+f2F2+f3F3)G(n)

=f1[F1G(n)

]+f2[F2G(n)

]+f3[F3G(n)]

=f1M1

G(n)

+f2M2

G(n)

+f3M3

G(n)

=(f1M1+f2M2+f3M3)

G(n)

=Mr

G(n)

第11页，讲稿共35页，2023年5月2日，星期三

30.自交群体的遗传模型.自交群体指每个个体自交所产生的后代群体.

F(n+1)=F(n)

F(n)

=f1(n)F

F1+f2(n)F2

F2+f3(n)F3

F3

=[F1F1/4+

F2/2+

F3/4

F3]

F(n)

=MS

F(n)第12页，讲稿共35页，2023年5月2日，星期三四.模型的分析

1.Hardy–Weinberg平衡(H–W平衡)

在前述假设下,一个随机交配的有性繁殖的群体中，等位基因频率保持不变，基因型频率至多经过一个世代也将保持不变。一个基因座，两个等位基因。

群体F

=(f1,f2,f3)T;f1,f2,f3

基因型AA,Aa,aa的频率;p=f1+f2/2,q=f3+f2/2基因A,a的频率.第13页，讲稿共35页，2023年5月2日，星期三

群体自身的随机交配,有模型:

F(n+1)

=Mr(n)F(n)其中

p(n)=f1(n)+f2(n)/2,q(n)=f3(n)+f2(n)/2.第14页，讲稿共35页，2023年5月2日，星期三

1.证明p(n+1)=p(n),q(n+1)=q(n)由模型可以得到f1(n+1)=p(n)f1(n)+p(n)f2(n)/2=p(n)[f1(n)+f2(n)/2]=

p(n)2

f3(n+1)=q(n)f3(n)+q(n)f2(n)/2=q(n)[f3(n)+f2(n)/2]=q(n)2

f2(n+1)=

q(n)

f1(n)+[p(n)+q(n)]f2(n)/2+p(n)f3(n)=q(n)[f1(n)+f2(n)/2]+p(n)[f3(n)+f2(n)/2]=2p(n)q(n)

p(n+1)=f1(n+1)+f2(n+1)/2=p(n),q(n+1)=f3(n+1)+f2(n+1)/2=q(n).

第15页，讲稿共35页，2023年5月2日，星期三

2.证明基因型频率不变(至多经过一个世代)

矩阵Mr的列和为1,因此它有全1的左特征向量1=(1,1,1)T,

且有1T

Mr=1T.

由此可知矩阵有特征值λ=1,并且只有这一个正特征根它所对应的右特征向量为F*

>0,(Perron–Frobenious定理)由MrF*=F*

不难得到F*=(p22pqq2)T就是模型的平衡态。对于任意的初始状态F(0)=(f1(0),f2(0),f3(0))T,

容易算出F(1)=MrF(0)=F*.

因此,至多经过一个世代,群体的基因型频率将处于平衡状态不变。第16页，讲稿共35页，2023年5月2日，星期三例1.M–N血型：血红细胞膜上的特异抗原类型，分别对抗体M和抗体N成阳性反应，这一性状分别由一对等位基因LM和LN控制，呈并显性。调查1279人三种基因型频数的比例为

LMLM:LMLN:LNLN=363:634:282由此可以算出群体的遗传结构为：

(f1f2f3)=(0.284，0.496，0.220)群体中基因LM和LN的频率分别为

p=0.532,q=0.468p2=0.283,2pq=0.498,q2=0.219与群体的遗传结构基本相同。群体为H–W群体。第17页，讲稿共35页，2023年5月2日，星期三例2.地中海贫血病，由一对等位基因Tm、Tn控制。调查意大利群体10000人三种基因型的频数为

TmTm:TmTn:TnTn=4:400:9596基因Tm和Tn的频率分别为

p=0.0204，q=0.9796

p2=0.00042，2pq=0.03997,q2=0.95961。与群体的遗传结构基本相同。群体为H–W群体第18页，讲稿共35页，2023年5月2日，星期三问题：P1758民族

O(%)

A(%）

B(%)

AB(%)汉

30.86

31.31

28.06

9.77维吾尔

27.50

29.22

31.92

11.36壮

47.28

21.25

27.57

3.90回

35.94

27.23

28.34

8.49哈萨克

37.97

22.82

29.83

9.38锡伯

24.42

25.00

40.12

10.46乌兹别克

25.58

25.58

38.76

10.08柯尔克孜

34.68

18.54

39.52

7.26塔克尔

21.62

40.54

35.14

2.70白

31.40

34.00

23.40

11.20傣

40.44

22.08

29.59

7.89景颇

37.81

34.83

20.39

6.97第19页，讲稿共35页，2023年5月2日，星期三

2.杂交育种的稳定性作物的品种：同型合子的基因型育种的过程：双亲杂型合子分离后代优良组合分离后代

优良组合…第20页，讲稿共35页，2023年5月2日，星期三问题：杂交后要自交多少世代就可以形成品种？

F(n)=MSF(n-1)=…=MSn-1F(1)=MSnF2.λ=1u1=(1,0,-1)Tλ=1u2=(0,0,1)Tλ=0.5u3=(1,-2,1)T.令则有分解式MS=U-1ΛU,MSn=U-1ΛnU.第21页，讲稿共35页，2023年5月2日，星期三由此不难算出可得

F(n)

=(0.5-0.5n+1,0.5n,0.5-0.5n+1)T.F(6)

=(0.4921875,0.015625,0.4921875)T.一般杂交后至少要经过六代的选择、自交才有可能作为一个稳定的品系参加品比试验。第22页，讲稿共35页，2023年5月2日，星期三五.状态转移的马尔可夫链模型

1.遗传模型的特征与马尔可夫链

10.研究的对象具有几个可能的不同状态。

20.各状态之间可能会发生随机转移的现象。

30.状态转移时将来的状态只与现在的状态有关，而独立与过去的状态。我们称描述这一类状态随机转移问题的模型为马尔可夫链模型。第23页，讲稿共35页，2023年5月2日，星期三2.状态转移的流程图：用图示来描述状态及其相互转移的关系。

F1

F21/2

1/2AAAaaa

F2F21/4

1/41/2AAAaaa第24页，讲稿共35页，2023年5月2日，星期三F1F1/2111/2AAAaaaF2F1/41/21/2¼1/21/21/2

AaaaAAF3F1/211/21AaaaAA第25页，讲稿共35页，2023年5月2日，星期三

随机交配群体

p12=qp23=q/2p21=p/2p32=p

p11=pp22=1/2p33=qF2:Aa

F1:AAF3:aa第26页，讲稿共35页，2023年5月2日，星期三其中，{fi}表示参与转移的状态的频率，pij

表示从状态j

向状态i

转移的概率。如果用F(n)=(f1(n),f2(n),f3(n))T表示第n代群体的遗传结构,它实际上描述了一个个体在群体中处于各种状态的概率,我们称之为状态向量.第27页，讲稿共35页，2023年5月2日，星期三矩阵M则描述了状态随机转移的状况,称之为转移矩阵.M非负且列和等于1。流程图所描述的状态转移的过程就可以用模型F(n+1)

=MF(n)

来描述.由此得到的状态向量序列{F(n)}称之为马尔可夫链.模型就称为马氏链模型.第28页，讲稿共35页，2023年5月2日，星期三

4.模型的分析

10.正则链.对于转移矩阵M,若存在正整数k,使得Mk>0(pij>0),则称M是正则的,M所描述的链为正则链.我们有:(1)M有唯一的不动点向量π

=(π1,…,πn)T,使得Mπ

=π.(2),其中1=(1,…,1)T.

(3).

正则链的状态向量一定趋向于一个稳定的分布π,并且这个分布与状态向量的初始取值无关.第29页，