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

Distributions of the HLA-A, HLA-B, HLA-C, HLA-DRB1, and HLA-DQB1 alleles and haplotype frequencies of 1763 stem cell donors in the Colombian Bone Marrow Registry typed by next-generation sequencing

The HLA compatibility continues to be the main limitation when finding compatible donors, especially if an identical match is not found within the patient's family group. The creation of bone marrow registries allowed a therapeutic option by identifying 10/10 compatible unrelated donors (URD). However, the availability and frequency of haplotypes and HLA alleles are different among ethnic groups and geographical areas, increasing the difficulty of finding identical matches in international registries. In this study, the HLA-A, -B, -C, -DRB1, and -DQB1 loci of 1763 donors registered in the Colombian Bone Marrow Registry were typed by next-generation sequencing. A total of 52 HLA-A, 111 HLA-B, 41 HLA-C, 47 HLA-DRB1, and 20 HLA-DQB1 alleles were identified. The 3 most frequent alleles for each loci were A*24:02g (20,8%), A*02:01g (16,1%), A*01:01g (7.06%); B*35:43g (7.69%), B*40:02g (7.18%), B*44:03g (6.07%); C*04:01g (15.40%), C*01:02g (10.49%), C*07:02g (10.44%); DRB1*04:07g (11.03%), DRB1*07:01g (9.78%), DRB1*08:02g (6.72%); DQB1*03:02g (20.96%), DQB1*03:01g (17.78%) and DQB1*02:01g (16.05%). A total of 497 HLA-A-C-B-DRB1-DQB1 haplotypes were observed with a frequency greater than or equal to 0.05% (> 0.05%); the haplotypes with the highest frequency were A*24:02g~B*35:43g~C*01:02g~DQB1*03:02g~DRB1*04:07g (3.34%), A*29:02g~B*44:03g~C*16:01g~DQB1*02:01g~DRB1*07:01g (2.04%), and A*01:01g~B*08:01g~C*07:01g~DQB1*02:01g~DRB1*03:01g (1.83%). This data will allow the new Colombian Bone Marrow Donor Registry to assess the genetic heterogeneity of the Colombian population and serve as a tool of interest for future searches of unrelated donors in the country.

HLA Diversity in Saudi Population: High Frequency of Homozygous HLA Alleles and Haplotypes

Human leukocyte antigens (HLA) diversity has a tremendous impact on shaping the transplantation practices, transfusion-associated graft versus host disease prevention strategies, and host-pathogen interactions. Here, we conducted a retrospective study of HLA class I and class II homozygosity at allelic and haplotype levels in unrelated individuals genotyped from 2012 to 2016 in a tertiary hospital in the capital of Saudi Arabia. Among 5,000 individuals, 2,773 individuals meet inclusion criteria and were retrospectively analyzed for HLA-A, -B, -C-DRB1, and -DQB1 homozygosity at allelic and haplotype levels. HLA molecular typing was performed using a commercial reverse sequence-specific oligonucleotide (rSSO) kit. We were able to identify 15 HLA-A, 20 HLA-B, 11 HLA-C, 13 HLA-DRB1, and five HLA-DQB1 homozygous alleles demonstrating a very low genetic diversity in the Saudi population. The highest homozygosity in HLA class I was found in locus C followed by A and B (20.3% > 16.1% > 15.5%; p < 0.001) where the most homozygote alleles were A*02 (9.2%), B*51 and B*50 (5.7% and 3.7%), and C*07, C*06, and C*15 (7.2%, 5.48%, and 3.3%) and in HLA class II, the highest homozygosity was found in locus DQB1 compared to DRB1 (31.71% > 19.2%; p < 0.001), with the most common homozygote alleles being DRB1*07 and DRB1*04 (5.33% and 4.2%) and DQB1*02, DQB1*06, and DQB1*03 (13.55%, 7.92%, and 7.64%). The frequency of finding an individual with one homozygote allele was (24.6%), two homozygote alleles (13.5%), three homozygote alleles (4.7%), four homozygote alleles (3.4%), and five alleles were (4.8%). The most frequent homozygote haplotypes are A*23∼C*06∼B*50∼DRB1*07∼DQB1*02 and A*02∼C*06∼B*50∼DRB1*07∼DQB1*02. This study shows low diversity of both class I and II alleles and haplotypes in the Saudi population, which would have a significant impact on shaping the transplantation practices, transfusion-associated graft versus host disease prevention strategies, and host-pathogen interactions.

