A combined DPA1~DPB1 amino acid epitope is the primary unit of selection on the HLA-DP heterodimer

Profiling the genetic diversity of the HLA system in Mexico using 9-locus allele and haplotype frequencies from donors in the NMDP Mexico donor center

Profiling the HLA diversity at the population level benefits multiple clinical and anthropological applications, such as tracing population migration, identifying genetic relationships between different groups, quantifying the added diversity in a global donor pool and matching for solid organ and stem cell transplantation. We calculated nine-locus HLA-A ∼ C ∼ B ∼ DRB1 ∼ DRB3/4/5 ∼ DQA1 ∼ DQB1 ∼ DPA1 ∼ DPB1 allele and haplotype frequencies in about 170,000 volunteer donor genotypes from the NMDP Mexico (NMDP MX, previously Be The Match Mexico) donor center. These donors are predominantly of Mexican ancestry recruited from multiple regions in Mexico. The goal of the study was to describe the HLA genetic profiles of the Mexican population and investigate the contribution of these donors' HLA in serving Mexican, US and international patients in need of hematopoietic cell transplants. Additionally, we estimated that almost all Mexican patients will have an available 5 of 8 or better matched donor in the NMDP MX donor center with matches also available for some of the Latino patients in the U.S. We demonstrate that Mexican populations clustered genetically and shared multiple frequent alleles and haplotypes with populations from the US Mexican or Chicano, US South/Central American Hispanic, and some Latino populations. Operationally, 78 % of NMDP Mexico donors contributed genotypes that were observed a total of three times or less on the registry, increasing the diversity of the overall NMDP registry. More than 300 donor collections were facilitated through the NMDP MX donor center serving mostly Hispanic/Latino patients in the US and abroad. This study highlights the importance of adding the NMDP MX donors to the worldwide donor pool and paves the way for a data-driven strategy for future planning and donor recruitment.

Unraveling the architecture of major histocompatibility complex class II haplotypes in rhesus macaques

The regions in the genome that encode components of the immune system are often featured by polymorphism, copy number variation, and segmental duplications. There is a need to thoroughly characterize these complex regions to gain insight into the impact of genomic diversity on health and disease. Here we resolve the organization of complete major histocompatibility complex (MHC) class II regions in rhesus macaques by using a long-read sequencing strategy (Oxford Nanopore Technologies) in concert with adaptive sampling. In particular, the expansion and contraction of the primate DRB-region appear to be a dynamic process that involves the rearrangement of different cassettes of paralogous genes. These chromosomal recombination events are propagated by a conserved pseudogene, DRB6, which features the integration of two retroviral elements. In contrast, the DRA locus appears to be protected from rearrangements, which may be owing to the presence of an adjacently located truncated gene segment, DRB9 With our sequencing strategy, the annotation, evolutionary conservation, and potential function of pseudogenes can be reassessed, an aspect that was neglected by most genome studies in primates. Furthermore, our approach facilitates the characterization and refinement of an animal model essential to study human biology and disease.

Natural Selection on HLA-DPB1 Amino Acids Operates Primarily on DP Serologic Categories

The DPB1 locus is notable among the classical HLA loci in that allele frequencies at this locus are consistent with genetic drift, whereas the frequencies of specific DPβ amino acids are consistent with the action of balancing selection. We investigated the influence of natural selection in shaping the diversity of three functional categories of DPB1 diversity defined by specific amino acid motifs, DPB1 T-cell epitopes, DPB1 supertypes and DP1-DP4 serologic categories (SCs), via Ewens-Watterson (EW) selective neutrality and asymmetric Linkage Disequilibrium (ALD) analyses in a worldwide sample of 136 populations. These EW analyses provide strong evidence for the operation of balancing selection on DP SCs, but no evidence for balancing selection on T-cell epitopes or supertypes. We further investigated the global distribution of SCs. Each SC is common in a different region of the world, with the DP1 SC most common in Southeast Asia and Oceania, the DP2 SC in North and South America, the DP3 SC in South America, and the DP4 SC in Europe. The DP2 SC is present in all populations, while 14% of populations are missing at least one DP1, DP3, or DP4 SC. We observed consistent DPA1∼DP SC haplotype associations across 10 populations from five global regions, and found that asymmetric linkage disequilibrium (LD) between the DPB1 locus and the four most-common DPA1 alleles (DPA1*01:03, *02:01, *02:02 and *03:01) is determined by variation at DPβ AA positions 85-87. These positions are in LD with both DPα positions 31 and 50. We conclude from these EW analyses that natural selection is primarily operating to maintain population-level diversity of DP SCs, rather than DPB1 alleles or other functional categories of DPB1 diversity.

Modulation of the pharmacokinetics of soluble ACE2 decoy receptors through glycosylation

The Spike of SARS-CoV-2 recognizes a transmembrane protease, angiotensin-converting enzyme 2 (ACE2), on host cells to initiate infection. Soluble derivatives of ACE2, in which Spike affinity is enhanced and the protein is fused to Fc of an immunoglobulin, are potent decoy receptors that reduce disease in animal models of COVID-19. Mutations were introduced into an ACE2 decoy receptor, including adding custom N-glycosylation sites and a cavity-filling substitution together with Fc modifications, which increased the decoy's catalytic activity and provided small to moderate enhancements of pharmacokinetics following intravenous and subcutaneous administration in humanized FcRn mice. Most prominently, sialylation of native glycans increases exposures by orders of magnitude, and the optimized decoy is therapeutically efficacious in a mouse COVID-19 model. Ultimately, an engineered and highly sialylated decoy receptor produced using methods suitable for manufacture of representative drug substance has high exposure with a 5- to 9-day half-life. Finally, peptide epitopes at mutated sites in the decoys generally have low binding to common HLA class II alleles and the predicted immunogenicity risk is low. Overall, glycosylation is a critical molecular attribute of ACE2 decoy receptors and modifications that combine tighter blocking of Spike with enhanced pharmacokinetics elevate this class of molecules as viable drug candidates.

HLA-DPB1 and DPA1 ~ DPB1 linkage mismatch affects the survival of recipients receiving HLA-14/14 matched unrelated donor HSCT

To analyse the effect of HLA-DPA1 and HLA-DPB1 allelic mismatches on the outcomes of unrelated donor haematopoietic stem cell transplantation (URD-HSCT), we collected 258 recipients with haematological disease who underwent HLA-10/10 matched URD-HSCT. HLA-A, -B, -C, -DRB1, -DQB1, -DRB3/4/5, -DQA1, -DPA1 and -DPB1 typing was performed for the donors and recipients using next-generation sequencing (NGS) technology. After excluding 8 cases with DQA1 or DRB3/4/5 mismatches, we included 250 cases with HLA-14/14 matching for further analysis. Our results showed that the proportion of matched DPA1 and DPB1 alleles was only 10.4% (26/250). The remaining 89.6% of donors and recipients demonstrated DPA1 or DPB1 mismatch. In the DPA1 matched and DPB1 mismatched group, accounting for 18.8% (47/250) of the cohort, DPB1*02:01/DPB1*03:01 allelic mismatches were associated with decreased 2-year OS and increased NRM. DPB1*02:02/DPB1*05:01 and DPB1*02:01/DPB1*05:01 mismatches showed no impact on outcomes. Moreover, the specific allelic mismatches observed were consistent with the DPB1 T-cell epitope (TCE) classification as permissive and non-permissive. We innovatively established an analysis method for DPA1 ~ DPB1 linkage mismatch for cases with both DPA1 and DPB1 mismatched, accounting for 70% (175/250) of the total. DPA1*02:02 ~ DPB1*05:01/DPA1*02:01 ~ DPB1*17:01 linkage mismatches were associated with lower 2-year OS, especially among AML/MDS recipients. DPA1*02:02 ~ DPB1*05:01/DPA1*01:03 ~ DPB1*02:01 linkage mismatches showed no impact on outcomes. In conclusion, applying the DPA1 ~ DPB1 linkage mismatch analysis approach can identify different types of mismatches affecting transplant outcomes and provide valuable insight for selecting optimal donors for AML/MDS and ALL recipients.

