A high-resolution HLA reference panel capturing global population diversity enables multi-ethnic fine-mapping in HIV host response
==================================================================================================================================

* Yang Luo
* Masahiro Kanai
* Wanson Choi
* Xinyi Li
* Kenichi Yamamoto
* Kotaro Ogawa
* Maria Gutierrez-Arcelus
* Peter K. Gregersen
* Philip E. Stuart
* James T. Elder
* Jacques Fellay
* Mary Carrington
* David W. Haas
* Xiuqing Guo
* Nicholette D. Palmer
* Yii-Der Ida Chen
* Jerome. I. Rotter
* Kent. D. Taylor
* Stephen. S. Rich
* Adolfo Correa
* James G. Wilson
* Sekar Kathiresan
* Michael H. Cho
* Andres Metspalu
* Tonu Esko
* Yukinori Okada
* Buhm Han
* NHLBI Trans-Omics for Precision Medicine (TOPMed) Consortium
* Paul J. McLaren
* Soumya Raychaudhuri

## Abstract

Defining causal variation by fine-mapping can be more effective in multi-ethnic genetic studies, particularly in regions such as the MHC with highly population-specific structure. To enable such studies, we constructed a large (N=21,546) high resolution HLA reference panel spanning five global populations based on whole-genome sequencing data. Expectedly, we observed unique long-range HLA haplotypes within each population group. Despite this, we demonstrated consistently accurate imputation at G-group resolution (94.2%, 93.7%, 97.8% and 93.7% in Admixed African (AA), East Asian (EAS), European (EUR) and Latino (LAT)). We jointly analyzed genome-wide association studies (GWAS) of HIV-1 viral load from EUR, AA and LAT populations. Our analysis pinpointed the MHC association to three amino acid positions (97, 67 and 156) marking three consecutive pockets (C, B and D) within the HLA-B peptide binding groove, explaining 12.9% of trait variance, and obviating effects of previously reported associations from population-specific HIV studies.

## Main

The HLA genes located within the MHC region encode proteins that play essential roles in immune responses including antigen presentation. They account for more heritability than all other variants together for many diseases1–4. It also has more reported GWAS trait associations than any other locus5. The extended MHC region spans 6Mb on chromosome 6p21.3 and contains more than 260 genes6. Due to population-specific positive selection it harbors unusually high sequence variation, longer haplotypes than most of the genome, and haplotypes that are specific to individual ancestral populations7,8. Consequently, the MHC is among the most challenging regions in the genome to analyze. Advances in HLA imputation have enabled population-specific association and fine-mapping studies of this locus2,9–12. But despite large effect sizes, fine-mapping in multiple populations simultaneously is challenging without a single large and high-resolution multi-ethnic reference panel. This has caused confusion in some instances. For example, defining the driving HLA alleles may inform the design of antigenic peptides for vaccines13,14 for HIV-1, which led to 770,000 deaths in 2018 alone15. However, multiple risk HLA risk alleles have been independently reported in different populations1,10,16, and it is not clear if they represent truly population-specific signals or are confounded by linkage.

## Results

### Performance evaluation of inferred classical HLA alleles

To build a large-scale multi-ethnic HLA imputation reference panel, we used high-coverage whole genome sequencing (WGS) datasets17–21 from the Japan Biological Informatics Consortium20, the BioBank Japan Project18, the Estonian Biobank22, the 1000 Genomes Project (1KG)21 and a subset of studies in the TOPMed program (**Supplementary Note, Supplementary Table 1-2**). To perform HLA typing using WGS data, we extracted reads mapped to the extended MHC region (chr6:25Mb-35Mb) and unmapped reads from 24,338 samples. We applied a population reference graph23–25, for the MHC region to infer classical alleles for three HLA class I genes (HLA*-A, -B* and -*C*) and five class II genes (HLA*-DQA1*, -*DQB1, -DRB1*, -*DPA1*, -*DPB1*) at G-group resolution, which determines the sequences of the exons encoding the peptide binding groove. We required samples to have >20x coverage across all HLA genes (**Supplementary Table 1, 3**). After quality control our panel included 21,546 individuals: 10,187 EUR, 7,849 AA, 2,069 EAS, 952 LAT and 489 SAS.

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/07/18/2020.07.16.20155606/T1)

Table 1. Effect estimates for the haplotypes defined by the three independent amino acids in HLA-B associated with HIV-1 viral load.
Only haplotypes with >1% frequency in the overall population are listed (**Supplementary Table 15**). Classical alleles of HLA-B are grouped based on the amino acid residues presented at position 97, 67 and 156 in HLA-B. For each haplotype, the multivariate effect is given as an effect size, taking the most frequent haplotype (97R-67S-156L) as the reference (effect size = 0). Heterogeneity p-value (P(het)) of each haplotype is calculated using a F-statistics with two degrees of freedom (**Methods**). Effect size and its standard error in each population are listed only for haplotypes that show evidence of heterogeneity (P-value < 0.05 /26, bolded). Unadjusted haplotype frequencies are given in each population.

To assess the accuracy of the WGS *HLA* allele calls, we compared the inferred *HLA* classical alleles to gold standard sequence-based typing (SBT) in 955 1KG subjects and 288 Japanese subjects and quantified concordance. In both cohorts we observed slightly higher average accuracy for class I genes, obtaining 99.0% (one-field, formally known as two-digit), 99.2% (amino acid) and 96.5% (G-group resolution), than class II genes, obtaining 98.7% (one-field), 99.7% (amino acid) and 96.7% (G-group resolution, **Methods, Supplementary Figure 1, Supplementary Tables 4-5, Extended Data 1**).

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/18/2020.07.16.20155606/F1/graphic-2.medium.gif)

[](http://medrxiv.org/content/early/2020/07/18/2020.07.16.20155606/F1/graphic-2)

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/18/2020.07.16.20155606/F1/graphic-3.medium.gif)

[](http://medrxiv.org/content/early/2020/07/18/2020.07.16.20155606/F1/graphic-3)

Figure 1. Global diversity of the MHC region.
(**a**) Principal component analysis of the pairwise IBD distance between 21,546 samples using MHC region markers. Allele diversity of (**b**) HLA*-B* and (**c**) HLA-*DQA1* among five continental populations (AA=Admixed African; EUR=European; LAT=Latino; EAS=East Asian; SAS=South Asian). The top two most common alleles within each population group are named, the remaining alleles are grouped as ‘others’.

### HLA diversity

To quantify MHC diversity, we calculated identity-by-descent (IBD) distances26 between all individuals using 38,398 MHC single nucleotide polymorphisms (SNPs) included in the multi-ethnic HLA reference panel (N=21,546) and applied principal component analysis (PCA, **Methods**). PCA distinguished EUR, EAS and AA as well as the admixed LAT and SAS samples (**Figure 1a, Supplementary Figure 2**). This reflected widespread *HLA* allele frequency differences between populations (**Figure 1b-c, Supplementary Figure 3**). Of 130 unique common (frequency > 1%) G-group alleles, 129 demonstrated significant differences of frequencies across populations (4 degree-of-freedom Chi-square test, p-value < 0.05/130, **Supplementary Figure 4**). The only exception was *DQA1*01:01:01G* which was nominally significant (unadjusted p-value = 0.047). These differences may be related to adaptive selection. For example, the *B*53:01:01G* allele is enriched in Admixed Africans (11.7% in AA versus 0.3% in others) and it has been previously associated with malaria protection27,28. Consistent with previous reports29,30, we observed that HLA*-B* had the highest allelic diversity (n=443) while HLA-*DQA1* had the least (n=17, **Supplementary Figure 5-6, Supplementary Table 6, Extended Data 1**).

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/18/2020.07.16.20155606/F2/graphic-4.medium.gif)

[](http://medrxiv.org/content/early/2020/07/18/2020.07.16.20155606/F2/graphic-4)

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/18/2020.07.16.20155606/F2/graphic-5.medium.gif)

