sgRNA脱靶位点检测

很多在线的gRNA脱靶检测软件,仅限选定的物种,所以对非模式物种来说并不友好。本文选取两个相对简单的脱靶检测软件,CHOPCHOP和CasOT进行脱靶位点的检测。前提是物种需要有公布的基因组序列或注释信息。

  1. CHOPCHOP

    首先在官网界面,输入基因名字或Fasta序列,选择相应的物种,然后点击运行即可搜寻候选的sgRNA。在结果文件页面,包含了10个候选的sgRNA,单击sgRNA进入单独界面,会详细罗列出在全基因组范围内潜在的脱靶位点。


    image.png
image.png

找寻脱靶位点基因组序列:
首先在NCBI上输入脱靶定位的序列,如NW_011876316,


F)NVUTM%_TG309{AO8YPK~5.png

在基因组中肉眼很难观测到,可以下载下来全序列然后MEGA寻找,当然这些都是最传统的办法,比较便捷的是可以通过NCBI提供的Genome Data Viewer直接查找而无需下载。在页面右侧单击Show in Genome Data Viewer(新建页面,保留原始页面)


image.png

image.png

然后选择Tools > Search,输入脱靶的核酸序列,然后找到具体对应基因组的位置,如下图所示:


image.png

image.png

本处比对到NW_011876316的945659——945681,然后回到NCBI比对的界面,在右侧的Change region shown,输入选择的区域,这里在我输入的序列在脱靶序列上下游个1000bp,便于后面进行基因扩增检测。


~3XU@A4$7%]``N7DUU23$X3.png

然后下载下来即可。

  1. CasOT

CasOT支持本地windows操作,且没有物种限制。软件是用Perl脚本,所以首先安装Perl,然后就是windows命令行本地化运行:

win+R
cmd
# 进入软件所在目录位置
cd G:\01_J.F.XIE\CasOT-1.0
# 检测软件是否可运行
.\CasOT-1.0\casot.pl -h

# 1st. 启动可执行perl脚本
.\CasOT-1.0\casot.pl
#2sn. 输入关键参数
-t='' -g='' -e='' -o=tab -s=2 -n=4 -p=A

注意输入的路径!!!
当前所在路径和输入文件路径(当前文件夹用:.\)

--- Options ---

        -h: show this help

        -m=: single, paired, or target (default: single)
        -t=
        -g=
        -e=
        -o=: csv or tab (defualt: csv)

        -s=: 0-6 (nt) (default: 2; >4 is not suggested)
        -n=: 0-255 (nt) (default: 255)
        -p=: A, B, C or N (default: A)
                Level A: NGG
                Level B: NGG + NAG
                Level C: NGG + NAG + NNGG
                Level N: (no limit)

        -d=: 0-1000 (default: 100)

        -r=: yes or no (default: yes)
        -l=: two number between 18 and 30 (default: 19-20)

参数详见文档说明,或-h查看。
运行大约5分钟即可出结果。
结果说明:


image.png

图片.png

你可能感兴趣的:(sgRNA脱靶位点检测)