HLAforest中如何使用RNA-seq数据进行HLA分型

这篇文章将为大家详细讲解有关HLAforest中如何使用RNA-seq数据进行HLA分型，文章内容质量较高，因此小编分享给大家做个参考，希望大家阅读完这篇文章后对相关知识有一定的了解。

1. 修改config.sh

主要是修改HLAFOREST_HOME和NUM_THREADS这两个选项，示例如下

#!/bin/bash
# User modifiable variables
HLAFOREST_HOME=/soft/hlaforest/
NUM_THREADS=10 # the number of threads bowtie should use

HLAFOREST_HOME指定软件的安装目录，NUM_THREADS指定bowtie比对的线程数。

2. 修改 CallSimulation.sh

调整CONFIG_PATH的设置

CONFIG_PATH=/soft/hlaforest/scripts/config.sh

3. 修改 CallHaplotypesPE.sh

调整CONFIG_PATH的设置

CONFIG_PATH=/soft/hlaforest/scripts/config.sh

4. 添加PATH环境变量

在环境变量的配置文件中,将scripts目录添加到PATH变量中，方便程序调用，写法如下

export PATH=/soft/hlaforest/scripts:$PATH

上述参数都调整好之后，软件就可以运行了。需要注意的是该软件依赖bioperl和bowtie,由于我的系统是已经装过了的，所以这里没给出这两个软件的安装过程。

软件自带测试数据集，用法如下

CallHaplotypesPE.sh test2/ test2/gm12878_short_1.fastq test2/gm12878_short_2.fastq

只需要指定双端测序的原始文件就可以了，默认分型结果保存在haplotypes.txt中，该文件内容示意如下

ROOT:A ROOT:A:11 ROOT:A:11:01 ROOT:A:11:01:04
ROOT:A ROOT:A:01 ROOT:A:01:01 ROOT:A:01:01:01 ROOT:A:01:01:01:01
ROOT:B ROOT:B:55 ROOT:B:55:01 ROOT:B:55:01:01
ROOT:B ROOT:B:08 ROOT:B:08:01 ROOT:B:08:01:01

关于HLAforest中如何使用RNA-seq数据进行HLA分型就分享到这里了，希望以上内容可以对大家有一定的帮助，可以学到更多知识。如果觉得文章不错，可以把它分享出去让更多的人看到。