Certain phosphorylated peptides are differentially presented by major histocompatibility complex (MHC) molecules on cancer cells characterized by aberrant phosphorylation. Phosphopeptides presented in complex with the human leukocyte antigen HLA-A*02:01 provide a stability advantage over their non-phosphorylated counterparts. This stability is thought to contribute to enhanced immunogenicity. Whether tumor-associated phosphopeptides presented by other common alleles exhibit immunogenicity and structural characteristics similar to those presented by A*02:01 is unclear. Therefore, we determined the identity, structural features, and immunogenicity of phosphopeptides presented by the prevalent alleles HLA-A*03:01, HLA-A*11:01, HLA-C*07:01, and HLA-C*07:02.

T-cell immunotherapies are often thwarted by the limited presentation of tumor-specific antigens abetted by the downregulation of human leukocyte antigen (HLA). We showed that drugs inhibiting ALK and RET produced dose-related increases in cell-surface HLA in tumor cells bearing these mutated kinases in vitro and in vivo, as well as elevated transcript and protein expression of HLA and other antigen-processing machinery. Subsequent analysis of HLA-presented peptides after ALK and RET inhibitor treatment identified large changes in the immunopeptidome with the appearance of hundreds of new antigens, including T-cell epitopes associated with impaired peptide processing (TEIPP) peptides. ALK inhibition additionally decreased PD-L1 levels by 75%. Therefore, these oncogenes may enhance cancer formation by allowing tumors to evade the immune system by downregulating HLA expression. Altogether, RET and ALK inhibitors could enhance T-cell-based immunotherapies by upregulating HLA, decreasing checkpoint blockade ligands, and revealing new, immunogenic, cancer-associated antigens.

HLA-E shares several peptide sequences with HLA-class Ia molecules. Therefore, anti-HLA-E antibodies that recognize the shared sequences may bind to HLA-class Ia alleles. This hypothesis was validated with a murine anti-HLA-E monoclonal antibody (mAb) MEM-E/02, which reacted with microbeads coated with several HLA-B and HLA-C antigens. In this report, the hypothesis was reexamined with another mAb 3D12, considered to be specific for HLA-E. The antibody binding is evaluated by measuring mean fluorescence index [MFI] with Luminex Multiplex Flow-Cytometric technology. The peptide-inhibition experiments are carried out with synthetic shared peptides, most prevalent to HLA-E and HLA-Ia alleles. The results showed that mAb 3D12 simulated MEM-E/02 in recognizing several HLA-B and HLA-C antigens. Both 3D12 and MEM-E/02 did not bind to HLA-A, HLA-F and HLA-G molecules. As observed with MEM-E/02, binding of 3D12 to HLA-E is inhibited by the peptides sequences (115)QFAYDGKDY(123) and (137)DTAAQI(142). Decrease in binding of mAb 3D12 to HLA class Ia, after heat treatment of antigen coated microbeads, supports the contention that the epitope may be located at the outside of the "thermodynamically stable" α-helix conformations of HLA-E. Several sequence and structure-based web-tools were employed to validate the discontinuous epitopes recognized by the mAbs. The scores obtained by these web-tools distinguished the shared peptide sequences that inhibited the mAb binding to HLA-E. Furthermore, ElliPro web tool points out that both mAbs recognize the conformational discontinuous epitopes (the shared inhibitory peptide sequences) in the secondary structure of the HLA-E molecule. The study favors the contention that the domain of the shared inhibitory peptide sequences may be the most immunogenic site of HLA-E molecule. It also postulates and clarifies that amino acid substitution on or near the binding domains may account for the lack of cross reactivity of 3D12 and MEM-E/02 with HLA-A, HLA-F and HLA-G molecules.

