HLA-VBSeq:对全基因组数据进行HLA分型

2023-10-12 02:50

本文主要是介绍HLA-VBSeq:对全基因组数据进行HLA分型,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

欢迎关注"生信修炼手册"

HLA-VBseq 利用全基因组测序的数据,可以提供8位的HLA分型结果,其文献链接如下

https://bmcgenomics.biomedcentral.com/articles/10.1186/1471-2164-16-S2-S7

在该文献中,利用30X的全基因组数据,对HLA-VBSeq, PHLAT, HLAminer这3款软件的分型结果进行了评估,准确率汇总如下

可以看到,只有HLA-VBSeq提供了8位的分型结果,准确率高达99.94%;对于2位到4位的分型结果,其准确率也高于另外两款软件。

同时还评估了不同测序量时,各种软件提供的4位分型结果的准确率,结果如下

在不同条件下,HLA-VBseq的准确率都是最高的。由此可见,该软件的分型效果还是相当不错的,官网如下

http://nagasakilab.csml.org/hla/

该软件采用java语言开发,直接下载HLAVBseq.jar就可以了,除了该文件之外,还需要下载以下几个文件

  1. bamNameIndex.jar

  2. SamToFastq.jar

  3. parse_result.pl

  4. hla_all.fasta

  5. Allelelist.txt

前三个程序在处理fastq文件时会用到;后两个文件是从IMGA/HLA数据库下载的,如果觉得官网提供的版本较老,可以从IMGA/HLA数据库下载最新版。

软件的步骤较多,首先将fastq序列与参考基因组进行比对,得到bam文件,然后对该bam文件进行操作。步骤如下:

1. 挑选位于HLA 基因区域的reads

利用samtools view 命令挑选出比对到HLA区域的reads , 命令如下

samtools view -hb align.bam  chr6:29907037-29915661 chr6:31319649-31326989 chr6:31234526-31241863 chr6:32914391-32922899 chr6:32900406-32910847 chr6:32969960-32979389 chr6:32778540-32786825 chr6:33030346-33050555 chr6:33041703-33059473 chr6:32603183-32613429 chr6:32707163-32716664 chr6:32625241-32636466 chr6:32721875-32733330 chr6:32405619-32414826 chr6:32544547-32559613 chr6:32518778-32554154 chr6:32483154-32559613 chr6:30455183-30463982 chr6:29689117-29699106 chr6:29792756-29800899 chr6:29793613-29978954 chr6:29855105-29979733 chr6:29892236-29899009 chr6:30225339-30236728 chr6:31369356-31385092 chr6:31460658-31480901 chr6:29766192-29772202 chr6:32810986-32823755 chr6:32779544-32808599 chr6:29756731-29767588 | samtools fastq - -1 R1.fq -2 R2.fq

需要注意的是,在使用view命令时,虽然也可以直接提供一个bed格式的文件来挑选特定区域的reads,但是这种用法不会利用到bam文件的索引,所以速度很慢。对于全基因组数据,bam文件很大,上述写法虽然冗长,但是执行效率高。

2. 挑选没比对上的reads

利用samtools view 命令挑选出没有比对上参考基因组的reads, 命令如下:

samtools view -hb  -f 12 /home/pub/output/WGS/18B0315D/6343/6343_final.bam | samtools fastq - -1 unmapped_R1.fq -2 unmapped_R2.fq
3. 合并reads

将比对到HLA区域的reads和没比对上参考基因组的reads合并,命令如下

cat R1.fq unmapped_R1.fq > R1.fastq
cat R2.fq unmapped_R2.fq > R2.fastq
4. 与HLA参考reads比对

利用bwa软件,将上一步得到的reads与HLA参考序列比对,命令如下

bwa index hla_all.fasta
bwa mem -t 8 -P -L 10000 -a hla_all.fasta R1.fastq R2.fastq > out.sam
5. 运行HLA-VBSeq

HLA-VBSeq支持双端或者单端测序的数据,这里以双端数据为例,用法如下

