「什么是哈溫平衡？」

?

哈迪-溫伯格(Hardy-Weinberg)法則 哈迪-溫伯格（Hardy-Weinberg）法則是群體遺傳中最重要的原理，它解釋了繁殖如何影響群體的基因和基因型頻率。這個(gè)法則是用Hardy,G.H (英國(guó)數(shù)學(xué)家) 和Weinberg，W.（德國(guó)醫(yī)生）兩位學(xué)者的姓來(lái)命名的，他們于同一年（1908年）各自發(fā)現(xiàn)了這一法則。他們提出在一個(gè)不發(fā)生突變、遷移和選擇的無(wú)限大的隨機(jī)交配的群體中，基因頻率和基因型頻率將逐代保持不變。---百度百科

?

「怎么做哈溫平衡檢驗(yàn)？」

?

「卡方適合性檢驗(yàn)！」 ，一個(gè)群體是否符合這種狀況，即達(dá)到了遺傳平衡，也就是一對(duì)等位基因的3種基因型的比例分布符合公式：p2+2pq+q2=1,p+q=1,(p+q)2=1.基因型MM的頻率為p2,NN的頻率為q2,MN的頻率為2pq。MN:MN：NN＝P2：2pq：q2。MN這對(duì)基因在群體中達(dá)此狀態(tài)，就是達(dá)到了遺傳平衡。如果沒(méi)有達(dá)到這個(gè)狀態(tài)，就是一個(gè)遺傳不平衡的群體。但隨著群體中的隨機(jī)交配，將會(huì)保持這個(gè)基因頻率和基因型分布比例，而較易達(dá)到遺傳平衡狀態(tài)。應(yīng)用Hardy-Weinberg遺傳平衡吻合度檢驗(yàn)方法，把計(jì)算得到的基因頻率代入，計(jì)算基因型平衡頻率，再乘以總?cè)藬?shù)，求得預(yù)期值（e）。把觀察數(shù)（O）與預(yù)期值（e）作比較，進(jìn)行χ2檢驗(yàn)。病例組和對(duì)照組的基因型分布的觀察值和預(yù)期值差異無(wú)顯著性（P>0.05），符合遺傳平衡定律.

?

「哈溫平衡過(guò)濾和MAF過(guò)濾的區(qū)別？」

?

之前，我對(duì)這兩個(gè)概念有點(diǎn)混淆，后來(lái)明白過(guò)來(lái)了。這兩個(gè)概念一個(gè)是對(duì)基因頻率進(jìn)行的篩選，一個(gè)是對(duì)基因型頻率進(jìn)行的篩選。對(duì)于一個(gè)位點(diǎn)“AA AT TT”，其中A的頻率為基因頻率，AA為基因型頻率。MAF直接是對(duì)基因頻率進(jìn)行篩選，而哈溫平衡檢驗(yàn)，則是根據(jù)基因型推斷出理想的（AA，AT，TT）的分布，然后和實(shí)際觀察的進(jìn)行適合性檢驗(yàn)，然后得到P值，根據(jù)P值進(jìn)行篩選。即P值越小，說(shuō)明該位點(diǎn)越不符合哈溫平衡。

?

「兩個(gè)目的：」

• 計(jì)算所有位點(diǎn)的哈溫檢測(cè)結(jié)果
• 刪除SNP中不符合哈溫平衡的位點(diǎn)

1. 計(jì)算所有位點(diǎn)的HWE的P值

plink.hwe的數(shù)據(jù)格式：

• CHR 染色體
• SNP SNP的ID
• TEST 類(lèi)型
• A1 minor 位點(diǎn)
• A2 major 位點(diǎn)
• GENO 基因型分布：A1A1, A1A2, A2A2
• O(HET) 觀測(cè)雜合度頻率
• E(HET) 期望雜合度頻率
• P 哈溫平衡的卡方檢驗(yàn)P-value值

結(jié)果預(yù)覽：

2. 提取哈溫p值小于0.0001的位點(diǎn)

這里我們使用awk：

共有123個(gè)位點(diǎn)，其中UNAFF為45個(gè)位點(diǎn)。

3.  設(shè)定過(guò)濾標(biāo)準(zhǔn)1e-4

日志：

可以看到，共有45個(gè)SNP根據(jù)哈溫的P值過(guò)濾掉了，和上面手動(dòng)計(jì)算的一樣。

4. 可視化

R代碼：

哈溫的P值直方圖：

過(guò)濾掉SNP位點(diǎn)的P值：

過(guò)濾后的結(jié)果文件