Clustering HLA class I superfamilies using structural interaction patterns

Prognostic impact of HLA supertype mismatch in single-unit cord blood transplantation

The "human leukocyte antigen (HLA) supertype" is a functional classification of HLA alleles, which was defined by structural features and peptide specificities, and has been reportedly associated with the clinical outcomes of viral infections and autoimmune diseases. Although the disparity in each HLA locus was reported to have no clinical significance in single-unit cord blood transplantation (sCBT), the clinical significance of the HLA supertype in sCBT remains unknown. Therefore, we retrospectively analyzed clinical data of 1603 patients who received sCBT in eight institutes in Japan between 2000 and 2017. Each HLA allele was categorized into 19 supertypes, and the prognostic effect of disparities was then assessed. An HLA-B supertype mismatch was identified as a poor prognostic factor (PFS: hazard ratio [HR] = 1.23, p = 0.00044) and was associated with a higher cumulative incidence (CI) of relapse (HR = 1.24, p = 0.013). However, an HLA-B supertype mismatch was not associated with the CI of acute and chronic graft-versus-host-disease. The multivariate analysis for relapse and PFS showed the significance of an HLA-B supertype mismatch independent of allelic mismatches, and other previously reported prognostic factors. HLA-B supertype-matched grafts should be selected in sCBT.

The Prognostic Significance of Selected HLA Alleles on Prostate Cancer Outcome

Recently, we have shown that HLA-A*02:01 and HLA-A*24:02 in de novo metastatic prostate cancer (MPCa) have an important role in disease progression. Since de novo MPCa represents a small group among patients diagnosed with prostate cancer (PCa), it was obvious to try to extend the validity of our results to larger cohorts of PCa patients. Herein, we analyzed patients irrespective of their disease status at diagnosis to include, besides patients with MPCa, those with localized PCa (LPCa). Our goal was to specify the prognostic value of HLA-A*02:01 and HLA-A*24:02 for overall survival (OS) prospectively and for early biochemical recurrence (BCR) and castrate resistance (CR) as additional clinical endpoints in a prospective/retrospective manner, to improve clinical decisions for patients covering all stages of PCa. On univariate analysis, HLA-A alleles were significantly associated as prognostic biomarkers with early BCR (p = 0.028; HR = 1.822), OS (p = 0.013; HR = 1.547) and showed a trend for CR (p = 0.150; HR = 1.239). On multivariate analysis, HLA-A alleles proved to be independent prognosticators for early BCR (p = 0.017; HR = 2.008), CR (p = 0.005; HR = 1.615), and OS (p = 0.002; HR = 2.063). Kaplan-Meier analyses revealed that patients belonging to the HLA-A*02:01+HLA-A*24:02- group progressed much faster to BCR and CR and had also shorter OS compared to HLA-A*24:02+ patients. Patients being HLA-A*02:01-HLA-A*24:02- exhibited varying clinical outcomes, pointing to the presence of additional HLA-A alleles with potential prognostic value. Our data underline the HLA-A alleles as valuable prognostic biomarkers for PCa that may assist with the appropriate treatment and follow-up schedule based on the risk for disease progression to avoid over-diagnosis and over-treatment.

Predicted impact of the viral mutational landscape on the cytotoxic response against SARS-CoV-2

The massive assessment of immune evasion due to viral mutations that increase COVID-19 susceptibility can be computationally facilitated. The adaptive cytotoxic T response is critical during primary infection and the generation of long-term protection. Here, potential HLA class I epitopes in the SARS-CoV-2 proteome were predicted for 2,915 human alleles of 71 families using the netMHCIpan EL algorithm. Allele families showed extreme epitopic differences, underscoring genetic variability of protective capacity between humans. Up to 1,222 epitopes were associated with any of the twelve supertypes, that is, allele clusters covering 90% population. Next, from all mutations identified in ~118,000 viral NCBI isolates, those causing significant epitope score reduction were considered epitope escape mutations. These mutations mainly involved non-conservative substitutions at the second and C-terminal position of the ligand core, or total ligand removal by large recurrent deletions. Escape mutations affected 47% of supertype epitopes, which in 21% of cases concerned isolates from two or more sub-continental areas. Some of these changes were coupled, but never surpassed 15% of evaded epitopes for the same supertype in the same isolate, except for B27. In contrast to most supertypes, eight allele families mostly contained alleles with few SARS-CoV-2 ligands. Isolates harboring cytotoxic escape mutations for these families co-existed geographically within sub-Saharan and Asian populations enriched in these alleles according to the Allele Frequency Net Database. Collectively, our findings indicate that escape mutation events have already occurred for half of HLA class I supertype epitopes. However, it is presently unlikely that, overall, it poses a threat to the global population. In contrast, single and double mutations for susceptible alleles may be associated with viral selective pressure and alarming local outbreaks. The integration of genomic, geographical and immunoinformatic information eases the surveillance of variants potentially affecting the global population, as well as minority subpopulations.