NKG2A/CD94 Is a New Immune Receptor for HLA-G and Distinguishes Amino Acid Differences in the HLA-G Heavy Chain

Natural killer (NK) cell therapies are a tool to antagonize a dysfunctional immune system. NK cells recognize malignant cells, traffic to a tumor location, and infiltrate the solid tumor. The immune checkpoint molecule human leukocyte antigen (HLA)-G is upregulated on malignant cells but not on healthy surrounding cells, the requirement of understanding the basis of receptor mediated events at the HLA-G/NK cell interface becomes obvious. The NK cell receptors ILT2 and KIR2DL4 have been described to bind to HLA-G; however, their differential function and expression levels on NK cell subsets suggest the existence of an unreported receptor. Here, we performed a ligand-based receptor capture on living cells utilizing sHLA-G*01:01 molecules coupled to TriCEPS and bound to NK cells followed by mass spectrometric analyses. We could define NKG2A/CD94 as a cognate receptor of HLA-G. To verify the results, we used the reciprocal method by expressing recombinant soluble heterodimeric NKG2A/CD94 molecules and used them to target HLA-G*01:01 expressing cells. NKG2A/CD94 could be confirmed as an immune receptor of HLA-G*01:01. Despite HLA-G is marginal polymorphic, we could previously demonstrate that the most common allelic subtypes HLA-G*01:01/01:03 and 01:04 differ in peptide repertoire, their engagement to NK cells, their catalyzation of dNK cell proliferation and their impact on NK cell development. Continuing these studies with regard to NKG2A/CD94 engagement we engineered recombinant single antigen presenting K562 cells and targeted the surface expressed HLA-G*01:01, 01:03 or 01:04 molecules with NKG2A/CD94. Specificity and sensitivity of HLA-G*01:04/NKG2A/CD94 engagement could be significantly verified. The binding affinity decreases when using K562-G*01:03 or K562-G*01:01 cells as targets. These results demonstrate that the ligand-receptor assignment between HLA-G and NKG2A/CD94 is dependent of the amino acid composition in the HLA-G heavy chain. Understanding the biophysical basis of receptor-mediated events that lead to NK cell inhibition would help to remove non-tumor reactive cells and support personalized mild autologous NK cell therapies.

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.

Diagnostic and prognostic biomarkers of Human Leukocyte Antigen complex for hepatitis B virus-related hepatocellular carcinoma

Background: Hepatitis B virus infection had been identified its relationship with liver diseases, including liver tumors. We aimed to explore diagnostic and prognostic values between the Human Leukocyte Antigen (HLA) complex and hepatocellular carcinoma (HCC).

Methods: We used the GSE14520 dataset to explore diagnostic and prognostic significance between HLA complex and HCC. A nomogram was constructed to predict survival probability of HCC prognosis. Gene set enrichment analysis was explored using gene ontologies and metabolic pathways. Validation of prognostic values of the HLA complex was performed in the Kaplan-Meier Plotter website.

Results: We found that HLA-C showed the diagnostic value (P <0.0001, area under curve: 0.784, sensitivity: 93.14%, specificity: 62.26%). In addition, HLA-DQA1 and HLA-F showed prognostic values for overall survival, and HLA-A, HLA-C, HLA-DPA1 and HLA-DQA1 showed prognostic values for recurrence-free survival (all P ≤ 0.05, elevated 0.927, 0.992, 1.023, 0.918, 0.937 multiples compared to non-tumor tissues, respectively). Gene set enrichment analysis found that they were involved in antigen processing and toll like receptor signalling pathway, etc. The nomogram was evaluated for survival probability of HCC prognosis. Validation analysis indicated that HLA-C, HLA-DPA1, HLA-E, HLA-F and HLA-G were associated with HCC prognosis of overall survival (all P ≤ 0.05, elevated 0.988 and 0.997 multiples compared to non-tumor tissues, respectively).

Conclusion: HLA-C might be a diagnostic and prognostic biomarker for HCC. HLA-DPA1 and HLA-F might be prognostic biomarkers for HCC.

