Physiochemical disparity of mismatched HLA class I alloantigens and risk of acute GVHD following HSCT

Maximizing matching, equity and survival in kidney transplantation using molecular HLA immunogenicity quantitation

HLA matching improves long-term outcomes of kidney transplantation, yet implementation challenges persist, particularly within the African American (Black) patient demographic due to donor scarcity. Consequently, kidney survival rates among Black patients significantly lag behind those of other racial groups. A refined matching scheme holds promise for improving kidney survival, with prioritized matching for Black patients potentially bolstering rates of HLA-matched transplants. To facilitate quantity, quality and equity in kidney transplants, we propose two matching algorithms based on quantification of HLA immunogenicity using the hydrophobic mismatch score (HMS) for prospective transplants. We mined the national transplant patient database (SRTR) for a diverse group of donors and recipients with known racial backgrounds. Additionally, we use novel methods to infer survival assessment in the simulated transplants generated by our matching algorithms, in the absence of actual target outcomes, utilizing modified unsupervised clustering techniques. Our allocation algorithms demonstrated the ability to match 87.7% of Black and 86.1% of White recipients under the HLA immunogenicity threshold of 10. Notably, at the lowest HMS threshold of 0, 4.4% of Black and 12.1% of White recipients were matched, a marked increase from the 1.8% and 6.6% matched under the prevailing allocation scheme. Furthermore, our allocation algorithms yielded similar or improved survival rates, as illustrated by Kaplan-Meier (KM) curves, and enhanced survival prediction accuracy, evidenced by C-indices and Integrated Brier Scores.

Improving Access to HLA-Matched Kidney Transplants for African American Patients

Introduction

Kidney transplants fail more often in Black than in non-Black (White, non-Black Hispanic, and Asian) recipients. We used the estimated physicochemical immunogenicity for polymorphic amino acids of donor/recipient HLAs to select weakly immunogenic kidney transplants for Black vs. White or non-Black patients.

Methods

OPTN data for 65,040 donor/recipient pairs over a 20-year period were used to calculate the individual physicochemical immunogenicity by hydrophobic, electrostatic and amino acid mismatch scores (HMS, EMS, AMS) and graft-survival outcomes for Black vs. White or vs. non-Black recipients, using Kaplan-Meier survival and Cox regression analyses. Simulations for re-matching recipients with donors were based on race-adjusted HMS thresholds with clinically achievable allocations.

Results

The retrospective median kidney graft survival was 12.0 years in Black vs. 18.6 years in White (6.6-year difference; p>0.001) and 18.4 years in non-Black (6.4-year difference; p>0.01) recipients. Only 0.7% of Blacks received transplants matched at HLA-A/B/DR/DQ (HMS=0) vs. 8.1% in Whites (p<0.001). Among fully matched Blacks (HMS=0), graft survival was 16.1-years and in well-matched Blacks (HMS ≤ 3.0) it was 14.0-years. Whites had 21.6-years survival at HMS ≤ 3.0 and 18.7-years at HMS ≤ 7.0 whereas non-Blacks had 22.0-year at HMS ≤ 3.0 and 18.7-year at HMS ≤ 7.0, confirming that higher HMS thresholds produced excellent survival. Simulation of ABO-compatible donor-recipient pairs using race-adjusted HMS thresholds identified weakly immunogenic matches at HMS=0 for 6.1% Blacks and 18.0% at HMS ≤ 3.0. Despite prioritizing Black patients, non-Black patients could be matched at the same level as in current allocation (47.0% vs 56.5%, at HMS ≤ 7.0).

Conclusions

Race-adjusted HMS (EMS, AMS)-based allocation increased the number of weakly immunogenic donors for Black patients, while still providing excellent options for non-Black recipients.

Low Hydrophobic Mismatch Scores Calculated for HLA-A/B/DR/DQ Loci Improve Kidney Allograft Survival

We evaluated the impact of human leukocyte antigen (HLA) disparity (immunogenicity; IM) on long-term kidney allograft survival. The IM was quantified based on physicochemical properties of the polymorphic linear donor/recipient HLA amino acids (the Cambridge algorithm) as a hydrophobic, electrostatic, amino acid mismatch scores (HMS\AMS\EMS) or eplet mismatch (EpMM) load. High-resolution HLA-A/B/DRB1/DQB1 types were imputed to calculate HMS for primary/re-transplant recipients of deceased donor transplants. The multiple Cox regression showed the association of HMS with graft survival and other confounders. The HMS integer 0–10 scale showed the most survival benefit between HMS 0 and 3. The Kaplan–Meier analysis showed that: the HMS=0 group had 18.1-year median graft survival, a 5-year benefit over HMS>0 group; HMS ≤ 3.0 had 16.7-year graft survival, a 3.8-year better than HMS>3.0 group; and, HMS ≤ 7.8 had 14.3-year grafts survival, a 1.8-year improvement over HMS>7.8 group. Stratification based on EMS, AMS or EpMM produced similar results. Additionally, the importance of HLA-DR with/without -DQ IM for graft survival was shown. In our simulation of 1,000 random donor/recipient pairs, 75% with HMS>3.0 were re-matched into HMS ≤ 3.0 and the remaining 25% into HMS≥7.8: after re-matching, the 13.5 years graft survival would increase to 16.3 years. This approach matches donors to recipients with low/medium IM donors thus preventing transplants with high IM donors.

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.

Quantity and Quality Reconstitution of NKG2A+ Natural Killer Cells Are Associated with Graft-versus-Host Disease after Allogeneic Hematopoietic Cell Transplantation

The immune mechanism underlying graft-versus-host disease (GVHD) after allogeneic stem cell transplantation (HSCT) remains unclear. Natural killer (NK) cells play a crucial role in mediating pathogen-specific immunity and are the first donor-derived lymphocytes reconstituted post-HSCT. However, NK cells vary at different stages after HSCT. Here, we found that the absolute NKG2A⁺ subset cell counts and the percentages of NKG2A⁺ among NK cells were significantly reduced in GVHD patients after HSCT compared with those from non-GVHD patients. Moreover, the reduction in NKG2A⁺ NK cells in post-HSCT GVHD patients was ascribed to increased apoptosis and a decreased proliferation capacity while retaining a strong graft-versus-leukemia effect. In vitro assays showed that co-culture of T cells with NKG2A⁺ NK cells significantly reduced IFN-γ secretion by T cells and increased IL-4 secretion. Moreover, the CD25 expression level was decreased, whereas the number of cells with the CD4⁺CD25⁺FOXP3⁺ phenotype was increased. In addition, the NKG2A⁺ NK cells induced T cell apoptosis and decreased T cell proliferation during the co-culture process. Importantly, NKG2A⁺ NK cells mainly regulated activated but not resting T cells. In vivo assays showed that the serologic IL-10 level was evidently lower in GVHD than in non-GVHD patients, whereas the IL-1β, IFN-γ, and tumor necrosis factor-α levels were higher in GVHD patients. Furthermore, the NKG2A⁺ NK cell ratio from GVHD patients was markedly increased by the presence of exogenous IL-10 but not by other cytokines. In contrast, the NKG2A⁺ cell ratio from non-GVHD patients was not increased by IL-10. Therefore, post-HSCT GVHD may be ascribed to the reduced induction of NKG2A⁺ NK cells by IL-10, which further overactivates T cells.

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.