The pockets guide to HLA class I molecules

Human leukocyte antigens (HLA) are cell-surface proteins that present peptides to T cells. These peptides are bound within the peptide binding cleft of HLA, and together as a complex, are recognised by T cells using their specialised T cell receptors. Within the cleft, the peptide residue side chains bind into distinct pockets. These pockets ultimately determine the specificity of peptide binding. As HLAs are the most polymorphic molecules in humans, amino acid variants in each binding pocket influences the peptide repertoire that can be presented on the cell surface. Here, we review each of the 6 HLA binding pockets of HLA class I (HLA-I) molecules. The binding specificity of pockets B and F are strong determinants of peptide binding and have been used to classify HLA into supertypes, a useful tool to predict peptide binding to a given HLA. Over the years, peptide binding prediction has also become more reliable by using binding affinity and mass spectrometry data. Crystal structures of peptide-bound HLA molecules provide a means to interrogate the interactions between binding pockets and peptide residue side chains. We find that most of the bound peptides from these structures conform to binding motifs determined from prediction software and examine outliers to learn how these HLAs are stabilised from a structural perspective.

Integrative Analysis of HTNV Glycoprotein Derived MHC II Epitopes by In Silico Prediction and Experimental Validation

Hantaan virus (HTNV), the causative pathogen of hemorrhagic fever with renal syndrome (HFRS), is a negative RNA virus belonging to the Orthohantaviridae family. HTNV envelope glycoprotein (GP), encoded by the genomic medium segment, is immunogenic and is therefore a promising vaccine candidate. Major histocompatibility complex class I (MHC-I) epitopes derived from HTNV has been extensively studied, but little is known of MHC-II epitopes. In silico predictions based on four databases indicated that the full-length HTNV GP has 1121 15-mer epitopes, of which 289 had a high score for binding to the human and murine MHC-II superfamily. It found that epitope ILTVLKFIANIFHTS could potentially bind most MHC-II molecules covering human and murine haplotypes. Dominant epitopes were validated by enzyme-linked immunospot assay of splenocytes from immunized mice; 6 of 10 epitopes supported the predictions including TATYSIVGPANAKVP, TKTLVIGQCIYTITS, FSLLPGVAHSIAVEL, CETYKELKAHGVSCP, CGLYLDRLKPVGSAY, and NLGENPCKIGLQTSS. Conservation analysis of dominant epitopes revealed host–virus interactions without geographic stratification, thus meeting the requirements of candidate vaccines for large-population prophylaxis. These findings provide insight into hantavirus antigenicity and suggest that vaccines targeting MHC-II could provide immune protection in large population to complement symptomatic therapies for the treatment of HFRS.

Neoantigen vaccination induces clinical and immunologic responses in non-small cell lung cancer patients harboring EGFR mutations

Background

Neoantigen (NeoAg) peptides displayed at the tumor cell surface by human leukocyte antigen molecules show exquisite tumor specificity and can elicit T cell mediated tumor rejection. However, few NeoAgs are predicted to be shared between patients, and none to date have demonstrated therapeutic value in the context of vaccination.

Methods

We report here a phase I trial of personalized NeoAg peptide vaccination (PPV) of 24 stage III/IV non-small cell lung cancer (NSCLC) patients who had previously progressed following multiple conventional therapies, including surgery, radiation, chemotherapy, and tyrosine kinase inhibitors (TKIs). Primary endpoints of the trial evaluated feasibility, tolerability, and safety of the personalized vaccination approach, and secondary trial endpoints assessed tumor-specific immune reactivity and clinical responses. Of the 16 patients with epidermal growth factor receptor (EGFR) mutations, nine continued TKI therapy concurrent with PPV and seven patients received PPV alone.

Results

