Pocketcheck: updating the HLA class I peptide specificity roadmap

Behçet's Disease: A Comprehensive Review on the Role of HLA-B*51, Antigen Presentation, and Inflammatory Cascade

Behçet's disease (BD) is a complex, recurring inflammatory disorder with autoinflammatory and autoimmune components. This comprehensive review aims to explore BD's pathogenesis, focusing on established genetic factors. Studies reveal that HLA-B*51 is the primary genetic risk factor, but non-HLA genes (ERAP1, IL-10, IL23R/IL-12RB2), as well as innate immunity genes (FUT2, MICA, TLRs), also contribute. Genome-wide studies emphasize the significance of ERAP1 and HLA-I epistasis. These variants influence antigen presentation, enzymatic activity, and HLA-I peptidomes, potentially leading to distinct autoimmune responses. We conducted a systematic review of the literature to identify studies exploring the association between HLA-B*51 and BD and further highlighted the roles of innate and adaptive immunity in BD. Dysregulations in Th1/Th2 and Th17/Th1 ratios, heightened clonal cytotoxic (CD8+) T cells, and reduced T regulatory cells characterize BD's complex immune responses. Various immune cell types (neutrophils, γδ T cells, natural killer cells) further contribute by releasing cytokines (IL-17, IL-8, GM-CSF) that enhance neutrophil activation and mediate interactions between innate and adaptive immunity. In summary, this review advances our understanding of BD pathogenesis while acknowledging the research limitations. Further exploration of genetic interactions, immune dysregulation, and immune cell roles is crucial. Future studies may unveil novel diagnostic and therapeutic strategies, offering improved management for this complex disease.

Contemplating immunopeptidomes to better predict them

The identification of T-cell epitopes is key for a complete molecular understanding of immune recognition mechanisms in infectious diseases, autoimmunity and cancer. T-cell epitopes further provide targets for personalized vaccines and T-cell therapy, with several therapeutic applications in cancer immunotherapy and elsewhere. T-cell epitopes consist of short peptides displayed on Major Histocompatibility Complex (MHC) molecules. The recent advances in mass spectrometry (MS) based technologies to profile the ensemble of peptides displayed on MHC molecules - the so-called immunopeptidome - had a major impact on our understanding of antigen presentation and MHC ligands. On the one hand, these techniques enabled researchers to directly identify hundreds of thousands of peptides presented on MHC molecules, including some that elicited T-cell recognition. On the other hand, the data collected in these experiments revealed fundamental properties of antigen presentation pathways and significantly improved our ability to predict naturally presented MHC ligands and T-cell epitopes across the wide spectrum of MHC alleles found in human and other organisms. Here we review recent computational developments to analyze experimentally determined immunopeptidomes and harness these data to improve our understanding of antigen presentation and MHC binding specificities, as well as our ability to predict MHC ligands. We further discuss the strengths and limitations of the latest approaches to move beyond predictions of antigen presentation and tackle the challenges of predicting TCR recognition and immunogenicity.

A plausible contributor to multiple sclerosis; presentation of antigenic myelin protein epitopes by major histocompatibility complexes

BACKGROUND

Multiple sclerosis (MS) can be induced upon successful presentation of myelin antigens by MHC I/II. Antigenic similarity between the myelin and viral proteins may worsen the immunological responses.

METHODOLOGY

Antigenic regions within myelin proteins; PLP1, MBP, MOG, and MAG were analyzed using SVMTrip and EMBOSS. Homology search identified sequence similarity between the predicted host epitopes and viral proteins. NetMHCpan predicted MHC I/II binding followed by peptide-protein docking through the HPEPDOCK server. Thereafter we analyzed conformational flexibility and stability of 15 protein-peptide complexes based on high docking scores. The binding free energy was calculated using conventional (MD) and Gaussian accelerated molecular dynamics simulation.

RESULTS