HLA-DPA1*02:01~B1*01:01 is a risk haplotype for primary sclerosing cholangitis mediating activation of NKp44+ NK cells

OBJECTIVE

Primary sclerosing cholangitis (PSC) is characterised by bile duct strictures and progressive liver disease, eventually requiring liver transplantation. Although the pathogenesis of PSC remains incompletely understood, strong associations with HLA-class II haplotypes have been described. As specific HLA-DP molecules can bind the activating NK-cell receptor NKp44, we investigated the role of HLA-DP/NKp44-interactions in PSC.

DESIGN

Liver tissue, intrahepatic and peripheral blood lymphocytes of individuals with PSC and control individuals were characterised using flow cytometry, immunohistochemical and immunofluorescence analyses. HLA-DPA1 and HLA-DPB1 imputation and association analyses were performed in 3408 individuals with PSC and 34 213 controls. NK cell activation on NKp44/HLA-DP interactions was assessed in vitro using plate-bound HLA-DP molecules and HLA-DPB wildtype versus knock-out human cholangiocyte organoids.

RESULTS

NKp44+NK cells were enriched in livers, and intrahepatic bile ducts of individuals with PSC showed higher expression of HLA-DP. HLA-DP haplotype analysis revealed a highly elevated PSC risk for HLA-DPA1*02:01~B1*01:01 (OR 1.99, p=6.7×10-50). Primary NKp44+NK cells exhibited significantly higher degranulation in response to plate-bound HLA-DPA1*02:01-DPB1*01:01 compared with control HLA-DP molecules, which were inhibited by anti-NKp44-blocking. Human cholangiocyte organoids expressing HLA-DPA1*02:01-DPB1*01:01 after IFN-γ-exposure demonstrated significantly increased binding to NKp44-Fc constructs compared with unstimulated controls. Importantly, HLA-DPA1*02:01-DPB1*01:01-expressing organoids increased degranulation of NKp44+NK cells compared with HLA-DPB1-KO organoids.

CONCLUSION

Our studies identify a novel PSC risk haplotype HLA-DP A1*02:01~DPB1*01:01 and provide clinical and functional data implicating NKp44+NK cells that recognise HLA-DPA1*02:01-DPB1*01:01 expressed on cholangiocytes in PSC pathogenesis.

HLA-DPB1*13:01 associates with enhanced, and KIR2DS4*001 with diminished protection from developing severe COVID-19

Extreme polymorphism of HLA and killer-cell immunoglobulin-like receptors (KIR) differentiates immune responses across individuals. Additional to T cell receptor interactions, subsets of HLA class I act as ligands for inhibitory and activating KIR, allowing natural killer (NK) cells to detect and kill infected cells. We investigated the impact of HLA and KIR polymorphism on the severity of COVID-19. High resolution HLA class I and II and KIR genotypes were determined from 403 non-hospitalized and 1575 hospitalized SARS-CoV-2 infected patients from Italy collected in 2020. We observed that possession of the activating KIR2DS4*001 allotype is associated with severe disease, requiring hospitalization (OR = 1.48, 95% CI 1.20-1.85, pc  = 0.017), and this effect is greater in individuals homozygous for KIR2DS4*001 (OR = 3.74, 95% CI 1.75-9.29, pc  = 0.003). We also observed the HLA class II allotype, HLA-DPB1*13:01 protects SARS-CoV-2 infected patients from severe disease (OR = 0.49, 95% CI 0.33-0.74, pc  = 0.019). These association analyses were replicated using logistic regression with sex and age as covariates. Autoantibodies against IFN-α associated with COVID-19 severity were detected in 26% of 156 hospitalized patients tested. HLA-C*08:02 was more frequent in patients with IFN-α autoantibodies than those without, and KIR3DL1*01502 was only present in patients lacking IFN-α antibodies. These findings suggest that KIR and HLA polymorphism is integral in determining the clinical outcome following SARS-CoV-2 infection, by influencing the course both of innate and adaptive immunity.

Report from the extended HLA-DPA1 ~ promoter ~ HLA-DPB1 haplotype of the 18th international HLA and immunogenetics workshop

The primary goal of the HLA-DPA1 ~ promoter ~ HLA-DPB1 haplotype component of the 18th IHIWS was to characterise the extended haplotypes within the HLA-DP region and survey the extent of genetic diversity in this region across human populations. In this report, we analysed single-nucleotide polymorphisms (SNPs) in 255 subjects from 6 different cohorts. The results from the HLA-DP haplotype component have validated findings from the initial pilot study. SNPs in this region were inherited in strong linkage, particularly HLA-DPA1, SNP-linked promoter haplotypes and motifs in exon 2 of HLA-DPB1. We reported 17 SNP-linked haplotypes in the promoter region. Together with HLA-DPA1 and HLA-DPB1 alleles, they formed 74 distinct extended HLA-DP haplotypes in 438 sequences. We also observed the presence of region-specific alleles and promoter haplotypes. Our approach involved phasing extended SNPs including promoter SNPs, HLA-DPA1 and HLA-DPB1 alleles, in a 22 kb region, GRCh38/hg38 (chr6:33,064,111-33,086,679), followed by clustering of these SNPs as one extended haplotype. This hierarchical clustering revealed four major clades, suggesting that haplotypes within each clade may have diverged from a common ancestral haplotype and undergone similar evolutionary processes. The correlation between HLA-DPA1 and the promoter region raises questions about the role of HLA-DPA1 antigen in the heterodimer. This finding requires validation on a larger sample size specifically designed for anthropological analysis. Nevertheless, the results from this study highlight the clinical potential of selecting better-matched donors for patients awaiting haematopoietic stem cell transplants from genetically overlapping groups that share common ancestral haplotypes.

Accurate prediction of HLA class II antigen presentation across all loci using tailored data acquisition and refined machine learning

Accurate prediction of antigen presentation by human leukocyte antigen (HLA) class II molecules is crucial for rational development of immunotherapies and vaccines targeting CD4+ T cell activation. So far, most prediction methods for HLA class II antigen presentation have focused on HLA-DR because of limited availability of immunopeptidomics data for HLA-DQ and HLA-DP while not taking into account alternative peptide binding modes. We present an update to the NetMHCIIpan prediction method, which closes the performance gap between all three HLA class II loci. We accomplish this by first integrating large immunopeptidomics datasets describing the HLA class II specificity space across all loci using a refined machine learning framework that accommodates inverted peptide binders. Next, we apply targeted immunopeptidomics assays to generate data that covers additional HLA-DP specificities. The final method, NetMHCIIpan-4.3, achieves high accuracy and molecular coverage across all HLA class II allotypes.

HLA-DP on Epithelial Cells Enables Tissue Damage by NKp44+ Natural Killer Cells in Ulcerative Colitis

BACKGROUND & AIMS

Ulcerative colitis (UC) is characterized by severe inflammation and destruction of the intestinal epithelium, and is associated with specific risk single nucleotide polymorphisms in HLA class II. Given the recently discovered interactions between subsets of HLA-DP molecules and the activating natural killer (NK) cell receptor NKp44, genetic associations of UC and HLA-DP haplotypes and their functional implications were investigated.

METHODS

HLA-DP haplotype and UC risk association analyses were performed (UC: n = 13,927; control: n = 26,764). Expression levels of HLA-DP on intestinal epithelial cells (IECs) in individuals with and without UC were quantified. Human intestinal 3-dimensional (3D) organoid cocultures with human NK cells were used to determine functional consequences of interactions between HLA-DP and NKp44.

RESULTS

These studies identified HLA-DPA1∗01:03-DPB1∗04:01 (HLA-DP401) as a risk haplotype and HLA-DPA1∗01:03-DPB1∗03:01 (HLA-DP301) as a protective haplotype for UC in European populations. HLA-DP expression was significantly higher on IECs of individuals with UC compared with controls. IECs in human intestinal 3D organoids derived from HLA-DP401pos individuals showed significantly stronger binding of NKp44 compared with HLA-DP301pos IECs. HLA-DP401pos IECs in organoids triggered increased degranulation and tumor necrosis factor production by NKp44+ NK cells in cocultures, resulting in enhanced epithelial cell death compared with HLA-DP301pos organoids. Blocking of HLA-DP401-NKp44 interactions (anti-NKp44) abrogated NK cell activity in cocultures.