Out of 29 patients enrolled in the trial, 24 were immunized with personalized NeoAg peptides. Aside from transient rash, fatigue and/or fever observed in three patients, no other treatment-related adverse events were observed. Median progression-free survival and overall survival of the 24 vaccinated patients were 6.0 and 8.9 months, respectively. Within 3–4 months following initiation of PPV, seven RECIST-based objective clinical responses including one complete response were observed. Notably, all seven clinical responders had EGFR-mutated tumors, including four patients that had continued TKI therapy concurrently with PPV. Immune monitoring showed that five of the seven responding patients demonstrated vaccine-induced T cell responses against EGFR NeoAg peptides. Furthermore, two highly shared EGFR mutations (L858R and T790M) were shown to be immunogenic in four of the responding patients, all of whom demonstrated increases in peripheral blood neoantigen-specific CD8+ T cell frequencies during the course of PPV.

Conclusions

These results show that personalized NeoAg vaccination is feasible and safe for advanced-stage NSCLC patients. The clinical and immune responses observed following PPV suggest that EGFR mutations constitute shared, immunogenic neoantigens with promising immunotherapeutic potential for large subsets of NSCLC patients. Furthermore, PPV with concurrent EGFR inhibitor therapy was well tolerated and may have contributed to the induction of PPV-induced T cell responses.

Global analysis of shared T cell specificities in human non-small cell lung cancer enables HLA inference and antigen discovery

To identify disease-relevant T cell receptors (TCRs) with shared antigen specificity, we analyzed 778,938 TCRβ chain sequences from 178 non-small cell lung cancer patients using the GLIPH2 (grouping of lymphocyte interactions with paratope hotspots 2) algorithm. We identified over 66,000 shared specificity groups, of which 435 were clonally expanded and enriched in tumors compared to adjacent lung. The antigenic epitopes of one such tumor-enriched specificity group were identified using a yeast peptide-HLA A^∗02:01 display library. These included a peptide from the epithelial protein TMEM161A, which is overexpressed in tumors and cross-reactive epitopes from Epstein-Barr virus and E. coli. Our findings suggest that this cross-reactivity may underlie the presence of virus-specific T cells in tumor infiltrates and that pathogen cross-reactivity may be a feature of multiple cancers. The approach and analytical pipelines generated in this work, as well as the specificity groups defined here, present a resource for understanding the T cell response in cancer.

Pharmacogenetic and Association Studies on the Influence of HLA Alleles and Rivastigmine on the Iranian Patients with Late-Onset Alzheimer's Disease

Alzheimer's disease (AD) is a progressive neurodegenerative disorder affecting cognitive function. A number of allelic genes from HLA complex have shown variable associations with AD in different populations. In this study, we investigated the association of DQB1*06:00/x, DRB1*04:00/x, DRB1*15:00/x, and B*07:00/x genotypes with AD and their relevance to the efficacy of rivastigmine treatment in the Iranian population. Our findings suggest that DQB1*06:00/x genotype offers strong protection against AD (P = 0.0074), while B*07:00/x genotype imposes a significant susceptibility for sporadic Alzheimer's disease (SAD) (P = 0.009). Interestingly, B*07:00/x genotype does not show any apparent associations with familial Alzheimer's disease (FAD). Our studies also suggest a pharmacogenetic relationship between drug treatment and presence of a particular genotype in the Iranian LOAD patient population. The Clinical Dementia Rating analysis showed that LOAD patients carrying DRB1*04:00/x genotype tend to display a downward trend in the disease severity and symptoms after 2-year follow-up with rivastigmine treatment. Moreover, in our total patient population, the carriers of DQB1*06:00/x and B*07:00/x alleles have better and worse responses to rivastigmine respectively. We also measured the clinical relevance of the testing for these genotypes employing prevalence-corrected positive predictive value (PcPPV) formula. The PcPPV of testing for DQB1*06:00/x in the Iranian LOAD patients was 1.17% which means that people carrying this genotype have half of the probability of the absolute risk for developing LOAD, whereas the PcPPV of testing for B*07:00/x was 4.45% for SAD, which can be interpreted as a doubling chance for developing LOAD among the Iranian population carrying this genotype. These results also suggest that DQβ1 peptide containing positively charged AAs histidine³⁰ and arginine⁵⁵ and HLA class I β chain containing negatively charges aspartic acid¹¹⁴ and glutamic acid^45,152 in their binding groove plays important roles in protection against and susceptibility for LOAD respectively.

