Y Zhang1,2,3, Y Song2,3, H Cao2,3, X Mo2,3,4, H Yang2,5, J Wang2,5, Z Lu1, T Zhang2,3. 1. State Key Lab of Bioelectronics, School of Biological Science and Medical Engineering, Southeast University, Nanjing, China. 2. Shenzen Key Laboratory of Neurogenomics, BGI-Shenzhen, Shenzhen, China. 3. China National GeneBank-Shenzhen, BGI-Shenzhen, Shenzhen, China. 4. BGI Education Center, University of Chinese Academy of Sciences, Shenzhen, China. 5. James D. Watson Institute of Genome Sciences, Hangzhou, China.
Abstract
BACKGROUND: HLA-DRB3, DRB4 and DRB5 (DRB3/4/5) are paralogues of HLA-DRB1. They have important roles in transplantation and have been reported to be related to many diseases. HLA typing methods for DRB3/4/5 based on NGS data have many limitations now, such as need of polymerase chain reaction (PCR) or low accuracy. MATERIALS AND METHODS: We present a HLA typing method for DRB3/4/5 based on read mapping and haplotype assembly from NGS data. Also, copy number of DRB3/4/5 is determined by a k-means clustering method according to ratio of sequencing depth between DRB3/4/5 and DRB1. RESULTS: We achieved 100%, 100%, 100% accuracy on simulated data and 95.88%, 98.89%, 99.34% accuracy on MHC capture Illumina sequencing data at 4-digit resolution with 30-fold coverage for DRB3/4/5 separately. We also explored the DRB3/4/5 profiles in five continental populations through low coverage WGS data generated by the 1000 Genome Project. We found that frequency of DRB4 in African were significantly lower than that in all other populations. CONCLUSION: Our method for DRB3/4/5 typing has high accuracy. It is a good supplement to regular HLA typing and could help in disease studies, medical applications and human population diversity studies.
BACKGROUND:HLA-DRB3, DRB4 and DRB5 (DRB3/4/5) are paralogues of HLA-DRB1. They have important roles in transplantation and have been reported to be related to many diseases. HLA typing methods for DRB3/4/5 based on NGS data have many limitations now, such as need of polymerase chain reaction (PCR) or low accuracy. MATERIALS AND METHODS: We present a HLA typing method for DRB3/4/5 based on read mapping and haplotype assembly from NGS data. Also, copy number of DRB3/4/5 is determined by a k-means clustering method according to ratio of sequencing depth between DRB3/4/5 and DRB1. RESULTS: We achieved 100%, 100%, 100% accuracy on simulated data and 95.88%, 98.89%, 99.34% accuracy on MHC capture Illumina sequencing data at 4-digit resolution with 30-fold coverage for DRB3/4/5 separately. We also explored the DRB3/4/5 profiles in five continental populations through low coverage WGS data generated by the 1000 Genome Project. We found that frequency of DRB4 in African were significantly lower than that in all other populations. CONCLUSION: Our method for DRB3/4/5 typing has high accuracy. It is a good supplement to regular HLA typing and could help in disease studies, medical applications and human population diversity studies.