SV变异检测软件(1)-breakdancer的安装及简单使用

小的Indel(Insertion 和 Deletion),指的是在基因组的某个位置上所发生的小片段序列的插入或者删除,其长度通常在50bp以下;
大的结构性变异,这种类型比较多,包括长度在50bp以上的长片段序列的插入或者删除、染色体倒位,染色体内部或染色体之间的序列易位,拷贝数变异,以及一些形式更为复杂的变异。

为了和SNP变异作区分,上述这两类变异通常也被称为基因组结构性变异(Structural variation,简称SV)。这里值得一提的是,研究人员对基因组的结构性变异发生兴趣,主要是由于这几年的研究发现:(1)虽然还未被广泛公认,但研究人员发现SV对基因组的影响比起SNP来说还要大;(2)基因组上的SV比起SNP而言,似乎更能用于解释人类群体多样性的特征;(3)稀有且相同的一些结构性变异往往和疾病(包括一些癌症)的发生相关联甚至还是其致病的诱因。

目前常用于检测SV变异的软件有很多,这里介绍下breakdancerCREST的安装和简单使用。
首先介绍下breakdancer。
这里先贴出软件GitHub的地址:https://github.com/genome/breakdancer
关于软件的说明可在GitHub上自行阅读。
sourceforge的下载地址:https://sourceforge.net/projects/breakdancer/files/
通过git clone命令或者wget命令都可以下载breakdancer的安装包,这里比较简单,略过不讲。

breakdancer的安装过程比较容易各种报错
首先先安装libgd库,具体可参考这里

然后安装以下的perl模块,安装方式同上:GD、GDGraph、GDGraph-histogram、GDTextUtil、Math-CDF、Statistics-Descriptive;也可自行通过cpanm方式自动安装

然后安装samtools,版本在0.1.7以下

最后一步才是安装breakdancer:

    cd breakdancer
    make && make install

安装成功后生成breakdancer-max的执行文件
以上各个安装过程若无root权限,请自行修改./config文件的安装路径

安装成功后进入perl文件夹下找到bam2cfg.pl文件

    perl bam2cfg.pl 

安装成功就可以看到该文件的使用说明,若失败根据提示安装对应的perl模块即可。

breakdancer检测SV变异的使用方法如下:

    perl bam2cfg.pl -v (float)sample.bam > sample.cfg

生成配置文件,配置文件如下所示:

    readgroup:sample
    platform:illumina
    map:sample.bam
    readlen:
    lib:
    num:
    lower:
    upper:
    mean:
    std:
    SWnormality:
    exe:samtools view

接着使用breakdancer-max鉴定结构变异:

    breakdancer_max  sample.cfg > sample.out

以上各步骤中请自行将安装路径加入环境变量中。breakdancer的检测量很大,检测时间会很久,建议加上nohup &命令放入后台运行。后面再进行CREST的介绍。

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

推荐阅读更多精彩内容

  • 什么是基因组变异 基因组变异是一个定义比较模糊的概念. 所谓的变异是相对于一个完美的“参考基因组”而言。但是其实完...
    xuzhougeng阅读 5,055评论 0 22
  • 首先,在开始之前我觉得有必要稍微科普缓冲一下,以便不使得不熟悉生物信息或基因组的客官们疑惑。O(∩_∩)O! 基因...
    王诗翔阅读 2,250评论 0 15
  • 刘小泽写于18.8.10,补充于18.8.14-15这之间经历了第一期授课结课,(回中山办购房手续)遥墙机场读了1...
    刘小泽阅读 9,345评论 6 41
  • 虽然用简书很久了,但是从没记录过文章,今天就记录下我的处女作吧。关于近期所追剧的几点记录。那就说说最近,一直...
    傻泡泡儿阅读 1,993评论 0 0
  • 柠檬~别名黎檬子!常吃能减肥! 柠檬~维生素含量丰富,能防止和清除皮肤色素沉着,具有美容增白的作用。 柠檬~含柠檬...
    珏瑀钰阅读 405评论 2 11