PLP1, MBP, MAG and MOG contained numerous antigenic epitopes. MBP and MOG epitopes had sequence similarity to HHV-6 BALF5; EBNA1 and CMV glycoprotein M (gM), and EBV LMP2B, gp350/220; HHV-8 ORFs respectively. Many herpes virus proteins like tegument, envelope glycoproteins, and ORFs of EBV, CMV, HHV-6, and HHV-8 demonstrated sequence similarity with MAG and PLP1. Some antigenic peptides were also linear B-cell epitopes and influenced cytokine production by T-cell. MHC I allele HLA-B*57:01 bound to PLP1 peptide and HLA-A*68:02 bound to a MAG peptide strongly. MHC II alleles HLA-DRB1*04:05 and HLA-DR1*01:01 associated with MAG- and MOG-derived peptides, respectively, demonstrating high HPEPDOCK scores. MD simulations established stable binding of certain peptides with the MHC namely HLA-B*51:01-MBP(DYKSAHKGFKGVDAQGTLSKIFKL), HLA-B*57:01-PLP1(PDKFVGITYALTVVWLLVFACSAVPVYIYF), HLA-DR1*01:01-MOG(VEDPFYWVSPGVLVLLAVLPVLLLQITVGLVFLCLQYR) and HLA-DRB1*04:05-MAG(TWVQVSLLHFVPTREA).

CONCLUSIONS

Cross-reactivity between self-antigens and pathogen derived immunodominant epitopes may induce MS. Our study supported the role of specific MHC alleles as a contributing MS risk factor.

Proteomic Profiling and T Cell Receptor Usage of Abacavir Susceptible Subjects

Type B adverse drug reactions (ADRs) represent a significant threat as their occurrence arises unpredictable and despite proper application of the drug. The severe immune reaction Abacavir Hypersensitivity Syndrome (AHS) that arises in HIV⁺ patients treated with the antiretroviral drug Abacavir (ABC) strongly correlates to the presence of the human leukocyte antigen (HLA) genotype HLA-B*57:01 and discriminates HLA-B*57:01⁺ HIV⁺ patients from ABC treatment. However, not all HLA-B*57:01⁺ HIV⁺ patients are affected by AHS, implying the involvement of further patient-specific factors in the development of AHS. The establishment of a reliable assay to classify HLA-B*57:01 carriers as ABC sensitive or ABC tolerant allowed to investigate the T cell receptor (TCR) Vβ chain repertoire of effector cells and revealed Vβ6 and Vβ24 as potential public TCRs in ABC sensitive HLA-B*57:01 carriers. Furthermore, distinct effects of ABC on the cellular proteome of ABC sensitive and tolerant volunteers were observed and suggest enhanced activation and maturation of dentritic cells (DC) in ABC sensitive volunteers. Analysis of ABC-naïve cellular proteomes identified the T cell immune regulator 1 (TCIRG1) as a potential prognostic biomarker for ABC susceptibility and the involvement of significantly upregulated proteins, particularly in peptide processing, antigen presentation, interferon (IFN), and cytokine regulation.

Unravelling the Proteomics of HLA-B*57:01+ Antigen Presenting Cells during Abacavir Medication

Type B adverse drug reactions (ADRs) are unpredictable based on the drug’s pharmacology and represent a key challenge in pharmacovigilance. For human leukocyte antigen (HLA)-mediated type B ADRs, it is assumed that the protein/small-molecule interaction alters the biophysical and mechanistic properties of the antigen presenting cells. Sophisticated methods enabled the molecular appreciation of HLA-mediated ADRs; in several instances, the drug molecule occupies part of the HLA peptide binding groove and modifies the recruited peptide repertoire thereby causing a strong T-cell-mediated immune response that is resolved upon withdrawal of medication. The severe ADR in HLA-B*57:01⁺ patients treated with the antiretroviral drug abacavir (ABC) in anti-HIV therapy is an example of HLA-drug-T cell cooperation. However, the long-term damages of the HLA-B*57:01-expressing immune cells following ABC treatment remain unexplained. Utilizing full proteome sequencing following ABC treatment of HLA-B*57:01⁺ cells, we demonstrate stringent proteomic alteration of the HLA/drug presenting cells. The proteomic content indisputably reflects the cellular condition; this knowledge directs towards individual pharmacovigilance for the development of personalized and safe medication.

The Loss of HLA-F/KIR3DS1 Ligation Is Mediated by Hemoglobin Peptides