CONCLUSIONS

We identified an UC risk HLA-DP haplotype that engages NKp44 and activates NKp44+ NK cells, mediating damage to intestinal epithelial cells in an HLA-DP haplotype-dependent manner. The molecular interaction between NKp44 and HLA-DP401 in UC can be targeted by therapeutic interventions to reduce NKp44+ NK cell-mediated destruction of the intestinal epithelium in UC.

Unrelated Stem Cell Donor HLA Match Likelihood in the US Registry Incorporating HLA-DPB1 Permissive Mismatching

Donor-recipient HLA matching at the DPB1 locus improves the outcomes of hematopoietic stem cell transplantation (HCT). Retrospective outcome studies found that in HCTs matched for all 8 alleles of the A, B, C, and DRB1 loci at high resolution (8/8 match), few transplantations were also allele-matched at the DPB1 locus. DPB1 allele matching was once thought to be logistically impractical; however, a DPB1-permissive mismatch model based on T cell epitope (TCE) reactivity expands the proportion of suitable donors. To understand the likelihood of finding a DPB1-permissive donor, we sought to expand population genetic match likelihood models for the US unrelated donor registry, the National Marrow Donor Program (NMDP). After extending HLA haplotype frequency estimates to include the DPB1 locus, our models found that the likelihood of having a DPB1-permissive donor was not much lower than likelihood of 8/8 matching. A maximum of 5 additional donors would need to be typed to find a more optimal DPB1-permissive donor at least 90% of the time. Linkage disequilibrium patterns between the DPB1 locus and other classical HLA loci varied markedly by haplotype and population, indicating that the known recombination hotspot between DQ and DP gene complexes has not had a uniform impact; thus, DPB1-permissive donors are easier to identify within minority populations. DPB1 TCE categories were highly predictable from HLA typing at other loci when imputed with extended haplotype frequency data. Our overall results indicate that registry search strategies that seek a more optimally matched HCT donor encompassing HLA-DPB1 permissibility are likely to be highly productive.

A universal MHCII technology platform to characterize antigen-specific CD4+ T cells

CD4+ T cells are critical to the immune system and perform multiple functions; therefore, their identification and characterization are crucial to better understanding the immune system in both health and disease states. However, current methods rarely preserve their ex vivo phenotype, thus limiting our understanding of their in vivo functions. Here we introduce a flexible, rapid, and robust platform for ex vivo CD4+ T cell identification. By combining MHCII allele purification, allele-independent peptide loading, and multiplexed flow cytometry technologies, we can enable high-throughput personalized CD4+ T cell identification, immunophenotyping, and sorting. Using this platform in combination with single-cell sorting and multimodal analyses, we identified and characterized antigen-specific CD4+ T cells relevant to COVID-19 and cancer neoantigen immunotherapy. Overall, our platform can be used to detect and characterize CD4+ T cells across multiple diseases, with potential to guide CD4+ T cell epitope design for any disease-specific immunization strategy.

Human T cells recognize HLA-DP-bound peptides in two orientations

Human leukocyte antigen (HLA) molecules present small peptide antigens to T cells, thereby allowing them to recognize pathogen-infected and cancer cells. A central dogma over the last 50+ y is that peptide binding to HLA molecules is mediated by the docking of side chains of particular amino acids in the peptide into pockets in the HLA molecules in a conserved N- to C-terminal orientation. Whether peptides can be presented in a reversed C- to N-terminal orientation remains unclear. Here, we performed large-scale identification of peptides bound to HLA-DP molecules and observed that in addition to peptide binding in an N- to C-terminal orientation, in 9 out of 14 HLA-DP allotypes, reverse motifs are found, compatible with C- to N-terminal peptide binding. Moreover, we isolated high-avidity human cytomegalovirus (CMV)-specific HLA-DP-restricted CD4⁺ T cells from the memory repertoire of healthy donors and demonstrate that such T cells recognized CMV-derived peptides bound to HLA-DPB1*01:01 or *05:01 in a reverse C- to N-terminal manner. Finally, we obtained a high-resolution HLA-DPB1*01:01-CMVpp65_(142-158) peptide crystal structure, which is the molecular basis for C- to N-terminal peptide binding to HLA-DP. Our results point to unique features of HLA-DP molecules that substantially broaden the HLA class II bound peptide repertoire to combat pathogens and eliminate cancer cells.

The changing landscape of HLA typing: Understanding how and when HLA typing data can be used with confidence from bench to bedside

Human leukocyte antigen (HLA) genes are extraordinary for their extreme diversity and widespread impact on human health and disease. More than 30,000 HLA alleles have been officially named and more alleles continue to be discovered at a rapid pace. HLA typing systems which have been developed to detect HLA diversity have advanced rapidly and are revolutionizing our understanding of HLA's clinical importance. However, continuous improvements in knowledge and technology have created challenges for clinicians and scientists. This review explains how differences in HLA typing systems can impact the HLA types that are assigned. The consequences of differences in laboratory testing methods and reference databases are described. The challenges of using HLA types that are not equivalent are illustrated. A fundamental understanding of the continual expansion of our understanding of HLA diversity and limitations in some of the typing data is essential for using typing data appropriately in clinical and research settings.

HLA-DPB1 allele frequencies in the West Midlands region of the United Kingdom: A critical evaluation against the common, intermediate and well-documented allele catalogues CWD 2.0.0, EFI CWD and CIWD 3.0.0

The Birmingham H&I laboratory performed HLA typing on 456 potential deceased solid organ donors in the UK between 2014 and 2016. Accurate DPB1 typing is essential for determining HLA compatibility in transplantation, thus we report HLA-DPB1 for potential deceased solid organ donors. To correctly interpret HLA typing data, laboratories must understand both international and local HLA allele frequencies. In this analysis, we determined HLA-DPB1 allele and genotype frequencies for these 456 donors. HLA-DPB1 diversity was evaluated against the common and well-documented (CWD) alleles 2.0.0 catalogue, the European Federation for Immunogenetics (EFI) CWD catalogue and the common, intermediate and well-documented (CIWD) 3.0.0 catalogue. Additionally, we determined which alleles are common in our local deceased donor population. We observed 27 HLA-DPB1 alleles with DPB1*04:01 being the most frequently observed (allele frequency = 0.4342). All alleles detected locally were present in CIWD 3.0.0, however, DPB1*124:01 and *135:01 were not present in CWD 2.0.0 and DPB1*104:01 and *135:01 were not present in EFI CWD. Twenty of 27 DPB1 alleles identified were defined as locally common and also listed as common in CIWD 3.0.0 representing 62.5% of common alleles in the subset of CIWD 3.0.0 from individuals of a European geographic, ancestral or ethnic background. The alleles HLA-DPB1*16:01 and *20:01 are locally common but not listed as common in EFI CWD and DPB1*104:01 is not listed as common in CWD 2.0.0 catalogue. Our analysis showed that the detected alleles and locally common alleles within our population were aligned with the CIWD 3.0.0 catalogue.

The combinatorial diversity of KIR and HLA class I allotypes in Peninsular Malaysia

Killer cell immunoglobulin-like receptors (KIRs) interact with polymorphic human leucocyte antigen (HLA) class I molecules, modulating natural killer (NK) cell functions and affecting both the susceptibility and outcome of immune-mediated diseases. The KIR locus is highly diverse in gene content, copy number and allelic polymorphism within individuals and across geographical populations. To analyse currently under-represented Asian and Pacific populations, we investigated the combinatorial diversity of KIR and HLA class I in 92 unrelated Malay and 75 Malaysian Chinese individuals from the Malay Peninsula. We identified substantial allelic and structural diversity of the KIR locus in both populations and characterized novel variations at each analysis level. The Malay population is more diverse than Malay Chinese, likely representing a unique history including admixture with immigrating populations spanning several thousand years. Characterizing the Malay population are KIR haplotypes with large structural variants present in 10% individuals, and KIR and HLA alleles previously identified in Austronesian populations. Despite the differences in ancestries, the proportion of HLA allotypes that serve as KIR ligands is similar in each population. The exception is a significantly reduced frequency of interactions of KIR2DL1 with C2+ HLA-C in the Malaysian Chinese group, caused by the low frequency of C2+ HLA. One likely implication is a greater protection from preeclampsia, a pregnancy disorder associated with KIR2DL1, which shows higher incidence in the Malay than in the Malaysian Chinese. This first complete, high-resolution, characterization of combinatorial diversity of KIR and HLA in Malaysians will form a valuable reference for future clinical and population studies.