High-resolution analysis of the HLA-A, -B, -C and -DRB1 alleles and national and regional haplotype frequencies based on 120 926 volunteers from the Italian Bone Marrow Donor Registry

HLA genes are highly polymorphic and structurally complex. They are located in the major histocompatibility complex (MHC) on chromosome 6, and the frequency of alleles and haplotypes varies widely among human populations. In this paper, we calculated the allele and haplotype frequencies using the HLA data of more than 120 000 Italian unrelated bone marrow donors enrolled in the national registry (IBMDR) and typed them with a high-resolution (HR) method for the HLA-A, -B, -C and -DRB1 alleles. The allele frequency data were obtained by manual counting; haplotype frequencies were calculated using the expectation maximisation (EM) algorithm. The total numbers of observed alleles were 226 for HLA-A, 343 for HLA-B, 201 for HLA-C and 210 for HLA-DRB1, which account for 5.4%, 6.7%, 5.2% and 8.5%, respectively, of each locus allele (IPD-IMGT/HLA Database Release 3.32, April 2018). The three most frequent Italian haplotypes were HLA-A*01:01~B*08:01~C*07:01~DRB1*03:01 (2.5%), A*02:01~B*18:01~C*07:01~DRB1*11:04 (1.1%) and A*30:01~B*13:02~C*06:02~DRB1*07:01 (1.1%). Moreover, for a relevant subset of the examined population (>100 000 individuals), the birthplace was available, and thus, we grouped the frequency data based on the corresponding Italian geographic areas, describing the HLA specificity of the Italian regional populations. The haplotype frequencies were also compared between national and regional data, and we observed remarkable differences in the regional haplotype frequencies, particularly in Sardinia. This study represents a valid tool to identify a more efficient haematopoietic stem cell unrelated donor recruitment and selection strategy, as well as for population genetic and HLA-disease association fields.

Releasing the concept of HLA-allele specific peptide anchors in viral infections: A non-canonical naturally presented human cytomegalovirus-derived HLA-A*24:02 restricted peptide drives exquisite immunogenicity

T‐cell receptors possess the unique ability to survey and respond to their permanently modified ligands, self HLA‐I molecules bound to non‐self peptides of various origin. This highly specific immune function is impaired following hematopoietic stem cell transplantation (HSCT) for a timespan of several months needed for the maturation of T‐cells. Especially, the progression of HCMV disease in immunocompromised patients induces life‐threatening situations. Therefore, the need for a new immune system that delivers vital and potent CD8+ T‐cells carrying TCRs that recognize even one human cytomegalovirus (HCMV) peptide/HLA molecule and clear the viral infection long term becomes obvious. The transcription and translation of HCMV proteins in the lytic cycle is a precisely regulated cascade of processes, therefore, it is a highly sensitive challenge to adjust the exact time point of HCMV‐peptide recruitment over self‐peptides. We utilized soluble HLA technology in HCMV‐infected fibroblasts and sequenced naturally sHLA‐A*24:02 presented HCMV‐derived peptides. One peptide of 14 AAs length derived from the IE2 antigen induced the strongest T‐cell responses; this peptide can be detected with a low ranking score in general peptide prediction databanks. These results highlight the need for elaborate and HLA‐allele specific peptide selection.

HLA-F*01:01 presents peptides with N-terminal flexibility and a preferred length of 16 residues

HLA-F belongs to the non-classical HLA-Ib molecules with a marginal polymorphic nature and tissue-restricted distribution. HLA-F is a ligand of the NK cell receptor KIR3DS1, whose activation initiates an antiviral downstream immune response and lead to delayed disease progression of HIV-1. During the time course of HIV infection, the expression of HLA-F is upregulated while its interaction with KIR3DS1 is diminished. Understanding HLA-F peptide selection and presentation is essential to a comprehensive understanding of this dynamic immune response and the molecules function. In this study, we were able to recover stable pHLA-F*01:01 complexes and analyze the characteristics of peptides naturally presented by HLA-F. These HLA-F-restricted peptides exhibit a non-canonical length without a defined N-terminal anchor. The peptide characteristics lead to a unique presentation profile and influence the stability of the protein. Furthermore, we demonstrate that almost all source proteins of HLA-F-restricted peptides are described to interact with HIV proteins. Understanding the balance switch between HLA-Ia and HLA-F expression and peptide selection will support to understand the role of HLA-F in viral pathogenesis.

