這篇文章主要為大家展示了“HLAreporter軟件有什么用”，內(nèi)容簡而易懂，條理清晰，希望能夠幫助大家解決疑惑，下面讓小編帶領(lǐng)大家一起研究并學(xué)習(xí)一下“HLAreporter軟件有什么用”這篇文章吧。

企業(yè)建站必須是能夠以充分展現(xiàn)企業(yè)形象為主要目的，是企業(yè)文化與產(chǎn)品對外擴展宣傳的重要窗口，一個合格的網(wǎng)站不僅僅能為公司帶來巨大的互聯(lián)網(wǎng)上的收集和信息發(fā)布平臺，創(chuàng)新互聯(lián)公司面向各種領(lǐng)域：成都自上料攪拌車等網(wǎng)站設(shè)計、成都全網(wǎng)營銷解決方案、網(wǎng)站設(shè)計等建站排名服務(wù)。

在前面的文章中，我們詳細(xì)介紹了HLA Allel的命名格式，示意圖如下

從示意圖可以看出，一個HLA Allel 可以分成四個字段，在加上最后的修飾后綴，共5個字段；在定義HLA 分型結(jié)果的分辨率時，會根據(jù)分型結(jié)果的最大位數(shù)來判斷，如果只給出了字段一，即血清學(xué)分類的信息，代表是2位的分型結(jié)果；如果最多給出了字段二，即對應(yīng)的蛋白信息，代表是4位的分型結(jié)果；如果最多給出了字段三，即CDS區(qū)信息，代表是8位的分型結(jié)果；如果分型結(jié)果給出了最后的后綴，代表是9位的分型結(jié)果。

HLA分型的技術(shù)手段很多，比如芯片，高通量測序等；不同手段識別到的HLA Allel 分辨率不同，如果只能給出2位的分型結(jié)果，則屬于低分辨率；如果給出4位分型結(jié)果，屬于中分辨率；能夠給出8位或以上分型結(jié)果，屬于高分辨率。

本篇文章主要介紹HLA  reporter 軟件，該軟件可以利用高通量測序的結(jié)果進(jìn)行HLA 分型。其分型結(jié)果分辨率高，最多可到9位。下圖是該軟件與其他同類軟件的比較結(jié)果

可以看到，在測試樣本中，HLAreporter 軟件分型結(jié)果由于HLAminer 的結(jié)果。該截圖來自于官方文獻(xiàn)，鏈接如下

https://genomemedicine.biomedcentral.com/articles/10.1186/s13073-015-0145-3

軟件的官網(wǎng)如下

http://paed.hku.hk/genome/software.html

安裝過程如下

wget http://paed.hku.hk/genome/software/HLAreporter.zip
unzip HLAreporter.zip
cd HLAreporter.v103

HLAreporter.sh就是該軟件的運行腳本。其原理如下

由于HLA的多態(tài)性，直接比對reads到參考基因組是不行的 ，HLAreporter 設(shè)計了一個CRP(comprehensive reference panel)參考基因組 , CRP 構(gòu)建時參考了IMGT/HLA 數(shù)據(jù)庫中已知的HLA  Allel信息，通過bwa將reads與CRP比對，提取比對到某個HLA基因的reads,然后進(jìn)行組裝，將組裝的contig與數(shù)據(jù)庫比較，確定最終的Allel。

該軟件的用法如下

bash HLAreporter.sh test HLA_B test.R1.fq test.R2.fq

第一個參數(shù)test代表樣本名稱；第二個參數(shù)代表檢測的HLA基因，第三個和第四個參數(shù)代表雙端測序的fastq序列。

整個pipeline 分成了3個部分

1. bwa 比對 CRP 數(shù)據(jù)庫

當(dāng)你提供了原始的fastq 數(shù)據(jù)時，會自動調(diào)用4digit_map_HLA.sh腳本進(jìn)行比對
該腳本核心內(nèi)容如下