The Human Leucocyte Antigens (HLA) work in concert with other immune factors to modulate immunity to viral infections. Extensive variation has been reported in the genetic sequences and functions of classical HLA class I genes in many (mostly Western) populations, and several HLA associations with infectious disease outcomes have been reported. Little is known about their role in the susceptibility or resistance to hepatitis viruses in Central African populations. The aim of this study was to determine variants of two HLA class I genes (HLA-A and -C) in adults infected with hepatitis B (HBV)- or -C (HCV) virus in Cameroon. In this case-control study, a total of 169 unrelated adults comprising 68 HCV-infected, 38 HBV-infected and 63 uninfected (controls) individuals participated. Each consented participant was screened for HBV, HCV, and HIV infections and willingly donated a single blood sample for genomic DNA isolation and some clinical laboratory tests. HLA-A and HLA-C were genotyped using previously described sequence-based techniques (SBT). A total of 54 HLA alleles were identified in the study population (27 HLA-A and 27 HLA-C). HLA-A∗23:01 and HLA-C∗07:01 were the most common alleles with genotype frequencies of 31.4% and 29.3%, respectively. Hepatitis individuals were six times more likely to be HLA-A∗30:01 carriers than uninfected controls (OR = 6.30, p = 0.020 (HBV); OR = 6.21, p = 0.010 (HCV), respectively). Similarly, carriers of HLA-C∗17:01 were over-represented in the HBV-infected compared to the uninfected control group (21.9% vs. 6.4%, respectively) suggesting that this allele could play a role in the susceptibility to HBV infection. These findings demonstrate that carriers of HLA-A∗30:01 were over-represented in the hepatitis group compared to uninfected controls while HLA-C∗17:01 was completely absent in the HCV + group.

The etiology of Posner-Schlossman syndrome (PSS) remains unknown. The association of human leukocyte antigens (HLA) allelic diversity with PSS has been poorly investigated. To evaluate the association of allelic polymorphisms of class I HLA-A, -B and -C and class II HLA-DRB1 and -DQB1 with PSS, 100 unrelated patients with PSS and 128 age- and ethnically matched control subjects were recruited from a southern Chinese Han population. Polymorphisms in exons 2-4 for HLA-A, -B, -C loci, exon 2 for HLA-DRB1 and exons 2,3 for HLA-DQB1 were analyzed for association with PSS at allele and haplotype levels. The allele frequency of HLA-C*1402 in PSS patients was significantly higher than that in controls (P = 0.002, OR = 4.12). This association survived the Bonferroni correction (Pc = 0.04). The allele frequency of HLA-B*1301 in PSS patients was lower than that in the control group (P = 0.003, OR = 0.21), although this association did not survive the Bonferroni correction (Pc = 0.16). In PSS patients, the haplotype frequencies of HLA-A*1101~C*1402 and B*5101~C*1402 were higher than that in controls (P = 0.03, OR = 4.44; P = 0.02, OR = 3.20; respectively), while the HLA-B*1301~C*0304 was lower than that in controls (P = 0.007, OR = 0.23), although these associations did not survive the Bonferroni correction (Pc > 0.16). This study for the first time demonstrated that polymorphisms at the HLA-B and HLA-C loci were nominally associated with PSS in the southern Chinese Han population. Our results suggest that HLA-C*1402, A*1101~C*1402 and B*5101~C*1402 might be risk factors for PSS, whereas HLA-B*1301 plus B*1301~C*0304 might be protective factors against PSS, but even larger datasets are required to confirm these findings. Findings from this study provide valuable new clues for investigating the mechanisms and development of new diagnosis and treatment for PSS.

Heterologous immunity of virus-specific T cells poses a potential barrier to transplantation tolerance. Cross-reactivity to HLA-A and -B molecules has broadly been described, whereas responses to allo-HLA-C have remained ill defined. In contrast to the transplant setting, HLA-C is the only polymorphic HLA molecule expressed by extravillous trophoblasts at the maternal-fetal interface during pregnancy. Uncontrolled placental viral infections, accompanied by a pro-inflammatory milieu, can alter the activation status and stability of effector T cells. Potential cross-reactivity of maternal decidual virus-specific T cells to fetal allo-HLA-C may thereby have detrimental consequences for the success of pregnancy. To explore the presence of cross-reactivity to HLA-C and the other non-classical HLA antigens expressed by trophoblasts, HLA-A and -B-restricted CD8+ T cells specific for Epstein-Barr virus, Cytomegalovirus, Varicella-Zoster virus, and Influenza virus were tested against target cells expressing HLA-C, -E, and -G molecules. An HLA-B*08:01-restricted EBV-specific T cell clone displayed cross-reactivity against HLA-C*01:02. Furthermore, cross-reactivity of HLA-C-restricted virus-specific CD8+ T cells was observed for HCMV HLA-C*06:02/TRA CD8+ T cell lines and clones against HLA-C*03:02. Collectively, these results demonstrate that cross-reactivity against HLA-C can occur and thereby may affect pregnancy outcome.