Human Leukocyte Antigen-A, B, C, DRB1, and DQB1 Allele and Haplotype Frequencies in a Subset of 237 Donors in the South African Bone Marrow Registry

Human leukocyte antigen- (HLA-) A, HLA-B, HLA-C, HLA-DRB1, and HLA-DQB1 allele and haplotype frequencies were studied in a subset of 237 volunteer bone marrow donors registered at the South African Bone Marrow Registry (SABMR). Hapl-o-Mat software was used to compute allele and haplotype frequencies from individuals typed at various resolutions, with some alleles in multiple allele code (MAC) format. Four hundred and thirty-eight HLA-A, 235 HLA-B, 234 HLA-DRB1, 41 HLA-DQB1, and 29 HLA-C alleles are reported. The most frequent alleles were A∗02:02g (0.096), B∗07:02g (0.082), C∗07:02g (0.180), DQB1∗06:02 (0.157), and DRB1∗15:01 (0.072). The most common haplotype was A∗03:01g~B∗07:02g~C∗07:02g~DQB1∗06:02~DRB1∗15:01 (0.067), which has also been reported in other populations. Deviations from Hardy-Weinberg equilibrium were observed in A, B, and DRB1 loci, with C~DQB1 being the only locus pair in linkage disequilibrium. This study describes allele and haplotype frequencies from a subset of donors registered at SABMR, the only active bone marrow donor registry in Africa. Although the sample size was small, our results form a key resource for future population studies, disease association studies, and donor recruitment strategies.

HLA-G peptide preferences change in transformed cells: impact on the binding motif

HLA-G is known for its strictly restricted tissue distribution. HLA-G expression could be detected in immune privileged organs and many tumor entities such as leukemia, multiple myeloma, and non-Hodgkin and Hodgkin’s lymphoma. This functional variability from mediation of immune tolerance to facilitation of tumor immune evasion strategies might translate to a differential NK cell inhibition between immune-privileged organs and tumor cells. The biophysical invariability of the HLA-G heavy chain and its contrary diversity in immunity implicates a strong influence of the bound peptides on the pHLA-G structure. The aim was to determine if HLA-G displays a tissue-specific peptide repertoire. Therefore, using soluble sHLA-G technology, we analyzed the K562 and HDLM-2 peptide repertoires. Although both cell lines possess a comparable proteome and recruit HLA-G-restricted peptides through the same peptide-loading pathway, the peptide features appear to be cell specific. HDLM-2 derived HLA-G peptides are anchored by an Arg at p1 and K562-derived peptides are anchored by a Lys. At p2, no anchor motif could be determined while peptides were anchored at pΩ with a Leu and showed an auxiliary anchor motif Pro at p3. To appreciate if the peptide anchor alterations are due to a cell-specific differential peptidome, we performed analysis of peptide availability within the different cell types. Yet, the comparison of the cell-specific proteome and HLA-G-restricted ligandome clearly demonstrates a tissue-specific peptide selection by HLA-G molecules. This exclusive and unexpected observation suggests an exquisite immune function of HLA-G.

The molecular basis for peptide repertoire selection in the human leucocyte antigen (HLA) C*06:02 molecule

Human leukocyte antigen (HLA)-C*06:02 is identified as the allele associated with the highest risk for the development of the autoimmune skin disease psoriasis. However, the diversity and mode of peptide presentation by the HLA-C*06:02 molecule remains unclear. Here, we describe the endogenous peptide repertoire of ∼3,000 sequences for HLA-C*06:02 that defines the peptide-binding motif for this HLA allomorph. We found that HLA-C*06:02 predominantly presents nonamer peptides with dominant arginine anchors at the P2 and P7 positions and a preference for small hydrophobic residues at the C terminus (PΩ). To determine the structural basis of this selectivity, we determined crystal structures of HLA-C*06:02 in complex with two self-peptides (ARTELYRSL and ARFNDLRFV) and an analogue of a melanocyte autoantigen (ADAMTSL5, VRSRR-abu-LRL) implicated in psoriasis. These structures revealed that HLA-C*06:02 possesses a deep peptide-binding groove comprising two electronegative B- and E-pockets that coincide with the preference for P2 and P7 arginine anchors. The ADAMTSL5 autoantigen possessed a P7-Leu instead of the P7-Arg residue, but nevertheless was accommodated within the HLA-C*06:02 antigen-binding cleft. Collectively, our results provide the structural basis for understanding peptide repertoire selection in HLA-C*06:02.