Sequence-structure-function relationships in class I MHC: A local frustration perspective

Class I Major Histocompatibility Complex (MHC) binds short antigenic peptides with the help of Peptide Loading Complex (PLC), and presents them to T-cell Receptors (TCRs) of cytotoxic T-cells and Killer-cell Immunglobulin-like Receptors (KIRs) of Natural Killer (NK) cells. With more than 10000 alleles, human MHC (Human Leukocyte Antigen, HLA) is the most polymorphic protein in humans. This allelic diversity provides a wide coverage of peptide sequence space, yet does not affect the three-dimensional structure of the complex. Moreover, TCRs mostly interact with HLA in a common diagonal binding mode, and KIR-HLA interaction is allele-dependent. With the aim of establishing a framework for understanding the relationships between polymorphism (sequence), structure (conserved fold) and function (protein interactions) of the human MHC, we performed here a local frustration analysis on pMHC homology models covering 1436 HLA I alleles. An analysis of local frustration profiles indicated that (1) variations in MHC fold are unlikely due to minimally-frustrated and relatively conserved residues within the HLA peptide-binding groove, (2) high frustration patches on HLA helices are either involved in or near interaction sites of MHC with the TCR, KIR, or tapasin of the PLC, and (3) peptide ligands mainly stabilize the F-pocket of HLA binding groove.

Scrutinizing human MHC polymorphism: Supertype analysis using Poisson-Boltzmann electrostatics and clustering

Peptide-binding MHC proteins are thought the most variable across the human population; the extreme MHC polymorphism observed is functionally important and results from constrained divergent evolution. MHCs have vital functions in immunology and homeostasis: cell surface MHC class I molecules report cell status to CD8+ T cells, NKT cells and NK cells, thus playing key roles in pathogen defence, as well as mediating smell recognition, mate choice, Adverse Drug Reactions, and transplantation rejection. MHC peptide specificity falls into several supertypes exhibiting commonality of binding. It seems likely that other supertypes exist relevant to other functions. Since comprehensive experimental characterization is intractable, structure-based bioinformatics is the only viable solution. We modelled functional MHC proteins by homology and used calculated Poisson-Boltzmann electrostatics projected from the top surface of the MHC as multi-dimensional descriptors, analysing them using state-of-the-art dimensionality reduction techniques and clustering algorithms. We were able to recover the 3 MHC loci as separate clusters and identify clear sub-groups within them, vindicating unequivocally our choice of both data representation and clustering strategy. We expect this approach to make a profound contribution to the study of MHC polymorphism and its functional consequences, and, by extension, other burgeoning structural systems, such as GPCRs.

How an alloreactive T-cell receptor achieves peptide and MHC specificity

T-cell receptor (TCR) allorecognition is often presumed to be relatively nonspecific, attributable to either a TCR focus on exposed major histocompatibility complex (MHC) polymorphisms or the degenerate recognition of allopeptides. However, paradoxically, alloreactivity can proceed with high peptide and MHC specificity. Although the underlying mechanisms remain unclear, the existence of highly specific alloreactive TCRs has led to their use as immunotherapeutics that can circumvent central tolerance and limit graft-versus-host disease. Here, we show how an alloreactive TCR achieves peptide and MHC specificity. The HCV1406 TCR was cloned from T cells that expanded when a hepatitis C virus (HCV)-infected HLA-A2^- individual received an HLA-A2⁺ liver allograft. HCV1406 was subsequently shown to recognize the HCV nonstructural protein 3 (NS3):1406-1415 epitope with high specificity when presented by HLA-A2. We show that NS3/HLA-A2 recognition by the HCV1406 TCR is critically dependent on features unique to both the allo-MHC and the NS3 epitope. We also find cooperativity between structural mimicry and a crucial peptide "hot spot" and demonstrate its role, along with the MHC, in directing the specificity of allorecognition. Our results help explain the paradox of specificity in alloreactive TCRs and have implications for their use in immunotherapy and related efforts to manipulate TCR recognition, as well as alloreactivity in general.

Human leukocyte antigen supertype matching after myeloablative hematopoietic cell transplantation with 7/8 matched unrelated donor allografts: a report from the Center for International Blood and Marrow Transplant Research