The present study investigates the human leucocyte antigen (HLA) allele and haplotype frequencies in Japanese population. We carried out the frequency analysis in 5824 families living across Japanese archipelago. The studied population has mainly been typed for the purpose of transplant, especially the hematopoietic stem cell transplantation (HSCT). We determined HLA class I (A, B, and C) and HLA class II (DRB1) using Luminex technology. The haplotypes were directly counted by segregation. A total of 44 HLA-A, 29 HLA-C, 75 HLA-B, and 42 HLA-DRB1 alleles were identified. In the HLA haplotypes of A-C-B-DRB1 and C-B, the pattern of linkage disequilibrium peculiar to Japanese population has been confirmed. Moreover, the haplotype frequencies based on family study was compared with the frequencies estimated by maximum likelihood estimation (MLE), and the equivalent results were obtained. The allele and haplotype frequencies obtained in this study could be useful for anthropology, transplantation therapy, and disease association studies.

HLA class I glycoproteins contain the functional sites that bind peptide antigens and engage lymphocyte receptors. Recently, clinical application of sequence-based HLA typing has uncovered an unprecedented number of novel HLA class I alleles. Here we define the nature and extent of the variation in 3,489 HLA-A, 4,356 HLA-B and 3,111 HLA-C alleles. This analysis required development of suites of methods, having general applicability, for comparing and analyzing large numbers of homologous sequences. At least three amino-acid substitutions are present at every position in the polymorphic α1 and α2 domains of HLA-A, -B and -C. A minority of positions have an incidence >1% for the 'second' most frequent nucleotide, comprising 70 positions in HLA-A, 85 in HLA-B and 54 in HLA-C. The majority of these positions have three or four alternative nucleotides. These positions were subject to positive selection and correspond to binding sites for peptides and receptors. Most alleles of HLA class I (>80%) are very rare, often identified in one person or family, and they differ by point mutation from older, more common alleles. These alleles with single nucleotide polymorphisms reflect the germ-line mutation rate. Their frequency predicts the human population harbors 8-9 million HLA class I variants. The common alleles of human populations comprise 42 core alleles, which represent all selected polymorphism, and recombinants that have assorted this polymorphism.

The individual cellular immune response to intracellular antigens is modeled by the highly polymorphic major histocompatibility complex (HLA) class I molecules. The epitopes presented and the T cell repertoire that recognizes them depend on the HLA constitution of the individual. Therefore, to monitor and to modify an individual's HLA class I-driven cellular immune response, it is necessary to know the HLA class I alleles of the person and the possible epitopes of the target antigen presented by those alleles. In particular, this is necessary in order to design peptide-based vaccines and immune therapies for the treatment of diseases caused by viruses, intracellular parasites or cancer, and to monitor the immune response during those treatments. We describe a new set of HLA-A, -B, and -C locus-specific primers for the polymerase chain reaction (PCR) amplification of the whole coding sequence of these genes from complementary DNA (cDNA). We describe their use for typing and for the production of a library of recombinant HLA class I genes. We discuss two downstream applications of this gene collection: production of soluble HLA molecules and discovery of new epitopes.

