基于GWAS結(jié)果的SNP篩選與坐標(biāo)檢索指南

1、關(guān)于潛關(guān)聯(lián)SNP的篩選

篩選疾病風(fēng)險(xiǎn)SNP時(shí),選擇GWAS p值小于0.05的結(jié)果 并不是一個(gè)最佳的策略,尤其是在全基因組范圍的關(guān)聯(lián)分析中讥脐。這是因?yàn)樵贕WAS中,通常要對(duì)數(shù)百萬(wàn)個(gè)SNP進(jìn)行檢驗(yàn)啼器,因此使用標(biāo)準(zhǔn)的p值閾值(例如0.05)會(huì)導(dǎo)致大量的假陽(yáng)性結(jié)果旬渠。為了解決這一問(wèn)題,GWAS分析中通常使用更嚴(yán)格的多重檢驗(yàn)校正閾值端壳。

為什么不建議使用 p < 0.05告丢?

  • 多重檢驗(yàn)問(wèn)題:在GWAS中,數(shù)百萬(wàn)個(gè)SNP同時(shí)被檢驗(yàn)损谦。如果僅使用p < 0.05的標(biāo)準(zhǔn)岖免,大量的SNP可能由于隨機(jī)波動(dòng)而被錯(cuò)誤地標(biāo)記為顯著。這會(huì)導(dǎo)致大量的假陽(yáng)性照捡。
  • 假陽(yáng)性率過(guò)高:考慮到GWAS分析中通常涉及數(shù)百萬(wàn)次比較颅湘,使用p < 0.05 的標(biāo)準(zhǔn)會(huì)產(chǎn)生大量的假陽(yáng)性。以p < 0.05篩選時(shí)栗精,大約會(huì)有5%的無(wú)效 SNP 被錯(cuò)誤地認(rèn)為是顯著的闯参。

推薦的 p 值閾值:Bonferroni 校正False Discovery Rate (FDR)

為了控制假陽(yáng)性,GWAS中通常會(huì)使用以下更嚴(yán)格的校正方法:

(1). Bonferroni 校正

  • Bonferroni校正是一種較為保守的多重比較校正方法悲立,它將標(biāo)準(zhǔn)的p值閾值除以檢驗(yàn)的總次數(shù)(SNP數(shù)目)鹿寨。
  • 假設(shè)你檢驗(yàn)了1,000,000個(gè)SNP,Bonferroni校正后的p值閾值會(huì)是 0.05 / 1,000,000 = 5 × 10??薪夕。
  • 在GWAS中脚草,p < 5 × 10?? 通常被認(rèn)為是全基因組顯著性閾值,用于識(shí)別有顯著關(guān)聯(lián)的SNP寥殖。

(2). False Discovery Rate (FDR) 控制

  • FDR 控制是一種稍微寬松的校正方法玩讳,可以控制錯(cuò)誤識(shí)別的比例涩蜘,而不是完全排除所有的假陽(yáng)性嚼贡。
  • 常用的FDR閾值是 q < 0.05,這意味著允許最高5%的假陽(yáng)性結(jié)果同诫。

結(jié)論

雖然p < 0.05在其他統(tǒng)計(jì)分析中是常見(jiàn)的顯著性標(biāo)準(zhǔn)粤策,但在GWAS中并不合適,因?yàn)樯婕暗酱罅康腟NP測(cè)試误窖。更合理的方法是使用p < 5 × 10??作為全基因組顯著性的標(biāo)準(zhǔn)叮盘,或者使用FDR校正來(lái)篩選有顯著關(guān)聯(lián)的SNP秩贰。

某些情況下,GWAS 樣本量相對(duì)較小柔吼,或者研究對(duì)象具有獨(dú)特的遺傳背景毒费,嚴(yán)格的閾值可能過(guò)于保守,導(dǎo)致更高的假陰性率愈魏,從而漏掉真正與性狀相關(guān)的SNP觅玻,無(wú)法捕獲潛在信號(hào)。

解決方案:

還有 80% 的精彩內(nèi)容
最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
禁止轉(zhuǎn)載培漏,如需轉(zhuǎn)載請(qǐng)通過(guò)簡(jiǎn)信或評(píng)論聯(lián)系作者溪厘。
支付 ¥6.00 繼續(xù)閱讀
  • 序言:七十年代末,一起剝皮案震驚了整個(gè)濱河市牌柄,隨后出現(xiàn)的幾起案子畸悬,更是在濱河造成了極大的恐慌,老刑警劉巖珊佣,帶你破解...
    沈念sama閱讀 206,126評(píng)論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件蹋宦,死亡現(xiàn)場(chǎng)離奇詭異,居然都是意外死亡咒锻,警方通過(guò)查閱死者的電腦和手機(jī)妆档,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,254評(píng)論 2 382
  • 文/潘曉璐 我一進(jìn)店門(mén),熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái)虫碉,“玉大人贾惦,你說(shuō)我怎么就攤上這事《嘏酰” “怎么了须板?”我有些...
    開(kāi)封第一講書(shū)人閱讀 152,445評(píng)論 0 341
  • 文/不壞的土叔 我叫張陵,是天一觀的道長(zhǎng)兢卵。 經(jīng)常有香客問(wèn)我习瑰,道長(zhǎng),這世上最難降的妖魔是什么秽荤? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 55,185評(píng)論 1 278
  • 正文 為了忘掉前任甜奄,我火速辦了婚禮,結(jié)果婚禮上窃款,老公的妹妹穿的比我還像新娘课兄。我一直安慰自己,他們只是感情好晨继,可當(dāng)我...
    茶點(diǎn)故事閱讀 64,178評(píng)論 5 371
  • 文/花漫 我一把揭開(kāi)白布烟阐。 她就那樣靜靜地躺著,像睡著了一般。 火紅的嫁衣襯著肌膚如雪蜒茄。 梳的紋絲不亂的頭發(fā)上唉擂,一...
    開(kāi)封第一講書(shū)人閱讀 48,970評(píng)論 1 284
  • 那天,我揣著相機(jī)與錄音檀葛,去河邊找鬼玩祟。 笑死,一個(gè)胖子當(dāng)著我的面吹牛屿聋,可吹牛的內(nèi)容都是我干的卵凑。 我是一名探鬼主播,決...
    沈念sama閱讀 38,276評(píng)論 3 399
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼胜臊,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼勺卢!你這毒婦竟也來(lái)了?” 一聲冷哼從身側(cè)響起象对,我...
    開(kāi)封第一講書(shū)人閱讀 36,927評(píng)論 0 259
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤黑忱,失蹤者是張志新(化名)和其女友劉穎,沒(méi)想到半個(gè)月后勒魔,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體甫煞,經(jīng)...
    沈念sama閱讀 43,400評(píng)論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 35,883評(píng)論 2 323
  • 正文 我和宋清朗相戀三年冠绢,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了抚吠。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 37,997評(píng)論 1 333
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡弟胀,死狀恐怖楷力,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情孵户,我是刑警寧澤萧朝,帶...
    沈念sama閱讀 33,646評(píng)論 4 322
  • 正文 年R本政府宣布,位于F島的核電站夏哭,受9級(jí)特大地震影響检柬,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜竖配,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,213評(píng)論 3 307
  • 文/蒙蒙 一何址、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧进胯,春花似錦用爪、人聲如沸。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 30,204評(píng)論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)班眯。三九已至希停,卻和暖如春烁巫,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背宠能。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 31,423評(píng)論 1 260
  • 我被黑心中介騙來(lái)泰國(guó)打工亚隙, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留,地道東北人违崇。 一個(gè)月前我還...
    沈念sama閱讀 45,423評(píng)論 2 352
  • 正文 我出身青樓阿弃,卻偏偏與公主長(zhǎng)得像,于是被迫代替她去往敵國(guó)和親羞延。 傳聞我的和親對(duì)象是個(gè)殘疾皇子渣淳,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 42,722評(píng)論 2 345

推薦閱讀更多精彩內(nèi)容