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

吓死我了!

前面我们讲到了可以把突变文件注释到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

图文编辑:吃瓜群众

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

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

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

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏跟着阿笨一起玩NET

[C#]工具类—FTP上传下载

  不错的文章:http://www.cnblogs.com/greatverve/archive/2012/03/03/csharp-ftp.html

1181
来自专栏成长道路

JDBC动态SQL语句连接orcale数据库的工具类

import java.sql.Connection; import java.sql.DriverManager; import java.sql.P...

2480
来自专栏积累沉淀

Hive2.0.0操作HBase 1.2.1报错解决

首先看错  org.apache.hive.service.cli.HiveSQLException: Failed to open new session: ...

2339
来自专栏谈补锅

记录C#常用的代码片段

using Newtonsoft.Json; using Newtonsoft.Json.Linq;

922
来自专栏菩提树下的杨过

MSDN官方的ASP.Net异步页面的经典示例代码

示例1.演示异步获取一个网址的内容,处理后显示在OutPut这一Label上 using System; using System.Web; using S...

1985
来自专栏张善友的专栏

Using sqlite with .NET

The other day I found that there is a .NET wrapper for sqlite. sqlite is a very ...

2278
来自专栏君赏技术博客

写了一个Swift版本的提示HUD 暂时没找到合适的可以用一下

灰常灰常感谢@地选之猿提出这篇文章代码已经不能使用。这个库是当时学 Swift 才写的已经属于 Swift2.0代码,现在已经修复完成。 因为之前的邮箱不能使...

623
来自专栏听雨堂

想修改CSS

      下载了一个“通用”的CSS文件,本来想偷懒的,结果发现有问题,就是它用的颜色是变量定义的,无法识别。我又找不到在哪里可以定义。 BODY{     ...

20410
来自专栏互联网开发者交流社区

WinForm之窗体应用程序

1853
来自专栏我和未来有约会

silverlight向服务器post数据类

using System; using System.Net; using System.Windows; using System.Windows.Co...

1965

扫码关注云+社区