位置: IT常识 - 正文

SNP注释神器——VEP(生信)(ssms注释)

编辑:rootadmin
SNP注释神器——VEP(生信) VEP使用简介1.背景知识2.vep功能2.1 可注释变异类型2.2 注释结果信息2.3 过滤筛选3.软件下载与安装3.1 安装vep的方式4.vep的使用4.1.使用基本的参数进行测试4.2 在docker下使用1.背景知识

推荐整理分享SNP注释神器——VEP(生信)(ssms注释),希望有所帮助,仅作参考,欢迎阅读内容。

文章相关热门搜索词:snpeff注释vcf,sprimg注解,snpeff注释,snpeff注释,snp和indel注释结果,snp注释结果怎么看,snp和indel注释结果,snpeff注释,内容如对您有帮助,希望把文章链接给更多的朋友!

二代测序发现的变异数量是巨大的,以人类基因组为例,1个人类基因组包含近3500000个SNV突变和1000个拷贝数变异,其中约20000-25000个变异是在编码区,10000个位点发生了氨基酸编码改变,仅有50-100个发生了蛋白截断或功能缺失,手动检查每个位点非常耗时且有些不切实际。 Variant Effect Predictor(VEP)是一个功能强大的工具集,用于编码和非编码区域的基因组变异的分析、注释和优先排序。它提供了对广泛的基因组注释集合的访问,具有各种接口以适应不同的需求,以及用于配置和扩展分析的简单参数。它是开源的,免费使用,并支持结果完全复现。 默认情况下,VEP 使用基于归一化的等位基因匹配算法来识别与输入变体匹配的已知变体。由于输入和已知变体可能具有多个替代(ALT)或变体等位基因,因此每对参考(REF)和ALT等位基因都经过归一化并独立比较以达到潜在的匹配。VCF 允许在同一行上编码多个等位基因类型,而 dbSNP 将单独的 rsID 标识符分配给同一位点的不同等位基因类型。这意味着来自同一输入变体的不同等位基因可能会被分配不同的已知变体标识符。 并且,VEP 支持使用 HGVS 作为输入。此功能目前正在开发中,并非所有 HGVS 表示法类型都受支持。完全支持与基因组 (g.) 或编码 (c.) 序列相关的符号;蛋白质(P.)符号以有限的方式得到支持,因为确定可能产生单个蛋白质变化的多种可能的潜在基因组序列变化涉及复杂性。如果无法解析特定表示法,将发出警告。

2.vep功能2.1 可注释变异类型

明确位置及序列改变的变异,全外显子组和全基因组检测出的变异,包含SNVs,insertions、deletion等。 VEP支持多种物种,也可自行构建数据库,不过使用ensemble数据库中所注释好的内容为最佳选择,可在https://ftp.ensembl.org/pub/release-108/variation/ 查看是否支持该物种。

2.2 注释结果信息SNP注释神器——VEP(生信)(ssms注释)

(1)功能注释信息,主要是把变异在基因组的位置和功能影响标出。包括基因信息和变异影响,信息来自RefSeq或Ensembl; (2)频率信息,标记群体数据库频率,包含数据库总的频率及子群体频率。这部分信息主要在疾病研究中根据疾病的流行病学研究,有针对性的筛选低频或罕见变异,为人类疾病研究提供基础。高通量测序检出的变异数目众多,对于遗传病来说,要筛选出致病变异,需要过滤掉一些常见变异,群体变异库的使用可以进一步缩小筛选范围; (3)软件预测结果信息,突变的致病性预测。可有效缩小变异范围,协助研究者更快找到致病基因; (4)疾病数据库信息,HGMD、Clinvar等数据库信息,协助用户快速定位已知致病位点;

2.3 过滤筛选

VEP提供注释结果的简单过滤,主要包含: (1)筛选功能突变,筛选出coding区和剪切区域的突变; (2)按照频率过滤,筛选出低于指定频率的突变; (3) 按照遗传模式过滤。

3.软件下载与安装3.1 安装vep的方式

(1)安装vep a. 直接使用conda来安装vep软件,省去了多余复杂的配置问题,但是由于vep涉及到的软件版本和依赖较多,因此建议新建一个虚拟环境来安装

conda create -n vep # 需要创建一个虚拟环境安装vepconda activate vep # 激活环境conda install -c bioconda ensembl-vep==108 ## 安装vpe,108是最新版本

b. 使用docker安装

docker pull ensemblorg/ensembl-vep

