Dynamical System Modeling to Simulate Donor T Cell Response to Whole Exome Sequencing-Derived Recipient Peptides Demonstrates Different Alloreactivity Potential in HLA-Matched and -Mismatched Donor-Recipient Pairs.

Immune reconstitution kinetics and subsequent clinical outcomes in HLA-matched recipients of allogeneic stem cell transplantation (SCT) are variable and difficult to predict. Considering SCT as a dynamical system may allow sequence differences across the exomes of the transplant donors and recipients to be used to simulate an alloreactive T cell response, which may allow better clinical outcome prediction. To accomplish this, whole exome sequencing was performed on 34 HLA-matched SCT donor-recipient pairs (DRPs) and the nucleotide sequence differences translated to peptides. The binding affinity of the peptides to the relevant HLA in each DRP was determined. The resulting array of peptide-HLA binding affinity values in each patient was considered as an operator modifying a hypothetical T cell repertoire vector, in which each T cell clone proliferates in accordance with the logistic equation of growth. Using an iterating system of matrices, each simulated T cell clone's growth was calculated with the steady-state population being proportional to the magnitude of the binding affinity of the driving HLA-peptide complex. Incorporating competition between T cell clones responding to different HLA-peptide complexes reproduces a number of features of clinically observed T cell clonal repertoire in the simulated repertoire, including sigmoidal growth kinetics of individual T cell clones and overall repertoire, Power Law clonal frequency distribution, increase in repertoire complexity over time with increasing clonal diversity, and alteration of clonal dominance when a different antigen array is encountered, such as in SCT. The simulated, alloreactive T cell repertoire was markedly different in HLA-matched DRPs. The patterns were differentiated by rate of growth and steady-state magnitude of the simulated T cell repertoire and demonstrate a possible correlation with survival. In conclusion, exome wide sequence differences in DRPs may allow simulation of donor alloreactive T cell response to recipient antigens and may provide a quantitative basis for refining donor selection and titration of immunosuppression after SCT.

The influence of lymphoid reconstitution kinetics on clinical outcomes in allogeneic stem cell transplantation.

Lymphoid recovery following myeloablative stem cell transplantation (SCT) displays a logistic pattern of exponential growth followed by a plateau. Within this logistic framework, lymphoid recovery is characterized by the parameters R (slope of ascent), a (time of maximal rate of ascent) and K (plateau), the 'steady-state' lymphocyte count. A retrospective analysis of allogeneic SCT performed from 2008 to 2013 was undertaken to compare lymphoid recovery and clinical outcomes in 131 patients with acute myelogenous leukemia, acute lymphocytic leukemia, and myelodysplastic syndromes. Using Prism software, a logistic curve was successfully fit to the absolute lymphocyte count recovery in all patients. Patients were classified according to the magnitude and rate of lymphoid recovery; pattern A achieved an absolute lymphocyte counts (ALC) of >1000/µL by day 45, pattern B an ALC 500 < x < 1000/µL, and pattern C an ALC <500/µL. Pattern A was characterized by a higher mean K (p < .0001) compared with patterns B and C. Patients with patterns B and C were more likely to have mixed T cell chimerism at 90 d following SCT (p = .01). There was a trend towards improved survival (and relapse-free survival) in those with pattern A and B at 1 year compared to pattern C (p = .073). There was no difference in cGVHD (p = .42) or relapse (p = .45) between pattern types. Cytomegalovirus (CMV), aGVHD, and all relapse were heralded by deviation from logistic behavior. Pattern C patients were more likely to require donor lymphocyte infusion (DLI) (p = .017). Weaning of tacrolimus post-transplant was associated with a second, separate logistic expansion in some patients. This study demonstrated that lymphoid reconstitution follows a prototypical logistic recovery and that pattern observed correlates with T cell chimerism and need for DLI, and may influence survival.

Sequence homology between HLA-bound cytomegalovirus and human peptides: A potential trigger for alloreactivity.

Human cytomegalovirus (hCMV) reactivation may often coincide with the development of graft-versus-host-disease (GVHD) in stem cell transplantation (SCT). Seventy seven SCT donor-recipient pairs (DRP) (HLA matched unrelated donor (MUD), n = 50; matched related donor (MRD), n = 27) underwent whole exome sequencing to identify single nucleotide polymorphisms (SNPs) generating alloreactive peptide libraries for each DRP (9-mer peptide-HLA complexes); Human CMV CROSS (Cross-Reactive Open Source Sequence) database was compiled from NCBI; HLA class I binding affinity for each DRPs HLA was calculated by NetMHCpan 2.8 and hCMV- derived 9-mers algorithmically compared to the alloreactive peptide-HLA complex libraries. Short consecutive (≥6) amino acid (AA) sequence homology matching hCMV to recipient peptides was considered for HLA-bound-peptide (IC50<500nM) cross reactivity. Of the 70,686 hCMV 9-mers contained within the hCMV CROSS database, an average of 29,658 matched the MRD DRP alloreactive peptides and 52,910 matched MUD DRP peptides (p<0.001). In silico analysis revealed multiple high affinity, immunogenic CMV-Human peptide matches (IC50<500 nM) expressed in GVHD-affected tissue-specific manner. hCMV+GVHD was found in 18 patients, 13 developing hCMV viremia before GVHD onset. Analysis of patients with GVHD identified potential cross reactive peptide expression within affected organs. We propose that hCMV peptide sequence homology with human alloreactive peptides may contribute to the pathophysiology of GVHD.