bwa aln exon23_high_resolution_multi_ref.fa $1 > $3_1_exon23_high_resolution_multi_ref.sai
bwa aln exon23_high_resolution_multi_ref.fa $2 > $3_2_exon23_high_resolution_multi_ref.sai
bwa sampe exon23_high_resolution_multi_ref.fa $3_1_exon23_high_resolution_multi_ref.sai $3_2_exon23_high_resolution_multi_ref.sai $1 $2 > $3_exon23_high_resolution_multi_ref.sam
samtools view -bS $3_exon23_high_resolution_multi_ref.sam > $3_exon23_high_resolution_multi_ref.bam
samtools view -b -F 4 $3_exon23_high_resolution_multi_ref.bam > $3_exon23_high_resolution_multi_ref_mappedreads.bam
samtools sort $3_exon23_high_resolution_multi_ref_mappedreads.bam $3_exon23_high_resolution_multi_ref_mappedreads_sorted
samtools index $3_exon23_high_resolution_multi_ref_mappedreads_sorted.bam

$1和$2分別對應(yīng)雙端測序的R1和R2端reads, $3表示樣本名稱，通過調(diào)用bwa sampe，將原始的雙端reads與exon23_high_resolution_multi_ref.fa比對，生成exon23_high_resolution_multi_ref_mappedreads_sorted.bam 文件。

2. 識別HLA Allel

對于HLA-A , HLA-B, HLA-C 這三個I型基因而言，調(diào)用4digit_predict_classI_ssake.sh腳本；對于HLA II 型基因，調(diào)用4digit_predict_classII_main.sh腳本。

3. 給出Allel 在不同人群中的頻率

利用ANFD數(shù)據(jù)庫中的信息，給出每個Allel在不同人群中的頻率，通過腳本HLAfreq.sh實現(xiàn)。

最終輸出結(jié)果保存在results目錄下，示意如下

------------------------------------------------------------------
HLAreporter version 1.03
Report created
Fri Apr 17 11:24:24 HKT 2015
------------------------------------------------------------------
Allele
    DRB1*01:01:01G
Allele
    DRB1*04:07:01G
Phase
    DRB1*04:08:01
Alleles detected
DRB1*01:01:01G
DRB1*01:17
DRB1*04:08:01
DRB1*04:07:01G
Gap location:
97    208
Non-polymorphic gap:
111 bp
HLA data quality profile:
10xcov%    100    20xcov%    100    30xcov%    97    50xcov%    89
------------------------------------------------------------------
HLA allele frequency
Four populations in Europe China Japan Africa are shown
By allele frequency net database (www.allelefrequencies.net)
------------------------------------------------------------------
[Allele]    [EUR]    [CHN]    [JPN]    [AFR]
DRB1*01:01    0.0860    0.0230    0.0582    0.0130
DRB1*04:07    0.0112    0.0030    0.0057    0.0030
DRB1*04:08    0.0039    0.0020    0.0000    0.0000
DRB1*01:17    0.0000    0.0000    0.0000    0.0000

HLAreporter目前只支持對以下11個基因的分型

1. HLA_A

2. HLA_B

3. HLA_C

4. HLA_DRB1

5. HLA_DRB2

6. HLA_DRB3

7. HLA_DRB4

8. HLA_DRB5

9. HLA_DQB1

10. HLA_DPB1

11. HLA_DQA1

以上是“HLAreporter軟件有什么用”這篇文章的所有內(nèi)容，感謝各位的閱讀！相信大家都有了一定的了解，希望分享的內(nèi)容對大家有所幫助，如果還想學(xué)習(xí)更多知識，歡迎關(guān)注創(chuàng)新互聯(lián)行業(yè)資訊頻道！

當(dāng)前標(biāo)題：HLAreporter軟件有什么用
文章URL：http://chinadenli.net/article14/gjecge.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián)，為您提供品牌網(wǎng)站制作、定制開發(fā)、網(wǎng)站導(dǎo)航、品牌網(wǎng)站設(shè)計、關(guān)鍵詞優(yōu)化、品牌網(wǎng)站建設(shè)

廣告

聲明：本網(wǎng)站發(fā)布的內(nèi)容（圖片、視頻和文字）以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主，如果涉及侵權(quán)請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場，如需處理請聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載，或轉(zhuǎn)載時需注明來源： 創(chuàng)新互聯(lián)