High-Accuracy HLA Type Inference from Whole-Genome Sequencing Data Using Population Reference Graphs - Inserm - Institut national de la santé et de la recherche médicale Accéder directement au contenu
Article Dans Une Revue PLoS Computational Biology Année : 2016

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

Résumé

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.
Fichier principal
Vignette du fichier
pcbi.1005151.pdf (1.47 Mo) Télécharger le fichier
Origine : Fichiers éditeurs autorisés sur une archive ouverte
Loading...

Dates et versions

inserm-02149588 , version 1 (06-06-2019)

Identifiants

Citer

Alexander T 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, 2016, 12 (10), pp.e1005151. ⟨10.1371/journal.pcbi.1005151⟩. ⟨inserm-02149588⟩
28 Consultations
100 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More