[](http://medrxiv.org/content/early/2020/07/18/2020.07.16.20155606/F2/graphic-5)

Figure 2. Pairwise LD and haplotype structure for six classical HLA genes in five population groups.
(**a**) shows the pairwise normalized entropy (ε) measuring the difference of the haplotype frequency distribution for linkage disequilibrium and linkage equilibrium among five population groups. It takes values between 0 (no LD) to 1 (perfect LD). (**b**) shows the haplotype structures of the eight classical HLA genes in each population. The tile in a bar represents an *HLA* allele, and its height corresponds to the frequencies of the *HLA* allele. The gray lines connecting between two alleles represent *HLA* haplotypes. The width of these lines corresponds to the frequencies of the haplotypes. The most frequent long-range HLA haplotypes within each population is bolded and highlighted in a color described by the key at the bottom.

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/18/2020.07.16.20155606/F3/graphic-6.medium.gif)

[](http://medrxiv.org/content/early/2020/07/18/2020.07.16.20155606/F3/graphic-6)

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/18/2020.07.16.20155606/F3/graphic-7.medium.gif)

[](http://medrxiv.org/content/early/2020/07/18/2020.07.16.20155606/F3/graphic-7)

Figure 3. The multi-ethnic HLA reference panel shows improvement in allele diversity and imputation accuracy.
(**a**). The number of HLA alleles at the two-field resolution included in the multi-ethnic HLA reference panel (N = 21,546) compared to the European only Type 1 Diabetes Genetics Consortium48 (T1DGC) panel (N = 5,225) as well as a subset of the multi-ethnic HLA panel down-sampled to the same size as T1DGC. (**b)**. The correlation between imputed and typed dosages of classical *HLA* alleles using the multi-ethnic HLA reference panel at one-filed (red), two-field (blue) and G-group resolution (black) of the 955 1000 Genomes subjects. (**c**). The imputation accuracy for five classical HLA genes at one-field, two-field and G-group resolution. (**d**). The imputation accuracy at G-group resolution of the 1000 Genomes subjects stratified by four diverse ancestries when using three different imputation reference panels as described in (**a**).

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/18/2020.07.16.20155606/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2020/07/18/2020.07.16.20155606/F4)

Figure 4. Stepwise conditional analysis of the allele and amino acid positions of classical HLA genes to HIV-1 viral load.
Each circle point represents the linear regression -log10(*P**binary*) for all classical *HLA* alleles. Each diamond point represents -log10(*P**omnibus*) for the tested amino acid positions in HLA (blue=HLA-*A*; yellow=HLA-*C*; red=HLA-*B*; lightblue=HLA-*DRB1*; green=HLA*-DQA1*; purple=HLA-*DQB1*, darkgreen=HLA-*DPA1*; lightgreen=HLA-*DPB1*). Association at amino acid positions with more than two alleles was calculated using a multi-degree-of-freedom omnibus test. The dashed blacked line represents the significance threshold of *P* = 5 × 10−8. Each panel shows the association plot in the process of stepwise conditional omnibus test. (**a**) One-field classical allele *B*57* (*P* = 9.84 × 10−138) and (**b**) amino acid position 97 in HLA-B (*P* *omnibus*= 2.86 × 10−184) showed the strongest association signal. Results conditioned on position 97 in HLA-B showed a secondary signal at (**c**) classical allele *B*81:0101:G* (*P* = 4.53 × 10−23) and (**d**) position 67 in HLA-B (*P* *omnibus*= 1.08 × 10−40). Results conditioned on position 97 and 67 in HLA-B showed the same classical allele (**e**) *B*81:0101G* (*P* = 2.70 × 10−23) and (**f**) third signal at position 156 in HLA-B (*P* *omnibus*= 1.92 × 10−30). Results conditioned on position 97, 67 and 156 int HLA-B showed a fourth signal at (**g**) HLA*-A*31* (*P* = 2.45 × 10−8) and (**h**) position 77 in HLA-A (*P**omnibus* = 5.35 × 10−7) outside HLA-B.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/07/18/2020.07.16.20155606/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2020/07/18/2020.07.16.20155606/F5)

