Human leukocyte antigen evolutionary divergence influences outcomes of paediatric patients and young adults affected by malignant disorders given allogeneic haematopoietic stem cell transplantation from unrelated donors.

High genetic heterogeneity in the human leukocyte antigen (HLA) increases the likelihood of efficient immune response to pathogens and tumours. As measure of HLA diversity, HLA evolutionary divergence (HED) has been shown to predict the response of tumours to immunotherapy and haematopoietic stem cell transplantation (HSCT) in adults. We retrospectively investigated the association of HED with outcomes of 153 paediatric/young adults patients, treated for malignant disorders with HSCT from 9-10/10 HLA-matched unrelated donors. HED was calculated as pairwise genetic distance between alleles in patient HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1, using the locus median to stratify patients with 'high' or 'low' HED. Patients with high HED-B and -DRB1 showed significantly improved disease-free survival (DFS), especially when combined (70.8% vs 53.7% p = 0.008). High HED-B + -DRB1 was also associated with improved overall survival (OS) (82.1 vs 66.4% p = 0.014), and concomitant reduction of non-relapse-mortality (5.1% vs 21.1% p = 0.006). The impact on OS and DFS of combined HED-B + -DRB1 was confirmed in multivariate analysis [hazard ratio (HR) 0.39, p = 0.009; and HR 0.45, p = 0.007 respectively]. Only high HED scores for HLA-DPB1 were associated, in univariate analysis, with reduced incidence of relapse (15.9% vs 31.1%, p = 0.03). These results support HED as prognostic marker in allogeneic HSCT and, if confirmed in larger cohorts, would allow its use to inform clinical risk and potentially influence clinical practice.

Permissive HLA-DPB1 mismatches in HCT depend on immunopeptidome divergence and editing by HLA-DM.

In hematopoietic cell transplantation (HCT), permissive HLA-DPB1 mismatches between patients and their unrelated donors are associated with improved outcomes compared with nonpermissive mismatches, but the underlying mechanism is incompletely understood. Here, we used mass spectrometry, T-cell receptor-ß (TCRß) deep sequencing, and cellular in vitro models of alloreactivity to interrogate the HLA-DP immunopeptidome and its role in alloreactive T-cell responses. We find that permissive HLA-DPB1 mismatches display significantly higher peptide repertoire overlaps compared with their nonpermissive counterparts, resulting in lower frequency and diversity of alloreactive TCRß clonotypes in healthy individuals and transplanted patients. Permissiveness can be reversed by the absence of the peptide editor HLA-DM or the presence of its antagonist, HLA-DO, through significant broadening of the peptide repertoire. Our data establish the degree of immunopeptidome divergence between donor and recipient as the mechanistic basis for the clinically relevant permissive HLA-DPB1 mismatches in HCT and show that permissiveness is dependent on HLA-DM-mediated peptide editing. Its key role for harnessing T-cell alloreactivity to HLA-DP highlights HLA-DM as a potential novel target for cellular and immunotherapy of leukemia.

Immunopeptidome Analysis of HLA-DPB1 Allelic Variants Reveals New Functional Hierarchies.

HLA-DP alleles can be classified into functional T cell epitope (TCE) groups. TCE-1 and TCE-2 are clearly defined, but TCE-3 still represents an heterogeneous group. Because polymorphisms in HLA-DP influence the presented peptidome, we investigated whether the composition of peptides binding in HLA-DP may be used to refine the HLA-DP group classification. Peptidomes of human HLA-DP-typed B cell lines were analyzed with mass spectrometry after immunoaffinity chromatography and peptide elution. Gibbs clustering was performed to identify motifs of binding peptides. HLA-DP peptide-binding motifs showed a clear association with the HLA-DP allele-specific sequences of the binding groove. Hierarchical clustering of HLA-DP immunopeptidomes was performed to investigate the similarities and differences in peptidomes of different HLA-DP molecules, and this clustering resulted in the categorization of HLA-DP alleles into 3-DP peptidome clusters (DPC). The peptidomes of HLA-DPB1*09:01, -10:01, and -17:01 (TCE-1 alleles) and HLA-DPB1*04:01, -04:02, and -02:01 (TCE-3 alleles) were separated in two maximal distinct clusters, DPC-1 and DPC-3, respectively, reflecting their previous TCE classification. HLA-DP alleles categorized in DPC-2 shared certain similar peptide-binding motifs with DPC-1 or DPC-3 alleles, but significant differences were observed for other positions. Within DPC-2, divergence between the alleles was observed based on the preference for different peptide residues at position 9. In summary, immunopeptidome analysis was used to unravel functional hierarchies among HLA-DP alleles, providing new molecular insights into HLA-DP classification.