Polymorphism clustering of the 21.5 kb DPA-promoter-DPB region reveals novel extended full-length haplotypes

DPB1 and DPA1 genes share the same promoter region. Single-nucleotide polymorphisms (SNPs) within the regulatory regions of DP have been reported. This study hypothesizes that by including the SNPs in the promoter region of DP, extended haplotypes are defined, and promoter polymorphism is more extensive than what is currently reported. To identify the SNPs in the region of interest, the DP region spanning 21.5 kb was amplified in three separate long-ranged polymerase chain reactions. A DNA panel consisting of 100 samples was selected to represent a broad range of DPB1 alleles. The panel was amplified and sequenced using a dual sequencing strategy. Binary alignment map (BAM) alignments were generated and the mapped sequence alignments were analyzed using Integrative Genomics Viewer. A total of 76 SNPs were identified, and SNPs were clustered into 12 SNP-linked haplotypes. Multiple sequence alignments of promoter sequences indicated four distinct lineages within the connective region (CR) between two genes. The relationship between DPA1, CR, DPB1, and amino acid motifs was found to be correlated with HV1 and HV6. Of the 12 promoter haplotypes, DPB1 alleles observed with ProDP-4 were in complete linkage with HV1/2/5/6, the rs9277534G SNP, and the highly immunogenic T-cell epitope group. Multiple extended haplotypes of different intronic subtypes of the same DPB1 alleles were also identified. This new view of the full DP haplotype shows the relation of polymorphism, genes, and alleles, and provides a basis for future functionality related nomenclature. The novel clustering of the DP-extended haplotype warrants future investigations of DP haplotype matching in the outcome of haematopoietic stem cell transplantation (HSCT).

HLA-DPB1 Reactive T Cell Receptors for Adoptive Immunotherapy in Allogeneic Stem Cell Transplantation

HLA-DPB1 antigens are mismatched in about 80% of allogeneic hematopoietic stem cell transplantations from HLA 10/10 matched unrelated donors and were shown to be associated with a decreased risk of leukemia relapse. We recently developed a reliable in vitro method to generate HLA-DPB1 mismatch-reactive CD4 T-cell clones from allogeneic donors. Here, we isolated HLA-DPB1 specific T cell receptors (TCR DP) and used them either as wild-type or genetically optimized receptors to analyze in detail the reactivity of transduced CD4 and CD8 T cells toward primary AML blasts. While both CD4 and CD8 T cells showed strong AML reactivity in vitro, only CD4 T cells were able to effectively eliminate leukemia blasts in AML engrafted NOD/SCID/IL2Rγc^−/− (NSG) mice. Further analysis showed that optimized TCR DP and under some conditions wild-type TCR DP also mediated reactivity to non-hematopoietic cells like fibroblasts or tumor cell lines after HLA-DP upregulation. In conclusion, T cells engineered with selected allo-HLA-DPB1 specific TCRs might be powerful off-the-shelf reagents in allogeneic T-cell therapy of leukemia. However, because of frequent (common) cross-reactivity to non-hematopoietic cells with optimized TCR DP T cells, safety mechanisms are mandatory.

Fine mapping of the MHC region identifies major independent variants associated with Han Chinese primary biliary cholangitis

The genetic association of primary biliary cholangitis with major histocompatibility complex (MHC) has been widely confirmed among different ethnicities. To map specific MHC region variants associated with PBC in a Han Chinese cohort, we imputed HLA antigens and amino acids (AA) in 1126 PBC cases and 1770 healthy control subjects using a Han-MHC reference database. We demonstrate that HLA-DRB1 and/or HLA-DQB1 contributed the strongest signals, and that HLA-DPB1 was a separate independent locus. Regression analyses with classical HLA alleles indicate that HLA-DQB1*03:01 or HLA-DQβ1-Pro55, HLA-DPB1*17:01 or HLA-DPβ1-Asp84 and HLA-DRB1*08:03 could largely explain MHC association with PBC. Forward stepwise regression analyses with HLA amino acid variants localize the major signals to HLA-DRβ1-Ala74, HLA-DQβ1-Pro55 and HLA-DPβ1-Asp84. Electrostatic potential calculations implicated AA variations at HLA-DQβ1 position 55 and HLA-DPβ1 position 84 as critical to peptide binding properties. Furthermore, although several critical Han Chinese AA variants differed from those shown in European populations, the predicted effects on antigen binding are likely to be very similar or identical and underlie the major component of MHC association with PBC.

Epitope-based human leukocyte antigen matching for transplantation: a personal perspective of its future

PURPOSE OF REVIEW

This study reflects my personal experience with the characterization of human leukocyte antigen (HLA) epitopes and their significance in HLA matching for transplantation. It offers a subjective assessment what further studies are needed to have this concept be applied in the clinical setting.

RECENT FINDINGS

This study addresses the structural characteristics of antibody-reactive HLA epitopes determined by different methods, eplet-associated antibody analysis and acceptable mismatching for sensitized patients and eplet immunogenicity and determination of mismatch permissibility. BASIC IMPLICATIONS: for clinical practice and research consider the need for further studies of the structural basis of antibody-verified HLA epitopes determined in different techniques and their clinical relevance, the biological basis of epitope immunogenicity and determinations of permissible mismatches and a computerized clinical transplant database with an Artificial Intelligence component that can generate evidence-based information for the practical application of epitope-based HLA matching.

Amino acid signatures of HLA Class-I and II molecules are strongly associated with SLE susceptibility and autoantibody production in Eastern Asians

Human leukocyte antigen (HLA) is a key genetic factor conferring risk of systemic lupus erythematosus (SLE), but precise independent localization of HLA effects is extremely challenging. As a result, the contribution of specific HLA alleles and amino-acid residues to the overall risk of SLE and to risk of specific autoantibodies are far from completely understood. Here, we dissected (a) overall SLE association signals across HLA, (b) HLA-peptide interaction, and (c) residue-autoantibody association. Classical alleles, SNPs, and amino-acid residues of eight HLA genes were imputed across 4,915 SLE cases and 13,513 controls from Eastern Asia. We performed association followed by conditional analysis across HLA, assessing both overall SLE risk and risk of autoantibody production. DR15 alleles HLA-DRB1*15:01 (P = 1.4x10^-27, odds ratio (OR) = 1.57) and HLA-DQB1*06:02 (P = 7.4x10^-23, OR = 1.55) formed the most significant haplotype (OR = 2.33). Conditioned protein-residue signals were stronger than allele signals and mapped predominantly to HLA-DRB1 residue 13 (P = 2.2x10^-75) and its proxy position 11 (P = 1.1x10^-67), followed by HLA-DRB1-37 (P = 4.5x10^-24). After conditioning on HLA-DRB1, novel associations at HLA-A-70 (P = 1.4x10^-8), HLA-DPB1-35 (P = 9.0x10^-16), HLA-DQB1-37 (P = 2.7x10^-14), and HLA-B-9 (P = 6.5x10^-15) emerged. Together, these seven residues increased the proportion of explained heritability due to HLA to 2.6%. Risk residues for both overall disease and hallmark autoantibodies (i.e., nRNP: DRB1-11, P = 2.0x10^-14; DRB1-13, P = 2.9x10^-13; DRB1-30, P = 3.9x10^-14) localized to the peptide-binding groove of HLA-DRB1. Enrichment for specific amino-acid characteristics in the peptide-binding groove correlated with overall SLE risk and with autoantibody presence. Risk residues were in primarily negatively charged side-chains, in contrast with rheumatoid arthritis. We identified novel SLE signals in HLA Class I loci (HLA-A, HLA-B), and localized primary Class II signals to five residues in HLA-DRB1, HLA-DPB1, and HLA-DQB1. These findings provide insights about the mechanisms by which the risk residues interact with each other to produce autoantibodies and are involved in SLE pathophysiology.