(2)下载vep相应的数据库 cache 数据集的下载 cache是一个可下载的文件,包含一个物种的所有转录模型,调控特征和变异数据。下载使用cache是使用VEP的最快和最有效的方法,因为在大多数情况下,只建立一个初始网络连接,大多数数据是从本地磁盘读取的。使用离线模式消除所有网络连接的速度和/或隐私。cache 只有两种类型,分别是refseq和merged。在此基础上,VEP运行时还可以使用plugin和custom参数自定义添加数据库。 注意:强烈建议下载/使用对应于Ensembl VEP安装的VEPcache版本,即VEP cache版本108应该与Ensembl VEP工具版本108一起使用。这主要是由于VEP cache(数据内容和结构)是在每个Ensembl版本中生成的,关于此版本的数据和API更新,因此cache数据格式可能在版本之间有所不同(并且与新版本的Ensembl VEP工具不兼容)。 cache 包含的信息:转录本位置,序列,外显子和其他属性;每个转录本的基因、蛋白质、 HGNC 和其他标识符(在适用情况下,仅适用于 RefSeq cache);现有变异的位置、等位基因和频率;调控区域;SIFT,PolyPhen 的预测和分数; cache下载方式:进入vep官网:https://ftp.ensembl.org/pub/release-108/variation/ 点击vep 进入新的界面,选择需要的的物种和相应版本的cache文件包,点击右键,复制链接,在服务器中wget下载; (3) 下载vep相关的插件 vep插件放在github公共仓库中(链接:https://github.com/Ensembl/VEP_plugins),根据需要下载相应的插件; vep官网有vep插件的详细介绍(链接:https://asia.ensembl.org/info/docs/tools/vep/script/vep_plugins.html);

4.vep的使用4.1.使用基本的参数进行测试/opt/vep/src/ensembl-vep/vep -i ${sample}.vcf \ # 输入的vcf文件--dir_cache ${dir_cache} \ # cache文件包路径--fork 4 \ # 线程数--refseq /--merged \ # cache文件包的种类,vep 的cache文件包只有两种,分别是refseq,merged--format vcf \ --output_file test.vcf \ --custom ${dir_plugins}/MT.mitomap.vcf.gz,Mitomap,vcf,exact,0,AC,AF,Disease,DiseaseStatus \--offline \--assembly GRCh37 \ --use_given_ref \ --vcf \--force_overwrite \ 4.2 在docker下使用docker run -t -i -v $HOME/vep_data:/opt/vep/.vep ensemblorg/ensembl-vep

vep官网中有关运行vep有很多参数,比如 Basic options 、Cache options 、Other annotation sources、 Output format options 等,其包含很多参数选项和简介,根据不同的需求选择合适的参数;vep参数的网站 https://asia.ensembl.org/info/docs/tools/vep/script/vep_options.html

本文链接地址:https://www.jiuchutong.com/zhishi/287138.html 转载请保留说明!

上一篇:最大的苹果商店(最大的apple商店)

