前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >【直播】我的基因组 36:这些可能是somatic突变的位点究竟是什么?

【直播】我的基因组 36:这些可能是somatic突变的位点究竟是什么?

作者头像
生信技能树
发布2018-03-08 10:11:47
9010
发布2018-03-08 10:11:47
举报
文章被收录于专栏:生信技能树

吓死我了!

前面我们讲到了可以把突变文件注释到dbsnp数据库,而dbsnp数据库信息非常丰富,比如我们可以看

##INFO=<ID=SAO,Number=1,Type=Integer,Description="Variant Allele Origin: 0 - unspecified, 1 - Germline, 2 - Somatic, 3 - Both">

这个SAO 信息很重要,请自行搜索了解变异为什么区分成germline和somatic。下面我用脚本统计一下SAO的信息:

cat realign.dbsnp.vcf |perl -alne '{/(SAO=\d)/;print $1}' |sort |uniq -c

我有近76万的位点没有在dbsnp数据库里面,剩余的在dbsnp里面的,有1610个可以肯定是germline的mutation,还有142个是可能是somaticmutation,这些位点很可怕,somatic mutation意思是在人体成长过程中的后发突变,很大可能是癌症或者疾病相关突变!

我写脚本看了看这142个可能是somatic突变的位点所在的基因:

cat realign.dbsnp.vcf |perl -alne '{print if /SAO=3/;}' >dbsnp_both_germ_somatic.vcf

cat dbsnp_both_germ_somatic.vcf |grep GENEINFO |perl -alne '{/(GENEINFO=.*?);/;print $1}' |sort -u

MUC开头的我大概知道是什么意思,因为我有中耳炎,它突变我不觉得稀奇!

BRCA1,2都突变了,我觉得可能是科研界过度研究,这个我暂时也不害怕!

CD分子的突变,预示着免疫机理可能是有点麻烦,但毕竟只是可能!

其余的基因,我都不认识, 因为我知识水平有限!

当然,我只定位到了基因,没有看这个位点的突变对这个基因具体的影响!但是有研究表明它们可能是somatic的突变,其实就有一点风险了。(这就是为什么大家做基因检测看到的报告,都是有什么风险,很狗血,因为连负责这个检测报告的生信工程师也没办法给出肯定答案啊!)

有点不开心,不想看下去了。

文:Jimmy

图文编辑:吃瓜群众

本文参与 腾讯云自媒体同步曝光计划,分享自微信公众号。
原始发表:2017-01-01,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 生信技能树 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
相关产品与服务
数据库
云数据库为企业提供了完善的关系型数据库、非关系型数据库、分析型数据库和数据库生态工具。您可以通过产品选择和组合搭建,轻松实现高可靠、高可用性、高性能等数据库需求。云数据库服务也可大幅减少您的运维工作量,更专注于业务发展,让企业一站式享受数据上云及分布式架构的技术红利!
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档