The human leukocyte antigen (HLA)-Ib molecule, HLA-F, is known as a CD4⁺ T-cell protein and mediator of HIV progression. While HLA-Ia molecules do not have the chance to select and present viral peptides for immune recognition due to protein downregulation, HLA-F is upregulated. Post HIV infection, HLA-F loses the affinity to its activating receptor KIR3DS1 on NK cells leading to progression of the HIV infection. Several studies aimed to solve the question of the biophysical interface between HLA ligands and their cognate receptors. It became clear that even an invariant HLA molecule can be structurally modified by the variability of the bound peptide. We recently discovered the ability of HLA-F to select and present peptides and the HLA-F allele-specific peptide selection from the proteomic content using soluble HLA (sHLA) technology and a sophisticated MS method. We established recombinant K562 cells that express membrane-bound HLA-F*01:01, 01:03 or 01:04 complexes. While a recombinant soluble form of KIR3DS1 did not bind to the peptide-HLA-F complexes, acid elution of the peptides resulted in the presentation of HLA-F open conformers, and the binding of the soluble KIR3DS1 receptor increased. We used CD4⁺/HIV⁻ and CD4⁺/HIV⁺ cells and performed an MS proteome analysis. We could detect hemoglobin as significantly upregulated in CD4⁺ T-cells post HIV infection. The expression of cellular hemoglobin in nonerythroid cells has been described, yet HLA-Ib presentation of hemoglobin-derived peptides is novel. Peptide sequence analysis from HLA-F allelic variants featured hemoglobin peptides as dominant and shared. The reciprocal experiment of binding hemoglobin peptide fractions to the HLA-F open conformers resulted in significantly diminished receptor recognition. These results underpin the molecular involvement of HLA-F and its designated peptide ligand in HIV immune escape.

HLA-F Allele-Specific Peptide Restriction Represents an Exceptional Proteomic Footprint

Peptide-dependent engagement between human leucocyte antigens class I (HLA-I) molecules and their cognate receptors has been extensively analyzed. HLA-F belongs to the non-classical HLA-Ib molecules with marginal polymorphic nature and tissue restricted distribution. The three common allelic variants HLA-F*01:01/01:03/01:04 are distinguished by polymorphism outside the peptide binding pockets (residue 50, α1 or residue 251, α3) and are therefore not considered relevant for attention. However, peptide selection and presentation undergoes a most elaborated extraction from the whole available proteome. It is known that HLA-F confers a beneficial effect on disease outcome during HIV-1 infections. The interaction with the NK cell receptor initiates an antiviral downstream immune response and lead to delayed disease progression. During the time of HIV infection, HLA-F expression is upregulated, while its interaction with KIR3DS1 is diminished. The non-polymorphic nature of HLA-F facilitates the conclusion that understanding HLA-F peptide selection and presentation is essential to a comprehensive understanding of this dynamic immune response. Utilizing soluble HLA technology we recovered stable pHLA-F*01:01, 01:03 and 01:04 complexes from K562 cells and analyzed the peptides presented. Utilizing a sophisticated LC-MS-method, we analyzed the complete K562 proteome and matched the peptides presented by the respective HLA-F subtypes with detected proteins. All peptides featured a length of 8 to 24 amino acids and are not N-terminally anchored; the C-terminus is preferably anchored by Lys. To comprehend the alteration of the pHLA-F surface we structurally compared HLA-F variants bound to selected peptides. The peptides were selected from the same cellular content; however, no overlap between the proteomic source of F*01:01, 01:03 or 01:04 selected peptides could be observed. Recognizing the balance between HLA-F expression, HLA-F polymorphism and peptide selection will support to understand the role of HLA-F in viral pathogenesis.

The Mechanistic Differences in HLA-Associated Carbamazepine Hypersensitivity

Drug hypersensitivity reactions that resemble acute immune reactions are linked to certain human leucocyte antigen (HLA) alleles. Severe and life-threatening Stevens Johnson Syndrome and Toxic Epidermal Necrolysis following treatment with the antiepileptic and psychotropic drug Carbamazepine are associated with HLA-B*15:02; whereas carriers of HLA-A*31:01 develop milder symptoms. It is not understood how these immunogenic differences emerge genotype-specific. For HLA-B*15:02 an altered peptide presentation has been described following exposure to the main metabolite of carbamazepine that is binding to certain amino acids in the F pocket of the HLA molecule. The difference in the molecular mechanism of these diseases has not been comprehensively analyzed, yet; and is addressed in this study. Soluble HLA-technology was utilized to examine peptide presentation of HLA-A*31:01 in presence and absence of carbamazepine and its main metabolite and to examine the mode of peptide loading. Proteome analysis of drug-treated and untreated cells was performed. Alterations in sA*31:01-presented peptides after treatment with carbamazepine revealed different half-life times of peptide-HLA- or peptide-drug-HLA complexes. Together with observed changes in the proteome elicited through carbamazepine or its metabolite these results illustrate the mechanistic differences in carbamazepine hypersensitivity for HLA-A*31:01 or B*15:02 patients and constitute the bridge between pharmacology and pharmacogenetics for personalized therapeutics.

