分析ABSOLUTE数据结果

解析ABSOLUTE软件一文中我简单介绍了ABSOLUTE软件,因为ABSOLUTE解的多样性,这篇文章主要聚焦于ABSOLUTE结果的解读,对象是结果图。

主要的参考资料来源于ABSOLUTE官方文档:http://software.broadinstitute.org/cancer/software/genepattern/analyzing-absolute-data

纯度/倍性图以及模型打分

image
image

第一个图形展示了符合要求的多个解,因为纯度和倍性是关联的,所以每一个解可以看对一个纯度/倍性对。在第二个图形中展示了3个打分最高的解,每个解通过不同的模型进行评估,主要是SCNA,核型模型以及它们的整合模型。如果用户提供了突变数据,还有一个基于SSNV模型的评估。

第一个图形大致理解是很轻松的,但细看的话很多都很费解,比如\sigma的意思,图中虚线又是什么意思。我现在也没搞懂,所以就不误人了。

备选的拷贝轮廓

image

这个图显示了前面3个解对应的拷贝轮廓,像纯度,倍性,异质性片段的比例以及还有一些看不懂为什么是负数的参数值。其实有个参数跟我们输入时密切相关,就是\sigma_h,它标定了样本允许的最大的方差。

等位基因比例图

image

这个图如果输入了突变数据应该就会有。等位基因比例可以解释为每个cancer细胞平均的等位基因拷贝数,也称为多样性(multiplicity),它可以揭示亚克隆突变。这其实跟利用VCF文件计算vaf分布差不多,不过这里用copy number进行了校正!

这个图的解释可以参照【直播】我的基因组81:看看我的vcf文件的vaf分布情况这篇文章进行理解。

正常人的二倍体基因组位点只有杂合或者纯合两种情况,对于纯合那么vaf必然是1,对于杂合,必然是0.5。但是现实测序得到的结果远比这要复杂,尤其是测序深度不够的时候。因为测序本身具有随机性,而且还有很多系统误差。理想情况也只能像是扔硬币。

image

而文献里面对TCGA里面的癌症样本的somatic mutation的vaf

image

可以看出tumor里面的vaf分布其实已经不再是扔硬币那样的概率了,对于杂合位点来说。

原因很多,首先tumor不一定是单纯的二倍体了,其次tumor样品一般来说本身异质性高,而我们测序是混合多个细胞的,有一些突变有一些并不突变。而且纯合的somatic mutation几乎没有,因为somatic mutation是tumor过滤了normal后留下来的变异位点,不是遗传多样性,突变这个过程既然是后天产生的,就很难保证取样部分的几百万个细胞全部突变了。

选择解的一些建议

分析ABSOLUTE的结果高度依赖你的数据和你对你数据的理解。这里列出一些指南帮助你选择好的解:

  • 简单至上!在考虑所有的模型信息后,选择最简单的解。

    • 如果两个模型在拟合效果和其他标准上几乎一致,那么选择更简单的那个——例如,大多数正常的基因组是2倍体或3倍体(这里指的是癌症)
  • 每个拷贝数的峰值应该在一个整数下面。

    • 整数反映了基因组的状态。1是单倍体,2是双倍体,3是三倍体等

      image
  • 峰值通常为一个更容易得到的数值(例如1,2,4)

    • 相比于3倍体,4倍体更容易,因为生物学上基因组翻倍比扩增部分更容易。这并不是说3倍体不可能发生,只是说从生物学上理解应该更少见。
  • 拷贝数底部(左下角)的峰值应该接近0

    • 0显示了拷贝的缺失,它比一些中间值更简单

      image
  • SSNV等位基因比例图的\alpha/2线应该将最高的SSNV等位基因比例划分为两半

这些规则有时候是相互冲突的,所以读者在实际选择解时应该从自己理解的角度去判断。

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 216,470评论 6 501
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,393评论 3 392
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 162,577评论 0 353
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,176评论 1 292
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,189评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,155评论 1 299
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,041评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,903评论 0 274
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,319评论 1 310
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,539评论 2 332
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,703评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,417评论 5 343
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,013评论 3 325
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,664评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,818评论 1 269
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,711评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,601评论 2 353

推荐阅读更多精彩内容