Figure 5. Location and effect of three independently associated amino acid positions in HLA-B.
(**a**) Allele frequency of six residues at position 97 in HLA-B among three populations. (**b**) Effect on spVL (i.e., change in log10 HIV-1 spVL per allele copy) of individual amino acid residues at position 97 in HLA-B. Results were calculated per allele using linear regression models, including gender and principal components within each ancestry as covariates. (**c**) HLA-B (PDB ID code 2bvp) proteins. Omnibus and stepwise conditional analysis identified three independent amino acid positions (positions 97 (red), 67 (orange), and 156 (green) in HLA-B. (**d**) Effect on spVL (i.e., change in log10 HIV-1 spVL per allele copy) of individual amino acid residues at each position reported in this and previous work10,16. Results were calculated per allele using linear regression models. The x-axis shows the effect size and its standard errors in the joint analysis, and the y-axis shows the effect size and its standard error in individual populations (purple = Admixed American; blue = European and orange = Latino). (**e**) Variance of spVL explained by the haplotypes formed by different amino acid positions.

To understand the haplotype structure of HLA between pairs of HLA genes we calculated a multiallelic linkage disequilibrium (LD) measurement index31–33, ε, which is 0 when there is no LD and 1 when there is perfect LD (**Figure 2a**). We observed higher ε between *DQA1, DQB1*, and *DRB1;* between *DPA1* and *DPB1*; and between *B* and *C* (**Supplementary Figure 7**). The heterogeneity between different populations was underscored by the presence of population-specific common (frequency >1%) high resolution long-range haplotypes (HLA-*A∼C∼B∼DRB1∼DQA1∼DQB1∼DPA1∼DPB1*, **Figure 2b, Supplementary Figure 8-12, Extended Data 2, Methods**). The most common within-population haplotype was A24::DP6 (HLA*-A\*24:02:01G∼C\*12:02:01G∼B\*52:01:01G∼DRB1\*15:02:01G∼DQA1\*01:03:01G∼DQB1\*06 :01:01G∼DPA1\*02:01:01G∼DPB1\*09:01:01G*) found at a frequency of 3.61% in EAS (**Supplementary Figure 8**). This haplotype is strongly associated with immune-mediated traits such as HIV34 and ulcerative colitis35 in Japanese individuals. The next most common haplotype was the well-described European-specific ancestral haplotype A1::DP1 or 8.136,37 (frequency=2.76%,

HLA-*A\*01:01:01G∼C\*07:01:01G∼B\*08:01:01G∼DRB1\*03:01:01G∼DQA1\*05:01:01G∼DQB1\*02: 01:01G∼DPA1\*02:01:02G∼DPB1\*01:01:01G*, **Supplementary Figure 9**). This haplotype is associated with diverse immunopathological phenotypes in the European population, including systemic lupus erythematosus38, myositis39 and several other conditions36. We observed long-range haplotypes in admixed populations including A1::DP4 in SAS (frequency=1.86%,

**Supplementary Figure 10**), A30::DP1 in AA (frequency=1.18%,HLA*-A\*30:01:01G∼C\*17:01:01G∼B\*42:01:01:G∼DRB1\*03:02:01G∼DQA1\*04:01:01G∼DQB1\*04 :02:01G∼DPA1\*02:02:02G∼DPB1\*01:01:01G*, **Supplementary Figure 11**), and A29::DP11 in LAT (frequency=0.74%,

HLA*-A\*29:02:01G∼C\*16:01:01G∼B\*44\*03:01:G∼DRB1\*07:01:01G∼DQA1\*02:01:01G∼DQB1*0 2:01:01G∼DPA1\*02:01:01G∼DPB1\*11:01:01G*, **Supplementary Figure 12**).

These haplotypes also have associations with multiple diseases: for example *C\*06:02∼B\*57:01* is associated with psoriasis40 and *A\*30:01∼C\*17:01∼B*42:01* is associated with HIV41.

### HLA selection signature

Previous studies have suggested that recent natural selection favors African ancestry in the HLA region in admixed populations42–45. To test this hypothesis in our data, we obtained WGS data from a subset of individuals within two admixed populations (1,832 AA and 594 LAT, determined by the first three global principal components, **Supplementary Figure 13, Supplementary Note**). Admixed individuals have genomes that are a mosaic of different ancestries. If genetic variations or haplotypes from an ancestral population are advantageous, then they are under selection and are expected to have higher frequency than by chance. Using ELAI46, we quantified how much the ancestry proportions differed within the MHC from the genome-wide average. In AA, we observed that the average genome-wide proportion of African ancestry was 74.5%, compared to 78.0% in the extended MHC region, corresponding to a 3.42 (95% CI: 3.35-3.49) standard deviation increase. In LAT, we observed 5.76% African ancestry genome-wide versus 16.0% in the extended MHC region, representing an increase of 4.23 (95% CI: 4.14-4.31) standard deviations (**Methods, Supplementary Figure 14**). To ensure our results are robust to different local ancestry inference methods, we applied an alternative method called RFMix47 and observed a similarly consistent MHC-specific excess of African ancestry in LAT, and also an excess in AA that was more modest (**Supplementary Figure 14)**.

### Construction of a multi-ethnic HLA reference panel and its performance evaluation

Next, we constructed a multi-ethnic HLA imputation reference panel based on classical HLA alleles and 38,398 genomic markers in the extended MHC region using a novel HLA-focused pipeline HLA-TAPAS (HLA-Typing At Protein for Association Studies). Briefly, HLA-TAPAS can handle HLA reference panel construction (*MakeReference*); HLA imputation (*SNP2HLA*) and HLA association (*HLAassoc*) (**Methods, URLs**). Compared to a widely used HLA reference panel with European-only individuals (The Type 1 Diabetes Genetics Consortium48, T1DGC), this new reference panel has a six-fold increase in the number of observed *HLA* alleles and non-HLA genomic markers (**Supplementary Table 7**). We noted the difference in observed classical *HLA* alleles is mainly due to the inclusion of diverse populations rather than its size; after downsampling the reference panel to be the same size as T1DGC (N=5,225), there was still a three-fold increase in observed alleles (**Figure 3a**).

To empirically assess imputation accuracy of our reference panel, we first used the publicly available gold-standard *HLA* types (HLA-*A*, -*B*,-*C*, -*DRB1* and -*DQB1*) of 1,267 diverse samples from AA, EAS, EUR and LAT included in 1KG. We removed 955 overlapping samples within the reference panel, and to ensure a representative analysis we kept 6,007 markers overlapping with the *Global Genotyping Array* SNPs. Across the five genes, the average G-group resolution accuracies were 94.2%, 93.7%, 97.8% and 93.7% in AA, EAS, EUR and LAT (**Figure 3b-c, Supplementary Table 8, Methods, Extended Data 3**). Compared to the T1DGC panel, our multi-ethnic reference panel showed the most improvement for individuals of non-European descent; we obtained 4.27%, 2.96%, 2.90% and 1.05% improvement at G-group resolution for AA, EAS, LAT, and EUR individuals, respectively (**Figure 3d**). Increased diversity was responsible for the improvement; downsampling the reference panel be the same size as the T1DGC panel still yielded superior performance (**Figure 3d**). To validate our panel further, we imputed *HLA* alleles into a multi-ethnic cohort of 2,291 individuals from the Genotype and Phenotype (GaP) registry genotyped on the ImmunoChip array. We obtained SBT *HLA* type information for six classical class I and class II loci (HLA-*A*, -*B*, -*C*, -*DQA1*, -*DQB1*, -*DRB1*) in 75 samples with diverse ancestral background (25 EUR, 25 EAS and 25 AA, **Supplementary Figure 15, Methods**). Average accuracies were 99.0%, 95.7% and 97.0% for EUR, EAS and AA respectively when comparing SBT *HLA* alleles at G-group resolution (**Methods, Extended Data 3**). Similar to the 1KG analysis, the multi-ethnic reference panel showed significant improvement for individuals with non-European descent (6.3% and 11.1% improvement for EAS and African individuals respectively at G-group resolution), and a more modest 2% improvement in EUR (**Supplementary Figure 16, Supplementary Table 9**).

### Fine-mapping causal variants of HIV jointly in three populations in the MHC region

Next we investigated MHC effects within human immunodeficiency virus type 1 (HIV-1) set point viral load. Upon primary infection with HIV-1, the set point viral load is reached after the immune system has developed specific cytotoxic T lymphocytes (CTL) that are able to partially control the virus. It has been well-established that the set point viral load (spVL) varies in the infected population and positively correlates with rate of disease progression49. Previous studies suggested that HIV-1 infection has a strong genetic component, and specific HLA class I alleles explain the majority of genetic risk10,50. The existence of multiple independent, ancestry-specific, risk-associated alleles has been reported in both European1,10 and African American16 populations. However, without a multi-ethnic reference panel it has not been possible to determine if these signals are consistent across different ancestral groups.

To define the MHC allelic effects shared across multiple populations, we applied our multi-ethnic MHC reference panel to 7,445 EUR, 3,901 AA and 677 LAT HIV-1 infected subjects (**Methods, Supplementary Table 10**). Imputation resulted in 640 classical HLA alleles, 4,513 amino acids in HLA proteins and 49,321 SNPs in the extended MHC region for association and fine-mapping analysis. We confirmed 96.6% imputation accuracy of two-field (or four-digit) resolution with a minor allele frequency > 0.5% in this cohort by comparing imputed classical alleles to the SBT alleles in a subset of 1,067 AA subjects16(**Supplementary Figure 17, Extended Data 3**).

We next tested SNPs, amino acid positions and classical *HLA* alleles across the MHC for association to spVL. We performed this jointly in EUR, AA and LAT population using a linear regression model with sex, principal components and ancestry as covariates (**Methods**). In agreement with previous studies, we found the strongest spVL-associated classical *HLA* allele is *B*57* (effect size = −0.84, *P**binary* = 8.68 × 10−144). This corresponded to a single residue Val97 in HLA-B that tracks almost perfectly with *B*57* (*r*2 = 0.995) and showed the strongest association of any single residue (effect size = −0.84, *P**binary* = 5.99 × 10−145, **Supplementary Figure 18**).

Then to determine which amino acid positions have independent association with spVL, we tested each of the amino acid positions by grouping haplotypes carrying a specific residue at each position in an additive model2,9 (**Methods**). We found the strongest spVL-associated amino acid variant in HLA-B is as previously reported1,10,16 at position 97 (**Figure 4a-b, Supplementary Table 11**) which strikingly explains 9.06% of the phenotypic variance. Position 97 in HLA-B was more significant (*P**omnibus* = 2.86 × 10−184) than any single SNP or classical *HLA* allele, including *B*57* (**Supplementary Figure 18, Extended Data 4**). Of the six allelic variants (Val/Asn/Trp/Thr/Arg/Ser) at this position, the Val residue conferred the strongest protective effect (effect size = −0.88, *P* = 9.32 × 10−152, **Supplementary Figure 19**) relative to the mos common residue Arg (frequency = 47.8%). All six amino acid alleles have consistent frequencies and effect sizes across the three population groups (**Figure 5a-b, Supplementary Figure 20**).

We next wanted to test whether there were other independent effects outside of position 97 in HLA-B. After accounting for the effects of amino acid 97 in HLA-B using a conditional haplotype analysis (**Methods**), we observed a significant independent association at position 67 in HLA-B (*P**omnibus* = 2.82 × 10−39, **Figure 4c-d, Supplementary Table 11**). Considering this might be an artifact of forward search, we exhaustively tested all possible pairs of polymorphic amino acid positions in HLA-B. Of 7,260 pairs of amino acid positions, none obtained a better goodness-of-fit than the pair of positions 97 and 67, which collectively explained 11.2% variance in spVL (**Figure 5e, Supplementary Table 12**). At position 67, Met67 residue shows the most protective effect (effect size = −0.44, *P* = 1.19 × 10−59) among the five possible amino acids (Cys/Phe/Met/Ser/Tyr) relative to the most common residue Ser (frequency =10.0%).

Conditioning on positions 97 and 67 revealed an additional association at position 156 in HLA-B (*P**omnibus* = 1.92 × 10−30, **Figure 4e-f, Supplementary Table 11**). In agreement with the stepwise conditional analysis, when we tested all 287,980 possible combinations of three amino acid positions in HLA-B, the most statistically significant combination of amino acids sites is 67, 97 and 156 (*P* = 5.68 × 10−244, **Supplementary Table 13**). These three positions explained 12.9% of the variance (**Figure 5e**). At position 156, residue Arg shows the largest risk effect (effect size = 0.180, *P* = 8.92 × 10−14) among the four possible allelic variants (Leu/Arg/Asp/Trp), relative to the most common residue Leu (frequency = 35.1%).

These amino acid positions mark three consecutive pockets within the HLA-B peptide-binding groove (**Figure 5c**). Position 97 is located in the C-pocket and has an important role in determining the specificity of the peptide-binding groove51,52. Position 67 is in the B-pocket, and Met67 side chains occupy the space where larger B-pocket anchors reside in other peptide-MHC structures; its presence limits the size of potential peptide position P2 side chains52. Amino acid position 156 is part of the D-pocket and influences the conformation of the peptide-binding region53. These results are consistent with the observation that in HLA-*B*57*, the single most protective spVL-associated one-field allele (a single change at position 156 from Leu → Arg or equivalently HLA-*B*57:03 →* HLA-*B*57:02*) leads to an increased repertoire of HIV-specific epitope41,54.

Despite differences in the power to detect associations due to differences in allele frequencies (**Supplementary Figure 21**), we observed generally consistent effects of individual residues across populations (**Figure 5d, Supplementary Figure 22-23, Supplementary Table 14**). There are 26 unique haplotypes defined by the amino acids at positions 67, 97 and 156 in HLA-B (**Table 1, Supplementary Table 15**). When we tested for effect size heterogeneity by ancestry for each of these haplotypes (**Methods**), we observed only 2 of 26 haplotypes showed heterogeneity (F-test P-value < 0.05/26), possibly due to different interplay between genetic and environmental variation at population-level. These results support the concept that these positions mediate HIV-1 viral load in diverse ancestries.

To assess whether there were other independent MHC associations outside HLA-B, we conditioned on all amino acid positions in HLA-B and observed associations at HLA-A, including at position 77 in HLA-A (*P* *omnibus*= 9.10 × 10−7, **Figure 4g-h, Supplementary Table 11**), the classical *HLA* allele *HLA-A*31* (*P* *binary*= 2.45 × 10−8) and the *rs2256919* promoter SNP (*P**binary* =3.10 × 10−16, **Supplementary Figure 18**). These associations argue for an effect at HLA-A, but larger studies and functional studies will be necessary to define the driving effects.

## Discussion

In our study we demonstrated accurate imputation with a single large reference panel for HLA imputation. We have shown how this reference panel can be used to impute genetic variation at eight *HLA* classical genes accurately across a wide range of populations. Accurate imputation in multi-ethnic studies is essential for fine-mapping.

We showed the utility of this approach by defining the alleles that best explain HIV-1 viral load in infected individuals. Our work implicates three amino acid positions (97, 67 and 156) in HLA-B in conferring the known protective effect of HLA class I variation on HIV-1 infection. Combining all alleles at these three positions explained 12.9% of the variance in spVL (**Figure 5e**). These positions all fall within the peptide-binding groove of the respective MHC protein (**Figure 5c**), indicating that variation in the amino acid content of the peptide-binding groove is the major genetic determinant of HIV control. Supported by experimental studies54–57, positions highlighted in our work indicated a structural basis for the HLA association with HIV disease progression that is mediated by the conformation of the peptide within the class I binding groove. This result highlights how a study with ancestrally diverse populations can potentially point to causal variation by leveraging linkage disequilibrium difference between ethnic groups.

We note that previous studies have shown position 97 in HLA-B has the strongest association with HIV-1 spVL or case-control in African American and European populations, but highlighted different additional signals via conditional analysis (position 45, 67 in HLA-B and position 77, 95 in HLA-A in Europeans1,10,16 and position 63, 116 and 245 in HLA-B in African Americans16). These signals do not explain the signals we report here; after conditioning on positions 45, 63, 116, 245 of HLA-B and 95 of HLA-A, the association of the four identified amino acids identified in this study remained significant (*P* < 5 × 10−8). In contrast, our binding groove alleles explain these other alleles; conditioning on the four amino acid positions identified in this study (positions 67, 97 and 156 in HLA-B), all previously reported positions did not pass the significance threshold (*P* > 5 × 10−8, **Supplementary Figure 24**).

Furthermore, defining the effect sizes for *HLA* alleles across different populations is essential for defining risk of a wide-range of diseases in the clinical setting. There is increasing application of genome-wide genotyping by patients both by healthcare providers and direct-to-consumer vendors. The large effects of the MHC region for a wide-range of immune and non-immune traits, makes it essential to define *HLA* allelic effect sizes essential in multi-ethnic studies in order to build generally applicable clinical polygenic risk scores for many diseases in diverse populations58–61. Resources like the one we present here will be an essential ingredient in such studies.

## Methods

### Individuals included in the reference panel

Study participants were from the Jackson Heart Study (JHS, N = 3,027), Multi-Ethnic Study of Atherosclerosis (MESA, N=4,620), Chronic Obstructive Pulmonary Disease Gene (COPDGene) study (N=10,623), Estonian Biobank (EST, N=2,244), Japan Biological Informatics Consortium (JPN, N=295), Biobank Japan (JPN, N=1,025) and 1000 Genomes Project (1KG, N=2,504). Each study was previously approved by respective institutional review boards (IRBs), including for the generation of WGS data and association with phenotypes. All participants provided written consent. Further details of cohort descriptions and phenotype definitions are described in the **Supplementary Note**.

### HLA-TAPAS

HLA-TAPAS (HLA-Typing At Protein for Association Studies) is an HLA-focused pipeline that can handle HLA reference panel construction (*MakeReference*), HLA imputation (*SNP2HLA*), and HLA association (*HLAassoc*). It is an updated version of the SNP2HLA48 to build an imputation reference panel, perform *HLA* classical allele, amino acid and SNP imputation within the extended MHC region. Briefly, major updates include (1) using PLINK1.9 (**URLs**) instead of v1.07; (2) using BEAGLE v4.1 (**URLs**) instead of v3 for phasing and imputation; and (3) including custom R scripts for performing association and fine-mapping analysis at amino acid level in multiple ancestries. The source code is available for download (**URLs**).

### Construction of a multi-ethnic HLA reference panel using whole-genome sequences

To construct a multi-ethnic HLA imputation reference panel, we used 24,338 whole-genome sequences at different depths (**Supplementary Table 1**). Details of the construction using deep-coverage whole-genome sequencing are described in the **Supplementary Note**. Briefly, alignment and variant-calling for genomes sequenced by each cohort were performed independently. We performed local realignment and quality recalibration with the Genome Analysis Toolkit62 (GATK; version 3.6) on Chromosome 6:25,000,000-35,000,000. We detected single nucleotide variants (SNV) and indels using GATK with HaplotypeCaller. To eliminate false-positive sites called in the MHC region, we restrict our panel to SNVs reported in 1000 Genomes Project21 only.

We next inferred classical HLA alleles at G-group resolution for eight classical HLA genes (*HLA-A*, -*B*, -*C*, -*DQA1*, -*DQB1*, -*DRB1, -DPA1* and -*DPB1*) using a population reference graph24,25. To extend the reference panel versatility, we inferred amino acid variation, one-field and two-field resolution alleles from the inferred G-group alleles. After removing samples with low-coverage and failed genome-wide quality control (**Supplementary Table 3**), we constructed a multi-ethnic HLA imputation reference panel (N=21,546) using the HLA-TAPAS *MakeReference* module (**URLs, Method**).

### Sequence-based typing of *HLA* alleles

Purified DNA from the 75 donors from the GaP registry (at the Feinstein Institute for Medical Research) was sent to NHS Blood and Transplant, UK, where *HLA* typing was performed. Next-generation sequencing was done for HLA*-A, -B, -C*, -*DQB1*, -*DPB1* and *-DRB1*.

PCR-sequence-specific oligonucleotide probe sequencing was performed for HLA*-DQA1* in all samples. These typing methods yielded classical allele calls for seven genes at three-field (HLA-*A*, -*B*, -*C* and -*DQB1*) or G-group resolution (HLA-*DQA1*, -*DPB1* and *-DRB1*).

Genomic DNA from the 288 unrelated samples of Japanese ancestry underwent high-resolution allele typing (three-field alleles) of six classical HLA genes (HLA-*A*, -*B* and -*C* for class I; and HLA-*DRB1*, -*DQA1* and -*DPB1* for class II)20.

The 1000 Genomes panel consists of 1,267 individuals with information on five HLA genes (HLA-*A, -B, -C, -DQB1*, and *-DRB1*) at G-group resolution among four major ancestral groups (AA, EAS, EUR and LAT)7.

We obtained HLA typing of the 1,067 African American subjects included in the HIV-1 viral load study as described previously16,63. Briefly, seven classical HLA genes (HLA-*A, -B, -C, -DQA1*, -*DQB1* -*DRB1* and *-DPB1*) were obtained by sequencing exons 2 and 3 and/or single-stranded conformation polymorphism PCR, and was provided at two-field resolution.

### Accuracy measure between inferred and sequence-based typing *HLA* genotypes

Allelic variants at HLA genes can be typed at different resolutions: one-field HLA types specify serological activity, two-field HLA types specify the amino acids encoded by the exons of the HLA gene, and three-field types determine the full exonic sequence including synonymous variants. G-group resolution determines the sequences of the exons encoding the peptide binding groove, that is, exons 2 and 3 for class I and exon 2 class II genes. Thus, any polymorphism occurring in exon 4 of class I gene or exon 3 of class II gene was not defined. This means many G-group alleles can map to multiple three-field and two-field *HLA* alleles.

We calculated the accuracy at each *HLA* gene by summing across the dosage of each correctly inferred *HLA* allele or amino acid across all individuals (N), and divided by the total number of observations (2*N). That is, ![Formula][1]</img>  where *Accuracy*(*g*) represents the accuracy at a classical HLA gene (e.g. HLA-*B*). *D**i*represents the inferred dosage of an allele in individual *i*, and alleles *A*1*i,g* and *A*2*i,g* represent the true (SBT) *HLA* types for an individual *i*.

To evaluate the accuracy between the inferred and validated *HLA* types obtained from SBT at G-group resolution, we translated the highest resolution specified by the validation data to its matching G-group resolution based IMGT/HLA database (e.g. HLA-*A*01:01* → HLA-*A*01:01:01G*), and compared it to the primary output from *HLA*LA* or *HLA-TAPAS*. We also translated all G-group alleles to their matching amino acid sequences, and compared them against the validation alleles, we referred to this as the amino acid level.

To evaluate imputation performance in individual classical *HLA* alleles and amino acids, we calculated the dosage *r*2 correlation between imputed and SBT dosage. ![Formula][2]</img>  where *x**i* and *y**i* represents the inferred and SBT dosage of an allele in individual *i*. *N* represents the number of individuals.

### Principal component analysis

We performed a principal component analysis of the MHC region based on the identity-by-descent (IBD) distances between all 21,809 individuals included in the multi-ethnic reference panel. We computed the IBD distance using Beagle (Version 4.1, **URLs**) and averaged over 100 runs with all variants (54,474) included in the HLA reference panel. Due to uneven representation of different ethnicity groups (**Supplementary Table 2**), we applied a weighted PCA approach, where mean and standard deviation of the IBD matrix within an ethnicity group are weighted inversely proportional to the sample size.

### HLA haplotype frequency estimation

We applied an expectation-maximization algorithm approach implemented in Hapl-o-Mat64 (**URLs**) to estimate HLA haplotype frequency based on eight classical HLA alleles inferred at G-group resolution. We estimated haplotype frequencies both overall and within five continental populations (**Extended Data 2**).

### Local ancestry inference

To detect local ancestry in admixed samples, we first applied ELAI46 to chromosome 6 with 1000 Genomes Project21 as the reference panel. We extracted 63,998 common HapMap3 SNPs between the WGS (MESA cohort) and the 1000 Genome reference panel. We used the same set of SNPs for ELAI and RFMix analysis. We applied ELAI46 to 1,832 African Americans and 594 Latinos. For 1,832 African American individuals included in the study, we used genotypes of 99 CEU and 108 YRI in the 1000 Genome Project as reference panel, assuming admixture generation to be seven generations ago. We used two upper-layer clusters and 10 lower-layer clusters in the model. For Latinos, we selected 65 Latinos with Native American (NAT) ancestry > 75% included in the 1000 Genomes Project identified using the ADMIXTURE analysis65 and used these individuals with high NAT, as well as CEU and YRI from 1000 Genomes as reference panels. We assumed that the admixture time was 20 generations ago. For ELAI, we used three upper-layer clusters and 15 lower-layer clusters in the model.

To address the technical concerns that local ancestry methods are biased by the high LD of MHC region66,67, we performed an alternative method, RFMix47, for local ancestry inference that accounts for high LD and lack of parental reference panels. Similar deviation from genome-wide ancestry was observed using RFMix (**Supplementary Figure 14)**, indicating that the selection signals we observed here are robust to different inference methods.

### HLA imputation in the HIV-1 viral load GWAS data in three populations

We used genome-wide genotyping data from 12,023 HIV-1 infected individuals aggregated across more than 10 different cohorts (**Supplementary Table 10**). The details of these samples and quality control procedures have been described previously10,68. Using the HIV-1 viral load GWAS data, we extracted the genotypes of SNPs located in the extended MHC region (chr6:28-34Mb, **Supplementary Table 10**). We conducted genotype imputation of one-field, two-field and G-group classical *HLA* alleles and amino acid polymorphisms of the eight class I and class II HLA genes using the constructed multi-ethnic HLA imputation reference panel and the HLA-TAPAS pipeline.

After imputation, we obtained the genotypes of 640 classical alleles, 4,513 amino acid positions of the eight classical HLA genes, and 49,321 SNPs located in the extended MHC region. We excluded variants with MAF < 0.5% and imputation *r*2 < 0.5 for all association studies. In total, we tested 51,358 variants in our association and fine-mapping study.

### HLA association analysis

For the HIV-1 viral loads of EUR, AA and LAT samples, we conducted a joint haplotype-based association analysis using a linear regression model under the assumption of additive effects of the number of HLA haplotypes for each individual. Phased haplotypes at a locus (i.e., HLA amino acid position) were constructed from the phased imputed genotypes of variants in the locus (i.e., amino acid change or SNP) and were converted to a haplotype matrix where each row is observed haplotypes (in the locus), not genotypes.

For each amino acid position, we applied a conditional haplotype analysis. We tested a multiallelic association between the HIV-1 viral load and a haplotype matrix (of the position) with covariates, including sex, study-specific PCs, and a categorical variable indicating a population. That is ![Formula][3]</img>  where *x**i* is the amino acid haplotype formed by each of the *m* amino acid residues that occur at that position, and *c**j* are the covariates included in the model.

To get an omnibus *P*-value for each position, we estimated the effect of each amino acid by assessing the significance of the improvement in fit by calculating the in-model fit, compared to a null model following an F-distribution with degrees of freedom. This is implemented using an ANOVA test in R as described previously32,69. The most frequent haplotype was excluded from a haplotype matrix as a reference haplotype for association.

**For the conditional analysis**, we assumed that the null model consisted of haplotypes as defined by residues at previously defined amino acid positions. The alternative model is in addition of another position with *m* residues. We tested whether the addition of those amino acid positions, and the creation of additional haplotypes groups, improved on the previousset. We then assessed the significance of the improvement in the delta deviance (sum of squares) over the previous model using an F-test. We performed stepwise conditional analysis to identify additional independent signals by adjusting for the most significant amino acid position in each step until none met the significance threshold (*P* = 5 × 10−8). We restricted analysis to haplotypes that have a minimum of 10 occurrences within HLA-B, and removed any individual with rare haplotypes for the conditional analysis.

**For the exhaustive search**, we tested all possible amino acid pairs and triplets for association. For each set of amino acid positions, we used the groups of residues occurring at these positions to estimate effect size and calculated for each of these models the delta deviance in risk prediction and its p-values compared to the null model.

### Heterogeneity testing of effect sizes

We used interaction analyses with models that included haplotype-by-ancestry (*Haplotype × Ancestry*) interaction terms. The fit of nested models was compared to a null model using the *F*-statistic with two degrees of freedom, for which the association interaction P-value indicated whether the inclusion of the *Haplotype × Ancestry* interaction terms improved the model fit compared to the null model that did not include the interaction terms. Interaction P-values for all haplotypes formed by positions 97, 67 and 156 in HLA-B are listed in **Supplementary Table 15**. Haplotypes that had a significant Bonferroni-corrected *Haplotype × Ancestry* interaction heterogeneity P-value (P < 0.05/26) were considered to show evidence of significant effect size heterogeneity between ancestries.

## Data Availability

The source code is available for download at https://github.com/immunogenomics

## URLs

HLA-TAPAS, [https://github.com/immunogenomics/HLA-TAPAS](https://github.com/immunogenomics/HLA-TAPAS)

IMGT/HLA, [https://www.ebi.ac.uk/ipd/imgt/hla/](https://www.ebi.ac.uk/ipd/imgt/hla/);

GATK version 3.6, [https://software.broadinstitute.org/gatk/download/archive](https://software.broadinstitute.org/gatk/download/archive);

HLA*LA, [https://github.com/DiltheyLab/HLA-PRG-LA](https://github.com/DiltheyLab/HLA-PRG-LA);

PLINK 1.90, [https://www.cog-genomics.org/plink2](https://www.cog-genomics.org/plink2);

Beagle 4.1, [https://faculty.washington.edu/browning/beagle/b4\_1.html](https://faculty.washington.edu/browning/beagle/b4_1.html);

Hapl-o-Mat, [https://github.com/DKMS/Hapl-o-Mat/](https://github.com/DKMS/Hapl-o-Mat/);

1000 Genomes gold-standard HLA types, [http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/data\_collections/HLA\_types/](http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/data_collections/HLA_types/)

## Author contributions

Y. L. and S.R. conceived, designed and performed analyses, wrote the manuscript and supervised the research. M.K. implemented the omnibus test for the HIV-1 fine-mapping study. Y.L., W.C., M.K., P.E.S., J.T.E., and B.H. contributed to the development of the HLA-TAPAS pipeline. X.L. performed the selection analysis. J.T.E, M.G.-A. and P.K.G helped with the GaP data acquisition. K.Y., K.O., D.W.H., X.G., N.D.P., Y.I.C., J.I.R., K.D.T., S.S.R., A.C., J.G.W., S.K., M.H.C., A.M., T.E., and Y.O. contributed to the WGS data acquisition. J.F., M.C. and P.J.M contributed to the HIV-1 data acquisition. All authors contributed to the writing of the manuscript.

## Competing interests

M.H.C. has received consulting or speaking fees from Illumina and AstraZeneca, and grant support from GSK and Bayer.

## Acknowledgements

The study was supported by the National Institutes of Health (NIH) TB Research Unit Network, Grant U19 AI111224-01.

The views expressed in this manuscript are those of the authors and do not necessarily represent the views of the National Heart, Lung, and Blood Institute; the National Institutes of Health; or the U.S. Department of Health and Human Services.

The Genotype and Phenotype (GaP) Registry at The Feinstein Institute for Medical Research provided fresh, de-identified human plasma; blood was collected from control subjects under an IRB-approved protocol (IRB# 09-081) and processed to isolate plasma. The GaP is a sub-protocol of the Tissue Donation Program (TDP) at Northwell Health and a national resource for genotype-phenotype studies.

[https://www.feinsteininstitute.org/robert-s-boas-center-for-genomics-and-human-genetics/gap-registry/](https://www.feinsteininstitute.org/robert-s-boas-center-for-genomics-and-human-genetics/gap-registry/) A.M. is supported by Gentransmed grant 2014-2020.4.01.15-0012.; D.W.H. is supported by NIH grants AI110527, AI077505, TR000445, AI069439, and AI110527. D.H.S. was supported by R01 HL92301, R01 HL67348, R01 NS058700, R01 AR48797, R01 DK071891, R01 AG058921, the General Clinical Research Center of the Wake Forest University School of Medicine (M01 RR07122, F32 HL085989), the American Diabetes Association, and a pilot grant from the Claude Pepper Older Americans Independence Center of Wake Forest University Health Sciences (P60 AG10484). J.T.E. and P.E.S. were supported by NIH/NIAMS R01 AR042742, R01 AR050511, and R01 AR063611.

For some HIV cohort participants, DNA and data collection was supported by NIH/NIAID AIDS Clinical Trial Group (ACTG) grants UM1 AI068634, UM1 AI068636 and UM1 AI106701, and ACTG clinical research site grants A1069412, A1069423, A1069424, A1069503, AI025859, AI025868, AI027658, AI027661, AI027666, AI027675, AI032782, AI034853, AI038858, AI045008, AI046370, AI046376, AI050409, AI050410, AI050410, AI058740, AI060354, AI068636, AI069412, AI069415, AI069418, AI069419, AI069423, AI069424, AI069428, AI069432, AI069432, AI069434, AI069439, AI069447, AI069450, AI069452, AI069465, AI069467, AI069470, AI069471, AI069472, AI069474, AI069477, AI069481, AI069484, AI069494, AI069495, AI069496, AI069501, AI069501, AI069502, AI069503, AI069511, AI069513, AI069532, AI069534, AI069556, AI072626, AI073961, RR000046, RR000425, RR023561, RR024156, RR024160, RR024996, RR025008, RR025747, RR025777, RR025780, TR000004, TR000058, TR000124, TR000170, TR000439, TR000445, TR000457, TR001079, TR001082, TR001111, and TR024160.

Molecular data for the Trans-Omics in Precision Medicine (TOPMed) program was supported by the National Heart, Lung and Blood Institute (NHLBI). See the TOPMed Omics Support Table (**Supplementary Table 16**) for study specific omics support information. Core support including centralized genomic read mapping and genotype calling, along with variant quality metrics and filtering were provided by the TOPMed Informatics Research Center (3R01HL-117626-02S1; contract HHSN268201800002I). Core support including phenotype harmonization, data management, sample-identity QC, and general program coordination were provided by the TOPMed Data Coordinating Center (R01HL-120393; U01HL-120393; contract HHSN268201800001I). We gratefully acknowledge the studies and participants who provided biological samples and data for TOPMed.

The COPDGene project was supported by Award Number U01 HL089897 and Award Number U01 HL089856 from the National Heart, Lung, and Blood Institute. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Heart, Lung, and Blood Institute or the National Institutes of Health. The COPDGene project is also supported by the COPD Foundation through contributions made to an Industry Advisory Board comprised of AstraZeneca, Boehringer Ingelheim, GlaxoSmithKline, Novartis, Pfizer, Siemens and Sunovion. A full listing of COPDGene investigators can be found at: [http://www.copdgene.org/directory](http://www.copdgene.org/directory)

The Jackson Heart Study (JHS) is supported and conducted in collaboration with Jackson State University (HHSN268201800013I), Tougaloo College (HHSN268201800014I), theMississippi State Department of Health (HHSN268201800015I) and the University of Mississippi Medical Center (HHSN268201800010I, HHSN268201800011I and HHSN268201800012I) contracts from the National Heart, Lung, and Blood Institute (NHLBI) and the National Institute on Minority Health and Health Disparities (NIMHD). The authors also wish to thank the staffs and participants of the JHS.

MESA and the MESA SHARe project are conducted and supported by the National Heart, Lung, and Blood Institute (NHLBI) in collaboration with MESA investigators. Support for MESA is provided by contracts 75N92020D00001, HHSN268201500003I, N01-HC-95159, 75N92020D00005, N01-HC-95160, 75N92020D00002, N01-HC-95161, 75N92020D00003, N01-HC-95162, 75N92020D00006, N01-HC-95163, 75N92020D00004, N01-HC-95164, 75N92020D00007, N01-HC-95165, N01-HC-95166, N01-HC-95167, N01-HC-95168, N01-HC-95169, UL1-TR-000040, UL1-TR-001079, UL1-TR-001420. MESA Family is conducted and supported by the National Heart, Lung, and Blood Institute (NHLBI) in collaboration with MESA investigators. Support is provided by grants and contracts R01HL071051, R01HL071205, R01HL071250, R01HL071251, R01HL071258, R01HL071259, by the National Center for Research Resources, Grant UL1RR033176. The provision of genotyping data was supported in part by the National Center for Advancing Translational Sciences, CTSI grant UL1TR001881, and the National Institute of Diabetes and Digestive and Kidney Disease Diabetes Research Center (DRC) grant DK063491 to the Southern California Diabetes Endocrinology Research Center. This project has been funded in whole or in part with federal funds from the Frederick National Laboratory for Cancer Research, under Contract No. HHSN261200800001E. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government. This Research was supported in part by the Intramural Research Program of the NIH, Frederick National Lab, Center for Cancer Research.

*   Received July 16, 2020.
*   Revision received July 16, 2020.
*   Accepted July 18, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/)

## References

1.  1.International HIV Controllers Study et al. The major genetic determinants of HIV-1 control affect HLA class I peptide presentation. Science 330, 1551–1557 (2010).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzMzAvNjAxMC8xNTUxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDcvMTgvMjAyMC4wNy4xNi4yMDE1NTYwNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

2.  2.Raychaudhuri, S. et al. Five amino acids in three HLA proteins explain most of the association between MHC and seropositive rheumatoid arthritis. Nat. Genet. 44, 291–296 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.1076&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22286218&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

3.  3.Evans, D. M. et al. Interaction between ERAP1 and HLA-B27 in ankylosing spondylitis implicates peptide handling in the mechanism for HLA-B27 in disease susceptibility. Nat. Genet. 43, 761–767 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.873&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21743469&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

4.  4.Snyder, A. et al. Genetic basis for clinical response to CTLA-4 blockade in melanoma. N. Engl. J. Med. 371, 2189–2199 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1406498&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25409260&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000345976700007&link_type=ISI) 

5.  5.Buniello, A. et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Research vol. 47 D1005–D1012 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gky1120&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30445434&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

6.  6.Horton, R. et al. Gene map of the extended human MHC. Nat. Rev. Genet. 5, 889–899 (2004).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrg1489&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15573121&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000225416800011&link_type=ISI) 

7.  7.Gourraud, P.-A. et al. HLA diversity in the 1000 genomes dataset. PLoS One 9, e97282 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0097282&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24988075&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

8.  8.Robinson, J. et al. IPD-IMGT/HLA Database. Nucleic Acids Res. 48, D948–D955 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkz950&link_type=DOI) 

9.  9.Hu, X. et al. Additive and interaction effects at three amino acid positions in HLA-DQ and HLA-DR molecules drive type 1 diabetes risk. Nat. Genet. 47, 898–905 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3353&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26168013&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

10. 10.McLaren, P. J. et al. Polymorphisms of large effect explain the majority of the host genetic contribution to variation of HIV-1 virus load. Proc. Natl. Acad. Sci. U. S. A. 112, 14658–14663 (2015).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTEyLzQ3LzE0NjU4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDcvMTgvMjAyMC4wNy4xNi4yMDE1NTYwNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

11. 11.Tian, C. et al. Genome-wide association and HLA region fine-mapping studies identify susceptibility loci for multiple common infections. Nat. Commun. 8, 599 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-017-00257-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28928442&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

12. 12.Onengut-Gumuscu, S. et al. Type 1 Diabetes Risk in African-Ancestry Participants and Utility of an Ancestry-Specific Genetic Risk Score. Diabetes Care 42, 406–415 (2019).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGlhY2FyZSI7czo1OiJyZXNpZCI7czo4OiI0Mi8zLzQwNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA3LzE4LzIwMjAuMDcuMTYuMjAxNTU2MDYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

13. 13.Matthews, P. C. et al. Central role of reverting mutations in HLA associations with human immunodeficiency virus set point. J. Virol. 82, 8548–8559 (2008).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoianZpIjtzOjU6InJlc2lkIjtzOjEwOiI4Mi8xNy84NTQ4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDcvMTgvMjAyMC4wNy4xNi4yMDE1NTYwNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

14. 14.Wang, Y. Development of a human leukocyte antigen-based HIV vaccine. F1000Res. 7, (2018).
    
    

15. 15.WHO | Progress reports on HIV. (2020).
    
    

16. 16.McLaren, P. J. et al. Fine-mapping classical HLA variation associated with durable host control of HIV-1 infection in African Americans. Hum. Mol. Genet. 21, 4334–4347 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/dds226&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22718199&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

17. 17.Taliun, D. et al. Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program. bioRxiv 563866 (2019) doi: 10.1101/563866.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czo4OiI1NjM4NjZ2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA3LzE4LzIwMjAuMDcuMTYuMjAxNTU2MDYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

18. 18.Okada, Y. et al. Deep whole-genome sequencing reveals recent selection signatures linked to evolution and disease risk of Japanese. Nat. Commun. 9, 1631 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-03274-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

19. 19.Mitt, M. et al. Improved imputation accuracy of rare and low-frequency variants using population-specific high-coverage WGS-based imputation reference panel. Eur. J. Hum. Genet. 25, 869–876 (2017).
    
    

20. 20.Hirata, J. et al. Genetic and phenotypic landscape of the major histocompatibilty complex region in the Japanese population. Nat. Genet. (2019) doi: 10.1038/s41588-018-0336-0.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0336-0&link_type=DOI) 

21. 21.1000 Genomes Project Consortium et al. A global reference for human genetic variation. Nature 526, 68–74 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature15393&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26432245&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

22. 22.Nelis, M. et al. Genetic structure of Europeans: a view from the north--east. PLoS One 4, (2009).
    
    

23. 23.Dilthey, A., Cox, C., Iqbal, Z., Nelson, M. R. & McVean, G. Improved genome inference in the MHC using a population reference graph. Nature Genetics vol. 47 682–688 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3257&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25915597&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

24. 24.Dilthey, A. T. et al. High-Accuracy HLA Type Inference from Whole-Genome Sequencing Data Using Population Reference Graphs. PLoS Comput. Biol. 12, e1005151 (2016).
    
    

25. 25.Dilthey, A. T. et al. HLA*LA-HLA typing from linearly projected graph alignments. Bioinformatics 35, 4394–4396 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btz235&link_type=DOI) 

26. 26.Browning, B. L. & Browning, S. R. A fast, powerful method for detecting identity by descent. Am. J. Hum. Genet. 88, 173–182 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2011.01.010&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21310274&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

27. 27.Hill, A. V. et al. Common west African HLA antigens are associated with protection from severe malaria. Nature 352, 595–600 (1991).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/352595a0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1865923&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1991GB21100046&link_type=ISI) 

28. 28.Sanchez-Mazas, A. et al. The HLA-B landscape of Africa: Signatures of pathogen-driven selection and molecular identification of candidate alleles to malaria protection. Mol. Ecol. 26, 6238–6252 (2017).
    
    

29. 29.Maiers, M., Gragert, L. & Klitz, W. High-resolution HLA alleles and haplotypes in the United States population. Hum. Immunol. 68, 779–788 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.humimm.2007.04.005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17869653&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249751300009&link_type=ISI) 

30. 30.Gonzalez-Galarza, F. F. et al. Allele frequency net database (AFND) 2020 update: gold-standard data classification, open access genotype data and new query tools. Nucleic Acids Res. 48, D783–D788 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkz1029&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

31. 31.Nothnagel, M., Fürst, R. & Rohde, K. Entropy as a measure for linkage disequilibrium over multilocus haplotype blocks. Hum. Hered. 54, 186–198 (2002).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1159/000070664&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12771551&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000183533100004&link_type=ISI) 

32. 32.Okada, Y. et al. Construction of a population-specific HLA imputation reference panel and its application to Graves’ disease risk in Japanese. Nat. Genet. 47, 798–802 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3310&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26029868&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

33. 33.Okada, Y. eLD: entropy-based linkage disequilibrium index between multiallelic sites. Hum Genome Var 5, 29 (2018).
    
    

34. 34.Chikata, T. et al. Host-specific adaptation of HIV-1 subtype B in the Japanese population. J. Virol. 88, 4764–4775 (2014).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoianZpIjtzOjU6InJlc2lkIjtzOjk6Ijg4LzkvNDc2NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA3LzE4LzIwMjAuMDcuMTYuMjAxNTU2MDYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

35. 35.Nomura, E. et al. Mapping of a disease susceptibility locus in chromosome 6p in Japanese patients with ulcerative colitis. Genes Immun. 5, 477–483 (2004).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15215890&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

36. 36.Price, P. et al. The genetic basis for the association of the 8.1 ancestral haplotype (A1, B8, DR3) with multiple immunopathological diseases. Immunol. Rev. 167, 257–274 (1999).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1600-065X.1999.tb01398.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10319267&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000079851100021&link_type=ISI) 

37. 37.Horton, R. et al. Variation analysis and gene annotation of eight MHC haplotypes: the MHC Haplotype Project. Immunogenetics 60, 1–18 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00251-007-0262-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18193213&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000252468000001&link_type=ISI) 

