vep怎么使用？新手快速上手指南3步搞定

昨天帮师弟处理基因组数据，折腾vep差点原地爆炸。这玩意儿对生信新手太不友好了，文档写得跟天书似的，今天必须把踩坑记录捋清楚。

第一步：装软件就给我当头一棒

打开官网下载页直接懵圈，光安装方式就有三种。想着docker总该省心，结果docker pull ensemblorg/vep敲下去，2GB的文件下载慢得能泡三壶茶。转头试conda安装，conda install -c bioconda ensembl-vep输完报错提示缺perl库，气得我差点砸键盘。老实走perl脚本路线：

• curl -O 官网地址（这步就卡了十分钟找下载链接）
• tar zxvf *解压时权限报错
• sudo chmod -R 777 vep_folders才搞定

第二步：准备数据文件更抓狂

文档说需要缓存文件，./vep --cache --dir_cache ./cache一运行就报错“Can’t find cache species”。翻了三页论坛才搞明白要指定物种名，改成./vep -i * --species homo_sapiens --cache，结果又弹“No fasta file”。连夜下载35G的参考基因组，硬盘直接红了警告。

更绝的是在线模式试了下./vep --port 3337，等五分钟弹出“连接被外星人劫持”（是防火墙拦截）。把整个命令行写成脚本才跑通：

bash

./vep -i my_* \

--species human \

--cache \

--dir_cache /home/ubuntu/cache \

--fasta /data/* \

--force_overwrite \

--vcf

第三步：结果解读才是终极考验

盯着输出文件里的CSV格式发呆，“Allele”、“Consequence”这些字段还能蒙个大概，到“PolyPhen”和“SIFT”预测分值彻底抓瞎。又翻出三篇教程对比着看：

• 红色警报是transcript_ablation（直接破坏基因）
• 错义突变看REVEL评分>0.75要警惕
• 内含子区域出现splice_region_variant得重点查

完事儿用grep ‘deleterious’ *筛高危突变，结果文件比师弟的头发还稀疏——合着熬夜配环境就分析出七个有效变异，其中三个还是数据库里标着“常见无害”的。

等跑完所有样本已经天亮，看着命令行里满屏的WARNING: 跳过xx条数据陷入沉思...这玩意儿说三步搞定纯属诈骗，光安装就能劝退八成新手，但搞熟之后确实能少写五百行python解析代码（自我安慰完毕）。