The Human leukocyte antigen (HLA) region is a key genetic factor conferring risk of systemic lupus erythematosus (SLE). In spite of multiple SLE association signals identified in the HLA region, only amino-acid residues within HLA-DRB1 have been specifically described previously. In this study, we performed an imputation-based analysis on individuals with East Asian ancestry, and characterized SLE risk within the HLA region for all involved independent genes (HLA-DRB1, HLA-DPB1, HLA-DQB1, HLA-A, and HLA-B). Furthermore, we identified a characteristic SLE risk residue signature as well as a pattern of specific nRNP and Ro/La autoantibody residues located in the peptide-binding grooves, suggesting their key involvement in autoantibody production.

In silico prediction of nonpermissive HLA-DPB1 mismatches in unrelated HCT by functional distance

In silico prediction of high-risk donor-recipient HLA mismatches after unrelated donor (UD) hematopoietic cell transplantation (HCT) is an attractive, yet elusive, objective. Nonpermissive T-cell epitope (TCE) group mismatches were defined by alloreactive T-cell cross-reactivity for 52/80 HLA-DPB1 alleles (TCE-X). More recently, a numerical functional distance (FD) scoring system for in silico prediction of TCE groups based on the median impact of exon 2-encoded amino acid polymorphism on T-cell alloreactivity was developed for all DPB1 alleles (TCE-FD), including the 28/80 common alleles not assigned by TCE-X. We compared clinical outcome associations of nonpermissive DPB1 mismatches defined by TCE-X or TCE-FD in 8/8 HLA-matched UD-HCT for acute leukemia, myelodysplastic syndrome, and chronic myelogenous leukemia between 1999 and 2011 (N = 2730). Concordance between the 2 models was 92.3%, with most differences arising from DPB1*06:01 and DPB1*19:01 being differently assigned by TCE-X and TCE-FD. In both models, nonpermissive mismatches were associated with reduced overall survival (hazard ratio [HR], 1.15, P < .006 and HR, 1.12, P < .03), increased transplant-related mortality (HR, 1.31, P < .001 and HR, 1.26, P < .001) as well as acute (HR, 1.16, P < .02 and HR, 1.22, P < .001) and chronic (HR, 1.20, P < .003 and HR, 1.22, P < .001) graft-versus-host disease (GVHD). We show that in silico prediction of nonpermissive DPB1 mismatches significantly associated with major transplant outcomes is feasible for any DPB1 allele with known exon 2 sequence based on experimentally elaborated FD scores. This proof-of-principle observation opens new avenues for developing HLA risk-prediction models in HCT and has practical implications for UD searches.

Using Whole Genome Sequencing in an African Subphenotype of Myasthenia Gravis to Generate a Pathogenetic Hypothesis

Myasthenia gravis (MG) is a rare, treatable antibody-mediated disease which is characterized by muscle weakness. The pathogenic antibodies are most frequently directed at the acetylcholine receptors (AChRs) at the skeletal muscle endplate. An ophthalmoplegic subphenotype of MG (OP-MG), which is characterized by treatment resistant weakness of the extraocular muscles (EOMs), occurs in a proportion of myasthenics with juvenile symptom onset and African genetic ancestry. Since the pathogenetic mechanism(s) underlying OP-MG is unknown, the aim of this study was to use a hypothesis-generating genome-wide analysis to identify candidate OP-MG susceptibility genes and pathways. Whole genome sequencing (WGS) was performed on 25 AChR-antibody positive myasthenic individuals of African genetic ancestry sampled from the phenotypic extremes: 15 with OP-MG and 10 individuals with control MG (EOM treatment-responsive). Variants were called according to the Genome Analysis Toolkit (GATK) best practice guidelines using the hg38 reference genome. In addition to single variant association analysis, variants were mapped to genes (±200 kb) using VEGAS2 to calculate gene-based test statistics and HLA allele group assignment was inferred through "best-match" alignment of reads against the IMGT/HLA database. While there were no single variant associations that reached genome-wide significance in this exploratory sample, several genes with significant gene-based test statistics and known to be expressed in skeletal muscle had biological functions which converge on muscle atrophy signaling and myosin II function. The closely linked HLA-DPA1 and HLA-DPB1 genes were associated with OP-MG subjects (gene-based p < 0.05) and the frequency of a functional A > G SNP (rs9277534) in the HLA-DPB1 3'UTR, which increases HLA-DPB1 expression, differed between the two groups (G-allele 0.30 in OP-MG vs. 0.60 in control MG; p = 0.04). Furthermore, we show that rs9277534 is an HLA-DBP1 expression quantitative trait locus in patient-derived myocytes (p < 1 × 10^-3). The application of a SNP to gene to pathway approach to this exploratory WGS dataset of African myasthenic individuals, and comparing dichotomous subphenotypes, resulted in the identification of candidate genes and pathways that may contribute to OP-MG susceptibility. Overall, the hypotheses generated by this work remain to be verified by interrogating candidate gene and pathway expression in patient-derived extraocular muscle.

High-resolution characterization of allelic and haplotypic HLA frequency distribution in a Spanish population using high-throughput next-generation sequencing

Next-generation sequencing (NGS) at the HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, -DRB1 and -DRB3/4/5 loci was performed on 282 healthy unrelated individuals from different major regions of Spain. High-resolution HLA genotypes defined by full sequencing of class I loci and extended coverage of class II loci were obtained to determine allele frequencies and also to estimate extended haplotype frequencies. HLA alleles were typed at the highest resolution level (4-field level, 4FL); with exception of a minor deviation in HLA-DPA1, no statistically significant deviations from expected Hardy Weinberg Equilibrium (HWE) proportions were observed for all other HLA loci. This study provides new 4FL-allele and -haplotype frequencies estimated for the first time in the Spanish population. Furthermore, our results describe extended haplotypes (including the less frequently typed HLA-DPA1 and HLA-DQA1 loci) and show distinctive haplotype associations found at 4FL-allele definition in this Spanish population study. The distinctive allelic and haplotypic diversity found at the 4FL reveals the high level of heterozygosity and specific haplotypic associations displayed that were not apparent at 2-field level (2FL). Overall, these results may contribute as a useful reference source for future population studies, for HLA-disease association studies as a healthy control group dataset and for improving donor recruitment strategies of bone marrow registries. HLA genotyping data of this Spanish population cohort was also included in the 17th International Histocompatibility and Immunogenetics Workshop (IHIW) as part of the study of HLA diversity in unrelated worldwide populations using NGS.

Dissecting Genetic Control of HLA-DPB1 Expression and Its Relation to Structural Mismatch Models in Hematopoietic Stem Cell Transplantation

HLA expression levels have been suggested to be genetically controlled by single nucleotide polymorphisms (SNP) in the untranslated regions (UTR), and expression variants have been associated with the outcome of chronic viral infection and hematopoietic stem cell transplantation (HSCT). In particular, the 3′UTR rs9277534-G/A SNP in HLA-DPB1 has been associated with graft-versus-host-disease after HSCT (Expression model); however its relevance in different immune cells and its mode of action have not been systematically addressed. In addition, there is a strong though not complete overlap between the rs9277534-G/A SNP and structural HLA-DPB1 T cell epitope (TCE) groups which have also been associated with HSCT outcome (TCE Structural model). Here we confirm and extend previous findings of significantly higher HLA-DPB1 expression in B cell lines, unstimulated primary B cells, and monocytes homozygous for rs9277534-G compared to those homozygous for rs9277534-A. However, these differences were abrogated by interferon-γ stimulation or differentiation into dendritic cells. We identify at least seven 3′UTR rs9277534-G/A haplotypes differing by a total of 37 SNP, also characterized by linkage to length variants of a short tandem repeat (STR) in intron 2 and TCE group assignment. 3′UTR mapping did not show any significant differences in post-transcriptional regulation assessed by luciferase assays between two representative rs9277534-G/A haplotypes for any of eight overlapping fragments. Moreover, no evidence for alternative splicing associated with the intron 2 STR was obtained by RT-PCR. In an exemplary cohort of 379 HLA-DPB1 mismatched donor-recipient pairs, risk prediction by the Expression model and the Structural TCE model was 36.7% concordant, with the majority of discordances due to non-applicability of the Expression model. HLA-DPB1 from different TCE groups expressed in the absence of the 3′UTR at similar levels by transfected HeLa cells elicited significantly different mean alloreactive CD4+ T-cell responses, as assessed by CD137 upregulation assays in 178 independent cultures. Taken together, our data provide new insights into the cell type-specific and mechanistic basis of the association between the rs9277534-G/A SNP and HLA-DPB1 expression, and show that, despite partial overlap between both models in HSCT risk-prediction, differential alloreactivity determined by the TCE structural model occurs independently from HLA-DPB1 differential expression.