38. 38.Graham, R. R. et al. Visualizing human leukocyte antigen class II risk haplotypes in human systemic lupus erythematosus. Am. J. Hum. Genet. 71, 543–553 (2002).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/342290&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12145745&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000177489600008&link_type=ISI) 

39. 39.Miller, F. W. et al. Genome-wide association study identifies HLA 8.1 ancestral haplotype alleles as major genetic risk factors for myositis phenotypes. Genes Immun. 16, 470–480 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/gene.2015.28&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26291516&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

40. 40.Haapasalo, K. et al. The Psoriasis Risk Allele HLA-C*06:02 Shows Evidence of Association with Chronic or Recurrent Streptococcal Tonsillitis. Infect. Immun. 86, (2018).
    
    

41. 41.Kløverpris, H. N. et al. HIV control through a single nucleotide on the HLA-B locus. J. Virol. 86, 11493–11500 (2012).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoianZpIjtzOjU6InJlc2lkIjtzOjExOiI4Ni8yMS8xMTQ5MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA3LzE4LzIwMjAuMDcuMTYuMjAxNTU2MDYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

42. 42.Salter-Townshend, M. & Myers, S. Fine-Scale Inference of Ancestry Segments Without Prior Knowledge of Admixing Groups. Genetics 212, 869–889 (2019).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXRpY3MiO3M6NToicmVzaWQiO3M6OToiMjEyLzMvODY5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDcvMTgvMjAyMC4wNy4xNi4yMDE1NTYwNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