Cytomegalovirus kinetics after hematopoietic cell transplantation reveal peak titers with differential impact on mortality, relapse and immune reconstitution.

Even in the era of PCR-based monitoring, prophylaxis, and preemptive therapy, Cytomegalovirus (CMV) viremia remains a relevant cause of non-relapse mortality (NRM) after allogeneic hematopoietic cell transplantation (HCT). However, studies using binary analysis (presence/absence of CMV) reported contradicting data for NRM, overall survival and leukemia relapse. Here, we analyzed CMV replication kinetics in 11 508 whole blood PCR samples of 705 patients with HCT between 2012 and 2017. Using two independent models based on CMV peak titers and on the time point of first CMV reactivation, we stratified patients into risk cohorts. Each cohort had distinct cellular immune reconstitution profiles and differentiated for relevant clinical outcomes. Patients with high CMV peak titers had significantly reduced overall survival (HR 2.13, 95% CI 1.53-2.96; p < .0001), due to high NRM. Early impaired T cell reconstitution was a risk factor for high CMV peak titers, however relevant CMV viremia also related to boosted T cell reconstitution. Importantly, intermediate CMV peak titers associated with a significantly reduced relapse probability (HR 0.53, 95% CI 0.31-0.91; p = .022). In short, CMV kinetics models distinguished relevant clinical outcome cohorts beyond the R+ serostatus with distinct immune reconstitution patterns and resolve in part contradicting results of previous studies exclusively focused on the presence or absence of CMV.

Multiple Knockout of Classical HLA Class II ß-Chains by CRISPR/Cas9 Genome Editing Driven by a Single Guide RNA.

Comprehensive knockout of HLA class II (HLA-II) ß-chain genes is complicated by their high polymorphism. In this study, we developed CRISPR/Cas9 genome editing to simultaneously target HLA-DRB, -DQB1, and -DPB1 through a single guide RNA recognizing a conserved region in exon 2. Abrogation of HLA-II surface expression was achieved in five different HLA-typed, human EBV-transformed B lymphoblastoid cell lines (BLCLs). Next-generation sequencing-based detection confirmed specific genomic insertion/deletion mutations with 99.5% penetrance in sorted cells for all three loci. No alterations were observed in HLA-I genes, the HLA-II peptide editor HLA-DMB, or its antagonist HLA-DOB, showing high on-target specificity. Transfection of full-length HLA-DPB1 mRNA into knockout BLCLs fully restored HLA-DP surface expression and recognition by alloreactive human CD4 T cells. The possibility to generate single HLA-II-expressing BLCLs by one-shot genome editing opens unprecedented opportunities for mechanistically dissecting the interaction of individual HLA variants with the immune system.

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.

Clinical Utility of Quantitative PCR for Chimerism and Engraftment Monitoring after Allogeneic Stem Cell Transplantation for Hematologic Malignancies.

Although quantitative PCR (qPCR) has been explored for chimerism monitoring after allogeneic stem cell transplantation (SCT), evidence regarding its clinical utility compared with standard short tandem repeat (STR) is still limited. We retrospectively studied commercial qPCR and STR chimerism with respective positivity thresholds of .1% and 1% in 359 peripheral blood (PB) and 95 bone marrow (BM) samples from 30 adult patients after first HLA-matched SCT for myeloid malignancies or acute lymphatic leukemia. Concordance between the 2 methods was 79.5%, with all discordant samples positive in qPCR but negative in STR. Of the latter, sporadic qPCR positivity without clinical correlates was seen mostly in BM samples early post-transplant. In 7 of 21 patients with available follow-up samples in the first months after transplantation, qPCR but not STR revealed low levels (<1%) of sustained host chimerism in PB, reflecting delayed engraftment or persistent mixed chimerism (PMC). These conditions were associated with donor-recipient cytomegalovirus (CMV) serostatus and early CMV reactivation but not with immunosuppressive regimens or clinical outcome. qPCR predicted all 8/8 relapses with samples in the 6 months before onset by sustained positivity in both PB and BM compared with 1/8 relapses predicted by STR mainly in BM. The response kinetics to donor lymphocyte infusions for the treatment of PMC or relapse was shown by qPCR but not STR to be protracted over several months in 3 patients. Our results demonstrate the superior clinical utility of qPCR compared with STR for monitoring subtle changes of host chimerism associated with different clinical conditions, making a case for its use in the clinical follow-up of transplant patients.