In Silico Prediction Analysis of Idiotope-Driven T-B Cell Collaboration in Multiple Sclerosis

Memory B cells acting as antigen-presenting cells are believed to be important in multiple sclerosis (MS), but the antigen they present remains unknown. We hypothesized that B cells may activate CD4⁺ T cells in the central nervous system of MS patients by presenting idiotopes from their own immunoglobulin variable regions on human leukocyte antigen (HLA) class II molecules. Here, we use bioinformatics prediction analysis of B cell immunoglobulin variable regions from 11 MS patients and 6 controls with other inflammatory neurological disorders (OINDs), to assess whether the prerequisites for such idiotope-driven T–B cell collaboration are present. Our findings indicate that idiotopes from the complementarity determining region (CDR) 3 of MS patients on average have high predicted affinities for disease associated HLA-DRB1*15:01 molecules and are predicted to be endosomally processed by cathepsin S and L in positions that allows such HLA binding to occur. Additionally, complementarity determining region 3 sequences from cerebrospinal fluid (CSF) B cells from MS patients contain on average more rare T cell-exposed motifs that could potentially escape tolerance and stimulate CD4⁺ T cells than CSF B cells from OIND patients. Many of these features were associated with preferential use of the IGHV4 gene family by CSF B cells from MS patients. This is the first study to combine high-throughput sequencing of patient immune repertoires with large-scale prediction analysis and provides key indicators for future in vitro and in vivo analyses.

HLA-DP in unrelated hematopoietic cell transplantation revisited: challenges and opportunities

When considering HLA-matched hematopoietic cell transplantation (HCT), sibling and unrelated donors (UDs) are biologically different because UD-HCT is typically performed across HLA-DP disparities absent in sibling HCT. Mismatched HLA-DP is targeted by direct alloreactive T cell responses with important implications for graft-versus-host disease and graft-versus-leukemia. This concise review details special features of HLA-DP as model antigens for clinically permissive mismatches mediating limited T-cell alloreactivity with minimal toxicity, and describes future avenues for their exploitation in cellular immunotherapy of malignant blood disorders.

Directionality of non-permissive HLA-DPB1 T-cell epitope group mismatches does not improve clinical risk stratification in 8/8 matched unrelated donor hematopoietic cell transplantation