43. 43.Zhou, Q., Zhao, L. & Guan, Y. Strong Selection at MHC in Mexicans since Admixture. PLoS Genet. 12, e1005847 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1005847&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26863142&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

44. 44.Meyer, D.C, Aguiar, V. R., Bitarello, B. D.C, Brandt, D. Y. & Nunes, K. A genomic perspective on HLA evolution. Immunogenetics 70, 5–27 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00251-017-1017-3&link_type=DOI) 

45. 45.Norris, E. T. et al. Admixture-enabledselection for rapid adaptive evolution in the Americas. bioRxiv 783845 (2019) doi: 10.1101/783845.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czo4OiI3ODM4NDV2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA3LzE4LzIwMjAuMDcuMTYuMjAxNTU2MDYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

46. 46.Guan, Y. Detecting structure of haplotypes and local ancestry. Genetics 196, 625–642 (2014).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXRpY3MiO3M6NToicmVzaWQiO3M6OToiMTk2LzMvNjI1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDcvMTgvMjAyMC4wNy4xNi4yMDE1NTYwNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

47. 47.Maples, B. K., Gravel, S., Kenny, E. E. & Bustamante, C. D. RFMix: a discriminative modeling approach for rapid and robust local-ancestry inference. Am. J. Hum. Genet. 93, 278–288 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2013.06.020&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23910464&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

