热点新闻
全基因组CNV分析4. 拷贝数的增减
2024-10-19 11:53  浏览:1764  搜索引擎搜索“早勤网”
温馨提示:信息一旦丢失不一定找得到,请务必收藏信息以备急用!本站所有信息均是注册会员发布如遇到侵权请联系文章中的联系方式或客服删除!
联系我时,请说明是在早勤网看到的信息,谢谢。
展会发布 展会网站大全 报名观展合作 软文发布

拷贝数的增减官方会被称作copy number gains and losses。拷贝率(copy ratio)通常和真实的拷贝数的关联度一般会受到肿瘤细胞纯度(pure),基因倍数体(ploidy),还有就是亚克隆群体的大小。因为有这些不确定的情况,所以CNVkit的最终结果只提供了log2拷贝率的计算值。于此同时提供了多种计算拷贝数绝对值或者相对增减数量的方法。

在2倍体里,假设纯度是100%,一个纯合样本的CN增加了1,那么拷贝率就变成了2/3,用log2计算的话就是log(3/2)=0.585。相反,如果CN减少了1,那么拷贝率就变成了1/2,log2(1/2)=-1.0。

在CNVkit的Diagram里(如下图),为了避免结果图过于杂乱,默认阈值为0.5来添加基因标签。





diagram

在使用genemetrics指令的时候,需要根据样本的纯度,倍数体来选择合适的阈值。可以从0.1,0.2来开始过滤基因数量。

要计算每个片段数据可信度的时候,可以使用指令segmetircs。还可以在指令里添加-ci来计算可信范围。

最后是call指令,call指令会把log2值转换成CN的绝对值。当样本纯度和倍数体的确定的时候就可以这么做。

.cns文件也可以通过export指令被转换成BED或者VCF格式。BED,VCF格式的结果会给出每一个片段的CN的绝对值。当然用export之前最好先根据纯度和倍数体用call来设置一下阈值。

发布人:0d57****    IP:124.223.189***     举报/删稿
展会推荐
让朕来说2句
评论
收藏
点赞
转发