ORIGINAL RESEARCH article
Exploratory Study of Predicted Indirectly ReCognizable HLA Epitopes in Mismatched Hematopoietic Cell Transplantations
- 1Laboratory of Translational Immunology, University Medical Center Utrecht, Utrecht, Netherlands
- 2Department of Theoretical Biology and Bioinformatics, Utrecht University, Utrecht, Netherlands
- 3Matchis Foundation, Leiden, Netherlands
- 4Department of Immunohematology and Blood Transfusion, Leiden University Medical Center, Leiden, Netherlands
- 5HOVON Data Center, Department of Hematology, Erasmus MC Cancer Institute, Rotterdam, Netherlands
- 6Department of Hematology, Cancer Center Amsterdam, Amsterdam UMC, VU Medical Center, Amsterdam, Netherlands
- 7Department of Hematology, Academic Medical Center, University of Amsterdam, Amsterdam, Netherlands
- 8Department of Haematology, University of Groningen, University Medical Centre Groningen, Groningen, Netherlands
- 9Department of Hematology, Leiden University Medical Center, Leiden, Netherlands
- 10Department of Hematology, Radboud University Medical Center, Nijmegen, Netherlands
- 11Department of Hematology, Erasmus Medical Center-Daniel Den Hoed Cancer Center, Rotterdam, Netherlands
- 12Department of Hematology, University Medical Center Utrecht, Utrecht, Netherlands
HLA-mismatches in hematopoietic stem-cell transplantation are associated with an impaired overall survival (OS). The aim of this study is to explore whether the Predicted Indirectly ReCognizable HLA-Epitopes (PIRCHE) algorithm can be used to identify HLA-mismatches that are related to an impaired transplant outcome. PIRCHE are computationally predicted peptides derived from the patient's mismatched-HLA molecules that can be presented by donor-patient shared HLA. We retrospectively scored PIRCHE numbers either presented on HLA class-I (PIRCHE-I) or class-II (PIRCHE-II) for a Dutch multicenter cohort of 103 patients who received a single HLA-mismatched (9/10) unrelated donor transplant in an early phase of their disease. These patients were divided into low and high PIRCHE-I and PIRCHE-II groups, based on their PIRCHE scores, and compared using multivariate statistical analysis methods. The high PIRCHE-II group had a significantly impaired OS compared to the low PIRCHE-II group and the 10/10 reference group (HR: 1.86, 95%-CI: 1.02–3.40; and HR: 2.65, 95%-CI: 1.53–4.60, respectively). Overall, PIRCHE-II seem to have a more prominent effect on OS than PIRCHE-I. This impaired OS is probably due to an increased risk for severe acute graft-vs.-host disease. These data suggest that high PIRCHE-II scores may be used to identify non-permissible HLA mismatches within single HLA-mismatched hematopoietic stem-cell transplantations.
HLA mismatching is an important factor in the outcome of allogeneic hematopoietic stem-cell transplantation (HSCT), leading to an impaired overall survival (OS) (1). Therefore, donor and recipient are currently preferably matched for HLA-A, -B, -C, -DRB1, -DQB1 loci (10/10) (2). The detrimental effect of HLA mismatching on OS after HSCT is most pronounced in patients who are transplanted in the early phase of their disease (3), whereas HLA mismatching plays a less prominent role in patients with an advanced disease status before HSCT (3).
HLA-matched (10/10) unrelated donors are available for 23–67% of the patients, depending on the ethnicity (4–6). This underscores the urgent medical need, also in the era of novel transplantation concepts, to discriminate between detrimental HLA mismatches (non-permissible) and well-tolerated HLA (permissible) mismatches. In order to estimate the permissibility of HLA mismatches before transplantation, in silico methods can be used to predict whether donor T cells are able to recognize these HLA mismatches (7). We developed recently one of these methods, the so-called “Predicted Indirectly ReCognizable HLA Epitopes” (PIRCHE) algorithm (7–9). The PIRCHE algorithm identifies mismatched HLA-derived epitopes that can potentially be presented on HLA class-I (designated as PIRCHE-I) or HLA class-II molecules (designated as PIRCHE-II). Theoretically, PIRCHE-I lead to CD8+ T-cell responses and PIRCHE-II to CD4+ T-cell responses. The PIRCHE model has shown correlations with transplant outcome in HLA-C and HLA-DPB1 mismatched unrelated donor (MUD) HSCT (8, 9), in HLA-mismatched cord blood transplantation (10), and with de novo HLA antibody formation in organ transplantation and pregnancy (11–14). In the current retrospective explorative study we multivariately investigate the role of the PIRCHE algorithm in identifying non-permissible HLA mismatches in 9/10-matched HSCT transplantations.
Materials and Methods
This study included 685 patients who were transplanted for malignant diseases with MUDs at 8 Dutch transplant centers between 1989 and 2011. A total of 249 patients (36%) were transplanted with a 9/10 match, and 436 (64%) were transplanted with a 10/10 match. Clinical data were collected according to EBMT guidelines (accessible via: https://www.ebmt.org/patient-privacy-statement). All subjects gave informed consent to use their clinical data according to the JACIE guidelines. Since data was collected using the EBMT/JACIE guidelines, additional local ethical approval for conducting the current study was not required in accordance with the institutional requirements and national legislation. The disease status was defined for each individual disease category. Early-stage disease was defined as acute leukemia in first complete remission, chronic myeloid leukemia (CML) in first chronic phase, and myelodysplastic syndrome (MDS), non-Hogdkin lymphoma (NHL), or multiple myeloma (MM) untreated or in first complete remission. Intermediate-stage disease was defined as acute leukemia in second complete remission, CML in second chronic or accelerated phase, MDS in second or partial remission, and NHL or MM in partial remission, second complete remission, or stable disease. Advance-stage disease was defined as acute leukemia, CML, lymphoma, MDS, or MM in later disease stages as defined in de early or intermediate-stage disease. Unambiguous high-resolution HLA-A, -B, -C, -DRB1, and -DQB1 typing data were available for all donor-recipient pairs. HLA typing was performed using sequence-based typing. Supplementary Table 1 shows the population characteristics according to the match status for the complete cohort.
For the 9/10-MUD group, PIRCHE were identified for each donor-recipient pair, as described previously (8, 9). Briefly, NetChop 3.1 (15, 16) (predicting the peptide generation for HLA class-I presentation), and NetMHCpan 2.4 (17, 18) (predicting peptide binding affinity to HLA class-I molecules) were used to identify PIRCHE-I. NetMHCIIPan 3.0 (19, 20) (predicting peptide binding affinity to HLA class-II molecules) was used to identify PIRCHE-II. Only peptides with high binding affinities, i.e., predicted IC50 ≤500 nM (for PIRCHE-I) or ≤1,000 nM (for PIRCHE-II), were accepted as relevant binders.
The primary endpoint used in this study was OS, defined as time from HSCT to death due to any cause. To explain the cause of differences in the primary endpoint, several secondary endpoints were evaluated: disease-free survival (DFS; defined as survival without recurrence of the primary malignancy), non-relapse mortality (NRM; defined as mortality without previous progression of the primary malignancy), acute graft-vs.-host-disease (GVHD), chronic GVHD, and progression.
The association of PIRCHE on OS or DFS was studied with Cox proportional hazard models. Stratification was used to account for heterogeneity of diagnosis. A gamma-frailty term was used to adjust for center effects. Competing risk analyses were performed for NRM (treating progression as a competing risk), acute GVHD, chronic GVHD, and progression. Grading of acute GVHD was defined according to international consensus criteria (21), treating progression and NRM as competing risks. For progression of the primary malignancy, NRM was treated as a competing risk. For all statistical models, the covariates listed in Supplementary Table 1 were evaluated for inclusion in the multivariate models. All statistical analyses were based on multivariate models, except for the Kaplan-Meier analysis (Figure 1A). The Kaplan-Meier analysis was performed univariately to visualize the effect of PIRCHE on OS. Statistical analyses were performed with SPSS 21 (IBM SPSS Software) and with R version 3.4.1 (R Foundation for Statistical Computing). A p-value below 0.05 was considered statistically significant.
Figure 1. Early-stage disease patients transplanted with a high PIRCHE-II mismatch have an impaired 5-year OS compared to early-stage disease 10/10-matched transplantations (A) and a higher risk for severe aGVHD (B). (A) High PIRCHE-II had an impaired OS (36%) compared to the 10/10-MUD group (56%) and the low PIRCHE-II group (52%). (B) Patients with a high PIRCHE-II mismatch had a significantly increased risk of NRM and acute GVHD compared to 10/10-matched transplantations, whereas the risk for disease progression was not affected. Multivariate models included for NRM: Time to HSCT, KIR ligand status, patient CMV status, conditioning regimen intensity, patient age at transplantation, and donor age; for acute GVHD: donor age (for II-IV) and patient age (for III-IV); for chronic GVHD: aGVHD II-IV, conditioning regimen intensity, patient age at transplantation, stem cell source, and ATG; for DFS: patient age at transplantation; for progression: patient CMV status, conditioning regimen intensity, and ATG. 10/10-MUD: n = 212 patients; low PIRCHE-II group: n = 70 patients; high PIRCHE-II group: n = 33 patients. The number of outcome events in different groups: 10/10: n = 94; PIRCHE-II low: n = 34; PIRCHE-II high: n = 23. OS, overall survival; HR, hazard ratio; 95%-CI, 95%-confidence interval of hazard ratio. PIRCHE-II low, 0–13 PIRCHE-II; PIRCHE-II high, >13 PIRCHE-II; NRM, non-relapse, mortality; GVHD, graft-vs.-host-disease; DFS, disease-free survival.
The total 9/10-MUD group had significantly impaired OS compared to the 10/10-MUD reference group (Table 1). Since the impact of HLA mismatching on OS highly depends on the disease status before HSCT (3), we subsequently stratified the 9/10-MUD and 10/10-MUD groups according their disease status. Within the early-stage disease patients, 9/10-MUD had significantly impaired OS compared to the 10/10-MUD group [Hazard Ratio (HR): 1.61, 95% confidence interval (95%-CI): 1.09–2.37, p = 0.02], whereas for the intermediate and late-stage disease patients no differences in OS were observed between 9/10-MUD and 10/10-MUD (Table 1).
Since our data confirm that the effect of HLA mismatching on OS is more prominent in early-stage disease patients, we next investigated whether PIRCHE scores may aid to identify non-permissible HLA mismatches within this patient group. The baseline characteristics of the early-stage disease patients according to the match grade are listed in Table 2. For the 9/10-MUD group, receiver operating characteristic (ROC) curves with regard to OS were generated for different PIRCHE-I and PIRCHE-II cutoffs (Supplementary Table 2); the cutoffs that yielded the highest area under the curve were used to define low PIRCHE-I/-II and high PIRCHE-I/-II [PIRCHE-I: 0–1 (low) vs. >1 (high); PIRCHE-II: 0–13 (low) vs. >13 (high)]. Statistical comparisons of the co-variates between the 10/10-MUD group and the low or high PIRCHE-I/-II groups showed that only ATG was significantly different between 10/10-MUD and high PIRCHE-I and 10/10-MUD and low PIRCHE-II (data not shown). Both the high PIRCHE-I and the high PIRCHE-II group had a significantly impaired OS compared to the 10/10-MUD early-stage disease group, while only the PIRCHE-II high group had an impaired OS also compared to the PIRCHE-II low group (Table 3). To visualize the effect of PIRCHE-II on OS, univariate Kaplan-Meier analyses were performed (Figure 1A). The high PIRCHE-II group had an impaired 5-year OS (36%) compared to the 10/10-MUD group (56%) and the low PIRCHE-II group (52%). Cumulatively, these data suggest that high PIRCHE-II mismatches are non-permissible HLA mismatches.
Table 3. Hazard ratios of OS for PIRCHE groups compared to 10/10 or to low PIRCHE within the early-stage disease patients.
Next, we investigated the potential underlying causes of the impaired OS within the high PIRCHE-II group by studying the secondary endpoints (Figure 1B). Patients transplanted with a non-permissible PIRCHE-II had an increased NRM risk compared to the 10/10-MUD group. Therefore, we investigated potential causes of NRM, such as acute GVHD (22). In our cohort, ~62% of the patients who died of NRM had mild to severe aGVHD. Patients transplanted with a high PIRCHE-II mismatch had a significantly increased risk of severe acute GVHD (grade III-IV) compared to the 10/10-MUD group. In contrast, the risk for chronic GVHD did not differ from the 10/10-MUD group. Since the OS is not only influenced by GVHD, but also by disease progression, we further analyzed the DFS and disease progression risk. Patients transplanted with high PIRCHE-II had an impaired DFS compared to the 10/10-MUD group, whereas no difference in disease progression risk was observed between both groups, suggesting that PIRCHE-II does not affect disease progression. No differences in NRM, GVHD, DFS, and progression risk were observed between the 10/10-MUD group and the low PIRCHE-II group (data not shown), suggesting that low PIRCHE-II scores do not impact alloreactivity and the anti-tumor reactivity.
Several studies have previously shown the association of PIRCHE with alloreactive responses after transplantation. A high number of PIRCHE-II was associated with HLA antibody formation after pregnancy (12) and in different organ transplantation settings (11, 13, 14). Additional studies in cord blood transplantation (10) and in HLA-C and HLA-DPB1 mismatched HSCT (8, 9) have shown that PIRCHE is related with transplant outcome. In the current explorative study we aimed to investigate whether the PIRCHE algorithm can identify non-permissible HLA mismatches in 9/10-matched HSCT transplantations. Since the impact of HLA mismatching on OS is the most prominent in patient transplanted within the early-stage of their disease (3), our analyses were performed within this patient group. Our data suggest that PIRCHE is associated with an impaired overall survival in 9/10-MUD HSCT.
Previous studies in cord blood transplantations have suggested that PIRCHE-II has a higher impact on alloreactivity than PIRCHE-I (10). In pediatric patients, high PIRCHE-I scores were associated with a lower relapse incidence in these patients, whereas a trend for high PIRCHE-II toward an increased GVHD risk was observed (10). Cumulatively, these data suggest that PIRCHE-II might play a more prominent role in GVHD and consequently NRM than PIRCHE-I. Also the current study suggests a more prominent role for PIRCHE-II than for PIRCHE-I. However, we also observed a high association between PIRCHE-I and PIRCHE-II; the majority of the patients with high PIRCHE-II also had high PIRCHE-I (data not shown). When analyzing patients who had both high PIRCHE-I and high PIRCHE-II, this group also had an impaired overall survival compared to the 10/10-MUD HSCT (HR: 2.78; 95%-CI: 1.59–4.86; p = 0.0003). This observation suggests that we cannot rule out the effect of PIRCHE-I numbers on the impaired overall survival in the high PIRCHE-II group. Therefore, further validation studies in larger cohorts are warranted to investigate whether PIRCHE-II indeed play a more prominent role in alloreactivity than PIRCHE-I in HLA-mismatched transplantation.
Although our study suggests a more prominent role for PIRCHE-II in alloreactivity than PIRCHE-II, the exact mechanism behind this observation is speculative. Theoretically, high PIRCHE-II numbers are associated with a high level of CD4+ T cell alloreactivity. PIRCHE-II specific CD4+ T cells may impact transplant outcome via two ways [reviewed in (23)]: CD4+ T-helper cells may, after indirect recognition of mismatched HLA, provide help to mismatched HLA-specific CD8+ cytotoxic T cells, and consequently induce alloreactive responses via these CD8+ T cells. On the other hand, PIRCHE-II specific CD4+ T cells may play a role in the formation of HLA-specific antibodies by providing B cell help. A previous study did not show a correlation between engraftment, acute GvHD, and survival after HSCT and the number of mismatched B-cell epitopes between donor and recipient, as predicted by the HLAMatchmaker algorithm (24). Since the previous study did not show a correlation between the number of mismatched B-cell epitopes and survival after mismatched HSCT, it is likely that PIRCHE-II impact overall survival after 9/10-MUD HSCT via providing help to mismatched HLA-specific CD8+ cytotoxic T cells rather than via providing help to B cells.
In the current explorative study we show that PIRCHE is associated with an impaired overall survival in 9/10-MUD HSCT, which was the primary endpoint. The heterogeneity and the small size of the cohort may have biased our results. Consequently, especially the results obtained in the comparison between the low and high PIRCHE group should be treated with caution, as these groups are rather small and the number of events limited. Although a difference in overall survival was observed between low PIRCHE-II and high PIRCHE-II, no differences in NRM, aGVHD, cGVHD, and relapse (the secondary endpoints) were observed between the low and high PIRCHE-II group (data not shown), which may be due to the small group size. Therefore, large (prospective) studies consisting of homogeneous study populations are required to further validate our observations. Moreover, these cohorts are also required to determine the most optimal PIRCHE-II threshold for clinical practice, as the current explorative study does not aim for identifying clinically relevant PIRCHE cutoffs.
In our study we investigated whether the PIRCHE algorithm might identify non-permissible HLA-A, -B, -C, -DRB1, or DQB1 mismatches. The locus mismatches of the early-stage disease 9/10-MUD group within our study cohort are listed in Supplementary Table 3. Individual HLA locus mismatches may have a differential impact on OS (3), which may be due to a differential number of PIRCHE for different locus mismatches. Since our study cohort is rather small, we were unable to investigate the impact of individual HLA locus mismatches on transplant outcome and its' relation with PIRCHE in our cohort.
Although we exclusively focused on HLA-A, -B, -C, -DRB1, and DQB1 mismatches in the current study, compelling evidence indicate that HLA-DPB1 mismatches may impact overall survival after HSCT as well (25, 26). Previously, we have shown in patients transplanted with a 10/10-MUD with a HLA-DPB1 mismatch that HLA-DPB1 mismatches that resulted in a positive PIRCHE score had a higher risk for acute GVHD than those HLA-DPB1 mismatches with a PIRCHE score of zero (8). These data indicate that ideally HLA-DPB1-derived PIRCHE should also be included our analyses. However, HLA-DPB1 typing is lacking for 11% of the study population. Moreover, as the aim of the current explorative study is to study the role of PIRCHE in 9/10 MUD-HSCT, we cannot confirm the previous studies performed in 10/10-MUD HSCT. Nevertheless, our study outcome may have been biased by the occurrence of HLA-DPB1-mismatches among the 9/10-MUD group. The occurrence of HLA-DPB1 mismatches did not significantly differ between the 10/10-MUD and 9/10-MUD group of the early-stage disease patients (p = 0.57; data not shown). Also between the low and high PIRCHE group, the occurrence of HLA-DPB1 mismatches did not differ (p = 0.31 and p = 0.33 for PIRCHE-I and PIRCHE-II, respectively; data not shown), suggesting that our study outcome is likely not biased by the presence of HLA-DPB1 mismatches. When analyzing the number of HLA-DPB1-derived PIRCHE-I and PIRCHE-II, the number of PIRCHE-I and PIRCHE-II derived from HLA-DPB1 did not differ between the low PIRCHE and high PIRCHE group (p = 0.56 and p = 0.66 for low vs. high PIRCHE-I and low vs. high PIRCHE-II, respectively; data not shown). Nevertheless, further studies in 10/10-MUD HSCT with a HLA-DPB1 mismatch are required to further validate the impact of HLA-DPB1-derived PIRCHE.
The current explorative study focusses on the role of PIRCHE in 9/10 MUD-HSCT. Recently, PIRCHE has been studied in the context of haploidentical HSCT (27). Huo et al. did not find a correlation between PIRCHE and clinical outcome in haploidentical HSCT. The differences in study outcome between this latter study and our current study may (partially) be explained by differences in transplant protocol, in the number of HLA mismatches, and in study design. Moreover, our study confirmed that the impact of HLA matching on OS after HSCT is more prominent in patients who are transplanted in the early phase of their disease. Although our study is restricted to 9/10 MUD-HSCT, the differential impact of HLA mismatching on OS determined by the disease status may also occur in haploidentical HSCT. The analyses in the study of Huo et al. were not stratified according to disease status, which may have impacted the study outcome. Therefore, additional studies are required to further investigate these potential differences between 9/10 MUD-HSCT and haploidentical HSCT.
The PIRCHE calculations are based on the amino acid differences between donor and recipient as well as the ability of shared HLA class I and class II to present epitopes. Therefore, one might argue that the PIRCHE score may just result from the number of amino acid differences between donor and recipient. Since these amino acid differences are one of the underlying factors in the PIRCHE calculations, it is difficult to separate them from the fact that PIRCHE calculations also involve the factor HLA restriction. Therefore, we recalculate the PIRCHE-II scores when scrambling the shared DRB1-background of the patients; i.e., these donor-recipient couples had the same number of amino acid differences, while their presenting HLA-DRB1 allele was scrambled. In these patients, no differences were observed in the OS between the low and high PIRCHE-II group (HR: 1.50; 95%-CI: 0.78–2.90; p = 0.29). Furthermore, the impaired OS of the high PIRCHE-II group compared to the 10/10-MUD group was diminished when scrambling the HLA-DRB1 presenting allele (5-years overall survival of 42% for the high PIRCHE-II group). These data indicate that the HLA-DRB1 background of the presenting allele(s) indeed significantly contributes to the observed correlation between PIRCHE-II and OS and that this correlation is not solely the effect of the number of amino acid differences between donor and recipient. These observations are also in line with previous observations in kidney transplantation (11) and in pregnancy (unpublished observations) showing that the correlation of PIRCHE-II with the immunogenicity of HLA mismatches was lost when scrambling the HLA-DRB1 background of the responder.
To identify PIRCHE, we used the NetChop, NetMHCpan, and NetMHCIIpan algorithms. Therefore, the number of PIRCHE highly depends on these algorithms. Although these algorithms are improving over time, prediction of epitope presentation can still be further improved. The prediction of epitope presentation on HLA class II molecules is especially challenging, as HLA class II molecules have a more open binding groove, resulting in a more liberate peptide binding (28, 29). Since epitopes with different lengths can bind to HLA class II molecules, resulting in a different alignment of epitopes in the HLA class II binding groove, the NetMHCIIpan algorithm predicts a core binding motif consisting of nine amino acids (30). Although the prediction of this core binding motif of nine amino acids has been becoming more and more reliable (30, 31), a better prediction of these binding motifs will further improve the PIRCHE algorithm, particularly when estimating specific risks for individual recipient-donor combinations.
The PIRCHE model is not the first model that aims for predicting alloreactivity. However, all these established models have been based on predicting direct T-cell recognition of mismatched HLA, whereas the PIRCHE model predicts indirect T-cell recognition. The Histocheck model (32) has been evaluated in various in vitro studies and a large clinical study, but did not correlate with clinical outcome (33–36). Additional studies have shown that the number of amino acid differences within the alpha helices and the beta sheet of HLA class I mismatches are to a certain extend predictive for CD8+ cytotoxic T cell alloreactivity in vitro as determined by the CTLp assay (37, 38). HLA class I mismatches that yield more than 5 amino acid differences in the alpha helices and more than 5 amino acid differences in the beta sheet appeared to less immunogenic, whereas HLA class I mismatches that yield less amino acid differences appeared to be highly immunogenic (37, 38). Furthermore, other studies have shown that specific amino acid polymorphisms influence the survival after HSCT (39–42). Additional models based on the physiochemical disparity between HLA class I antigens have shown correlation with transplant outcome; HLA class I mismatches with a high physiochemical disparity are associated with a higher risk of aGVHD after HSCT (43). Moreover, the T-cell epitope (TCE) model for HLA-DP has also shown its effect on HSCT outcome in various studies (44, 45). More recently also an expression model has been evaluated, which showed that the level of cell-surface expression of HLA-DPB1 mismatches is highly predictive for GVHD after HSCT (46, 47). The expression model and the TCE model are highly correlated, suggesting that both models may partially overlap or act synergistically in the prediction of non-permissible HLA-DPB1 mismatches (47). The TCE model may also partly overlap or complement with the PIRCHE model for HLA-DP (8), suggesting that inclusion of both a direct recognition model and an indirect recognition model into one single prediction algorithm may significantly enhance the reliability of the risk classification of specific HLA mismatches.
In conclusion, our data suggest that high PIRCHE-II scores can be classified as non-permissive as these score are associated with an impaired OS after HSCT. These effects on OS are likely due to an increased risk of severe acute GVHD, without affecting progression risk. Our current data suggest that analyzing donors before transplantation for non-permissible PIRCHE-II as part of a donor selection algorithm could be a valuable tool for HLA-mismatched HSCT.
This retrospective study was conducted following the JACIE guidelines. Patients gave informed consent to used their data according JACIE.
Contribution: KG, KT, HvD, JC, CK, JK, and ES were involved in design of the work and interpretation of the data. MO, EB, BvdH, EM, SZ, MdG, PvdB, NS, and JC were involved in acquisition of the data. All authors were involved in drafting or revising the manuscript and approved the final version. All authors agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
We further acknowledge support from our funding agency: The KWF grant (UU 2015-7601) to JK and ES, The Netherlands Organization for Scientific Research, Computational Life Sciences Program (grant number 635.100.025) to JC, and NWO Aard en Levenswetenschappen grant (823.02.014) to CK.
Conflict of Interest Statement
The UMC Utrecht has filed a patent application on the prediction of an alloimmune response against mismatched HLA. ES is listed as inventor on this patent.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Grateful acknowledgment is made to all the clinicians and HLA laboratories that provided clinical data.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2019.00880/full#supplementary-material
1. Bray RA, Hurley CK, Kamani NR, Woolfrey A, Muller C, Spellman S, et al. National marrow donor program HLA matching guidelines for unrelated adult donor hematopoietic cell transplants. Biol Blood Marrow Transplant. (2008) 14:45–53. doi: 10.1016/j.bbmt.2008.06.014
3. Lee SJ, Klein J, Haagenson M, Baxter-Lowe LA, Confer DL, Eapen M, et al. High-resolution donor-recipient HLA matching contributes to the success of unrelated donor marrow transplantation. Blood. (2007) 110:4576–83. doi: 10.1182/blood-2007-06-097386
4. van Walraven SM, Brand A, Bakker JN, Heemskerk MB, Nillesen S, Bierings MB, et al. The increase of the global donor inventory is of limited benefit to patients of non-Northwestern European descent. Haematologica. (2017) 102:176–83. doi: 10.3324/haematol.2016.145730
6. Dehn J, Buck K, Maiers M, Confer D, Hartzman R, Kollman C, et al. 8/8 and 10/10 High-resolution match rate for the be the match unrelated donor registry. Biol Blood Marrow Transplant. (2015) 21:137–41. doi: 10.1016/j.bbmt.2014.10.002
8. Thus KA, Ruizendaal MT, de Hoop TA, Borst E, van Deutekom HW, Te Boome L, et al. Refinement of the definition of permissible HLA-DPB1 mismatches with predicted indirectly recognizable HLA-DPB1 epitopes. Biol Blood Marrow Transplant. (2014) 20:1705–10. doi: 10.1016/j.bbmt.2014.06.026
10. Thus KA, de Hoop TA, de Weger RA, Bierings MB, Boelens JJ, Spierings E. Predicted indirectly recognizable HLA epitopes class I promote antileukemia responses after cord blood transplantation: indications for a potential novel donor selection tool. Biol Blood Marrow Transplant. (2016) 22:170–3. doi: 10.1016/j.bbmt.2015.08.014
11. Otten HG, Calis JJ, Kesmir C, van Zuilen AD, Spierings E. Predicted indirectly recognizable HLA epitopes presented by HLA-DR correlate with the de novo development of donor-specific HLA IgG antibodies after kidney transplantation. Hum Immunol. (2013) 74:290–96. doi: 10.1016/j.humimm.2012.12.004
12. Geneugelijk K, Honger G, van Deutekom HW, Thus KA, Kesmir C, Hosli I, et al. Predicted indirectly recognizable HLA epitopes presented by HLA-DRB1 are related to HLA antibody formation during pregnancy. Am J Transplant. (2015) 15:3112–22. doi: 10.1111/ajt.13508
13. Chaigne B, Geneugelijk K, Bedat B, Ahmed MA, Honger G, de Seigneux S, et al. Immunogenicity of anti-HLA antibodies in pancreas and islet transplantation. Cell Transplant. (2016) 25:2041–50. doi: 10.3727/096368916X691673
14. Lachmann N, Niemann M, Reinke P, Budde K, Schmidt D, Halleck F, et al. Donor recipient matching based on predicted recognizable HLA epitopes predicts the incidence of de novo donor-specific HLA antibodies following renal transplantation. Am J Transplant. (2017) 17:3076–86. doi: 10.1111/ajt.14393
16. Nielsen M, Lundegaard C, Lund O, Kesmir C. The role of the proteasome in generating cytotoxic T-cell epitopes: insights obtained from improved predictions of proteasomal cleavage. Immunogenetics. (2005) 57:33–41. doi: 10.1007/s00251-005-0781-7
17. Nielsen M, Lundegaard C, Blicher T, Lamberth K, Harndahl M, Justesen S, et al. NetMHCpan, a method for quantitative predictions of peptide binding to any HLA-A and -B locus protein of known sequence. PLoS ONE. (2007) 2:e796. doi: 10.1371/journal.pone.0000796
18. Hoof I, Peters B, Sidney J, Pedersen LE, Sette A, Lund O, et al. NetMHCpan, a method for MHC class I binding prediction beyond humans. Immunogenetics. (2009) 61:1–13. doi: 10.1007/s00251-008-0341-z
19. Nielsen M, Lundegaard C, Lund O. Prediction of MHC class II binding affinity using SMM-align, a novel stabilization matrix alignment method. BMC Bioinformatics. (2007) 8:238. doi: 10.1186/1471-2105-8-238
22. Rowlings PA, Przepiorka D, Klein JP, Gale RP, Passweg JR, Henslee-Downey PJ, et al. IBMTR Severity Index for grading acute graft-versus-host disease: retrospective comparison with Glucksberg grade. Br J Haematol. (1997) 97:855–64. doi: 10.1046/j.1365-2141.1997.1112925.x
24. Duquesnoy RJ, Spellman S, Haagenson M, Wang T, Horowitz MM, Oudshoorn M. HLAMatchmaker-defined triplet matching is not associated with better survival rates of patients with class I HLA allele mismatched hematopoietic cell transplants from unrelated donors. Biology of Blood and Marrow Transplantation. (2008) 14:1064-1071. doi: 10.1016/j.bbmt.2008.07.001
25. Shaw BE, Potter MN, Mayor NP, Pay AL, Smith C, Goldman JM, et al. The degree of matching at HLA-DPB1 predicts for acute graft-versus-host disease and disease relapse following haematopoietic stem cell transplantation. Bone Marrow Transplant. (2003) 31:1001–8. doi: 10.1038/sj.bmt.1704029
26. Shaw BE, Gooley TA, Malkki M, Madrigal JA, Begovich AB, Horowitz MM, et al. The importance of HLA-DPB1 in unrelated donor hematopoietic cell transplantation. Blood. (2007) 110:4560–66. doi: 10.1182/blood-2007-06-095265
27. Huo MR, Li D, Chang YJ, Xu LP, Zhang XH, Liu KY, et al. Predicted indirectly recognizable HLA epitopes are not associated with clinical outcomes after haploidentical hematopoietic stem cell transplantation. Hum Immunol. (2018) 79:117–21. doi: 10.1016/j.humimm.2017.11.004
28. Castellino F, Zhong G, Germain RN. Antigen presentation by MHC class II molecules: invariant chain function, protein trafficking, and the molecular basis of diverse determinant capture. Hum Immunol. (1997) 54:159–169. doi: 10.1016/S0198-8859(97)00078-5
30. Karosiene E, Rasmussen M, Blicher T, Lund O, Buus S, Nielsen M. NetMHCIIpan-3.0, a common pan-specific MHC class II prediction method including all three human MHC class II isotypes, HLA-DR, HLA-DP and HLA-DQ. Immunogenetics. (2013) 65:711–24. doi: 10.1007/s00251-013-0720-y
31. Andreatta M, Karosiene E, Rasmussen M, Stryhn A, Buus S, Nielsen M. Accurate pan-specific prediction of peptide-MHC class II binding affinity with improved binding core identification. Immunogenetics. (2015) 67:641–50. doi: 10.1007/s00251-015-0873-y
32. Elsner HA, DeLuca D, Strub J, Blasczyk R. HistoCheck: rating of HLA class I and II mismatches by an internet-based software tool. Bone Marrow Transplant. (2004) 33:165–9. doi: 10.1038/sj.bmt.1704301
34. Heemskerk MB, Doxiadis II, Roelen DL, Claas FH, Oudshoorn M. The HistoCheck algorithm does not predict T-cell alloreactivity in vitro. Bone Marrow Transplant. (2005) 36:927–8. doi: 10.1038/sj.bmt.1705154
35. Askar M, Sobecks R, Morishima Y, Kawase T, Nowacki A, Makishima H, et al. Predictions in the face of clinical reality: Histocheck versus high-risk HLA allele mismatch combinations responsible for severe acute graft-versus-host disease. Biol Blood Marrow Transplant. (2011) 17:1409–15. doi: 10.1016/j.bbmt.2011.01.015
36. Spellman S, Klein J, Haagenson M, Askar M, Baxter-Lowe LA, He J, et al. Scoring HLA class I mismatches by HistoCheck does not predict clinical outcome in unrelated hematopoietic stem cell transplantation. Biol Blood Marrow Transplant. (2012) 18:739–46. doi: 10.1016/j.bbmt.2011.09.008
37. Heemskerk MB, Cornelissen JJ, Roelen DL, van Rood JJ, Claas FH, Doxiadis II, et al. Highly diverged MHC class I mismatches are acceptable for haematopoietic stem cell transplantation. Bone Marrow Transplant. (2007) 40:193–200. doi: 10.1038/sj.bmt.1705721
38. Heemskerk MB, Roelen DL, Dankers MK, van Rood JJ, Claas FH, Doxiadis II, et al. Allogeneic MHC class I molecules with numerous sequence differences do not elicit a CTL response. Hum Immunol. (2005) 66:969–76. doi: 10.1016/j.humimm.2005.06.007
39. Crivello P, Heinold A, Rebmann V, Ottinger HD, Horn PA, Beelen DW, et al. Functional distance between recipient and donor HLA-DPB1 determines nonpermissive mismatches in unrelated HCT. Blood. (2016) 128:120-129. doi: 10.1182/blood-2015-12-686238
40. Marino SR, Lin S, Maiers M, Haagenson M, Spellman S, Klein JP, et al. Identification by random forest method of HLA class I amino acid substitutions associated with lower survival at day 100 in unrelated donor hematopoietic cell transplantation. Bone Marrow Transplant. (2012) 47:217–26. doi: 10.1038/bmt.2011.56
41. Ferrara GB, Bacigalupo A, Lamparelli T, Lanino E, Delfino L, Morabito A, et al. Bone marrow transplantation from unrelated donors: the impact of mismatches with substitutions at position 116 of the human leukocyte antigen class I heavy chain. Blood. (2001) 98:3150–5. doi: 10.1182/blood.V98.10.3150
42. Pidala J, Wang T, Haagenson M, Spellman S, Askar M, Battiwalla M, et al. Amino acid substitution at peptide-binding pockets of HLA class I molecules increases risk of severe acute GVHD and mortality. Blood. (2013) 122:3651–8. doi: 10.1182/blood-2013-05-501510
43. Kosmoliaptsis V, Chaudhry AN, Sharples LD, Halsall DJ, Dafforn TR, Bradley JA, et al. Predicting HLA class I alloantigen immunogenicity from the number and physiochemical properties of amino acid polymorphisms. Transplantation. (2009) 88:791–8. doi: 10.1097/TP.0b013e3181b4a9ff
44. Fleischhauer K, Shaw BE, Gooley T, Malkki M, Bardy P, Bignon JD, et al. Effect of T-cell-epitope matching at HLA-DPB1 in recipients of unrelated-donor haemopoietic-cell transplantation: a retrospective study. Lancet Oncol. (2012) 13:366–74. doi: 10.1016/S1470-2045(12)70004-9
45. Zino E, Frumento G, Marktel S, Sormani MP, Ficara F, Di Terlizzi S, et al. A T-cell epitope encoded by a subset of HLA-DPB1 alleles determines nonpermissive mismatches for hematologic stem cell transplantation. Blood. (2004) 103:1417–24. doi: 10.1182/blood-2003-04-1279
46. Petersdorf EW, Malkki M, O'hUigin C, Carrington M, Gooley T, Haagenson MD, et al. High HLA-DP expression and graft-versus-host disease. N Engl J Med. (2015) 373:599–609. doi: 10.1056/NEJMoa1500140
Keywords: HLA, PIRCHE, Non-permissible mismatch, HSCT—hematopoietic stem cell transplant, HLA mismatch
Citation: Geneugelijk K, Thus KA, van Deutekom HWM, Calis JJA, Borst E, Keşmir C, Oudshoorn M, van der Holt B, Meijer E, Zeerleder S, de Groot MR, von dem Borne PA, Schaap N, Cornelissen J, Kuball J and Spierings E (2019) Exploratory Study of Predicted Indirectly ReCognizable HLA Epitopes in Mismatched Hematopoietic Cell Transplantations. Front. Immunol. 10:880. doi: 10.3389/fimmu.2019.00880
Received: 21 December 2018; Accepted: 05 April 2019;
Published: 24 April 2019.
Edited by:Michael Uhlin, Karolinska Institute (KI), Sweden
Reviewed by:Marco Andreani, Bambino Gesù Ospedale Pediatrico (IRCCS), Italy
Bjarne Kuno Møller, Aarhus University Hospital, Denmark
Copyright © 2019 Geneugelijk, Thus, van Deutekom, Calis, Borst, Keşmir, Oudshoorn, van der Holt, Meijer, Zeerleder, de Groot, von dem Borne, Schaap, Cornelissen, Kuball and Spierings. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Eric Spierings, email@example.com
†These authors have contributed equally to this work