Carbamazepine-Mediated Adverse Drug Reactions: CBZ-10,11-epoxide but Not Carbamazepine Induces the Alteration of Peptides Presented by HLA-B∗15:02

Among patients treated with the anticonvulsive and psychotropic drug carbamazepine (CBZ), approximately 10% develop severe and life-threatening adverse drug reactions. These immunological conditions are resolved upon withdrawal of the medicament, suggesting that the drug does not manifest in the body in long term. The HLA allele B∗15:02 has been described to be a genomic biomarker for CBZ-mediated immune reactions. It is not well understood if the immune reactions are triggered by the original drug or by its metabolite carbamazepine-10,11-epoxide (EPX) and how the interaction between the drug and the distinct HLA molecule occurs. Genetically engineered human B-lymphoblastoid cells expressing soluble HLA-B∗15:02 molecules were treated with the drug or its metabolite. Functional pHLA complexes were purified; peptides were eluted and sequenced. Applying mass spectrometric analysis, CBZ and EPX were monitored by analyzing the heavy chain and peptide fractions separately for the presence of the drug. This method enabled the detection of the drug in a biological situation post-pHLA assembly. Both drugs were bound to the HLA-B∗15:02 heavy chain; however, solely EPX altered the peptide-binding motif of B∗15:02-restricted peptides. This observation could be explained through structural insight; EPX binds to the peptide-binding region and alters the biochemical features of the F pocket and thus the peptide motif. Understanding the nature of immunogenic interactions between CBZ and EPX with the HLA immune complex will guide towards effective and safe medications.

The polymorphism at residue 156 determines the HLA-B*35 restricted peptide repertoire during HCMV infection

Peptide selection in infected cells is not fully understood yet, but several indications point to the fact that there are differences to uninfected cells, especially in productive HCMV infection, since HCMV evolved various strategies to disable the hosts immune system, including presentation of peptide-HLA complexes to immune effector cells. Therefore, peptide predictions for specific HLA alleles are limited in these cases and the naturally presented peptide repertoire of HCMV-infected cells is of major interest to optimize adoptive T cell therapies. The allotypes HLA-B*35:01 and B*35:08 differ at a single amino acid at position 156 and have been described to differ in their peptide features and in their association with the peptide loading complex. Virus specific T cells recognizing the allelic pHLA-B*35 complexes could be detected, indicating a significant role of this HLA subtypes in viral immunity. However, naturally selected and presented viral peptides have not been described so far. In this study, we analyzed the peptide binding repertoire for HLA-B*35:01 and HLA-B*35:08 in HCMV-infected cells. The isolated peptides from both allelic subtypes were of extraordinary length, however differed in their features, origin, and sequence. For these HCMV-originated peptides, no overlap in the peptide repertoire could be observed between the two allelic subtypes. These findings reveal the discrepancies between predicted and naturally presented immunogenic epitopes and support the need of comprehensive peptide recruitment data for personalized and effective cellular therapies.

Unveiling the Peptide Motifs of HLA-C and HLA-G from Naturally Presented Peptides and Generation of Binding Prediction Matrices

The classical HLA-C and the nonclassical HLA-E and HLA-G molecules play important roles both in the innate and adaptive immune system. Starting already during embryogenesis and continuing throughout our lives, these three Ags exert major functions in immune tolerance, defense against infections, and anticancer immune responses. Despite these important roles, identification and characterization of the peptides presented by these molecules has been lacking behind the more abundant HLA-A and HLA-B gene products. In this study, we elucidated the peptide specificities of these HLA molecules using a comprehensive analysis of naturally presented peptides. To that end, the 15 most frequently expressed HLA-C alleles as well as HLA-E*01:01 and HLA-G*01:01 were transfected into lymphoblastoid C1R cells expressing low endogenous HLA. Identification of naturally presented peptides was performed by immunoprecipitation of HLA and subsequent analysis of HLA-bound peptides by liquid chromatographic tandem mass spectrometry. Peptide motifs of HLA-C unveil anchors in position 2 or 3 with high variances between allotypes, and a less variable anchor at the C-terminal end. The previously reported small ligand repertoire of HLA-E was confirmed within our analysis, and we could show that HLA-G combines a large ligand repertoire with distinct features anchoring peptides at positions 3 and 9, supported by an auxiliary anchor in position 1 and preferred residues in positions 2 and 7. The wealth of HLA ligands resulted in prediction matrices for octa-, nona-, and decamers. Matrices were validated in terms of their binding prediction and compared with the latest NetMHC prediction algorithm NetMHCpan-3.0, which demonstrated their predictive power.