Human leucocyte antigen (HLA) class I molecules present endogenously processed antigens to T-cells and have been linked to differences in HIV-1 disease progression. HLA allelotypes show considerable geographical and inter-individual variation, as does the rate of progression of HIV-1 disease, with long-term non-progression (LTNP) of disease having most evidence of an underlying genetic contribution. However, most genetic analyses of LTNP have occurred in adults of European ancestry, limiting the potential transferability of observed associations to diverse populations who carry the burden of disease. This is particularly true of HIV-1 infected children. Here, using exome sequencing (ES) to infer HLA allelotypes, we determine associations with HIV-1 LTNP in two diverse African pediatric populations. We performed a case-control association study of 394 LTNPs and 420 rapid progressors retrospectively identified from electronic medical records of pediatric HIV-1 populations in Uganda and Botswana. We utilized high-depth ES to perform high-resolution HLA allelotyping and assessed evidence of association between HLA class I alleles and LTNP. Sixteen HLA alleles and haplotypes had significantly different frequencies between Uganda and Botswana, with allelic differences being more prominent in HLA-A compared to HLA-B and C allelotypes. Three HLA allelotypes showed association with LTNP, including a novel association in HLA-C (HLA-B∗57:03, aOR 3.21, Pc = 0.0259; B∗58:01, aOR 1.89, Pc = 0.033; C∗03:02, aOR 4.74, Pc = 0.033). Together, these alleles convey an estimated population attributable risk (PAR) of non-progression of 16.5%. We also observed novel haplotype associations with HLA-B∗57:03-C∗07:01 (aOR 5.40, Pc = 0.025) and HLA-B∗58:01-C∗03:02 (aOR 4.88, Pc = 0.011) with a PAR of 9.8%, as well as a previously unreported independent additive effect and heterozygote advantage of HLA-C∗03:02 with B∗58:01 (aOR 4.15, Pc = 0.005) that appears to limit disease progression, despite weak LD (r 2 = 0.18) between these alleles. These associations remained irrespective of gender or country. In one of the largest studies of HIV in Africa, we find evidence of a protective effect of canonical HLA-B alleles and a novel HLA-C association that appears to augment existing HIV-1 control alleles in pediatric populations. Our findings outline the value of using multi-ethnic populations in genetic studies and offer a novel HIV-1 association of relevance to ongoing vaccine studies.

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.

T-cell responses to minor histocompatibility antigens (mHAs) mediate graft-versus-leukemia (GVL) effects and graft-versus-host disease (GVHD) in allogeneic hematopoietic cell transplantation. Therapies that boost T-cell responses improve allogeneic hematopoietic cell transplant (alloHCT) efficacy but are limited by concurrent increases in the incidence and severity of GVHD. mHAs with expression restricted to hematopoietic tissue (GVL mHAs) are attractive targets for driving GVL without causing GVHD. Prior work to identify mHAs has focused on a small set of mHAs or population-level single-nucleotide polymorphism-association studies. We report the discovery of a large set of novel GVL mHAs based on predicted immunogenicity, tissue expression, and degree of sharing among donor-recipient pairs (DRPs) in the DISCOVeRY-BMT data set of 3231 alloHCT DRPs. The total number of predicted mHAs varied by HLA allele, and the total number and number of each class of mHA significantly differed by recipient genomic ancestry group. From the pool of predicted mHAs, we identified the smallest sets of GVL mHAs needed to cover 100% of DRPs with a given HLA allele. We used mass spectrometry to search for high-population frequency mHAs for 3 common HLA alleles. We validated 24 predicted novel GVL mHAs that are found cumulatively within 98.8%, 60.7%, and 78.9% of DRPs within DISCOVeRY-BMT that express HLA-A∗02:01, HLA-B∗35:01, and HLA-C∗07:02, respectively. We confirmed the immunogenicity of an example novel mHA via T-cell coculture with peptide-pulsed dendritic cells. This work demonstrates that the identification of shared mHAs is a feasible and promising technique for expanding mHA-targeting immunotherapeutics.

The frequency of recurrence of subacute thyroiditis (SAT) is rather high, reaching 20⁻30%. The reason for SAT relapse is still unknown. Recently, we have demonstrated the association between SAT and the presence of HLA-B*18:01, DRB1*01, and C*04:01, apart from the previously known HLA-B*35. The aim of the present study was to evaluate the correlation between SAT-associated HLA haplotypes and the risk of SAT recurrence. HLA-A, -B, -C, -DQB1 and -DRB1 were genotyped using a next-generation sequencing method in 49 SAT patients. The patients were divided into the following HLA groups: 1. HLA-B*35 and/or HLA-C*04, but without any other of the analyzed antigens; 2. HLA-DRB1*01, regardless of the co-presence of HLA-B*35 or -C*04:01, but without HLA-B*18:01; 3. HLA-B18 only, without any other antigen; 4. HLA-B*18:01 plus -B*35, regardless of the presence of any other analyzed antigens. The recurrence rate was compared between the groups. The recurrence rate was significantly increased in patients with HLA-B*18:01 plus HLA-B*35. In conclusion, the risk of SAT recurrence was HLA-dependent and the determining factor was the co-presence of HLA-B*18:01 and -B*35. In such high-risk patients, the steroid treatment regimen should be intensified with a slower dose reduction.

