跳到主要内容

東京大学 情報理工学系研究科 コンピュータ科学専攻 2015年8月実施 専門科目I 問題3

Author

祭音Myyura (co-authored with GPT 5.6 SOL)

Description

某基因有等位基因 A,aA,a,个体基因型为 AA,Aa,aaAA,Aa,aa 之一;子代从父母各自的两条常染色体中等概率继承一条。

(1)父母基因型均为 AaAa 时,求子代为 AAAA 的概率。

(2)已知父亲为 AaAa、子代为 aaaa。母亲两条染色体上的等位基因相互独立地以概率 2/3,1/32/3,1/3A,aA,a。求母亲最可能的基因型。

以下考虑足够大的随机交配种群。无世代重叠、迁移与突变,且同一对常染色体上的两个等位基因独立出现。第 nn 代出生后 A,aA,a 的频率分别为 pn,qnp_n,q_n,其中 pn+qn=1p_n+q_n=1,故 AA,Aa,aaAA,Aa,aa 的频率为 pn2,2pnqn,qn2p_n^2,2p_nq_n,q_n^2

(3)若 AA,Aa,aaAA,Aa,aa 存活至成年的概率(适合度)分别为 1,1s,12s1,1-s,1-2s,用 qn,sq_n,s 表示 qnqn+1q_n-q_{n+1}

(4)若三种基因型适合度分别为 1,1,01,1,0,且 q0=1/100q_0=1/100,求使 qn1/10000q_n\le1/10000 的最小 nn

Kai

(1)

父母各自传递 AA 的概率均为 1/21/2,故

Pr(AA)=14.\boxed{\Pr(AA)=\frac14}.

(2)

母亲的先验概率为

Pr(AA)=49,Pr(Aa)=49,Pr(aa)=19.\Pr(AA)=\frac49,\qquad \Pr(Aa)=\frac49,\qquad \Pr(aa)=\frac19.

在父亲已知为 AaAa 时,观测到子代 aaaa 的似然对母亲 AA,Aa,aaAA,Aa,aa 分别为 0,1/4,1/20,1/4,1/2。因而后验权重为

0:19:118=0:2:1.0:\frac19:\frac1{18}=0:2:1.

所以母亲最可能的基因型是 Aa\boxed{Aa}(后验概率为 2/32/3)。

(3)

选择后的平均适合度为

wˉ=pn2+2pnqn(1s)+qn2(12s)=12sqn.\bar w=p_n^2+2p_nq_n(1-s)+q_n^2(1-2s)=1-2sq_n.

成年个体产生的配子中,等位基因 aa 的频率为

qn+1=pnqn(1s)+qn2(12s)wˉ=qn(1ssqn)12sqn.q_{n+1} =\frac{p_nq_n(1-s)+q_n^2(1-2s)}{\bar w} =\frac{q_n(1-s-sq_n)}{1-2sq_n}.

因此

qnqn+1=sqn(1qn)12sqn.\boxed{q_n-q_{n+1}=\frac{s q_n(1-q_n)}{1-2s q_n}}.

(4)

只有 AA,AaAA,Aa 能繁殖,故

qn+1=pnqnpn2+2pnqn=qn1+qn.q_{n+1}=\frac{p_nq_n}{p_n^2+2p_nq_n} =\frac{q_n}{1+q_n}.

于是

1qn+1=1qn+1,qn=1100+n.\frac1{q_{n+1}}=\frac1{q_n}+1,\qquad q_n=\frac1{100+n}.

qn1/10000q_n\le1/10000 等价于 100+n10000100+n\ge10000,故最小值为

n=9900.\boxed{n=9900}.