The diversity of the human leukocyte antigen (HLA) class I and II alleles can be simplified by consolidating them into fewer supertypes based on functional or predicted structural similarities in epitope-binding grooves of HLA molecules. We studied the impact of matched and mismatched HLA-A (265 versus 429), -B (230 versus 92), -C (365 versus 349), and -DRB1 (153 versus 51) supertypes on clinical outcomes of 1934 patients with acute leukemias or myelodysplasia/myeloproliferative disorders. All patients were reported to the Center for International Blood and Marrow Transplant Research following single-allele mismatched unrelated donor myeloablative conditioning hematopoietic cell transplantation. Single mismatched alleles were categorized into six HLA-A (A01, A01A03, A01A24, A02, A03, A24), six HLA-B (B07, B08, B27, B44, B58, B62), two HLA-C (C1, C2), and five HLA-DRB1 (DR1, DR3, DR4, DR5, DR9) supertypes. Supertype B mismatch was associated with increased risk of grade II-IV acute graft-versus-host disease (hazard ratio =1.78, P=0.0025) compared to supertype B match. Supertype B07-B44 mismatch was associated with a higher incidence of both grade II-IV (hazard ratio=3.11, P=0.002) and III-IV (hazard ratio=3.15, P=0.01) acute graft-versus-host disease. No significant associations were detected between supertype-matched versus -mismatched groups at other HLA loci. These data suggest that avoiding HLA-B supertype mismatches can mitigate the risk of grade II-IV acute graft-versus-host disease in 7/8-mismatched unrelated donor hematopoietic cell transplantation when multiple HLA-B supertype-matched donors are available. Future studies are needed to define the mechanisms by which supertype mismatching affects outcomes after alternative donor hematopoietic cell transplantation.

HLA superfamily assignment is a predictor of immune response to cancer testis antigens and survival in ovarian cancer

OBJECTIVES

To characterize the association between major histocompatibility complex (MHC) types and spontaneous antibody development to the cancer testis (CT) antigen NY-ESO-1.

METHODS

Tumor expression of NY-ESO-1 and serum antibodies to NY-ESO-1 were characterized in addition to human leukocyte antigen (HLA) type for patients with epithelial ovarian cancer. HLA types were assigned to structure-based superfamilies and statistical associations were examined. HLA types were compared to existing reference libraries of HLA frequencies in a European-Caucasian American population.

RESULTS

Out of 126 patients identified, 81% were expression positive and 48% had spontaneous antibody responses to NY-ESO-1. There was an association between HLA-B superfamily and seropositivity among patients with tumors expressing NY-ESO-1 (p<0.001). The differences in HLA-B superfamily assignment were driven by HLA-B44. Among all patients, the B27 superfamily was over-represented compared with the general population (p<0.001).

CONCLUSIONS

HLA type appears to be associated with spontaneous anti-CT antigen antibodies, as well as with the overall risk of ovarian cancer.

Deciphering complex patterns of class-I HLA-peptide cross-reactivity via hierarchical grouping

T-cell responses in humans are initiated by the binding of a peptide antigen to a human leukocyte antigen (HLA) molecule. The peptide-HLA complex then recruits an appropriate T cell, leading to cell-mediated immunity. More than 2000 HLA class-I alleles are known in humans, and they vary only in their peptide-binding grooves. The polymorphism they exhibit enables them to bind a wide range of peptide antigens from diverse sources. HLA molecules and peptides present a complex molecular recognition pattern, as many peptides bind to a given allele and a given peptide can be recognized by many alleles. A powerful grouping scheme that not only provides an insightful classification, but is also capable of dissecting the physicochemical basis of recognition specificity is necessary to address this complexity. We present a hierarchical classification of 2010 class-I alleles by using a systematic divisive clustering method. All-pair distances of alleles were obtained by comparing binding pockets in the structural models. By varying the similarity thresholds, a multilevel classification was obtained, with 7 supergroups, each further subclassifying to yield 72 groups. An independent clustering performed based only on similarities in their epitope pools correlated highly with pocket-based clustering. Physicochemical feature combinations that best explain the basis of clustering are identified. Mutual information calculated for the set of peptide ligands enables identification of binding site residues contributing to peptide specificity. The grouping of HLA molecules achieved here will be useful for rational vaccine design, understanding disease susceptibilities and predicting risk of organ transplants.