High Frequencies of Anti-Host Reactive CD8+ T Cells Ignore Non-Hematopoietic Antigen after Bone Marrow Transplantation in a Murine Model.

BACKGROUND: Graft versus host disease (GvHD) occurs in 20% of cases with patients having an MHC I matched bone marrow transplantation (BMT). Mechanisms causing this disease remain to be studied. METHODS: Here we used a CD8+ T cell transgenic mouse line (P14/CD45.1+) and transgenic DEE mice bearing ubiquitously the glycoprotein 33-41 (GP33) antigen derived from the major lymphocytic choriomeningitis virus (LCMV) epitope to study mechanisms of tolerance in anti-host reactive CD8+ T cells after BMT. RESULTS: We found that anti-host reactive CD8+ T cells (P14 T cells) were not negatively selected in the thymus and that they were present in wild type (WT) recipient mice as well as in DEE recipient mice. Anti-host reactive CD8+ T cells ignored the GP33 antigen expressed ubiquitously by host cells but they could be activated ex vivo via LCMV-infection. Lipopolysaccharides (LPS) induced transient cell damage in DEE mice bearing anti-host reactive CD8+ T cells after BMT, suggesting that induction of host inflammatory response could break antigen ignorance. Introducing the GP33 antigen into BM cells led to deletion of anti-host reactive CD8+ T cells. CONCLUSION: We found that after BMT anti-host reactive CD8+ T cells ignored host antigen in recipients and that they were only deleted when host antigen was present in hematopoietic cells. Moreover, LPS-induced immune activation contributed to induction of alloreactivity of anti-host reactive CD8+ T cells after BMT.

The impact of amino acid variability on alloreactivity defines a functional distance predictive of permissive HLA-DPB1 mismatches in hematopoietic stem cell transplantation.

A major challenge in unrelated hematopoietic stem cell transplantation (HSCT) is the prediction of permissive HLA mismatches, ie, those associated with lower clinical risks compared to their nonpermissive counterparts. For HLA-DPB1, a clinically prognostic model has been shown to be matching for T cell epitope (TCE) groups assigned by cross reactivity of T cells alloreactive to HLA-DPB1∗09:01; however, the molecular basis of this observation is not fully understood. Here, we have mutated amino acids (aa) in 10 positions of HLA-DPB1∗09:01 to other naturally occurring variants, expressed them by lentiviral vectors in B cell lines, and quantitatively measured allorecognition by 17 CD4(+) T cell effectors from 6 unrelated individuals. A significant impact on the median alloresponse was observed for peptide contact positions 9, 11, 35, 55, 69, 76, and 84, but not for positions 8, 56, and 57 pointing away from the groove. A score for the "functional distance" (FD) from HLA-DPB1∗09:01 was defined as the sum of the median impact of polymorphic aa in a given HLA-DPB1 allele on T cell alloreactivity. Established TCE group assignment of 23 alleles correlated with FD scores of ≤0.5, 0.6 to 1.9 and ≥2 for TCE groups 1, 2, and 3, respectively. Based on this, prediction of TCE group assignment will be possible for any given HLA-DPB1 allele, including currently 367 alleles encoding distinct proteins for which T cell cross reactivity patterns are unknown. Experimental confirmation of the in silico TCE group classification was successfully performed for 7 of 7 of these alleles. Our findings have practical implications for the applicability of TCE group matching in unrelated HSCT and provide new insights into the molecular mechanisms underlying this model. The innovative concept of FD opens new potential avenues for risk prediction in unrelated HSCT.

Qualitative, rather than quantitative, differences between HLA-DQ alleles affect HLA-DQ immunogenicity in organ transplantation.