48. 48.Jia, X. et al. Imputing amino acid polymorphisms in human leukocyte antigens. PLoS One 8, e64683 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0064683&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23762245&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

49. 49.Mellors, J. W. et al. Quantitation of HIV-1 RNA in plasma predicts outcome after seroconversion. Ann. Intern. Med. 122, 573–579 (1995).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/0003-4819-122-8-199504150-00003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7887550&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995QR52500003&link_type=ISI) 

50. 50.Bartha, I. et al. Estimating the Respective Contributions of Human and Viral GeneticVariation to HIV Control. PLoS Comput. Biol. 13, e1005339 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1005339&link_type=DOI) 

51. 51.Blanco-Gelaz, M. A. et al. The amino acid at position 97 is involved in folding and surface expression of HLA-B27. Int. Immunol. 18, 211–220 (2006).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/intimm/dxh364&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16361312&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234436500020&link_type=ISI) 

52. 52.Stewart-Jones, G. B. E. et al. Structures of Three HIV-1 HLA-B*5703-Peptide Complexes and Identification of Related HLAs Potentially Associated with Long-Term Nonprogression. The Journal of Immunology vol. 175 2459–2468 (2005).
    
    

53. 53.Archbold, J. K. et al. Natural micropolymorphism in human leukocyte antigens provides a basis for genetic control of antigen recognition. J. Exp. Med. 206, 209–219 (2009).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamVtIjtzOjU6InJlc2lkIjtzOjk6IjIwNi8xLzIwOSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA3LzE4LzIwMjAuMDcuMTYuMjAxNTU2MDYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

