Sentieon | 应用教程: 使用DNAscope对HiFi长读长数据进行胚系变异检测分析

本文主要是介绍Sentieon | 应用教程: 使用DNAscope对HiFi长读长数据进行胚系变异检测分析,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

介绍

本文描述了使用Sentieon® DNAscope进行PacBio® HiFi数据胚系突变检测。PacBio® HiFi技术产⽣质量值超过Q20的高质量长读段,平均长度在10-25kb之间。准确的长读段可以对短读段和高噪音长读段方法无法检测的基因组重复区域进行精准的变异检测。
Sentieon® DNAscope能利用PacBio® HiFi数据高质量、长读长的优势,使用经过校准的机器学习模型进行快速、精准的变异检测。针对HiFi数据的DNAscope流程输入比对后的HiFi数据,并输出VCF格式变异检测结果。

本流程的使用需要使用202010.03或更新版本的Sentieon软件和可从Sentieon®获取的相关脚本。本流程需要Python版本>2.7或>3.3,bcftools版本为1.10及以上。Pythonbcftoolsbedtools需要在用户环境变量中。

输入数据要求

比对后数据

本流程的输入使用pbmm2minimap2比对后的PacBio® HiFi数据。使⽤pbmm2比对推荐使用-c 0 -y 70 --preset HIFI参数。这些设置关闭了pbmm2传统的比对⼀致性过滤,转而使用间隔压缩序列过滤并使用PacBio®推荐的HiFi数据比对设置。使用minimap2比对推荐使用-x map-hifi参数。该设置为minimap2针对HiFi数据的推荐设置。

参考基因组

DNAscope将根据高质量参考基因组对样本进行变异检测。除了参考基因组外,由samtools生成的基因组索引文件(.fai)也是必须的。我们建议使用没有补丁序列的参考基因组。

PacBio® HiFi数据Sentieon® DNAscope流程

流程概览

本流程会进行两轮变异调用,然后将两轮结果合并到⼀起生成最终的结果。具体步骤如下:

  • 本流程在第⼀轮调用中会检测样本中的变异位点;
  • 利⽤第⼀轮检出的SNV和长读长信息进行定相;
  • 第⼆轮调用:
    在定相区,从每个单倍型中分别进行变异调用;
    在非定相区,使用更准确的⼆倍体模型进行变异调用;
  • 将第⼀轮和第⼆轮的变异位点合并生成最终结果;
  • 通过提供的MHC bed⽂件,对MHC区域进行特殊处理,进⼀步提高变异检测准确性;本流程所需的DNAscope机器学习模型可从https://github.com/Sentieon/sentieon-models获取。

运行流程

通过⼀个包含多个单⼀Sentieon命令的脚本即可运行HiFi数据DNAscope流程。⼀行命令即可完成变异检测并运用机器学习模型。HiFi数据比对文件可以是用pbmm2比对和建立索引后的bam或cram⽂件。

dnascope_HiFi.sh [-h] -r REFERENCE -i HIFI_BAM -m MODEL [-d dbSNP] [-B MHC_INTERVAL] [-b
INTERVAL] [-t NUMBER_THREADS] [-h] [--] VARIANT_VCF

HiFi数据Sentieon® DNAscope流程必选参数如下:
-r REFERENCE:参考基因组fasta文件路径。请确保使用的参考基因组⽂件与比对阶段使用的⼀致。
-i HIFI_BAM:比对后的BAM文件路径。
-m MODEL:DNAscope HiFi模型⽂件。
HiFi数据Sentieon® DNAscope流程可选参数如下:
-d dbSNP:dbSNP数据库VCF文件路径。仅需⼀个文件。该文件将用于变异检测结果中位点refSNP ID的注释。
-B MHC_INTERVAL:MHC区间文件,BED格式。该文件将用于MHC区域变异检测的特殊处理。
-b INTERVAL:区间文件,BED格式。该文件将限制变异检测在该区间。
-t NUMBER_THREADS:并行线程数。该参数可选,默认使用计算机所有线程。
-h:打印帮助信息。
HiFi数据Sentieon® DNAscope流程位置参数如下:
VARIANT_VCF:变异检测输出文件名。本流程会输出bgzip压缩的VCF文件及其索引文件。

流程输出文件

本流程会输出⼀个以VCF4.2格式标准的bgzip压缩文件(.vcf.gz)及其索引文件(.vcf.gz.tbi)。

其他注意事项

目前,该流程仅推荐应用于⼆倍体样本。对于包含⼆倍体和单倍体的样本,应使⽤-b INTERVAL参数将变异检测限定到⼆倍体染⾊体。

Sentieon软件介绍