In 8/8 HLA-matched unrelated donor (UD) hematopoietic cell transplants (HCT), HLA-DPB1 mismatches between alleles from different T-cell epitope (TCE) groups (non-permissive mismatches) are associated with significantly higher risks of mortality compared with those between alleles from the same TCE group (permissive mismatches); however, the relevance of mismatch directionality, that is (host vs graft (uni-directional HvG), graft vs host (uni-directional GvH) or both (bi-directional) in the non-permissive setting is unknown. We show here significantly higher in vitro relative responses (RR) to bi-directional mismatches compared with uni-directional HvG or GvH mismatches in a total of 420 one-way mixed lymphocyte reactions between 10/10 matched pairs (RR 27.5 vs 7.5 vs 15.5, respectively, P<0.001). However, in 3281 8/8 matched UD HCT for leukemia or myelodysplastic syndrome, the hazards of transplant-related mortality (TRM) were similar for uni-directional HvG or GvH mismatches and bi-directional mismatches (hazard ratio (HR) 1.32, P=0.001 vs HR 1.28, P=0.005 and HR 1.34, P=0.046), compared with permissive mismatches. Similar results were observed for overall survival. No statistical differences between the uni- and the bi-directional non-permissive groups were detected in pairwise comparisons for any of the outcomes tested. We conclude that consideration of directionality does not improve risk stratification by non-permissive HLA-DPB1 TCE mismatches in UD searches.

Identification of amino acids in antigen-binding site of class II HLA proteins independently associated with hepatitis B vaccine response

BACKGROUND & AIMS

Genetic factors in class II human leukocyte antigen (HLA) have been reported to be associated with inter-individual variation in hepatitis B virus (HBV) vaccine response. However, the mechanism underlying the associations remains elusive. In particular, the broad linkage disequilibrium in HLA region complicates the localization of the independent effects of genetic variants. Thus, the present study aimed to identify the most probable causal variations in class II HLA loci involved in the immune response to HBV vaccine.

METHODS

We performed a case-control study to assess whether HLA-DRB1, -DQB1, and -DPB1 4-digit alleles were associated with the response to primary HBV vaccination in 574 healthy Japanese students. To identify causative variants, we next assessed independently associated amino acid variants in these loci using conditional logistic regression analysis. Furthermore, to clarify the functional effects of these variants on HLA proteins, we performed computational structural studies.

RESULTS

HLA-DRB1^∗01:01, HLA-DRB1^∗08:03, HLA-DQB1^∗05:01, and HLA-DPB1^∗04:02 were significantly associated with sufficient response, whereas HLA-DPB1^∗05:01 was associated with poor response. We then identified amino acids independently associated with sufficient response, namely, leucine at position 26 of HLA-DRβ1 and glycine-glycine-proline-methionine at positions 84-87 of HLA-DPβ1. These amino acids were located in antigen-binding pocket 4 of HLA-DR and pocket 1 of HLA-DP, respectively, which are important structures for selective binding of antigenic peptides. In addition, the detected variations in HLA-DP protein were responsible for the differences in the electrostatic potentials of the pocket, which can explain in part the sufficient/poor vaccine responses.

CONCLUSION

HLA-DRβ1 position 26 and HLA-DPβ1 positions 84-87 are independently associated with anti-HBs production against HBV vaccine. Our results suggest that HBsAg presentation through these HLA pocket structures plays an important role in the inter-individual variability of HBV vaccination.

Shared HLA Class I and II Alleles and Clonally Restricted Public and Private Brain-Infiltrating αβ T Cells in a Cohort of Rasmussen Encephalitis Surgery Patients

Rasmussen encephalitis (RE) is a rare pediatric neuroinflammatory disease characterized by intractable seizures and unilateral brain atrophy. T cell infiltrates in affected brain tissue and the presence of circulating autoantibodies in some RE patients have indicated that RE may be an autoimmune disease. The strongest genetic links to autoimmunity reside in the MHC locus, therefore, we determined the human leukocyte antigen (HLA) class I and class II alleles carried by a cohort of 24 RE surgery cases by targeted in-depth genomic sequencing. Compared with a reference population the allelic frequency of three alleles, DQA1*04:01:01, DQB1*04:02:01, and HLA-C*07:02:01:01 indicated that they might confer susceptibility to the disease. It has been reported that HLA-C*07:02 is a risk factor for Graves disease. Further, eight patients in the study cohort carried HLA-A*03:01:01:01, which has been linked to susceptibility to multiple sclerosis. Four patients carried a combination of three HLA class II alleles that has been linked to type 1 diabetes (DQA1*05:01:01:01~DQB1*02:01:01~DRB1*03:01:01:01), and five patients carried a combination of HLA class II alleles that has been linked to the risk of contracting multiple sclerosis (DQA1*01:02:01:01, DQB1*06:02:01, DRB1*15:01:01:01). We also analyzed the diversity of αβ T cells in brain and blood specimens from 14 of these RE surgery cases by sequencing the third complementarity regions (CDR3s) of rearranged T cell receptor β genes. A total of 31 unique CDR3 sequences accounted for the top 5% of all CDR3 sequences in the 14 brain specimens. Thirteen of these sequences were found in sequencing data from healthy blood donors; the remaining 18 sequences were patient specific. These observations provide evidence for the clonal expansion of public and private T cells in the brain, which might be influenced by the RE patient’s HLA haplotype.

Identification of DPB1 Permissive Unrelated Donors Is Highly Likely

HLA-DPB1 permissive matching based on T cell epitope (TCE) groups should be considered when selecting among equally matched HLA-A, -B, -C, -DRB1 unrelated hematopoietic stem cell donors to improve patient survival. Previous studies have defined 3 TCE groups based on functional assays of alloreactivity. Combinations of donor and recipient DPB1 alleles with low immunogenic potential identify permissive donors, who provide no increased risk of mortality compared with DPB1-matched donors. To determine the likelihood of identifying a DPB1 permissive-matched (includes both allele-matched and DPB1-permissive mismatched) unrelated donor for patients with high-resolution matches at 10/10 HLA-A, -B,- C, -DRB1, and -DQB1 in the Be The Match Registry, preliminary search requests from United States' transplant centers for 595 DPB1-typed patients were evaluated for existence of a DPB1 permissive-matched donor, identified either among already typed donors or by prospective DPB1 typing. The baseline DPB1 permissive match rate was 69% and improved to 80% after additional donor DPB1 typing (median, 4 donors per patient). When seeking a 10/10-matched, young (18- to 32-year-old) donor in the registry, the probability of finding a DPB1 permissive-matched donor started lower at 59% and improved to 70% after additional DPB1 testing. Our results show that most patients with a 10/10 match can find a DPB1 permissive-matched donor.

Antibodies against HLA-DP recognize broadly expressed epitopes

HLA matching and avoidance of pre-transplant donor-specific antibodies are important in selection of donors for solid organ transplant. Solid phase testing with single antigen beads allows resolution of antibody reactivity to the level of the allele. Single antigen bead testing results at a large transplant center were reviewed to identify selective reactivity patterns of anti-HLA antibodies. Many HLA-DP antibodies were identified in the context of other HLA antibodies, but some sera had antibodies against only HLA-DP. B cell flow crossmatch testing was positive for 2 out of 9 sera with HLA-DP antibodies. Many patterns of reactivity corresponded to epitopes in hypervariable regions C and F of DPB1, but some matched epitopes in other regions or DPA1. Through analysis of single antigen bead testing from a large number of patients, we report that anti-HLA-DP antibodies predominantly recognize broadly cross-reactive epitopes. The United Network for Organ Sharing has mandated HLA-DP typing on all deceased kidney donors, and HLA-DP epitopes should be considered as the major antigens for avoidance of pre-transplant donor-specific antibodies.

Genetic risk factors for sclerotic graft-versus-host disease

Sclerotic graft-versus-host disease (GVHD) is a distinctive phenotype of chronic GVHD after allogeneic hematopoietic cell transplantation, characterized by fibrosis of skin or fascia. Sclerotic GVHD has clinical and histopathological similarities with systemic sclerosis, an autoimmune disease whose risk is influenced by genetic polymorphisms. We examined 13 candidate single-nucleotide polymorphisms (SNPs) that have a well-documented association with systemic sclerosis to determine whether these SNPs are also associated with the risk of sclerotic GVHD. The study cohort included 847 consecutive patients who were diagnosed with chronic GVHD. Genotyping was performed using microarrays, followed by imputation of unobserved SNPs. The donor rs10516487 (BANK1: B-cell scaffold protein with ankyrin repeats 1) TT genotype was associated with lower risk of sclerotic GVHD (hazard ratio [HR], 0.43; 95% confidence interval [CI], 0.21-0.87; P = .02). Donor and recipient rs2056626 (CD247: T-cell receptor ζ subunit) GG or GT genotypes were associated with higher risk of sclerotic GVHD (HR, 1.57; 95% CI, 1.13-2.18; P = .007 and HR, 1.66; 95% CI, 1.19-2.32; P = .003, respectively). Donor and recipient rs987870 (5'-flanking region of HLA-DPA1) CC genotypes were associated with higher risk of sclerotic GVHD (HR, 2.50; 95% CI, 1.22-5.11; P = .01 and HR, 2.13; 95% CI, 1.00-4.54; P = .05, respectively). In further analyses, the recipient DPA1*01:03∼DPB1*04:01 haplotype and certain amino acid substitutions in the recipient P1 peptide-binding pocket of the HLA-DP heterodimer were associated with risk of sclerotic GVHD. Genetic components associated with systemic sclerosis are also associated with sclerotic GVHD. HLA-DP-mediated antigen presentation, T-cell response, and B-cell activation have important roles in the pathogenic mechanisms of both diseases.

Functional distance between recipient and donor HLA-DPB1 determines nonpermissive mismatches in unrelated HCT

The role of HLA amino acid (AA) polymorphism for the outcome of hematopoietic cell transplantation (HCT) is controversial, in particular for HLA class II. Here, we investigated this question in nonpermissive HLA-DPB1 T-cell epitope (TCE) mismatches reflected by numerical functional distance (FD) scores, assignable to all HLA-DPB1 alleles based on the combined impact of 12 polymorphic AAs. We calculated the difference in FD scores (ΔFD) of mismatched HLA-DPB1 alleles in patients and their 10/10 HLA-matched unrelated donors of 379 HCTs performed at our center for acute leukemia or myelodysplastic syndrome. Receiver-operator curve-based stratification into 2 ΔFD subgroups showed a significantly higher percentage of nonpermissive TCE mismatches for ΔFD >2.665, compared with ΔFD ≤2.665 (88% vs 25%, P < .0001). In multivariate analysis, ΔFD >2.665 was significantly associated with overall survival (hazard ratio [HR], 1.40; 95% confidence interval [CI], 1.05-1.87; P < .021) and event-free survival (HR, 1.39; 95% CI, 1.05-1.82; P < .021), compared with ΔFD ≤2.665. These associations were stronger than those observed for TCE mismatches. There was a marked but not statistically significant increase in the hazards of relapse and nonrelapse mortality in the high ΔFD subgroup, whereas no differences were observed for acute and chronic graft-versus-host disease. Seven nonconservative AA substitutions in peptide-binding positions had a significantly stronger impact on ΔFD compared with 5 others (P = .0025), demonstrating qualitative differences in the relative impact of AA polymorphism in HLA-DPB1. The novel concept of ΔFD sheds new light onto nonpermissive HLA-DPB1 mismatches in unrelated HCT.

First report on HLA-DPA1 gene allelic distribution in the general Lebanese population

Aims

HLA-DPA1 is an important marker in bone marrow and organ transplantation and a highly emerging screening parameter in histocompatibility laboratories. Being highly polymorphic, it has another significant value in detecting population origins and migrations. This is the first study to assess DPA1 allele frequencies in an Arab population.

Methods

The HLA DPA1 alleles were identified using the One-Lambda assays on a Luminex reverse SSO DNA typing system. Our study included 101 individuals coming from different Lebanese geographical areas representing the different communities and religious sects of the country.

Results

We compared the results of this study to 16 different populations and found very interesting similarities and differences between Lebanese people and individuals of European ancestry.

Conclusion

This study is the first to describe the different allelic frequencies of HLA-DPA1 in the Lebanese population and will serve as a template that can be later used for disease association studies both at the level of the country and internationally.

A Diverse Repertoire of CD4 T Cells Targets the Immediate-Early 1 Protein of Human Cytomegalovirus

T-cell responses to the immediate-early 1 (IE-1) protein of human cytomegalovirus (HCMV) are associated with protection from viral disease. Thus, IE-1 is a promising target for immunotherapy. CD8 T-cell responses to IE-1 are generally strong. In contrast, CD4 T-cell responses to IE-1 were described to be comparatively infrequent or undetectable in HCMV carriers, and information on their target epitopes and their function has been limited. To analyze the repertoire of IE-1-specific CD4 T cells, we expanded them from healthy donors with autologous IE-1-expressing mini-Epstein–Barr virus-transformed B-cell lines and established IE-1-specific CD4 T-cell clones. Clones from seven out of seven HCMV-positive donors recognized endogenously processed IE-1 epitopes restricted through HLA-DR, DQ, or DP. Three to seven IE-1 epitopes were recognized per donor. Cumulatively, about 27 different HLA/peptide class II complexes were recognized by 117 IE-1-specific clones. Our results suggest that a highly diversified repertoire of IE-1-specific CD4 T cells targeting multiple epitopes is usually present in healthy HCMV carriers. Therefore, multiepitope approaches to immunomonitoring and immunotherapy will make optimal use of this potentially important class of HCMV-specific effector cells.

HLA-DPB1 and HLA class I confer risk of and protection from narcolepsy

Type 1 narcolepsy, a disorder caused by a lack of hypocretin (orexin), is so strongly associated with human leukocyte antigen (HLA) class II HLA-DQA1(∗)01:02-DQB1(∗)06:02 (DQ0602) that very few non-DQ0602 cases have been reported. A known triggering factor for narcolepsy is pandemic 2009 influenza H1N1, suggesting autoimmunity triggered by upper-airway infections. Additional effects of other HLA-DQ alleles have been reported consistently across multiple ethnic groups. Using over 3,000 case and 10,000 control individuals of European and Chinese background, we examined the effects of other HLA loci. After careful matching of HLA-DR and HLA-DQ in case and control individuals, we found strong protective effects of HLA-DPA1(∗)01:03-DPB1(∗)04:02 (DP0402; odds ratio [OR] = 0.51 [0.38-0.67], p = 1.01 × 10(-6)) and HLA-DPA1(∗)01:03-DPB1(∗)04:01 (DP0401; OR = 0.61 [0.47-0.80], p = 2.07 × 10(-4)) and predisposing effects of HLA-DPB1(∗)05:01 in Asians (OR = 1.76 [1.34-2.31], p = 4.71 × 10(-05)). Similar effects were found by conditional analysis controlling for HLA-DR and HLA-DQ with DP0402 (OR = 0.45 [0.38-0.55] p = 8.99 × 10(-17)) and DP0501 (OR = 1.38 [1.18-1.61], p = 7.11 × 10(-5)). HLA-class-II-independent associations with HLA-A(∗)11:01 (OR = 1.32 [1.13-1.54], p = 4.92 × 10(-4)), HLA-B(∗)35:03 (OR = 1.96 [1.41-2.70], p = 5.14 × 10(-5)), and HLA-B(∗)51:01 (OR = 1.49 [1.25-1.78], p = 1.09 × 10(-5)) were also seen across ethnic groups in the HLA class I region. These effects might reflect modulation of autoimmunity or indirect effects of HLA class I and HLA-DP alleles on response to viral infections such as that of influenza.

Risk associations between HLA-DPB1 T-cell epitope matching and outcome of unrelated hematopoietic cell transplantation are independent of HLA-DPA1

HLA-DP antigens are beta-alpha heterodimers encoded by polymorphic HLA-DPB1 and -DPA1 alleles, respectively, in strong linkage disequilibrium (LD) with each other. Non-permissive unrelated donor (UD)-recipient HLA-DPB1 mismatches across three different T-cell epitope (TCE) groups are associated with increased mortality after hematopoietic SCT (HCT), but the role of HLA-DPA1 is unclear. We studied 1281 onco-hematologic patients after 10/10 HLA-matched UD-HCT facilitated by the National Marrow Donor Program. Non-permissive mismatches defined solely by HLA-DPB1 TCE groups were associated with significantly higher risks of TRM compared to permissive mismatches (hazard ratio (HR) 1.30, confidence interval (CI) 1.06-1.53; P=0.009) or allele matches. Moreover, non-permissive HLA-DPB1 TCE group mismatches in the graft versus host (GvH) direction significantly decreased the risk of relapse compared to permissive mismatches (HR 0.55, CI 0.37-0.80; P=0.002) or allele matches. Splitting each group into HLA-DPA1*02:01 positive or negative, in frequent LD with HLA-DPB1 alleles from two of the three TCE groups, or into HLA-DPA1 matched or mismatched, did not significantly alter the observed risk associations. Our findings suggest that the effects of clinically non-permissive HLA-DPB1 TCE group mismatches are independent of HLA-DPA1, and that selection of donors with non-permissive DPB1 TCE mismatches in GvH direction might provide some protection from disease recurrence.

Sequence-based typing of HLA: an improved group-specific full-length gene sequencing approach

Matching for HLA at the allele level is crucial for stem cell transplantation. The golden standard approach for allele definition of full gene polymorphism, the so-called high-resolution HLA typing, is sequence-based typing (SBT). Although the majority of the polymorphism for class I is located in exons 2 and 3 and for class II in exon 2, for allele definition it is necessary to unravel the complete coding and intron sequences leading to an ultrahigh HLA typing resolution at the allele level, i.e., a full-length gene polymorphism identification.This chapter describes our recently developed SBT method for HLA-A, -B, -C, and -DQB1, that is based on full-length hemizygous Sanger sequencing of the alleles, separated by group-specific amplification using the low-resolution typing result as reference starting point. Group-specific amplification has already been established for DRB. This method enables a cost-efficient, user-friendly SBT approach resulting in a timely unambiguous HLA typing to an ultrahigh resolution level with minimal hands-on time.

Effects of transmembrane region variability on cell surface expression and allorecognition of HLA-DP3

The functional relevance of polymorphisms outside the peptide binding groove of HLA molecules is poorly understood. Here we have addressed this issue by studying HLA-DP3, a common antigen relevant for functional matching algorithms of unrelated hematopoietic stem cell transplantation (HSCT) encoded by two transmembrane (TM) region variants, DPB1(*)03:01 and DPB1(*)104:01. The two HLA-DP3 variants were found at a overall allelic frequency of 10.4% in 201 volunteer stem cell donors, at a ratio of 4.2:1. No significant differences were observed in cell surface expression levels of the two variants on B lymphoblastoid cell lines (BLCL), primary B cells or monocytes. Three different alloreactive T cell lines or clones showed similar levels of activation marker CD107a and/or CD137 upregulation in response to HLA-DP3 encoded by DPB1(*)03:01 and DPB1(*)104:01, either endogenously on BLCL or after lentiveral-vector mediated transfer into the same cellular background. These data provide, for the first time, direct evidence for a limited functional role of a TM region polymorphism on expression and allorecognition of HLA-DP3 and are compatible with the notion that the two variants can be considered as a single functional entity for unrelated stem cell donor selection.

Multiple mismatches at the low expression HLA loci DP, DQ, and DRB3/4/5 associate with adverse outcomes in hematopoietic stem cell transplantation

A single mismatch in highly expressed HLA-A, -B, -C, and -DRB1 loci (HEL) is associated with worse outcomes in hematopoietic stem cell transplantation, while less is known about the cumulative impact of mismatches in the lesser expressed HLA loci DRB3/4/5, DQ, and DP (LEL). We studied whether accumulation of LEL mismatches is associated with deleterious effects in 3853 unrelated donor transplants stratified according to number of matches in the HEL. In the 8/8 matched HEL group, LEL mismatches were not associated with any adverse outcome. Mismatches at HLA-DRB1 were associated with occurrence of multiple LEL mismatches. In the 7/8 HEL group, patients with 3 or more LEL mismatches scored in the graft-versus-host vector had a significantly higher risk of mortality (1.45 and 1.43) and transplant-related mortality (1.68 and 1.54) than the subgroups with 0 or 1 LEL mismatches. No single LEL locus had a more pronounced effect on clinical outcome. Three or more LEL mismatches are associated with lower survival after 7/8 HEL matched transplantation. Prospective evaluation of matching for HLA-DRB3/4/5, -DQ, and -DP loci is warranted to reduce posttransplant risks in donor-recipient pairs matched for 7/8 HEL.