title: Imputing amino acid polymorphisms in human leukocyte antigens.
abstract: DNA sequence variation within human leukocyte antigen (HLA) genes mediate susceptibility to a wide range of human diseases. The complex genetic structure of the major histocompatibility complex (MHC) makes it difficult, however, to collect genotyping data in large cohorts. Long-range linkage disequilibrium between HLA loci and SNP markers across the major histocompatibility complex (MHC) region offers an alternative approach through imputation to interrogate HLA variation in existing GWAS data sets. Here we describe a computational strategy, SNP2HLA, to impute classical alleles and amino acid polymorphisms at class I (HLA-A, -B, -C) and class II (-DPA1, -DPB1, -DQA1, -DQB1, and -DRB1) loci. To characterize performance of SNP2HLA, we constructed two European ancestry reference panels, one based on data collected in HapMap-CEPH pedigrees (90 individuals) and another based on data collected by the Type 1 Diabetes Genetics Consortium (T1DGC, 5,225 individuals). We imputed HLA alleles in an independent data set from the British 1958 Birth Cohort (N = 918) with gold standard four-digit HLA types and SNPs genotyped using the Affymetrix GeneChip 500 K and Illumina Immunochip microarrays. We demonstrate that the sample size of the reference panel, rather than SNP density of the genotyping platform, is critical to achieve high imputation accuracy. Using the larger T1DGC reference panel, the average accuracy at four-digit resolution is 94.7% using the low-density Affymetrix GeneChip 500 K, and 96.7% using the high-density Illumina Immunochip. For amino acid polymorphisms within HLA genes, we achieve 98.6% and 99.3% accuracy using the Affymetrix GeneChip 500 K and Illumina Immunochip, respectively. Finally, we demonstrate how imputation and association testing at amino acid resolution can facilitate fine-mapping of primary MHC association signals, giving a specific example from type 1 diabetes.
authors.affiliation: Harvard-MIT (Massachusetts Institute of Technology) Division of Health Sciences and Technology, Boston, Massachusetts, United States of America.
authors.forename: Xiaoming
authors.initials: X
authors.lastname: Jia
authors.forename: Buhm
authors.initials: B
authors.lastname: Han
authors.forename: Suna
authors.initials: S
authors.lastname: Onengut-Gumuscu
authors.forename: Wei-Min
authors.initials: WM
authors.lastname: Chen
authors.forename: Patrick J
authors.initials: PJ
authors.lastname: Concannon
authors.forename: Stephen S
authors.initials: SS
authors.lastname: Rich
authors.forename: Soumya
authors.initials: S
authors.lastname: Raychaudhuri
authors.forename: Paul I W
authors.initials: PI
authors.lastname: de Bakker
mesh_terms: D000483:Alleles; D000596:Amino Acids / Q000235:genetics* / Q000276:immunology; D002874:Chromosome Mapping; D003922:Diabetes Mellitus, Type 1 / Q000208:ethnology / Q000235:genetics* / Q000276:immunology; D020022:Genetic Predisposition to Disease; D055106:Genome-Wide Association Study; D006680:HLA Antigens / Q000235:genetics* / Q000276:immunology; D060148:HapMap Project; D006239:Haplotypes; D015395:Histocompatibility Antigens Class I / Q000235:genetics* / Q000276:immunology; D000949:Histocompatibility Antigens Class II / Q000235:genetics* / Q000276:immunology; D006801:Humans; D015810:Linkage Disequilibrium; D020641:Polymorphism, Single Nucleotide; D044465:White People
chemical_list: D000596:Amino Acids; D006680:HLA Antigens; D015395:Histocompatibility Antigens Class I; D000949:Histocompatibility Antigens Class II