Sentieon为完整的纯软件基因变异检测二级分析方案,其分析流程完全忠于BWA、GATK、MuTect2、STAR、Minimap2、Fgbio、picard等金标准的数学模型。在匹配开源流程分析结果的前提下,大幅提升WGS、WES、Panel、UMI、ctDNA、RNA等测序数据的分析效率和检出精度,并匹配目前全部第二代、三代测序平台。

在这里插入图片描述
Sentieon软件团队拥有丰富的软件开发及算法优化工程经验,致力于解决生物数据分析中的速度与准确度瓶颈,为来自于分子诊断、药物研发、临床医疗、人群队列、动植物等多个领域的合作伙伴提供高效精准的软件解决方案,共同推动基因技术的发展。

截至2023年3月份,Sentieon已经在全球范围内为1300+用户提供服务,被世界一级影响因子刊物如NEJM、Cell、Nature等广泛引用,引用次数超过700篇。此外,Sentieon连续数年摘得了Precision FDA、Dream Challenges等多个权威评比的桂冠,在业内获得广泛认可。

软件试用:https://www.insvast.com/sentieon

这篇关于Sentieon | 应用教程: 使用DNAscope对HiFi长读长数据进行胚系变异检测分析的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/169059

相关文章

springboot security使用jwt认证方式

《springbootsecurity使用jwt认证方式》:本文主要介绍springbootsecurity使用jwt认证方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地... 目录前言代码示例依赖定义mapper定义用户信息的实体beansecurity相关的类提供登录接口测试提供一

go中空接口的具体使用

《go中空接口的具体使用》空接口是一种特殊的接口类型,它不包含任何方法,本文主要介绍了go中空接口的具体使用,具有一定的参考价值,感兴趣的可以了解一下... 目录接口-空接口1. 什么是空接口?2. 如何使用空接口?第一,第二,第三,3. 空接口几个要注意的坑坑1:坑2:坑3:接口-空接口1. 什么是空接

Java利用JSONPath操作JSON数据的技术指南

《Java利用JSONPath操作JSON数据的技术指南》JSONPath是一种强大的工具,用于查询和操作JSON数据,类似于SQL的语法,它为处理复杂的JSON数据结构提供了简单且高效... 目录1、简述2、什么是 jsONPath?3、Java 示例3.1 基本查询3.2 过滤查询3.3 递归搜索3.4

springboot security快速使用示例详解

《springbootsecurity快速使用示例详解》:本文主要介绍springbootsecurity快速使用示例,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝... 目录创www.chinasem.cn建spring boot项目生成脚手架配置依赖接口示例代码项目结构启用s

Python如何使用__slots__实现节省内存和性能优化

《Python如何使用__slots__实现节省内存和性能优化》你有想过,一个小小的__slots__能让你的Python类内存消耗直接减半吗,没错,今天咱们要聊的就是这个让人眼前一亮的技巧,感兴趣的... 目录背景:内存吃得满满的类__slots__:你的内存管理小助手举个大概的例子:看看效果如何?1.

Python中随机休眠技术原理与应用详解

《Python中随机休眠技术原理与应用详解》在编程中,让程序暂停执行特定时间是常见需求,当需要引入不确定性时,随机休眠就成为关键技巧,下面我们就来看看Python中随机休眠技术的具体实现与应用吧... 目录引言一、实现原理与基础方法1.1 核心函数解析1.2 基础实现模板1.3 整数版实现二、典型应用场景2

java中使用POI生成Excel并导出过程

《java中使用POI生成Excel并导出过程》:本文主要介绍java中使用POI生成Excel并导出过程,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录需求说明及实现方式需求完成通用代码版本1版本2结果展示type参数为atype参数为b总结注:本文章中代码均为

springboot简单集成Security配置的教程

《springboot简单集成Security配置的教程》:本文主要介绍springboot简单集成Security配置的教程,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,... 目录集成Security安全框架引入依赖编写配置类WebSecurityConfig(自定义资源权限规则

Spring事务中@Transactional注解不生效的原因分析与解决

《Spring事务中@Transactional注解不生效的原因分析与解决》在Spring框架中,@Transactional注解是管理数据库事务的核心方式,本文将深入分析事务自调用的底层原理,解释为... 目录1. 引言2. 事务自调用问题重现2.1 示例代码2.2 问题现象3. 为什么事务自调用会失效3

MySQL大表数据的分区与分库分表的实现

《MySQL大表数据的分区与分库分表的实现》数据库的分区和分库分表是两种常用的技术方案,本文主要介绍了MySQL大表数据的分区与分库分表的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有... 目录1. mysql大表数据的分区1.1 什么是分区?1.2 分区的类型1.3 分区的优点1.4 分