Detection of alloreactive anti-HLA antibodies is a frequent and mandatory test before and after organ transplantation to determine the antigenic targets of the antibodies. Nowadays, this test involves the measurement of fluorescent signals generated through antibody-antigen reactions on multi-beads flow cytometers. In this study, in a cohort of 1,066 patients from one country, anti-HLA class I responses were analyzed on a panel of 98 different antigens. Knowing that the immune system responds typically to "shared" antigenic targets, we studied the clustering patterns of antibody responses against HLA class I antigens without any a priori hypothesis, applying two unsupervised machine learning approaches. At first, the principal component analysis (PCA) projections of intra-locus specific responses showed that anti-HLA-A and anti-HLA-C were the most distantly projected responses in the population with the anti-HLA-B responses to be projected between them. When PCA was applied on the responses against antigens belonging to a single locus, some already known groupings were confirmed while several new cross-reactive patterns of alloreactivity were detected. Anti-HLA-A responses projected through PCA suggested that three cross-reactive groups accounted for about 70% of the variance observed in the population, while anti-HLA-B responses were mainly characterized by a distinction between previously described Bw4 and Bw6 cross-reactive groups followed by several yet undocumented or poorly described ones. Furthermore, anti-HLA-C responses could be explained by two major cross-reactive groups completely overlapping with previously described C1 and C2 allelic groups. A second feature-based analysis of all antigenic specificities, projected as a dendrogram, generated a robust measure of allelic antigenic distances depicting bead-array defined cross reactive groups. Finally, amino acid combinations explaining major population specific cross-reactive groups were described. The interpretation of the results was based on the current knowledge of the antigenic targets of the antibodies as they have been characterized either experimentally or computationally and appear at the HLA epitope registry.

Induced pluripotent stem cells (iPSCs) have strong potential in regenerative medicine applications; however, immune rejection caused by HLA mismatching is a concern. B2M gene knockout and HLA-homozygous iPSC stocks can address this issue, but the former approach may induce NK cell activity and fail to present antigens, and it is challenging to recruit rare donors for the latter method. Here, we show two genome-editing strategies for making immunocompatible donor iPSCs. First, we generated HLA pseudo-homozygous iPSCs with allele-specific editing of HLA heterozygous iPSCs. Second, we generated HLA-C-retained iPSCs by disrupting both HLA-A and -B alleles to suppress the NK cell response while maintaining antigen presentation. HLA-C-retained iPSCs could evade T cells and NK cells in vitro and in vivo. We estimated that 12 lines of HLA-C-retained iPSCs combined with HLA-class II knockout are immunologically compatible with >90% of the world's population, greatly facilitating iPSC-based regenerative medicine applications.

The HLA (Human Leukocyte Antigens) genes are well-documented targets of balancing selection, and variation at these loci is associated with many disease phenotypes. Variation in expression levels also influences disease susceptibility and resistance, but little information exists about the regulation and population-level patterns of expression. This results from the difficulty in mapping short reads originated from these highly polymorphic loci, and in accounting for the existence of several paralogues. We developed a computational pipeline to accurately estimate expression for HLA genes based on RNA-seq, improving both locus-level and allele-level estimates. First, reads are aligned to all known HLA sequences in order to infer HLA genotypes, then quantification of expression is carried out using a personalized index. We use simulations to show that expression estimates obtained in this way are not biased due to divergence from the reference genome. We applied our pipeline to the GEUVADIS dataset, and compared the quantifications to those obtained with reference transcriptome. Although the personalized pipeline recovers more reads, we found that using the reference transcriptome produces estimates similar to the personalized pipeline (r ≥ 0.87) with the exception of HLA-DQA1. We describe the impact of the HLA-personalized approach on downstream analyses for nine classical HLA loci (HLA-A, HLA-C, HLA-B, HLA-DRA, HLA-DRB1, HLA-DQA1, HLA-DQB1, HLA-DPA1, HLA-DPB1). Although the influence of the HLA-personalized approach is modest for eQTL mapping, the p-values and the causality of the eQTLs obtained are better than when the reference transcriptome is used. We investigate how the eQTLs we identified explain variation in expression among lineages of HLA alleles. Finally, we discuss possible causes underlying differences between expression estimates obtained using RNA-seq, antibody-based approaches and qPCR.

