PennCNV官网:https://penncnv.openbioinformatics.org/en/latest/
PennCNV软件的文献: DOI: 10.1101/gr.6861907
本次基于这个2007年发布的文献进行学习,后续根据官网的说明进行进一步学习。(所以文中的部分言论是以2007年的视角进行论述的)
杂交芯片(array-CGH)检测CNV的技术缺陷是仅限于检测数十或数千bp的CNV,而SNP芯片能够更加准确的检测SNP(更灵敏,并且SNP芯片不需要进行PCR扩增,不会造成基因组偏倚,所以获得的SNP检测信息更加真实),由于SNP芯片中SNP位点的距离可达几至几十Kb(不同密度的SNP芯片其距离范围不同),所以SNP芯片相比于杂交芯片更有利于发现大范围的CNV。
PennCNV是专门用于Illumina平台芯片数据进行CNV分析的,对于Illumina平台自带的CNV分析软件------BeadStudio软件的LOH-plus模块,LOH-plus模块通过在沿着染色体的滑动窗口中计算SNP的B等位基因频率的模式来鉴定拷贝数变化。虽然实现简单,但滑动窗口方法对于检测到的CNV具有有限且相对粗糙的边界分辨率。而PennCNV通过使用隐马尔可夫模型(HMM),纳入了log R比率和B等位基因频率,再加上基于贝叶斯模型将系谱等信息纳入矫正,能够显著提高CNV检测的分辨率(更加清晰的CNV编辑的鉴定)。
常规的算法将CNV的类型分为了三类:“loss,”“normal,” and “gain”,而PennCNV将CNV分为了五种类型(由六种可能的情况造成)(下表)。


为了充分利用每个SNP的所有可用信息,PennCNV将几个组件一起合并到隐马尔可夫模型(HMM)中,包括LRR,BAF,相邻SNP之间的距离和B等位基因的群体频率。LRR和BAF值都可以从BeadStudio显示和导出,前提是存在具有每个SNP的规范聚类位置的适当聚类文件。相邻SNP之间的距离决定了它们之间具有拷贝数状态变化的概率。每个SNP具有被称为A和B等位基因的两个等位基因,因此我们使用术语“B等位基因的群体频率”将其与测量等位基因强度比的BAF术语区分开。所有SNPs的B等位基因的群体频率值是从具有混合种族背景和正常表型的大量个体汇编的;然后确定每个拷贝数状态的拷贝数基因型的可能性。
我们认为后代中的绝大多数CNV是从父母遗传的(Locke et al.2006),因此在后代中推断但在亲本中未检测到的CNV的分数(CNV-NDP)可用作假阳性和假阴性率的复合量度。因此,我们检查了由三个不同平台产生的CNV调用中HapMap CEU + YRI后代中CNV-NDP的分数。

CNV学习2(illumina芯片分析CNV的主流软件------PennCNV)相关推荐

  1. RT-1052学习笔记 - GPIO架构分析

    RT-1052学习笔记 - GPIO架构分析 一.RT1052之GPIO与IOMUX 二.GPIO架构分析 分析GPIO的外设框架 库函数操作 1.初始化MUX 2.初始化PAD 3.初始化GPIO外 ...

  2. ARM学习笔记--LPC1788引脚分析

    http://zgmgypb.blog.163.com/blog/static/962028192012730104640980/ RM学习笔记--LPC1788引脚分析 一.电源分析 要弄清楚单片机 ...

  3. 借由ARM CORTEX-M芯片分析C程序加载和存储模型

    https://zhuanlan.zhihu.com/p/22048373 写文章 借由ARM CORTEX-M芯片分析C程序加载和存储模型 王小军 1 年前 阿军最近在忙着血氧手环嵌入式系统的技术预 ...

  4. 免费教材丨第56期:《深度学习导论及案例分析》、《谷歌黑板报-数学之美》

    小编说  离春节更近了!  本期教材        本期为大家发放的教材为:<深度学习导论及案例分析>.<谷歌黑板报-数学之美>两本书,大家可以根据自己的需要阅读哦! < ...

  5. 2016年大数据Spark“蘑菇云”行动代码学习之AdClickedStreamingStats模块分析

    2016年大数据Spark"蘑菇云"行动代码学习之AdClickedStreamingStats模块分析     系统背景:用户使用终端设备(IPAD.手机.浏览器)等登录系统,系 ...

  6. 《深度学习导论及案例分析》一2.11概率图模型的推理

    本节书摘来自华章出版社<深度学习导论及案例分析>一书中的第2章,第2.11节,作者李玉鑑 张婷,更多章节内容可以访问云栖社区"华章计算机"公众号查看. 2.11概率图模 ...

  7. LIVE555再学习 -- testOnDemandRTSPServer 源码分析

    一.简介 先看一下官网上的介绍: testOnDemandRTSPServercreates a RTSP server that can stream, via RTP unicast, from ...

  8. Qt学习笔记,再次分析EVA源码之后得出的结论-QListView,QListViewItem(Qt3);Q3ListView,Q3ListViewItem(Qt4)...

    Qt学习笔记,再次分析EVA源码之后得出的结论-QListView,QListViewItem(Qt3);Q3ListView,Q3ListViewItem(Qt4) 今天再次分析了Eva的源码,也看 ...

  9. OpenCV与图像处理学习九——连通区域分析算法(含代码)

    OpenCV与图像处理学习九--连通区域分析算法(含代码) 一.连通区域概要 二.Two-Pass算法 三.代码实现 一.连通区域概要 连通区域(Connected Component)一般是指图像中 ...

最新文章

  1. Java任务调度框架Quartz
  2. oracle xe gentoo,Oracle在gentoo下安装
  3. 退火算法 贪婪算法_算法贪婪策略简介
  4. 【计算机组成原理】内存
  5. 解决ssh登录后闲置时间过长而断开连接
  6. androidStudio快捷键概览
  7. JAVA之outofmemory
  8. 华为机试:提取不重复的整数
  9. AI CC呼叫中心源码
  10. 少年Pi的奇幻漂流-我们的后台自动化发布方案
  11. 面试中常见的问题总结
  12. blastn 输出结果每列啥意思_本地blast的详细用法
  13. RDS数据库cpu过高分析
  14. 液晶面板里面有些什么配件_液晶显示器参数和部件详细介绍
  15. 国标28181:什么是RTP协议与RTCP协议
  16. 主页劫持该怎么办?解决主页劫持的方法
  17. It is time you did something (注意这里要用过去式)
  18. sql复合索引使用和注意事项
  19. 销售管理之调拨销售管理
  20. 软件项目管理 7.5.项目进度模型(SPSP)

热门文章

  1. 金蝶专业版怎么反过账当月_金蝶专业版怎么反过账
  2. hiredis异步操作模型
  3. 短视频用户陷入视觉疲劳,AI能否救场
  4. linux版teamview15.5.3
  5. 【转】Java8 函数式编程详解
  6. 用友开发者社区 - http://udn.yonyou.com
  7. web前端面试十大 难点问题
  8. AvsP介绍——非常好用的AviSynth辅助工具
  9. 在阿里上班,是种什么体验?
  10. MATLAB 2008a 激活成功后还要激活(许可证过期)