Molecular disparity in human leukocyte antigens is associated with outcomes in haploidentical stem cell transplantation

Allogeneic hematopoietic stem cell transplantation (HSCT) is the only curative therapy for many advanced hematologic malignancies and nonmalignant hematologic disorders. With the success in the prophylaxis of graft-versus-host disease (GVHD) and graft rejection,^1,2  haploidentical related donors, who share 1 haplotype with recipients, have become widely accepted stem cell sources in clinical practice with comparable clinical outcomes with HLA-matched donor transplants.^3-5  The use of haploidentical donors significantly expanded the likelihood of finding a source of hematopoietic stem cells, especially in certain ethnic minority groups.^6-8 

It has been shown that HLA disparity has a negative effect on transplant outcomes of patients with HLA-mismatched unrelated donors. Compared with 8/8 matched unrelated donor transplant (MUD), a single mismatch at HLA-A, HLA-B, HLA-C, or HLA-DRB1 was associated with ∼10% decrease in 1-year survival and higher incidence of GVHD.⁹  In the setting of haploidentical HSCT (haplo-HSCT) performed with conventional GVHD prophylaxis, hyperacute reactions in graft-versus-host (GVH) and host-versus-graft (HVG) reactions occurred as a result of strong bidirectional alloreactivity.^10,11  However, with the use of posttransplant cyclophosphamide (PTCy) as GVHD prophylaxis, the influence of HLA mismatch in haplo-HSCT appears to be less significant and distinct from the impact seen in the transplantation with unrelated donors. Raiola et al recently studied a relatively large haplo-HSCT cohort and concluded that there is no correlation between the number of mismatched HLA antigens and clinical outcomes.¹²  However, although the immediate hyperacute reactions may be attenuated, higher disparities at particular HLA loci could perpetuate different alloreactive immune responses. The European Society for Blood and Marrow Transplantation reported recently that there is no influence of a cumulative number of mismatched HLA antigens on clinical outcomes in their haplo-HSCT cohort, yet an association between mismatched HLA-DRB1 and a higher risk of grade 2-4 GVHD was observed.¹³  Additionally, the molecular mechanisms behind T-cell alloreactivity are complex and presumably determined by permissibility and structural homology of HLA/peptide complexes. A recent study categorized the HLA mismatches into several of supertype groups, which were defined by anchor specificity of the presented peptide on HLA molecules. Compared with the supertype matched group, HLA-B supertype mismatch was associated with an increased risk of grade 2-4 acute GVHD (aGVHD) in mismatched unrelated donor transplants.¹⁴  Therefore, a comprehensive study of the structural and functional disparities at individual HLA loci, in both GVH and HVG directions, might shed light on minimizing risks and maximizing the benefit of alloreactive reactions in HLA-haploidentical transplantation.

HLAMatchmaker is a molecular matching algorithm that considers the functional components of epitopes exposed on HLA molecules. As the key determinants, eplets represent distinct configurations of amino acid polymorphisms that could elicit the immune response (Figure 1). HLAMatchmaker program is used to quantitatively determine the degree of mismatch by comparing the eplet repertories between donor and recipient.^15,16  It has been demonstrated that less ME load is associated with longer graft survival¹⁷  and reduced humoral sensitization in solid organ transplantation.¹⁸  Additionally, alloreactive T-cell clones that are specific to certain eplets were identified,^19-21  suggesting that HLAMatchmaker could likewise reflect compatibility at many polymorphic residues involved in T-cell receptor interaction.

Solomon et al recently reported a positive survival advantage associated with HLA-DPB1 nonpermissive mismatch defined by T-cell epitope (TCE) grouping in PTCy-based haplo-HSCT, which was attributed to a significant reduction in relapse rate without a parallel increase in the development of GVHD.²²  A recent study in unrelated donor transplantation demonstrated that HLA-DPB1 nonpermissive mismatches in GVH direction were associated with an increased risk of aGVHD.²³  Concurrently, a reduced risk for disease progression was observed with HLA-DPB1 mismatches.^24,25  Consequently, whether the directionality of nonpermissive mismatches affects the risk of relapse or GVHD in haplo-HSCT needs to be clarified.

