High-Accuracy HLA Type Inference from Whole-Genome Sequencing Data Using Population Reference Graphs

Abstract : Genetic variation at the Human Leucocyte Antigen (HLA) genes is associated with many autoimmune and infectious disease phenotypes, is an important element of the immunolog-ical distinction between self and non-self, and shapes immune epitope repertoires. Determining the allelic state of the HLA genes (HLA typing) as a by-product of standard whole-genome sequencing data would therefore be highly desirable and enable the immunoge-netic characterization of samples in currently ongoing population sequencing projects. Extensive hyperpolymorphism and sequence similarity between the HLA genes, however, pose problems for accurate read mapping and make HLA type inference from whole-genome sequencing data a challenging problem. We describe how to address these challenges in a Population Reference Graph (PRG) framework. First, we construct a PRG for 46 (mostly HLA) genes and pseudogenes, their genomic context and their characterized sequence variants, integrating a database of over 10,000 known allele sequences. Second, we present a sequence-to-PRG paired-end read mapping algorithm that enables accurate read mapping for the HLA genes. Third, we infer the most likely pair of underlying alleles at G group resolution from the IMGT/HLA database at each locus, employing a simple likelihood framework. We show that HLA*PRG, our algorithm, outperforms existing methods by a wide margin. We evaluate HLA*PRG on six classical class I and class II HLA genes (HLA-A,-B,-C,-DQA1,-DQB1,-DRB1) and on a set of 14 samples (3 samples with 2 x 100bp, 11 samples with 2 x 250bp Illumina HiSeq data). Of 158 alleles tested, we correctly infer 157 alleles (99.4%). We also identify and re-type two erroneous alleles in the original validation data. We conclude that HLA*PRG for the first time achieves accuracies comparable to gold-standard reference methods from standard whole-genome sequencing data, though high computational demands (currently~30-250 CPU hours per sample) remain a significant challenge to practical application.
Document type :
Journal articles
Complete list of metadatas

Cited literature [34 references]  Display  Hide  Download

https://www.hal.inserm.fr/inserm-02149588
Contributor : Ana Paula Dutra Azevedo <>
Submitted on : Thursday, June 6, 2019 - 2:54:25 PM
Last modification on : Thursday, August 8, 2019 - 2:10:37 PM

File

pcbi.1005151.pdf
Publisher files allowed on an open archive

Identifiers

Collections

Citation

Alexander Dilthey, Pierre-Antoine Gourraud, Alexander Mentzer, Nezih Cereb, Zamin Iqbal, et al.. High-Accuracy HLA Type Inference from Whole-Genome Sequencing Data Using Population Reference Graphs. PLoS Computational Biology, Public Library of Science, 2016, 12 (10), pp.e1005151. ⟨10.1371/journal.pcbi.1005151⟩. ⟨inserm-02149588⟩

Share

Metrics

Record views

30

Files downloads

72