54. 54.Gaiha, G. D. et al. Structural topology defines protective CD8+ T cell epitopes in the HIV proteome. Science 364, 480–484 (2019).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjQvNjQzOS80ODAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNy8xOC8yMDIwLjA3LjE2LjIwMTU1NjA2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

55. 55.Macdonald, W. A. et al. A naturally selected dimorphism within the HLA-B44 supertype alters class I structure, peptide repertoire, and T cell recognition. J. Exp. Med. 198, 679–691 (2003).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamVtIjtzOjU6InJlc2lkIjtzOjk6IjE5OC81LzY3OSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA3LzE4LzIwMjAuMDcuMTYuMjAxNTU2MDYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

56. 56.Kloverpris, H. N. et al. HLA-B*57 Micropolymorphism Shapes HLA Allele-Specific Epitope Immunogenicity, Selection Pressure, and HIV Immune Control. Journal of Virology vol. 86 919–929 (2012).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoianZpIjtzOjU6InJlc2lkIjtzOjg6Ijg2LzIvOTE5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDcvMTgvMjAyMC4wNy4xNi4yMDE1NTYwNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

57. 57.Carrington, M. & Walker, B. D. Immunogenetics of spontaneous control of HIV. Annu. Rev. Med. 63, 131–145 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-med-062909-130018&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22248321&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000301838400008&link_type=ISI) 

