其实前面我们已经分享了MiXCR,还有igblast,这两个免疫组库上游分析软件已经够用,如下:
理论上不应该再介绍过多软件和流程,避免增加大家的认知负担,但是看到一个很新的文章发表在NC杂志,时间是11 April 2019,标题是 Multi-region sequencing unveils novel actionable targets and spatial heterogeneity in esophageal squamous cell carcinoma,该研究对33个食管癌病人进行了多位点取样,其中10个病人的64个样品进行了TCR的免疫组库测序,就使用了IMonitor作为分析工具。
测序数据在:PRJNA511368 ,也是可以下载的,但是工具IMonitor其实引用并不多,是深圳华大基因员工开发。
IMonitor文章
首先下载和安装IMonitor
该软件的GitHub链接是:https://github.com/zhangwei2015/IMonitor
里面有一个压缩包,是打包好的软件,依赖于Linux系统的perl语言环境,
教程很简单:
FQ_run.sh:
perl ../IMonitor.pl -a data/XHS_1.fq.gz -b data/XHS_2.fq.gz -A1 data/1.adapter.list.gz -A2 data/2.adapter.list.gz -o . -n XHS -T TRB -k 100 -r ../Ref/TRB -d -m -Rs /opt/blc/genome/biosoft/R/bin/Rscript
FA_run.sh:
perl ../IMonitor.pl -i data/XHS.merged.fa.gz -o . -n XHS -T TRB -k 100 -r ../Ref/TRB -d -m -Rs /opt/blc/genome/biosoft/R/bin/Rscript
下载和安装该软件的代码是:
mkdir -p ~/biosoft/IMonitor
cd ~/biosoft/IMonitor
perl -v
# This is perl 5, version 18, subversion 4 (v5.18.4)
git clone https://github.com/zhangwei2015/IMonitor
我们这里仍然是使用在前面我们认识的免疫组库测序数据,是人类的,MiSeq测序仪,PE300测序策略,TRB,DNA测序,进行示范。