下一篇:RuoYi-Vue——图标使用(vue调整图标大小)

  • windows7系统如何重置电脑

    windows7系统如何重置电脑

  • 快手如何添加一键拨号(快手上面怎么添加好友)

    快手如何添加一键拨号(快手上面怎么添加好友)

  • 小米滑板车充电多久(小米滑板车充电口)

    小米滑板车充电多久(小米滑板车充电口)

  • 华为freebuds3怎么调音量(华为freebuds3怎么重新配对)

    华为freebuds3怎么调音量(华为freebuds3怎么重新配对)

  • 红米k30充电慢怎么回事(红米k30充电太慢)

    红米k30充电慢怎么回事(红米k30充电太慢)

  • 苹果手机安全码是什么东西(使用苹果手机安全吗)

    苹果手机安全码是什么东西(使用苹果手机安全吗)

  • 淘宝现金红包提现是什么意思(淘宝现金红包提现怎么做账?)

    淘宝现金红包提现是什么意思(淘宝现金红包提现怎么做账?)

  • 路由器灯变黄怎么修复(路由器灯变黄了)

    路由器灯变黄怎么修复(路由器灯变黄了)

  • 荣耀10怎么设置24小时(荣耀10怎么设置灭屏显示时间)

    荣耀10怎么设置24小时(荣耀10怎么设置灭屏显示时间)

  • 哔哩哔哩33卡可以拒收吗(哔哩哔哩33卡可以退吗)

    哔哩哔哩33卡可以拒收吗(哔哩哔哩33卡可以退吗)

  • oppoa57支持5gwifi吗(oppoa57支持高清通话吗)

    oppoa57支持5gwifi吗(oppoa57支持高清通话吗)

  • 爱奇艺会员别人怎么登录(爱奇艺会员别人登录怎么取消)

    爱奇艺会员别人怎么登录(爱奇艺会员别人登录怎么取消)

  • 腾讯视频老是卡顿怎么回事(腾讯视频老是卡住)

    腾讯视频老是卡顿怎么回事(腾讯视频老是卡住)

  • cpu后面带x是什么意思(cpu后面的x)

    cpu后面带x是什么意思(cpu后面的x)

  • 抖音下载的视频在哪个文件(抖音下载的视频如何去掉抖音水印)

    抖音下载的视频在哪个文件(抖音下载的视频如何去掉抖音水印)

  • 小米手机恢复出厂设置系统会还原吗(小米手机恢复出厂设置后小米账号还会有吗)

    小米手机恢复出厂设置系统会还原吗(小米手机恢复出厂设置后小米账号还会有吗)

  • 喵喵机扫码连接失败(喵喵机扫码连接没反应)

    喵喵机扫码连接失败(喵喵机扫码连接没反应)

  • 荣耀v20支持防水吗(荣耀v20防水等级是多少)

    荣耀v20支持防水吗(荣耀v20防水等级是多少)

  • 淘宝怎么关闭农村淘宝(新版淘宝怎样关闭农村淘宝)

    淘宝怎么关闭农村淘宝(新版淘宝怎样关闭农村淘宝)

  • 如何修改唯品会登录名(如何修改唯品会绑定的手机号)

    如何修改唯品会登录名(如何修改唯品会绑定的手机号)

  • iqoo是什么屏幕(iqoo是什么屏幕LCD吗)

    iqoo是什么屏幕(iqoo是什么屏幕LCD吗)

  • 小米手机微信扫一扫扫不出来(小米手机微信扫码黑屏怎么回事)

    小米手机微信扫一扫扫不出来(小米手机微信扫码黑屏怎么回事)

  • 全球购骑士卡怎么退订单(全球购骑士卡怎么用)

    全球购骑士卡怎么退订单(全球购骑士卡怎么用)

  • 西数蓝盘和紫盘的区别(西数蓝盘和紫盘绿盘哪个好)

    西数蓝盘和紫盘的区别(西数蓝盘和紫盘绿盘哪个好)

  • iphone5s有没有录屏功能(苹果5s有没有屏幕录制功能)

    iphone5s有没有录屏功能(苹果5s有没有屏幕录制功能)

  • 首字下沉两行怎么设置(首字下沉两行怎么做word)

    首字下沉两行怎么设置(首字下沉两行怎么做word)

  • 无需命令行,这款小工具可以帮你在 Win11 上傻瓜式安装安卓 App(不用命令打开组策略)

    无需命令行,这款小工具可以帮你在 Win11 上傻瓜式安装安卓 App(不用命令打开组策略)

  • 企业所得税亏损弥补的规定
  • 报税扣款锁定如何解除
  • 哪些产品享受教育优惠
  • 土地增值税与增值税
  • 多余的进项税该怎么处理
  • 资产负债表里面的应收账款怎么算
  • 上月发票错误退回怎么做账
  • 固定资产融资租出计入什么科目
  • 旅游开发行业的会计怎么做
  • 控股合并的会计准则
  • 增值税10%怎么算
  • 政府扶持资金科目怎么做?
  • 专票已认证但又没有发票
  • 报关单金额大于合同金额
  • 一般纳税人购进货物用于免征增值税项目
  • 普通发票负数发票可以作废吗
  • 境内个人去境外开设公司
  • 合伙企业累进税率
  • 上年工资计提多了才发现
  • 发票失联企业不处理的后果
  • 培训机构的收入有哪几方面
  • 月末汇兑损益怎么算
  • 资产利润率和资本利润率监管标准
  • 票据背书怎么做账
  • 品牌使用费分录
  • thinkphp5.0自定义路由
  • 写一篇推荐书
  • ahc软件进不去
  • 蓝牙有哪些功能
  • linux运维有前景吗
  • PHP:pg_close()的用法_PostgreSQL函数
  • cesium加载s3m
  • php文件怎么打开运行
  • 财务管理中成本控制的原则
  • reactjs路由跳转
  • 激光器原理及应用
  • blkid命令详解
  • php对接第三方支付
  • 普通发票收款人填管理员可以吗
  • 洗浴报销票子怎么写
  • 只有发票没有合同可以做账吗
  • 出口退税转内销还是转免税好
  • 个体工商户经营所得税税率表2023年
  • 无形资产一般包括哪些权利
  • 小企业如何核算成本费用
  • mysql两个数据库连接查询
  • 公司控股另一家公司好与坏
  • 公司减少注册资本的条件和程序
  • 租房交付确认书
  • 我国会计准则规定企业的会计核算基础是
  • 合并报表两种方法
  • 材料盘盈如何入账
  • 租赁设备如何做会计分录
  • 成本结转怎么做分录
  • 公户转私户的钱会退回吗
  • 废料销售的会计处理
  • 红字发票是干什么的
  • 应收账款和坏账准备的解题思路
  • 净资产增长率增长说明什么
  • 新成立公司怎么办理社保开户
  • sql server的 update from 语句的深究
  • YunDetectService.exe是什么进程?YunDetectService.exe禁止使用的方法
  • win7系统对拷的方法
  • 微软修复系统
  • centos安装软件教程
  • WinCtlAdAlt.exe - WinCtlAdAlt是什么进程
  • ubuntu系统电脑开机键盘亮 屏幕不亮
  • windows缓存写入失败,数据怎么找回
  • linux图形界面与命令行
  • linux常用性能监控命令
  • linux就该这么些
  • html中渐变
  • unity 动画播放
  • 注释标签的内容会显示在浏览器窗口中
  • jquery.form.js实现文件上传
  • 重庆国税电子税务局手机版
  • 如何查询企业税收情况
  • 广东国税局发票查询平台
  • 选矿比怎么算
  • 融资租赁 合同
  • 免责声明:网站部分图片文字素材来源于网络,如有侵权,请及时告知,我们会第一时间删除,谢谢! 邮箱:opceo@qq.com

    鄂ICP备2023003026号

    网站地图: 企业信息 工商信息 财税知识 网络常识 编程技术

    友情链接: 武汉网站建设