Understanding the obstacle of incompatibility at residue 156 within HLA-B*35 subtypes

Defining permissive and non-permissive mismatches for transplantation is a demanding challenge. Single mismatches at amino acid (AA) position 156 of human leucocyte antigen (HLA) class I have been described to alter the peptide motif, repertoire, or mode of peptide loading through differential interaction with the peptide-loading complex. Hence, a single mismatch can tip the balance and trigger an immunological reaction. HLA-B*35 subtypes have been described to evade the loading complex, 156 mismatch distinguishing B*35:01 and B*35:08 changes the binding groove sufficiently to alter the sequence features of the selected peptide repertoire. To understand the functional influences of residue 156 in B*35 variants, we analyzed the peptide binding profiles of HLA-B*35:01^156Leu, B*35:08^156Arg and B*35:62^156Trp. The glycoprotein tapasin represents a target for immune evasions and functions within the multimeric peptide-loading complex to stabilize empty class I molecules and promote acquisition of high-affinity peptides. All three B*35 subtypes showed a tapasin-independent mode of peptide acquisition. HLA-B*35-restricted peptides of low- and high-binding affinities were recovered in the presence and absence of tapasin and subsequently sequenced utilizing mass spectrometry. The peptides derived from B*35 variants differ substantially in their features dependent on their mode of recruitment; all peptides were preferentially anchored by Pro at p2 and Tyr, Phe, Leu, or Lys at pΩ. However, the Trp at residue 156 altered the p2 motif to an Ala and restricted the pΩ to a Trp. Our results highlight the importance of understanding the impact of key micropolymorphism and how a single AA mismatch orchestrates the neighboring AAs.

The diversity of the HLA-E-restricted peptide repertoire explains the immunological impact of the Arg107Gly mismatch

Human leukocyte antigen (HLA)-E molecules are potent inhibitors of NK cell-mediated killing. Low in polymorphisms, two alleles are widely expressed among diverse populations: HLA-E*01:01 and HLA-E*01:03. Both alleles are distinguished by one SNP resulting in the substitution Arg107Gly. Both alleles present a limited set of peptides derived from class I leader sequences physiologically; however, HLA-E*01:01 presents non-canonical peptides in the absence of HLA class I molecules. To further assess the functional differences between both alleles, we analyzed the peptide repertoire of HLA-E*01:03 by applying soluble HLA technology followed by mass-spectrometric peptide sequencing. HLA-E*01:03 restricted peptides showed a length of 9-17 amino acids and differed in their biophysical properties, no overlap in the peptide repertoire of both allelic variants could be observed; however, both alleles shared marginal peptides from the same proteomic content. Artificial APCs expressing empty HLA-E*01:01 or E*01:03 molecules were generated and stabilized using cognate HLA class I-derived peptide ligands to analyze the impact of residue 107 within the HLA-E heavy chain on the NKG2/CD94 receptor engagement. Differences in peptide stabilization could be translated to the density and half-life time of peptide-HLA-E molecules on the cell surface that subsequently impacted NK cell inhibition as verified by cytotoxicity assays. Taken together, these data illustrate functional differences of HLA-E allelic variants induced by a single amino acid. Furthermore, the function of HLA-E in pathophysiologic situations when the HLA processing machinery is interrupted seems to be more emphasized than previously described, implying a crucial role for HLA-E in tumor or viral immune episodes.

Human Leukocyte Antigen Diversity: A Southern African Perspective

Despite the increasingly well-documented evidence of high genetic, ethnic, and linguistic diversity amongst African populations, there is limited data on human leukocyte antigen (HLA) diversity in these populations. HLA is part of the host defense mechanism mediated through antigen presentation to effector cells of the immune system. With the high disease burden in southern Africa, HLA diversity data is increasingly important in the design of population-specific vaccines and the improvement of transplantation therapeutic interventions. This review highlights the paucity of HLA diversity data amongst southern African populations and defines a need for information of this kind. This information will support disease association studies, provide guidance in vaccine design, and improve transplantation outcomes.

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.