Extracellular vesicles can modulate diverse processes ranging from proliferation and tissue repair, to chemo-resistance and cellular differentiation. With the advent of tissue and immunological targeting, extracellular vesicles are also increasingly viewed as promising vectors to deliver peptide-based cancer antigens to the human immune system. Despite the clinical relevance and therapeutic potential of such 'cell-free' approaches, the natural antigen presentation landscape exported in extracellular vesicles is still largely uncharted, due to the challenging nature of such preparations and analyses. In the context of therapeutic vesicle production, a critical evaluation of the similarity in vesicular antigen presentation is also urgently needed. In this work, we compared the HLA-I peptide ligandomes of extracellular vesicles against that of whole-cells of the same cell line. We found that extracellular vesicles not only over-represent HLA-B complexes and peptide ligands, but also cysteinylated peptides that may modulate immune responses. Collectively, these findings describe the pre-existing provision of vesicular HLA complexes that may be utilized to carry peptide vaccines, as well as the propensity for different peptide and post-translationally modified ligands to be presented, and will outline critical considerations in devising novel EV vaccination strategies.

Human HLA-E can, in addition to self-antigens, also present pathogen-derived sequences, which elicit specific T-cell responses. T-cells recognize their antigen presented by HLA-E highly specifically and have unique functional and phenotypical properties. Pathogen specific HLA-E restricted CD8+ T-cells are an interesting new player in the field of immunology. Future work should address their exact roles and relative contributions in the immune response against infectious diseases.

The physiological functions of natural killer (NK) cells in human immunity and reproduction depend upon diverse interactions between killer cell immunoglobulin-like receptors (KIRs) and their HLA class I ligands: HLA-A, HLA-B, and HLA-C. The genomic regions containing the KIR and HLA class I genes are unlinked, structurally complex, and highly polymorphic. They are also strongly associated with a wide spectrum of diseases, including infections, autoimmune disorders, cancers, and pregnancy disorders, as well as the efficacy of transplantation and other immunotherapies. To facilitate study of these extraordinary genes, we developed a method that captures, sequences, and analyzes the 13 KIR genes and HLA-A, HLA-B, and HLA-C from genomic DNA. We also devised a bioinformatics pipeline that attributes sequencing reads to specific KIR genes, determines copy number by read depth, and calls high-resolution genotypes for each KIR gene. We validated this method by using DNA from well-characterized cell lines, comparing it to established methods of HLA and KIR genotyping, and determining KIR genotypes from 1000 Genomes sequence data. This identified 116 previously uncharacterized KIR alleles, which were all demonstrated to be authentic by sequencing from source DNA via standard methods. Analysis of just two KIR genes showed that 22% of the 1000 Genomes individuals have a previously uncharacterized allele or a structural variant. The method we describe is suited to the large-scale analyses that are needed for characterizing human populations and defining the precise HLA and KIR factors associated with disease. The methods are applicable to other highly polymorphic genes.

HLA-B∗46:01 was formed by an intergenic mini-conversion, between HLA-B∗15:01 and HLA-C∗01:02, in Southeast Asia during the last 50,000 years, and it has since become the most common HLA-B allele in the region. A functional effect of the mini-conversion was introduction of the C1 epitope into HLA-B∗46:01, making it an exceptional HLA-B allotype that is recognized by the C1-specific natural killer (NK) cell receptor KIR2DL3. High-resolution mass spectrometry showed that HLA-B∗46:01 has a low-diversity peptidome that is distinct from those of its parents. A minority (21%) of HLA-B∗46:01 peptides, with common C-terminal characteristics, form ligands for KIR2DL3. The HLA-B∗46:01 peptidome is predicted to be enriched for peptide antigens derived from Mycobacterium leprae. Overall, the results indicate that the distinctive peptidome and functions of HLA-B∗46:01 provide carriers with resistance to leprosy, which drove its rapid rise in frequency in Southeast Asia.