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

相关文章

Spring 请求之传递 JSON 数据的操作方法

《Spring请求之传递JSON数据的操作方法》JSON就是一种数据格式,有自己的格式和语法,使用文本表示一个对象或数组的信息,因此JSON本质是字符串,主要负责在不同的语言中数据传递和交换,这... 目录jsON 概念JSON 语法JSON 的语法JSON 的两种结构JSON 字符串和 Java 对象互转

C++如何通过Qt反射机制实现数据类序列化

《C++如何通过Qt反射机制实现数据类序列化》在C++工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作,所以本文就来聊聊C++如何通过Qt反射机制实现数据类序列化吧... 目录设计预期设计思路代码实现使用方法在 C++ 工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作。由于数据类

SpringBoot使用GZIP压缩反回数据问题

《SpringBoot使用GZIP压缩反回数据问题》:本文主要介绍SpringBoot使用GZIP压缩反回数据问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录SpringBoot使用GZIP压缩反回数据1、初识gzip2、gzip是什么,可以干什么?3、Spr

QT进行CSV文件初始化与读写操作

《QT进行CSV文件初始化与读写操作》这篇文章主要为大家详细介绍了在QT环境中如何进行CSV文件的初始化、写入和读取操作,本文为大家整理了相关的操作的多种方法,希望对大家有所帮助... 目录前言一、CSV文件初始化二、CSV写入三、CSV读取四、QT 逐行读取csv文件五、Qt如何将数据保存成CSV文件前言

SpringBoot集成Milvus实现数据增删改查功能

《SpringBoot集成Milvus实现数据增删改查功能》milvus支持的语言比较多,支持python,Java,Go,node等开发语言,本文主要介绍如何使用Java语言,采用springboo... 目录1、Milvus基本概念2、添加maven依赖3、配置yml文件4、创建MilvusClient

SpringValidation数据校验之约束注解与分组校验方式

《SpringValidation数据校验之约束注解与分组校验方式》本文将深入探讨SpringValidation的核心功能,帮助开发者掌握约束注解的使用技巧和分组校验的高级应用,从而构建更加健壮和可... 目录引言一、Spring Validation基础架构1.1 jsR-380标准与Spring整合1

通过Spring层面进行事务回滚的实现

《通过Spring层面进行事务回滚的实现》本文主要介绍了通过Spring层面进行事务回滚的实现,包括声明式事务和编程式事务,具有一定的参考价值,感兴趣的可以了解一下... 目录声明式事务回滚:1. 基础注解配置2. 指定回滚异常类型3. ​不回滚特殊场景编程式事务回滚:1. ​使用 TransactionT

MySQL 中查询 VARCHAR 类型 JSON 数据的问题记录

《MySQL中查询VARCHAR类型JSON数据的问题记录》在数据库设计中,有时我们会将JSON数据存储在VARCHAR或TEXT类型字段中,本文将详细介绍如何在MySQL中有效查询存储为V... 目录一、问题背景二、mysql jsON 函数2.1 常用 JSON 函数三、查询示例3.1 基本查询3.2

SpringBatch数据写入实现

《SpringBatch数据写入实现》SpringBatch通过ItemWriter接口及其丰富的实现,提供了强大的数据写入能力,本文主要介绍了SpringBatch数据写入实现,具有一定的参考价值,... 目录python引言一、ItemWriter核心概念二、数据库写入实现三、文件写入实现四、多目标写入

Java中使用Hutool进行AES加密解密的方法举例

《Java中使用Hutool进行AES加密解密的方法举例》AES是一种对称加密,所谓对称加密就是加密与解密使用的秘钥是一个,下面:本文主要介绍Java中使用Hutool进行AES加密解密的相关资料... 目录前言一、Hutool简介与引入1.1 Hutool简介1.2 引入Hutool二、AES加密解密基础