结构变异SV分析(WGS/WES)和融合基因检测(RNA)

69种工具检测SV对比 https://www.jianshu.com/p/85ff211b4b4e

https://zhuanlan.zhihu.com/p/87408242?utm_campaign=rss&utm_medium=rss&utm_source=rss&utm_content=title
最终决定使用wham和Manta组合,其实最优组合应该是wham,inGAP-SV组合或者inGAP-sv,PBHoney-NGM组合。

安装wham(现在实际是whamg) github

git clone --recursive  https://github.com/zeeev/wham.git
cd wham
make

安装Manta github

#强烈推荐使用二进制版本,不用编译,少了很多麻烦(虽然是centos6但是其他版本的linux也能使用)
wget https://github.com/Illumina/manta/releases/download/v1.6.0/manta-1.6.0.centos6_x86_64.tar.bz2
tar -jxvf manta-1.6.0.centos6_x86_64.tar.bz2
cd manta-1.6.0.centos6_x86_64/bin
./configManta.py -h
#Manta是基于python2,如果你的默认python是python3,可以这样调用。
python2 configManta.py

##我直接把上述命令和路径加进path
alias Manta="python2 /share/home/software/manta/manta-1.6.0.centos6_x86_64/bin/configManta.py"
#直接运行Manta 就行了。

运行Manta

Manta \
--bam NA12878_S1.bam \
--bam NA12891_S1.bam \
--bam NA12892_S1.bam \
--callRegions hg38.bed \
--referenceFasta hg38.fa \
--runDir ${MANTA_ANALYSIS_PATH}
参数说明:

要求输入的bam和fa必须有samtool的fai索引文件,bam还需要sort过。
--callRegions 适用于指定染色体,通过这种方式来只分析chr上的,而不分析scaffold。
hg38.bed内容如下

chr1    0   248956422
chr2    0   242193529
chr3    0   198295559
chr4    0   190214555
chr5    0   181538259
chr6    0   170805979
chr7    0   159345973
chr8    0   145138636
chr9    0   138394717
chr10   0   133797422
chr11   0   135086622
chr12   0   133275309
chr13   0   114364328
chr14   0   107043718
chr15   0   101991189
chr16   0   90338345
chr17   0   83257441
chr18   0   80373285
chr19   0   58617616
chr20   0   64444167
chr21   0   46709983
chr22   0   50818468
chrX    0   156040895
chrY    0   57227415
chrM    0   16569
运行whamg

目前无法成功安装wham

融合基因检测
http://www.bio-info-trainee.com/7258.html
Comparative assessment of methods for the fusion transcripts detection from RNA-Seq data 2016年 测评12个工具
Accuracy assessment of fusion transcript detection via read-mapping and de novo fusion transcript assembly-based methods 2019年测评 23个工具论文
决定使用STAR-fusion检测融合基因。
STAR-fusion是基于STAR比对的结果进行检测融合基因的。

#安装START-fusion
conda install -c bioconda star-fusion

注意:STAR和STAR-fusion的不同版本的兼容性有较大问题,开始分析前,一定先确定使用的版本之间产出的文件是否兼容。
参考:STAR-fusion的使用
STAR-fusion检出融合基因的后续分析

WGS和RNA-seq联合分析参考文章 2019 NC https://www.nature.com/articles/s41467-019-10680-5#Sec10

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

推荐阅读更多精彩内容