Prolonging the lifespan of transplanted organs is critical to combat the shortage of this life-saving resource. Chronic rejection, with irreversible demise of the allograft, is often caused by the development of donor-specific HLA antibodies. Currently, enumerating molecular (amino acid) mismatches between recipient and donor is promoted to identify patients at higher risk of developing HLA antibodies, for use in organ allocation, and immunosuppression-minimization strategies. We have counseled against the incorporation of such approaches into clinical use and hypothesized that not all molecular mismatches equally contribute to generation of donor-specific immune responses. Herein, we document statistical shortcomings in previous study design: for example, use of individuals who lack the ability to generate donor-specific-antibodies (HLA identical) as part of the negative cohort. We provide experimental evidence, using CRISPR-Cas9-edited cells, to rebut the claim that the HLAMatchmaker eplets represent "functional epitopes." We further used unique sub-cohorts of patients, those receiving an allograft with two HLA-DQ mismatches yet developing antibodies only to one mismatch (2MM1DSA), to interrogate differential immunogenicity. Our results demonstrate that mismatches of DQα05-heterodimers exhibit the highest immunogenicity. Additionally, we demonstrate that the DQα chain critically contributes to the overall qualities of DQ molecules. Lastly, our data proposes that an augmented risk to develop donor-specific HLA-DQ antibodies is dependent on qualitative (evolutionary and functional) divergence between recipient and donor, rather than the mere number of molecular mismatches. Overall, we propose an immunological mechanistic rationale to explain differential HLA-DQ immunogenicity, with potential ramifications for other pathological processes such as autoimmunity and infections.

Impact of the HLA Immunopeptidome on Survival of Leukemia Patients After Unrelated Donor Transplantation.

PURPOSE: Immunopeptidome divergence between mismatched HLA-DP is a determinant of T-cell alloreactivity and clinical tolerability after fully HLA-A, -B, -C, -DRB1, -DQB1 matched unrelated donor hematopoietic cell transplantation (UD-HCT). Here, we tested this concept in HLA-A, -B, and -C disparities after single class I HLA-mismatched UD-HCT. PATIENTS AND METHODS: We studied 2,391 single class I HLA-mismatched and 14,426 fully HLA-matched UD-HCT performed between 2008 and 2018 for acute leukemia or myelodysplastic syndromes. Hierarchical clustering of experimentally determined peptide-binding motifs (PBM) was used as a proxy for immunopeptidome divergence of HLA-A, -B, or -C disparities, allowing us to classify 1,629/2,391 (68.1%) of the HLA-mismatched UD-HCT as PBM-matched or PBM-mismatched. Risks associated with PBM-matching status were assessed by Cox proportional hazards models, with overall survival (OS) as the primary end point. RESULTS: Relative to full matches, bidirectional or unidirectional PBM mismatches in graft-versus-host (GVH) direction (PBM-GVH mismatches, 60.7%) were associated with significantly lower OS (hazard ratio [HR], 1.48; P < .0001), while unidirectional PBM mismatches in host-versus-graft direction or PBM matches (PBM-GVH matches, 39.3%) were not (HR, 1.13; P = .1017). PBM-GVH mismatches also had significantly lower OS than PBM-GVH matches in direct comparison (HR, 1.32; P = .0036). The hazards for transplant-related mortality and acute and chronic graft-versus-host disease but not relapse increased stepwise from full HLA matches to single PBM-GVH matches, and single PBM-GVH mismatches. A webtool for PBM-matching of single class I HLA-mismatched donor-recipient pairs was developed. CONCLUSION: PBM-GVH mismatches inform mortality risks after single class I HLA-mismatched UD-HCT, suggesting that prospective consideration of directional PBM-matching status might improve outcome. These findings highlight immunopeptidome divergence between mismatched HLA as a driver of clinical tolerability in UD-HCT.

Impact of CMV reactivation on relapse of acute myeloid leukemia after HCT is dependent on disease stage and ATG.

Cytomegalovirus (CMV) reactivation is a frequent complication after allogeneic hematopoietic cell transplantation (HCT), whose impact on clinical outcome, in particular on leukemic relapse, is controversial. We retrospectively analyzed 687 HCT recipients with acute myeloid leukemia (AML) and ciclosporin-based immunosuppression to better understand the differential impact of CMV on transplant outcomes depending on AML disease stage and in vivo T cell depletion with antithymocyte globulin (ATG). Without ATG, CMV reactivation associated with significantly reduced relapse, yet its effect was more pronounced for advanced disease AML (P = .0002) than for patients in first complete remission (CR1, P = .0169). Depending on the disease stage, ATG exposure abrogated relapse protection following CMV reactivation in advanced stages (P = .796), while it inverted its effect into increased relapse for CR1 patients (P = .0428). CMV reactivation was associated with significantly increased nonrelapse mortality in CR1 patients without ATG (P = .0187) but not in those with advanced disease and ATG. Following CMV reactivation, only patients with advanced disease had significantly higher event-free survival rates as compared with patients without CMV. Overall, our data suggest that both ATG and disease stage modulate the impact of post-HCT CMV reactivation in opposite directions, revealing a level of complexity that warrants future studies regarding the interplay between antivirus and antitumor immunity.

