【直播】我的基因组 37:gwas研究结果在我身上得到了验证

前面讲到了我的480万变异里面里面有一些在dbSNP数据库里面记录着可能是somatic的变异,让我有点担心,尤其是我还看到了好几个MUC系列基因,主要是黏液素家族基因,而Mucin gene expression in the effusions of otitis media with effusion.让我想起了儿时的种种不愉快,反正这个也不是什么隐私了,曝就曝吧!

中耳炎是一系列中耳发炎疾病之统称。其中以急性中耳炎(AOM)和中耳积水两型最为常见(OME)(Otitis Media With Effusion)。

如果要了解关于它的医学知识,建议自行搜索。

http://emedicine.medscape.com/article/994656-overview

http://emedicine.medscape.com/article/994656-treatment

http://emedicine.medscape.com/article/994656-medication

本文的重点是探究如何把现有的研究结果的我的全基因组测序数据结合起来

我的搜索关键词是:Otitis Media With Effusion gene ,谷歌可以得到一系列的结果,当然需要你英文水平还行才能看下去。

有些文章说的是某些SNP位点与该OME的相关程度,其它的是在探究OME患者跟正常人有哪些基因的表达是显著变化的,当然也可以是研究甲基化情况,或者什么miRNA的调控呀,转录因子的结合情况呀,蛋白或者代谢物的表达情况呀!

值得我关心的有3个;

https://www.ncbi.nlm.nih.gov/pubmed/27677580

https://www.ncbi.nlm.nih.gov/pubmed/23929584

https://www.ncbi.nlm.nih.gov/pubmed/16847180

先看第一个GWAS分析:

a total of 825 cases and 7,936 controls of European descent retained for analysis at CHOP together with 1219 cases and 1,067 controls at Generation R. Logistic regression analysis under an additive model was performed to assess the association between SNP genotype and AOM at each individual site, including the first three principal components and sex as covariates.

也就是这个科学团体研究了近万人才得到了这个可疑位点,跟OME显著相关的。

前面我们已经把我的vcf文件添加了dbSNP的标签,简单搜索一下发现,我的确有这个变异,唉。(虽然文章里面强调了这个研究来源是欧洲人,但是人种差异在这里并没有 体现出来。)

grep rs2932989 realign.dbsnp.vcf

6 159699284 rs2932989 T G 228.0 . DP=57;VDB=0.94206;SGB=-0.693147;MQSB=1;MQ0F=0;AC=2;AN=2;DP4=0,0,36,11;MQ=60;ASP;CAF=0.1034,0.8966;COMMON=1;G5;GNO;HD;KGPhase1;KGPhase3;RS=2932989;RSPOS=159699284;RV;SAO=0;SLO;SSR=0;VC=SNV;VLD;VP=0x05010000000515053f000100;WGT=1;dbSNPBuildID=101 GT:PL 1/1:255,141,0

我在IGV看了这个位点,唉,很靠谱的变异!

这次真的是身体力行的证明了GWAS研究是靠谱的,还是不开心,就到处为止吧。

参考:http://www.nature.com/articles/ncomms12792/figures/1

https://www.ncbi.nlm.nih.gov/projects/SNP/snp_ref.cgi?rs=2932989

文:Jimmy

图文编辑:吃瓜群众

原文发布于微信公众号 - 生信技能树(biotrainee)

原文发表时间:2017-01-01

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏PPV课数据科学社区

【译文】关于数据科学难以忽视的11大真相

1、Data is never clean. 2、You will spend most of your time cleaning and preparing...

3447
来自专栏Java学习网

游戏编程十年总结

游戏编程十年总结 自敲第一行代码起,已经十年多了,今天既不是十年整的日子,也不是一个有特定意义的日子,本来像这种大总结的文章,当择良辰吉日,斋戒沐浴三日,方可...

3087
来自专栏机器人网

想学程序猿拿高薪?你得先买一个教你写代码的机器人

都说程序猿是一个高薪水的职业,这是真的吗?小编没调查过,但是小编身边的程序猿朋友,虽说大框架眼镜下迷离的眼神+日渐闪亮的发际线,也挡不住写代码的热情,说不准厚厚...

2506
来自专栏生信技能树

【生信菜鸟经】漫谈如何跨越摆在生信入门路上的三大障碍

如何跨越摆在生信入门菜鸟面前的三大障碍的。 第一大障碍:透析数据背后的生物学知识,完成从测序数据到生物问题的连接; 我的做法: 通过谷歌百度了解fasta、f...

3247
来自专栏java一日一条

游戏编程十年总结

自敲第一行代码起,已经十年多了,今天既不是十年整的日子,也不是一个有特定意义的日子,本来像这种大总结的文章,当择良辰吉日,斋戒沐浴三日,方可动笔。一开始计划是写...

591
来自专栏BestSDK

【七夕佳节】今晚我们还是谈谈Bug好了

最近看了「简书」上的一篇文章《女生如何优雅的上下车》,可是,额,我腿短,下车时依然是艰难得一点点‘挪’出去,正如懂得很多道理,却依然过不好这一生,啧啧,又如程序...

2285
来自专栏知晓程序

王者荣耀有罪吗?用这款小程序,亮出你的态度

除了亲儿子腾讯投票,以及它的挑战者投票帮。今天,知晓程序(微信号 zxcx0101)还要给大家介绍一款专门的投票小程序。

932
来自专栏Crossin的编程教室

一行代码扫出“敬业福”

好吧,我承认有那么一点标题党。不过说起标题党这事儿,咱先来看看支付BAO,最近几天搞得全国人民都不安心工作的“集五福”: ? 好(shua)好(hou)的“集五...

3458
来自专栏程序员互动联盟

能让程序员心花怒放的七种礼物

一般偏感性的人都容易用小礼物哄的开心,比如从路边摘一朵送给一个小姑娘,她会开心好几天。但对一些偏理性的人,就很难用一个简单的礼物打动他了。以开发软件为职业的程...

3377
来自专栏玉树芝兰

知识工作者的便携利器

工业时代,工作时间工作,休息时间休息。现在你如果能25分钟不看手机、不刷朋友圈专注读书(哪怕是小说),那就叫完成了一个番茄钟,是非常难得的成就了(保守地说,也足...

702

扫描关注云+社区