Is Behçet's disease a 'class 1-opathy'? The role of HLA-B*51 in the pathogenesis of Behçet's disease

The association between carriage of the human leucocyte antigen (HLA)-B*51 allele and development of Behçet's disease (BD) has been known since the early 1970s, but the exact mechanisms responsible for its role in pathogenesis remain much-debated. In an effort to explain the disease process, it has been suggested that BD constitutes one of a newly termed group of diseases, the 'MHC-I-opathies'. Other MHC-I-opathies include ankylosing spondylitis and HLA-B*27-associated spondyloarthropathies and HLA-C*0602-associated skin psoriasis. Recent work analysing the peptidome of HLA-B*51 suggests that altered peptide presentation by HLA-B*51 is vital to the disease process. In this review, we argue that immune receptor interactions with HLA-B*51 or the HLA-B*51-peptide complex could lead to development of inflammation in BD. The evidence for CD8⁺ T cell involvement is weak, and based on emerging studies it seems more likely that natural killer (NK) or other cell interactions, perhaps mediated by leucocyte immunoglobulin-like receptor (LILR) or killer immunoglobulin-like receptor (KIR) receptors, are culpable in pathogenesis. HLA misfolding leading directly to inflammation is another hypothesis for BD pathogenesis that deserves greater investigation. Ultimately, greater understanding of HLA-B*51's unique role in BD will probably lead to improved development of therapeutic strategies.

HLA-E: Presentation of a Broader Peptide Repertoire Impacts the Cellular Immune Response-Implications on HSCT Outcome

The HLA-E locus encodes a nonclassical class Ib molecule that serves many immune functions from inhibiting NK cells to activating CTLs. Structural analysis of HLA-E/NKG2A complexes visualized fine-tuning of protective immune responses through AA interactions between HLA-E, the bound peptide, and NKG2A/CD94. A loss of cellular protection through abrogation of the HLA-E/NKG2A engagement is dependent on the HLA-E bound peptide. The role of HLA-E in posttransplant outcomes is not well understood but might be attributed to its peptide repertoire. To investigate the self-peptide repertoire of HLA-E^∗01:01 in the absence of protective HLA class I signal peptides, we utilized soluble HLA technology in class I negative LCL cells in order to characterize HLA-E^∗01:01-bound ligands by mass-spectrometry. To understand the immunological impact of these analyzed ligands on NK cell reactivity, we performed cellular assays. Synthesized peptides were loaded onto recombinant T2 cells expressing HLA-E^∗01:01 molecules and applied in cytotoxicity assays using the leukemia derived NK cell line (NKL) as effector. HLA-E in complex with the self-peptides demonstrated a shift towards cytotoxicity and a loss of cell protection. Our data highlights the fact that the HLA-E-peptidome is not as restricted as previously thought and support the suggestion of a posttransplant role for HLA-E.

A Micropolymorphism Altering the Residue Triad 97/114/156 Determines the Relative Levels of Tapasin Independence and Distinct Peptide Profiles for HLA-A(*)24 Allotypes