java -jar HLAVBSeq.jar hla_all.fasta out.sam result.txt --alpha_zero 0.01 --is_paired
6. 格式化结果

上一步就已经生成结果了,这一步只是格式化,下面的代码会筛选出HLA-A基因的分型结果

perl parse_result.pl Allelelist.txt result.txt | grep "^A\*" | sort -k2 -n -r > HLA.txt

格式化之后的结果,内容如下

A*01:01:01:01 17.4022266628604
A*11:01:01 12.0376819868684

共两列,第一列为Allel, 第二列为该Allel区域的平均测序深度。

扫描关注微信号,更多精彩内容等着你!

这篇关于HLA-VBSeq:对全基因组数据进行HLA分型的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

Java使用Curator进行ZooKeeper操作的详细教程

《Java使用Curator进行ZooKeeper操作的详细教程》ApacheCurator是一个基于ZooKeeper的Java客户端库,它极大地简化了使用ZooKeeper的开发工作,在分布式系统... 目录1、简述2、核心功能2.1 CuratorFramework2.2 Recipes3、示例实践3

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

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

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

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

Mysql删除几亿条数据表中的部分数据的方法实现

《Mysql删除几亿条数据表中的部分数据的方法实现》在MySQL中删除一个大表中的数据时,需要特别注意操作的性能和对系统的影响,本文主要介绍了Mysql删除几亿条数据表中的部分数据的方法实现,具有一定... 目录1、需求2、方案1. 使用 DELETE 语句分批删除2. 使用 INPLACE ALTER T

Python Dash框架在数据可视化仪表板中的应用与实践记录

《PythonDash框架在数据可视化仪表板中的应用与实践记录》Python的PlotlyDash库提供了一种简便且强大的方式来构建和展示互动式数据仪表板,本篇文章将深入探讨如何使用Dash设计一... 目录python Dash框架在数据可视化仪表板中的应用与实践1. 什么是Plotly Dash?1.1

基于Flask框架添加多个AI模型的API并进行交互

《基于Flask框架添加多个AI模型的API并进行交互》:本文主要介绍如何基于Flask框架开发AI模型API管理系统,允许用户添加、删除不同AI模型的API密钥,感兴趣的可以了解下... 目录1. 概述2. 后端代码说明2.1 依赖库导入2.2 应用初始化2.3 API 存储字典2.4 路由函数2.5 应

Redis 中的热点键和数据倾斜示例详解

《Redis中的热点键和数据倾斜示例详解》热点键是指在Redis中被频繁访问的特定键,这些键由于其高访问频率,可能导致Redis服务器的性能问题,尤其是在高并发场景下,本文给大家介绍Redis中的热... 目录Redis 中的热点键和数据倾斜热点键(Hot Key)定义特点应对策略示例数据倾斜(Data S

Python使用date模块进行日期处理的终极指南

《Python使用date模块进行日期处理的终极指南》在处理与时间相关的数据时,Python的date模块是开发者最趁手的工具之一,本文将用通俗的语言,结合真实案例,带您掌握date模块的六大核心功能... 目录引言一、date模块的核心功能1.1 日期表示1.2 日期计算1.3 日期比较二、六大常用方法详

Python使用DrissionPage中ChromiumPage进行自动化网页操作

《Python使用DrissionPage中ChromiumPage进行自动化网页操作》DrissionPage作为一款轻量级且功能强大的浏览器自动化库,为开发者提供了丰富的功能支持,本文将使用Dri... 目录前言一、ChromiumPage基础操作1.初始化Drission 和 ChromiumPage

Python实现将MySQL中所有表的数据都导出为CSV文件并压缩

《Python实现将MySQL中所有表的数据都导出为CSV文件并压缩》这篇文章主要为大家详细介绍了如何使用Python将MySQL数据库中所有表的数据都导出为CSV文件到一个目录,并压缩为zip文件到... python将mysql数据库中所有表的数据都导出为CSV文件到一个目录,并压缩为zip文件到另一个