58. 58.Khera, A. V. et al. Genome-wide polygenic scores for common diseases identify individuals with risk equivalent to monogenic mutations. Nat. Genet. 50, 1219–1224 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018&ndash;0183-z&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30104762&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

59. 59.Khera, A. V. et al. Polygenic Prediction of Weight and Obesity Trajectories from Birth to Adulthood. Cell 177, 587–596.e9 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2019.03.028&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31002795&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

60. 60.Torkamani, A. & Topol, E. Polygenic Risk Scores Expand to Obesity. Cell vol. 177 518–520 (2019).
    
    

61. 61.Martin, A. R. et al. Clinical use of current polygenic risk scores may exacerbate health disparities. Nat. Genet. 51, 584–591 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0379-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30926966&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

62. 62.Van der Auwera, G. A. et al. From FastQ data to high confidence variant calls: the Genome Analysis Toolkit best practices pipeline. Curr. Protoc. Bioinformatics 43, 11.10.1–33 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=doi:10.1002/0471250953.bi1110s43&link_type=DOI) 

63. 63.Julg, B. et al. Possession of HLA class II DRB1*1303 associates with reduced viral loads in chronic HIV-1 clade C and B infection. J. Infect. Dis. 203, 803–809 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiq122&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21257739&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

64. 64.Schäfer, C., Schmidt, A. H. & Sauter, J. Hapl-o-Mat: open-source software for HLA haplotype frequency estimation from ambiguous and heterogeneous data. BMC Bioinformatics 18, 284 (2017).
    
    

65. 65.Alexander, D. H., Novembre, J. & Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19, 1655–1664 (2009).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjk6IjE5LzkvMTY1NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA3LzE4LzIwMjAuMDcuMTYuMjAxNTU2MDYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

66. 66.Price, A. L. et al. Long-range LD can confound genome scans in admixed populations. American journal of human genetics vol. 83 132–5; author reply 135–9 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2008.06.009&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000257784000021&link_type=ISI) 

67. 67.Pasaniuc, B. et al. Analysis of Latino populations from GALA and MEC studies reveals genomic loci with biased local ancestry estimation. Bioinformatics 29, 1407–1415 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btt166&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23572411&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000319428600007&link_type=ISI) 

68. 68.McLaren, P. J. et al. Association study of common genetic variants and HIV-1 acquisition in 6,300 infected cases and 7,200 controls. PLoS Pathog. 9, e1003515 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1003515&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23935489&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom) 

69. 69.Okada, Y. et al. Contribution of a Non-classical HLA Gene, HLA-DOA, to the Risk of Rheumatoid Arthritis. Am. J. Hum. Genet. 99, 366–374 (2016).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F07%2F18%2F2020.07.16.20155606.atom)

 [1]: /embed/graphic-10.gif
 [2]: /embed/graphic-11.gif
 [3]: /embed/graphic-12.gif