While many HLA class I molecules interact directly with the peptide loading complex (PLC) for conventional loading of peptides certain class I molecules are able to present peptides in a way that circumvents the PLC components. We investigated micropolymorphisms at position 156 of HLA-A(*)24 allotypes and their effects on PLC dependence for assembly and peptide binding specificities. HLA-A(*)24:06(156Trp) and HLA-A(*)24:13(156Leu) showed high levels of cell surface expression while HLA-A(*)24:02(156Gln) was expressed at low levels in tapasin deficient cells. Peptides presented by these allelic variants showed distinct differences in features and repertoire. Immunoprecipitation experiments demonstrated all the HLA-A(*)24/156 variants to associate at similar levels with tapasin when present. Structurally, HLA-A(*)24:02 contains the residue triad Met97/His114/Gln156 and a Trp156 or Leu156 polymorphism provides tapasin independence by stabilizing these triad residues, thus generating an energetically stable and a more peptide receptive environment. Micropolymorphisms at position 156 can influence the generic peptide loading pathway for HLA-A(*)24 by altering their tapasin dependence for peptide selection. The trade-off for this tapasin independence could be the presentation of unusual ligands by these alleles, imposing significant risk following hematopoietic stem cell transplantation (HSCT).

HLA-E: a novel player for histocompatibility

The classical class I human leukocyte antigens (HLA-A, -B, and -C) present allele-specific self- or pathogenic peptides originated by intracellular processing to CD8(+) immune effector cells. Even a single mismatch in the heavy chain (hc) of an HLA class I molecule can impact on the peptide binding profile. Since HLA class I molecules are highly polymorphic and most of their polymorphisms affect the peptide binding region (PBR), it becomes obvious that systematic HLA matching is crucial in determining the outcome of transplantation. The opposite holds true for the nonclassical HLA class I molecule HLA-E. HLA-E polymorphism is restricted to two functional versions and is thought to present a limited set of highly conserved peptides derived from class I leader sequences. However, HLA-E appears to be a ligand for the innate and adaptive immune system, where the immunological response to peptide-HLA-E complexes is dictated through the sequence of the bound peptide. Structural investigations clearly demonstrate how subtle amino acid differences impact the strength and response of the cognate CD94/NKG2 or T cell receptor.

Soluble HLA technology as a strategy to evaluate the impact of HLA mismatches

HLA class I incompatibilities still remain one of the main barriers for unrelated bone marrow transplantation (BMT); hence the molecular understanding of how to mismatch patients and donors and still have successful clinical outcomes will guide towards the future of unrelated BMT. One way to estimate the magnitude of polymorphisms within the PBR is to determine which peptides can be selected by individual HLA alleles and subsequently presented for recognition by T cells. The features (structure, length, and sequence) of different peptides each confer an individual pHLA landscape and thus directly shape the individual immune response. The elution and sequencing of peptides by mass spectrometric analysis enable determining the bona fide repertoire of presented peptides for a given allele. This is an effective and simple way to compare the functions of allelic variants and make a first assessment of their degree of permissivity. We describe the methodology used for peptide sequencing and the limitations of peptide prediction tools compared to experimental methods. We highlight the altered peptide features that are observed between allelic variants and the need to discover the altered peptide repertoire in situations of "artificial" graft versus host disease (GvHD) that occur in HLA-specific hypersensitive immune responses to drugs.

Behçet disease-associated MHC class I residues implicate antigen binding and regulation of cell-mediated cytotoxicity

The HLA protein, HLA-B*51, encoded by HLA-B in MHC, is the strongest known genetic risk factor for Behçet disease (BD). Associations between BD and other factors within the MHC have been reported also, although strong regional linkage disequilibrium complicates their confident disentanglement from HLA-B*51. In the current study, we examined a combination of directly obtained and imputed MHC-region SNPs, directly obtained HLA-B locus types, and imputed classical HLA types with their corresponding polymorphic amino acid residues for association with BD in 1,190 cases and 1,257 controls. SNP mapping with logistic regression of the MHC identified the HLA-B/MICA region and the region between HLA-F and HLA-A as independently associated with BD (P < 1.7 × 10(-8)). HLA-B*51, -A*03, -B*15, -B*27, -B*49, -B*57, and -A*26 each contributed independently to BD risk. We directly examined rs116799036, a noncoding SNP upstream of HLA-B that was recently suggested to underlie the association of HLA-B*51 with BD, but we were unable to replicate that finding in our collection. Instead, we mapped the BD association to seven MHC class I (MHC-I) amino acid residues, including anchor residues that critically define the selection and binding of peptides to MHC-I molecules, residues known to influence MHC-I-killer immunoglobulin-like receptor interactions, and a residue located in the signal peptide of HLA-B. The locations of these variants collectively implicate MHC-I peptide binding in the pathophysiology of BD. Furthermore, several lines of evidence suggest a role for altered regulation of cellular cytotoxicity in BD pathogenesis.