【直播】我的基因组 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 条评论
登录 后参与评论

相关文章

来自专栏企鹅号快讯

程序猿必须知道的一些有用的网站

来自:开源中国 原文:https://github.com/sdmg15/Best-websites-a-programmer-should-visit 在学习...

20810
来自专栏知晓程序

本周最新 10 款小程序,快来选出你最爱的那个 | 晓榜 #21

如果感到无聊,不妨用「微博热门」小程序,看看最流行的微博,热门资讯、搞笑视频,总有一款是你爱的。

461
来自专栏腾讯大讲堂的专栏

自由测试人 Jarod 的一天

上午10:05 五道口漫咖啡,Jarod摆弄着新淘来的Nexus5手机,时而饮一口桌上的焦糖拿铁,间或偷眼瞄一下邻桌的长腿妹子。 上午10:30 Alliso...

1765
来自专栏云计算D1net

SDN基础分析浅谈

SDN(Software Defined Network)是个有意思的概念。ONF(Open Network Foundation)这样定义SDN: In th...

2435
来自专栏Bug生活2048

程序猿的无所不能

我是一名程序员,长相一般,水平一般,什么都一般的普通程序员。可在朋友眼中,我是无所不能的,只因为我是程序猿……

711
来自专栏腾讯移动品质中心TMQ的专栏

陪你度过漫长岁月:WiFi管家测试一纸芳华诉流年

就在两小时前,WiFi管家刚刚发布了3.0版本。做完最后一个专项汇报,坐在屏幕前揉了揉肩膀松了口气。一种特别的情绪突然蔓延开来。

750
来自专栏程序人生

那些年,我追过的语言

程序君也年轻过,年轻的代价就是盲目追随。 从MS-DOS6.0开始,程序君就是微软的狂热拥趸。 这种狂热自win95走上高潮(有谁还记得win95光盘里带的Go...

3295
来自专栏IT 指南者专栏

推荐 6 款实用好玩的小程序

想必很多人都知道,现在小程序非常火,特别是在微信超过 10 亿日活的超强生态之下,各种各样的小程序目不暇接。

1032
来自专栏城域网光通信

基于反向复用技术的100G光传输系统

基于反向复用技术的100G光传输系统

1340
来自专栏EAWorld

DevOps与合规性:鱼和熊掌兼得指南

编者按:很多行业身处强力监管领域,因而格外强调合规性。反映在IT上就是开发、部署和运维等规范(比如开发团队不能碰生产日志)的不可或缺。本文中提到的一些方法(如自...

2994

扫码关注云+社区