In this study, we aimed to comprehensively assess the ability of HLAMatchmaker algorithm in predicting the risks or benefits associated with HLA molecular disparity in haplo-HSCT. In addition, we investigated the effect of HLA-PDB1 nonpermissive mismatches, either in the HVG or GVH direction, on clinical outcomes in a relatively large cohort of patients receiving haplo-HSCT at our institution.

Our cohort included 278 consecutively treated patients, 18 years of age or older, with hematologic malignancies who underwent first or second haplo-HSCT at The University of Texas MD Anderson Cancer Center between July 2009 and January 2015. The majority of patients received T cell–replete bone marrow graft (n = 247, 90%) and all received PTCy-based GVHD prophylaxis with tacrolimus and mycophenolate, as previously described by our group.²⁶  Donors were first-degree relatives who share 1 haplotype match with the recipients. All patients with a significant level (mean fluorescence intensity >2000) of donor-specific anti-HLA antibodies (DSA) received desensitization treatment before transplant per institutional protocol.²⁷  Comorbidities before transplant were evaluated using the Hematopoietic Cell Transplant-Comorbidity Index (HCT-CI).²⁸  Hematological malignancies were risk-stratified using the refined Disease Risk Index (DRI).²⁹  Clinical and laboratory data were collected via electronic medical records. All patients provided written informed consent for transplant following the Declaration of Helsinki. The study protocol was approved by the institutional review board of The University of Texas MD Anderson Cancer Center.

Patients eligible for this study had donor and recipient HLA typing performed at the HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DRB3/4/5, HLA-DQA, HLA-DQB1, and HLA-DPB1 loci using sequenced-based typing methods at high resolution.²³  For patients or donors who had no DQA typing in our record (∼15% of the cohort), inferred typing was assigned based on HLA-DR-DQ association at http://www.allelefrequencies.net/. Alleles reported with National Marrow Donor Program codes were deduced to the most common allele in the coded string. ME load at each HLA locus and total ME of class I or class II loci were quantified by using HLAMatchmaker module incorporated in HLA Fusion software v4.4. The analysis was performed in both GVH and HVG directions. The software identified antibody-verified eplets (ME-Ab) with informative HLA antibodies and theoretically predicted eplets (ME-ALL) based on the crystalized HLA molecule models.¹⁶  Eplet repertoires are listed in HLA Epitope Registry (http://www.epitopes.net/downloads.html). Inter-locus eplets were assigned as a given eplet identified in 2 or more loci within class I or class II molecules. While comparing eplets between donor and recipient to identify the mismatched eplets, the inter-locus eplets mismatched on 1 locus but matched on another locus were removed to truly reflect the molecular disparity.

