【直播】我的基因组 38:我得了艾滋病?我是暴躁狂?

请原谅我这次当了标题党,本题应该是,把vcf格式的变异文件注释上gwas的研究结果

这个同样是为了丰富我们的variation文件,加强我们队每个变异位点的变异情况的理解。前面我们也刚刚提到一个GWAS研究rs2932989位点在我身上得到了验证,跟OME显著相关的。

GWAS数据里面收集这这些年来, 所有已经发表的GWAS研究的结果,最重要的结果,就是哪些位点与哪些性状是有联系的,而且联系的概率是多少。当然,这个主要是基于人群的,如果发表的研究主要是基于西方人种,对我这个黄皮肤可能参考价值不大。但是可以先做着试试看。

比如前面我们就提到了大批量的中国人群的GWAS研究表明,STAT4上的rs7574865和HLA-DQ的 rs9275319 --》中国人群中乙型肝炎病毒(HBV)相关肝细胞癌(HCC)遗传易感基因,我就自己分析了一下我的测序数据里面有没有这个变异,但是那样的分析只能一个个来,效率毕竟太低了,我得把GWAS研究全部肉眼看一遍,然后一个个去我的数据里面寻找,是否出现了该变异。

这时候,我们作为一个高效的数据分析人员,就得拿出不一样的能力,

http://snpeff.sourceforge.net/SnpSift.html#gwasCatalog