Melanoma Differentiation Trajectories Determine Sensitivity toward Pre-Existing CD8+ Tumor-Infiltrating Lymphocytes.

The highly plastic nature of melanoma enables its transition among diverse cell states to survive hostile conditions. However, the interplay between specific tumor cell states and intratumoral T cells remains poorly defined. With MAPK inhibitor‒treated BRAFV600-mutant tumors as models, we linked human melanoma state transition to CD8+ T cell responses. Repeatedly, we observed that isogenic melanoma cells could evolve along distinct differentiation trajectories on single BRAF inhibitor (BRAFi) treatment or dual BRAFi/MEKi treatment, resulting in BRAFi‒induced hyperdifferentiated and BRAFi/MEKi‒induced dedifferentiated resistant subtypes. Taking advantage of patient-derived autologous CD8+ tumor-infiltrating lymphocytes (TILs), we demonstrate that progressive melanoma cell state transition profoundly affects TIL function. Tumor cells along the hyperdifferentiation trajectory continuously gained sensitivity toward tumor-reactive CD8+ TILs, whereas those in the dedifferentiation trajectory acquired T cell resistance in part owing to the loss of differentiation antigens. Overall, our data reveal the tight connection of MAPKi‒induced temporary (drug-tolerant transition state) and stable (resistant state) phenotype alterations with T cell function and further broaden the current knowledge on melanoma plasticity in terms of sculpting local antitumor immune responses, with implications for guiding the optimal combination of targeted therapy and immunotherapy.

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

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

Alloreactive T Cell Receptor Diversity against Structurally Similar or Dissimilar HLA-DP Antigens Assessed by Deep Sequencing.

T cell alloreactivity is mediated by a self-human leukocyte antigen (HLA)-restricted T cell receptor (TCR) repertoire able to recognize both structurally similar and dissimilar allogeneic HLA molecules (i.e., differing by a single or several amino acids in their peptide-binding groove). We hypothesized that thymic selection on self-HLA molecules could have an indirect impact on the size and diversity of the alloreactive response. To test this possibility, we used TCR Vß immunophenotyping and immunosequencing technology in a model of alloreactivity between self-HLA selected T cells and allogeneic HLA-DPB1 (DPB1) differing from self-DPB1*04:02 by a single (DPB1*02:01) or several (DPB1*09:01) amino acids in the peptide-binding groove. CD4+ T cells from three different self-DPB1*04:01,*04:02 individuals were stimulated with HeLa cells stably transduced with the relevant peptide processing machinery, co-stimulatory molecules, and HLA-DP. Flow cytometric quantification of the DPB1-specific T cell response measured as upregulation of the activation marker CD137 revealed significantly lower levels of alloreactivity against DPB1*02:01 compared with DPB1*09:01 (mean CD4+CD137+ frequency 35.2 ± 9.9 vs. 61.5 ± 7.7%, respectively, p < 0.0001). These quantitative differences were, however, not reflected by differences in the breadth of the alloreactive response at the Vß level, with both alloantigens eliciting specific responses from all TCR-Vß specificities tested by flow cytometry, albeit with higher levels of reactivity from most Vß specificities against DPB1*09:01. In line with these observations, TCRB-CDR3 immunosequencing showed no significant differences in mean clonality of sorted CD137+CD4+ cells alloreactive against DPB1*02:01 or DPB1*09:01 [0.39 (0.36-0.45) and 0.39 (0.30-0.46), respectively], or in the cumulative frequencies of the 10 most frequent responding clones (55-67 and 58-62%, respectively). Most of the clones alloreactive against DPB1*02:01 (68.3%) or DPB1*09:01 (75.3%) were characterized by low-abundance (i.e., they were not appreciable among the pre-culture T cells). Interestingly, however, their cumulative frequency was lower against DPB1*02:01 compared with DPB1*09:01 (mean cumulative frequency 35.3 vs. 50.6%, respectively). Our data show that, despite lower levels of alloreactivity, a similar clonal diversity can be elicited by structurally similar compared with structurally dissimilar HLA-DPB1 alloantigens and demonstrate the power of TCRB immunosequencing in unraveling subtle qualitative changes not appreciable by conventional methods.

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.