HLA-DPB1 mismatches were classified into permissive and nonpermissive mismatches based on TCE algorithms (version 2.0) at the IPD-IMGT/HLA Web site (https://www.ebi.ac.uk/ipd/imgt/hla/dpb.html).³⁰  As previously described,²³  the vectors of HLA-DPB1 mismatches, either in the GVH or HVG direction, were assigned. Transplantations were therefore classified into 4 groups: (1) HLA-DPB1 matched; (2) HLA-DPB1 permissive mismatched; (3) nonpermissive mismatches in HVG direction; and (4) nonpermissive mismatches in GVH direction. Ligand-ligand model³¹  was used to predict natural killer (NK) cell alloreactivity based on the high-resolution HLA typing of the donor and recipient. Briefly, killer immunoglobulin receptor (KIR) ligand HLA-C and HLA-B molecules were grouped into 3 major categories (C1, C2, Bw4) based on the specific amino acid sequence that defines specific KIR ligand binding (https://www.ebi.ac.uk/ipd/kir/ligand.html). The NK cell alloreactivity in the GVH direction was assigned when the recipient lacked at least 1 of the HLA ligands that were present in the donor. This theoretical prediction of KIR benefit has been commonly used in haplo-HSCT studies.^32,33 

Clinical endpoints included overall survival (OS), progression-free survival (PFS), relapse, nonrelapse mortality (NRM), acute and chronic GVHD, neutrophil engraftment, and platelet engraftment. All outcomes were measured from the time of stem cell infusion. OS was based on death from any cause. PFS events included death or relapse. NRM was defined as death without a previous relapse. Patients without an event were censored at the time of last contact. Time to neutrophil engraftment was defined as the first of 3 consecutive days with an absolute neutrophil count exceeding 0.5 × 10³/µL. Time to platelet engraftment was defined as the first day of a platelet count >20 × 10³/µL without transfusion support for 7 consecutive days. Cumulative incidence with competing risks method was used to calculate aGVHD, chronic GVHD, NRM, relapse, and engraftment. For NRM, relapse was the competing risk, and for relapse, the competing risk was NRM. For acute and chronic GVHD, death without the event and relapse were the competing risks, whereas death without engraftment was considered a competing risk for engraftment. Associations between survival outcomes (OS and PFS) and each HLA locus ME number were determined using univariable and multivariable Cox proportional hazards regression models, whereas the impact of the number of HLA locus ME on relapse, NRM, and acute and chronic GVHD was determined using univariable and multivariable proportional subdistribution hazards regression model. HLA locus ME associated with P < .10 by the univariable analysis were included in the multivariable analysis along with other potential prognostic factors including sex, age, HCT-CI, DRI, transplant protocol, conditioning regimen intensity, donor relation, stem cell source, and ABO match. Given that ME-Ab and ME-ALL at the same locus are highly associated and overlapped, only 1 factor with the lowest P value was chosen and included in multivariable analysis. All variables of interest were tested for the proportional hazard assumption and interaction terms. All tests were 2-sided. The type 1 error rate was fixed at 0.05. No adjustments for multiple testing were made. Stata statistical software (SE 13, Stata Corp LP, College Station, TX) was used for statistical analyses.

The median age of this cohort was 48 years (range, 18-72 years). The clinical characteristics are summarized in Table 1. Acute myeloid leukemia/myelodysplastic syndrome was the main indication for haplo-HSCT (163 patients, 58.6%). One hundred and thirty-four patients had high- (40.7%) or very high-risk diseases (7.6%) based on the DRI. Myeloablative conditioning was used in approximately one-half of the patients. Fifty-one percent of the patients had HCT-CI score ≥3. CMV serostatus was nonreactive in 13% of the recipients. NK cell alloreactivity was present in 28% of transplants. HLA-DPB1 nonpermissive mismatches, either in GVH or HVG directions, were noted in one-quarter of the transplantations. In our cohort, the NRM at 1 and 3 years was 25.9% and 29.7% and relapse rate at 1 and 3 years was 24.1% and 27.4%, respectively. PFS and OS at 3 years were 42.9% and 47.4%.

Both ME-Ab and ME-ALL load were quantified on each locus in the HVG or GVH direction (supplemental Table 1). The median of the total number of class I ME and class II ME in HVG direction was 24 (range, 16-32) and 21 (range, 11-31), respectively. The median of the total number of class I ME and class II ME in the GVH direction was 26 (range, 17-35) and 23 (range, 14-32), respectively. No correlation was identified between the number of ME in HVG direction vs that in GVH direction (supplemental Figure 1A,C), suggesting that the ME load recognized from the donor perspective is different from that recognized from the recipient perspective. In contrast, the number of theoretical ME-ALL and ME-Ab at 1 specific locus was highly and positively correlated. For instance, the linear regression between ME-ALL and ME-Ab at HLA-A locus and class I were analyzed, the r² value was 0.78 and 0.80, respectively (supplemental Figure 1B,D). Similar to a recent study in solid organ transplantation,³⁴  no significant differences in clinical impact was recognized between antibody-verified ME and theoretical ME in the present study.

We observed a significant association between higher ME at HLA-DRB3/4/5 in HVG direction and an inferior OS on the univariable analysis (hazard ratio [HR], 1.03; 95% confidence interval [CI], 1.00-1.05; P = .025) (Table 2). After adjusting for potential confounders, the association remained significantly in multivariable analysis (HR, 1.03; 95% CI, 1.00-1.05; P = .038) (Table 2). ME load at HLA-A locus and class I in the HVG direction was associated with a better survival on the univariable analysis (Table 2). After adjusting for other patients, disease, and transplant characteristics, an increased HLA-A ME load in the HVG direction remained significantly associated with better OS with a 3% mortality reduction for each ME (HR, 0.97; 95% CI, 0.95-0.99; P = .003). For transplants with lower ME load at HLA-A locus, a 3-year OS rate was only 39%, significantly lower than 54.3% seen in transplants with a higher number of ME at HLA-A locus (Figure 2A). Consistently, the risk of mortality decreased significantly when the total number of ME from class I increased on the multivariable analysis (HR, 0.99; 95% CI, 0.97-0.99; P = .047) (Table 2). Because no association was found between OS and ME from HLA-B or HLA-C, the survival advantage from class I ME is likely attributed to ME from HLA-A locus. Other variables associated with worse OS were poorer HCT-Cl score, high- and very-high DRI, and older age at transplantation (Table 2). No significant influence on OS from ME at other individual locus and total class II loci was detected, indicating the alloimmunity elicited by the mismatches from each HLA locus could be different.

Further studies identified the association between better PFS and a higher load of ME at HLA-A and class I in the HVG direction as presented in Table 3 and Figure 2B. As a continuous variable, each class I ME was significantly associated with better PFS with 2% mortality or relapse reduction (HR, 0.98; 95% CI, 0.97-0.99; P = .039), and each ME in HLA-A was significantly associated better PFS with 4% mortality or relapse reduction on the multivariable analysis (HR, 0.96; 95% CI, 0.05-0.98; P < .001). The survival advantage of HLA-A ME load was primarily attributed to a reduced risk in relapse, which progressively decreased with an increase in the number of A-ME (subdistribution hazard ratio [SHR], 0.95; 95% CI, 0.92-0.98; P = .004) (Figure 2C; Table 4). The 3-year cumulative incidence of relapse was 21.6% in the higher A-ME group compared with 34.3% seen in the lower A-ME group (Figure 2C). As described previously, the antibody-verified ME was included in the theoretical version of ME and they were highly correlated. Both ME-Ab and ME-ALL at the HLA-A locus were significantly correlated with reduced risk of relapse in our analysis (Table 4). Neither HLA-A ME nor class I ME was significantly associated with NRM (Table 5), whereas an inferior NRM was significantly associated with increased ME at HLA-DRB1 (SHR, 1.04; 95% CI, 1.00-1.09; P = .048) but not with increased ME from class II (SHR, 1.01; 95% CI, 0.99-1.04; P = .053).

A previous study suggested that HLA-DPB1 disparity might be associated with relapse protection in haplo-HSCT with PTCy.²²  We further characterized the relapse protection by separately assessing the nonpermissive mismatches in HVG or GVH directions predicted by TCE grouping. A nonsignificant trend of relapse protection was found with HLA-DP nonpermissive in HVG direction (SHR, 0.54; 95% CI, 0.22-1.33; P = .181) on univariable analysis. In contrast, we observed a significant reduction in the incidence of relapse associated with HLA-DPB1 nonpermissive mismatch in the GVH direction (SHR, 0.25; 95% CI, 0.09-0.76; P = .014). The 3-year cumulative incidence of relapse was 15.3% in the nonpermissive GVH group, compared with 31.3% in HLA-DPB1 matched group (Figure 3A). We further investigated the role of nonpermissive mismatches in either HVG or GVH vectors on multivariable analysis and confirmed that only GVH nonpermissive mismatch was significantly associated with relapse protection (SHR, 0.19; 95% CI, 0.06-0.59; P = .004) (Table 4) without impact on OS (Table 2) or risk of GVHD (Figure 3B). Unlike observed in transplant with MUD donor,²³  no significantly higher risk of grade 3-4 aGVHD was observed in nonpermissive GVH group (SHR, 1.00; 95% CI, 0.27-3.74; P = .995) (Table 6) with univariable analysis. On multivariable analysis, the HLA-DPB1 ME in GVH direction was not significantly correlated with a higher risk of aGVHD or grade 3-4 aGVHD (SHR, 0.78; 95% CI, 0.65-1.02; P = .112).

On univariable analysis, no significant association was found between a higher load of ME from an individual locus and risk of aGVHD (Table 7). The impact on grade 3-4 aGVHD was further investigated and the significance was noted on multivariable analysis with a significantly higher incidence of grade 3-4 aGVHD in the group with higher HLA-B ME (SHR, 1.05; 95% CI, 1.01-1.12; P = .036) and class I ME in GVH direction (SHR, 1.09; 95% CI, 1.00-1.18; P = .041) (Table 6). As shown in Figure 4A, the 100-day cumulative incidence of grade 3-4 aGVHD in a higher B-ME group was 10% compared with 7% in the group with lower B-ME.

The effect of HLA disparity on engraftment and chronic GVHD was also evaluated. Different from a recent study on HLA mismatch in haplo-HSCT,³⁵  we did not find any association between ME load and neutrophil and platelet engraftment. The association between ME load and chronic GVHD was not identified. No significant impact from NK alloreactivity was noted on OS (HR, 1.13; 95% CI, 0.79-1.61; P = .499), PFS (HR, 1.11; 95% CI, 0.79-1.57; P = .520), or relapse (SHR, 1.20; 95% CI, 0.74-1.94; P = .441).

It has been well accepted that a greater HLA disparity is associated with worse clinical outcomes and a higher risk of GVHD in HSCT with mismatched donors receiving conventional GVHD prophylaxis.^36,37  However, the immune system seems to respond differently in the setting of haplo-HSCT with GVHD prophylaxis using PTCy regimen. As haplo-HSCT is increasingly used for patients requiring transplantation, it is essential to understand the immunogenicity to alloantigen on the mismatched graft from haploidentical donors. In the present study, with patient cohort receiving PTCy-based GVHD prevention and bone marrow graft, we demonstrated that HLA disparity quantified at molecular mismatch level, and not by the cumulative number of mismatched antigens, could be relevant to the clinical outcome of patients receiving haploidentical hematopoietic cell transplantation.

With the rapid progress in molecular typing and protein modeling, the degree of HLA disparity has recently been assessed at the molecular level in solid organ transplantation to precisely determine HLA alloimmune risk.^38,39  Different algorithms of molecular mismatching in predicting immunogenicity have been developed with different focuses, such as the number of mismatched amino acids or physiochemical properties of the amino acid substitution.⁴⁰  Initially, HLAMatchmaker quantified mismatched epitopes and predicted humoral response along with antibody-mediated rejection. The revised HLAMatchmaker algorithm considers the structural proximity that is accessible by antibodies and includes all 3-dimensional patch of linear or discontinuous polymorphic amino acids within 3 Å radium. It defines both antibody-verified eplets and theoretical eplets based on stereochemical modeling of HLA molecule and antigen-antibody complex. In the present study, we demonstrated that the revised HLAMatchmaker could be used for the assessment of histocompatibility in haplo-HSCT at the molecular level.

Here, we have found that ME load at different loci has different effects on the clinical outcomes of haplo-HSCT. Patients with higher ME load at class I loci and HLA-A locus in the HVG direction was found to be associated with improved OS. Our study also suggested that this survival benefit of A-ME may primarily result from a decreased risk of relapse. Kasamon et al studied a cohort of haplo-HSCT patients receiving bone marrow graft and high-dose PTCy and demonstrated that presence of higher allelic level mismatches (3/4 vs 1/2) was associated with a nonsignificant trend to less relapse, and the higher number of class I mismatches (2/3 vs 0/1) was significantly associated with lower relapse rate.⁴¹  Interestingly, further analysis indicated that the significant protective effect was mainly associated with mismatches in the HVG direction. The authors were unable to disentangle the effect from HVG versus that from GVH direction because their mismatches at the antigen level were tightly linked from both vectors. In our study, using ME instead of mismatched antigens/alleles, we showed that the ME at the same locus but from different directions are not correlated, which allowed us to separately assess their effect (supplemental Figure 1). Surprisingly, we found that a lower risk of relapse was significantly associated with HLA disparity in the HVG but not in the GVH direction, indicating that alloreactive eplets present in the HVG vector could be involved in a subsequent immune response. Several clinical studies demonstrated the association between the rejection of donor bone marrow and the regression of hematological tumor without GVHD.^42-44  In addition, several studies in animal models (reviewed by Li and Sykes⁴⁵ ) suggested that HVG alloresponse possibly induces an antitumor response by enhancing the pro-inflammatory environment and leads to subsequent activation of particular cell subsets.^46,47  Fleischhauer et al studied the effect of HLA-DPB1 nonpermissive mismatches on clinical outcomes in bidirectional vs unidirectional groups.⁴⁸  The proliferative response in mixed lymphocyte reaction assay was significantly higher in the presence of nonpermissive mismatches. Interestingly, there was no significant difference in the levels of response invoked by alloreactivity in the HVG and GVH direction, whereas a significantly higher response was observed in the bidirectional group. They postulated that in response to allogeneic HLA on the responder cells (graft), the stimulator cells (host) release cytokines that would enhance the proliferation of responder cells.⁴⁸  In this study, the exact mechanism underlying relapse protection by the HVG alloreactivity is unclear and remains to be explored in future studies. Different from what we observed here, a recent study using peripheral blood graft reported that class II ME in the GVH direction was associated with decreased relapse rate in haplo-HSCT.³⁵  This inconsistency may be attributed to different stem cell sources, transplant protocol, disease status, DSA, and different ethnic cohorts. Recent studies in haplo-HSCT with PTCy-based prophylaxis suggested that HLA disparity is associated with reduced relapse without concurrent increased risk of GVHD, likely because of the distinctive ability of PTCy to manipulate allocative T cells and regulatory T cells.⁴⁹  In our study, the administration of PTCy appears not to completely mitigate the GVHD in response to allogeneic HLA disparity, as the ME load at HLA-B locus and class I in the GVH direction were correlated with a higher risk of developing grade 3-4 aGVHD. The distinctive allorecognition responses to HLA disparity from different loci, may be dictated by the different expression levels of HLA molecules, or specific reactivity of potential effector cells.^41,50  The modulation of alloreactive T cells by PTCy adds an extra layer of complexity. Our current observational study is unable to clarify this complex association and selecting the optimal donor could be difficult when both favorable and undesirable factors exist concomitantly. However, our results suggest that donor selection could be optimized based on the specific impact from ME at 1 HLA locus, in particular, the donor with a higher A-ME load might be preferred over other donors for the patients with a higher risk of relapse.

The role of HLA-DPB1 disparity in HSCT with MUD has been reported by several studies.^23,24,51,52  HLA-DPB1 mismatches were classified into nonpermissive and permissive mismatches, the latter was defined as the net result of elicited alloreactivity is toward to graft versus leukemia (GVL) with clinically tolerable GVHD effect.^53,54  One model predicted the permissiveness using the TCE grouping⁵⁵  and another model was based on specific allele expression levels of the HLA-DPB1 allele.⁵⁶  In HSCT with unrelated donors, nonpermissive mismatches are generally associated with increased risk of GVHD that overweighs the GVL benefit and results in an inferior survival.^24,57  The study on HSCT with MUD advocated that HLA-DPB1 disparity in HVG vector, in addition to the disparity in GVH direction, might trigger donor effector cells to exert GVHD and GVL effect.⁵³  In our haplo-HSCT cohort, only GVL but not GVHD effect was observed with GVH nonpermissive mismatches of HLA-DPB1, perhaps because of the use of PTCy. Neither HLA-DPB1 nonpermissive mismatch in HVG direction nor ME in HLA-DPB1 was associated with decreased risk of relapse. The HLA-DPB1 disparity predicted by TCE model seems to be more clinically relevant compared with ME load, possibly because TCE classification considers permissiveness based on the functional “thymic selection” with proven T-cell alloreactivity patterns.⁵³ 

Our study has several limitations including its retrospective nature and a large number of analyses on a relatively small number of patients, which may result in the false-positive potential resulting from lower statistical power. Although it is well accepted that each eplet does not contribute to immunogenicity equally, we were able to study the disparity only by ME load but not by individual mismatched eplets. Studies with a larger cohort of patients are warranted to fully characterize the eplets with a specific impact in haplo-HSCT. Additionally, inferred DQA was used in about 15% of our cohort, which may not be completely accurate. Although the prevalence of DSA in our cohort was relatively low (9.6%), we were unable to exclude the influence of DSAs on transplant outcomes because the majority of patients with DSAs were desensitized before the transplant. In line with a haplo-HSCT study on patients with acute leukemia,⁵⁸  we found no association between clinical outcomes and KIR mismatches assessed by the ligand–ligand mismatched model. Compared with the gene–gene model that requires KIR gene typing on both recipient and donor, it has been shown that the model we used predicts a considerably lower number of KIR alloreactivity,³¹  which might underpower the KIR benefit given the low number of occurrences. Despite these limitations, to the best of our knowledge, this study is the largest study of the clinical impact of HLA disparity measured at the molecular level in haplo-HSCT.

In conclusion, molecular HLA disparity at different loci appears to have a different effect on clinical outcomes in haplo-HSCT. Higher ME load at HLA-A locus was associated with improved OS primarily because of a reduction in relapse rate posttransplant. HLA-DPB1 nonpermissive mismatch in GVH direction demonstrated also significant relapse protection without increasing the risk of GVHD. Molecular mismatch analysis could be performed routinely in the future in haploidentical donor transplants to assist in risk stratification and optimal donor selection.

The authors thank Rene Duquesnoy for helpful discussion and preparation of the manuscript and Kevin Harrell and Jar-How Lee from Thermo Fisher Scientific for their help in eplet analysis.

Contribution: J.Z., S.O.C., and K.C. designed the study, contributed to data collection and interpretation, and manuscript writing; J.Z. wrote the initial draft of the manuscript; P.K. and M.Y. contributed statistical analysis, interpretation of statistical results, and reviewed and approved the manuscript; Y.C. contributed to data collection and data analysis and reviewed and approved the manuscript; G.R. contributed to data collection and reviewed and approved the manuscript; S.S., P.K., and D.P. contributed to data interpretation and reviewed and approved the manuscript; and S.O.C. and R.E.C. contributed to treatment of patients, reviewed, edited, and approved the final version of the manuscript.

Conflict-of-interest disclosure: The authors declare no competing financial interests.

The current affiliation for S.O.C. is Division of Hematology/Oncology, Department of Medicine, Chao Family Comprehensive Cancer Center, University of California, Irvine, Orange, CA.

Correspondence: Jun Zou, Department of Laboratory Medicine, HLA Laboratory, University of Texas MD Anderson Cancer Center, 6565 MD Anderson Blvd, Houston, TX 77030; e-mail: jzou@mdanderson.org; or Stefan O. Ciurea, Division of Hematology/Oncology, Department of Medicine, Chao Family Comprehensive Cancer Center, University of California, Irvine, 101 The City Dr South, ZOT 4061, Orange, CA 92868; e-mail: sciurea@hs.uci.edu.