GWAS研究现状及资源下载 (http://www.bio-info-trainee.com/719.html)

但是,这个软件总是在这个分析要点上面报错,我很郁闷呀!其实这个我自己写一个脚本也可以做,但是重复造轮子不是很划算,而且保证正确性更耗时。但是作为一个一流的生信工程师,我当然不会在这样的阴沟里翻船呀,大不了换一个工具呗。反正snpEFF和ANNOVAR几乎在功能上面没有本质的区别。

生信必备技能之首,就是Google搜索啦,结果如下:

http://annovar.openbioinformatics.org/en/latest/user-guide/region/#identify-variants-reported-in-previously-published-gwas

自己看前面的帖子下载并且按照ANNOVAR咯,然后

~/biosoft/ANNOVAR/annovar/annotate_variation.pl -build hg19 -downdb gwasCatalog humandb/ ## 下载gwas数据库

~/biosoft/ANNOVAR/annovar/convert2annovar.pl -format vcf4old realign.vcf >realign.annovar

~/biosoft/ANNOVAR/annovar/annotate_variation.pl -regionanno -build hg19 -out ex1 -dbtype gwasCatalog realign.annovar humandb/

不知道为什么,annovar下载的是UCSC的gwas总结信息,比snpEFF下载的NCBI要多很多。

## 然后执行annovar的两个标准步骤就好了;注释好的文件就是 ex1.hg19_gwasCatalog,共11041行记录,这次真的吓到我了!

随便打开浏览了一下,我勒个去,扑面而来的都是各种疾病呀,居然还让我看到了这货:

gwasCatalog Name=AIDS progression 1 146560564 146560564 G A hom

纳尼!!!艾滋病????还是纯合突变!!!

这我果断不信呀!!!

这破工具注释结果页太简陋了,我打开gwas文件看了看:

1703 chr1 146560563 146560564 rs11239930 21502085 Troyer JL 2011-05-01 J Infect Dis Genome-wide association study implicates PARD3B-based AIDS restriction. AIDS progression 755 European ancestry seroconverter cases NA 1q21.1 NBPF14 rs11239930-? 0.58 2E-6 1.47 [NR] Affymetrix [700022] N

还好我英文不错,上面写的是这个突变有利于限制艾滋病的发展!某种意义上的超人!

感兴趣的可以查看原文:http://www.nature.com/ng/journal/v36/n6/full/ng1369.html

或去看看这个变异位点的详细信息:http://www.snpedia.com/index.php/Rs11239930

我总结了一下我的gwas注释结果:

cut -f 2 ex1.hg19_gwasCatalog |sort |uniq -c |sort -k1 -n -r |less

排在第一的赫然是,思覺失調症(Schizophrenia),精神疾病的一種。其特徵為患者出現異常的社會行為!(我表示一脸懵逼!)

我已经尽可能的不暴露自己的真实信息了,大家对这个注释结果有没有什么想法呢?

其实GWAS并不是以疾病风险因子为主,很多研究就是纯粹的探索一些表型,比如身高体重呀,跳高跳远呀,各种鬼扯的天赋基因都来自于此~

里面居然有58个位点是躁狂抑郁症Bipolar Disorder -(科学网—我们都是潜在的bipolar患者- 廖晓琳的博文http://blog.sciencenet.cn/blog-349936-584626.html )

克隆氏症(又名克羅恩氏症, Crohn's Disease ),慢性腸道炎的一種!

其余的我都不想翻译了,蓝瘦香菇!

而且,极端的不靠谱,我不否认这些研究结果里面有真实的信息,但是假阳性太多了,而且你们看我博客里面对于这类研究发表的杂志探索,就明白了。大部分是在plos ONE杂志上面发表的,那个神刊!!!

文:Jimmy

图文编辑:吃瓜群众

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

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

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

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏知晓程序

深夜爆料!微信小程序之间可以互相跳转了

果不其然,今天,知晓程序(微信号 zxcx0101)无意中发现,小程序可以跳转其他小程序了!

411
来自专栏FreeBuf

OpenATS续篇:搭建自己的卫星地球站

? 继上篇文章:DIY天线自动追踪系统OpenATS 的发表已经一年多了,这一年多的时间都没有更新这个小项目,因为我本人真的是太懒了。 其实OpenATS存在...

2756
来自专栏FreeBuf

「闭塞网络」攻击思路总结

如何从不联网的电脑中盗取1KB名叫《新建文本文档》内容为 test 的txt文件?实际上这个是FreeBuf小酒馆里面的某个问题,当然,回答也很喜感….. ?...

2115
来自专栏用户2442861的专栏

Linux网络编程书籍

原文转载:http://blog.csdn.net/21aspnet/article/details/6694485

963
来自专栏ThoughtWorks

程序员比谁都浪漫。

【2014.2.14 情人节快乐!】 以下不同的示爱方式,都来自程序员。 谁说程序员不浪漫!你最喜欢哪一款? (图片:网络) ? ? ? ? ? ? ? ? ?...

34711
来自专栏Jerry的SAP技术分享

金庸和古龙,Netweaver和微服务,以及SAP Hybris Revenue Cloud

这周Jerry在长沙客户现场待了几天,感谢易总和彩亮的款待。终于有机会和关注这个公众号的一些CRM顾问们进行线下互动,感觉很不错。得知公众号里某些文章帮助顾问们...

1483
来自专栏AI科技评论

学界 | FPL 2017最佳论文:如何对FPGA云发动DoS攻击?

AI科技评论按:第27届现场可编程逻辑与应用国际会议(The International Conference on Field-Programmable Lo...

36810
来自专栏企鹅号快讯

程序员一年写百万行代码是什么体验?这肯定是个Bug

2017 骄傲和遗憾都将成为过去,2018 元旦将至,春节还会远吗? 年年岁岁没对象,岁岁年年撸代码,或许这一年你还在new对象,但是依旧没有对象,让我们情场失...

1918
来自专栏何俊林

Android TV开发总结(六)构建一个TV app的直播节目实例

近年来,Android TV的迅速发展,传统的有线电视受到较大的冲击,在TV上用户同样也可以看到各个有线电视的直播频道,相对于手机,这种直播节目,体验效果更佳,...

26710
来自专栏大数据文摘

可视化:覆盖全球的网络攻击如何展现?

1656

扫码关注云+社区