Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
HLA and KIR typing Lj Lee Bioinformatics 1 06-03-2014 11:01 AM
4 Approaches for Preparing DNA for HLA Typing Genohub General 3 06-03-2014 05:43 AM
HLA typing vasvale Literature Watch 0 11-07-2008 12:53 PM
PubMed: Multiplex HLA-Typing by Pyrosequencing. Newsbot! Literature Watch 0 10-08-2008 05:00 AM

Thread Tools
Old 05-08-2017, 03:16 AM   #1
Location: Italia

Join Date: Jul 2012
Posts: 27
Default HLA typing with bwakit

Dear all,

has somebody attempted HLA typing with bwakit and can share some experience?
I have whole exome sequencing data of human individuals for who I would like to do HLA typing. I ran the bwakit as instructed in the Readme, but have problems interpreting the output. Here are the first few lines from the out.hla.all for HLA-B for one individual.
HLA-B*50:01:01  HLA-B*55:01:01  0       0       2
HLA-B*50:01:02  HLA-B*55:01:01  0       0       2
HLA-B*50:01:04  HLA-B*55:01:01  0       0       2
HLA-B*50:01:01  HLA-B*55:01:03  0       0       2
HLA-B*50:01:02  HLA-B*55:01:03  0       0       2
HLA-B*50:01:04  HLA-B*55:01:03  0       0       2
HLA-B*50:01:01  HLA-B*55:01:07  0       0       2
HLA-B*50:01:02  HLA-B*55:01:07  0       0       2
HLA-B*50:01:04  HLA-B*55:01:07  0       0       2
Here Heng Li explains the meaning of columns 3 to 5:
column 3: #mismatches on the primary exon(s)
column 4: #mismatches on all considered exons
column 5: #exons used in typing

My questions are:
1. What are columns 1 and 2? I would have expected results for one allele, not two.
2. What is the "best" or the "true" allele? The file simply gives the first line of the out.hla.all for each gene, but from the evidence in the out.hla.all file, the first line does not seem to be better than the second or third.
3. Would you say that given high coverage exome sequencing data, one should be able to clearly identify (with bwakit) the true allele an individual has for each of the HLA genes?

Thank you

Last edited by evakoe; 05-08-2017 at 06:08 AM.
evakoe is offline   Reply With Quote
Old 05-17-2017, 02:45 AM   #2
Location: Italia

Join Date: Jul 2012
Posts: 27

For the people interested in this post:
1. I assume the first two columns are the HLA alleles for the two human alleles, so the HLA-B allele is heterozygous in this case.
2. The coverage was probably not high enough to determine a six digit resolution. After all, the first four digits are allthe same, and the last two are variable. I ran OptiType on the same sample which gave the output
B1	B2
B*50:01	B*55:01
3. The coverage on this sample is about 25X, which apparently is enough for four digit resolution with bwakit, but not for six digit resolution.

Last edited by evakoe; 05-17-2017 at 02:51 AM.
evakoe is offline   Reply With Quote
Old 06-21-2017, 01:16 AM   #3
Location: Heidelberg

Join Date: Jun 2010
Posts: 14

hi evakoe,

which command did you use? cuz i'm trying to do the same but it's not working.
i'm using this:
run-bwamem -o NA12878-hs38 -t 20 -H -s /gpfs/data_jrnas1/ref_data/Hsapiens/hs38DH/hs38DH.fa NA12878_2.fastq.gz NA12878_1.fastq.gz |sh
and it doesn't give me sorted bam or HLA typing.
afadda is offline   Reply With Quote
Old 07-19-2017, 07:21 AM   #4
Location: Italia

Join Date: Jul 2012
Posts: 27

The command I used was
run-bwamem -t 4 -R "@RG\tID:myID\tSM:mySM\tPL:myPL\tLB:myLB\tPU:myPU" -H -d -s -o sample BWAindex sample_1.fastq sample_2.fastq | sh
and it ran as expected
evakoe is offline   Reply With Quote

bwakit, hla, hla typing, typing

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 05:56 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO