GWAS+自然選擇:62個樣本的GWAS分析,沒信號,如何巧妙的發文章


歡迎來到"bio生物信息"的世界

6天前,BMC Genomics 推了一篇文獻“Population history and genetic adaptation of the Fulani nomads: inferences from genome-wide data and the lactase persistence trait”。

要不是這個標題起的太大,又是Population history,又是 genetic adaptation ,我可能都不會點開來看。

但是既然看了,還是要寫點總結,畢竟,每一篇能發出來的文章,總是有那么一點值得我們借鑒的。

這篇文章很聰明的一點是,他避開了GWAS沒信號這個問題,轉而將重點放在:結合表型和人群歷史以及遺傳適應性來講。

下面講講整篇文章的研究思路。

0. 研究背景

Fulani nomads是在非洲西部的一群游牧民族。

游牧民族最大的特點是什么。當然是(牛/羊)奶多。

所以乳糖耐性就是一個明顯受到選擇的表型。

這也是他們本次研究采用乳糖耐性作為表型的原因。

講完背景,我們接下來看看他們做了什么工作。

1. 研究的人群的歷史。

62個樣本,講人群歷史當然很干巴巴。

那么,就需要增大樣本量了。

增大樣本量,首選就是千人基因組公共數據庫了。 全世界各個群體的基因型數據都有,不嫖白不嫖。

QUQsbR.png

這張圖用的是EIGENSOFT 做PCA分析,ADMIXTURE推斷人群結構。

2.祖先特異性推斷

對之前文獻發表的乳糖耐性相關變異位點進行祖先特異性推斷。

乳糖耐性相關的基因周圍的位點多樣性降低了,這表明乳糖耐性在這個群體是受到強烈的選擇。

QU8esg.md.png

3.Fulani人群與乳糖耐性相關的顯著位點分析

這一步就是用的GWAS分析。

當然,62個樣本,找不到信號也很正常。

因此,作者還結合了Fulani人群的iHS值和XPEHH值,放在一起和GWAS找出來的最高信號一起討論。

這么一結合,文章的檔次就提升了。

QUxCF0.md.png

4.結束語

各位找不到GWAS信號的同行,也不用太灰心。

畢竟GWAS只是找遺傳位點的一個方法而已,可我們的研究內容不只找遺傳位點啊。

項目推不下去的時候,就借鑒別人的研究思路,為自己的項目錦上添花。


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM