Dynamical System Modeling of Immune Reconstitution after Allogeneic Stem Cell Transplantation Identifies Patients at Risk for Adverse Outcomes

Joint models quantify associations between immune cell kinetics and allo-immunological events after allogeneic stem cell transplantation and subsequent donor lymphocyte infusion

Alloreactive donor-derived T-cells play a pivotal role in alloimmune responses after allogeneic hematopoietic stem cell transplantation (alloSCT); both in the relapse-preventing Graft-versus-Leukemia (GvL) effect and the potentially lethal complication Graft-versus-Host-Disease (GvHD). The balance between GvL and GvHD can be shifted by removing T-cells via T-cell depletion (TCD) to reduce the risk of GvHD, and by introducing additional donor T-cells (donor lymphocyte infusions [DLI]) to boost the GvL effect. However, the association between T-cell kinetics and the occurrence of allo-immunological events has not been clearly demonstrated yet. Therefore, we investigated the complex associations between the T-cell kinetics and alloimmune responses in a cohort of 166 acute leukemia patients receiving alemtuzumab-based TCD alloSCT. Of these patients, 62 with an anticipated high risk of relapse were scheduled to receive a prophylactic DLI at 3 months after transplant. In this setting, we applied joint modelling which allowed us to better capture the complex interplay between DLI, T-cell kinetics, GvHD and relapse than traditional statistical methods. We demonstrate that DLI can induce detectable T-cell expansion, leading to an increase in total, CD4+ and CD8+ T-cell counts starting at 3 months after alloSCT. CD4+ T-cells showed the strongest association with the development of alloimmune responses: higher CD4 counts increased the risk of GvHD (hazard ratio 2.44, 95% confidence interval 1.45-4.12) and decreased the risk of relapse (hazard ratio 0.65, 95% confidence interval 0.45-0.92). Similar models showed that natural killer cells recovered rapidly after alloSCT and were associated with a lower risk of relapse (HR 0.62, 95%-CI 0.41-0.93). The results of this study advocate the use of joint models to further study immune cell kinetics in different settings.

Allogeneic haematopoietic cell transplants as dynamical systems: influence of early-term immune milieu on long-term T-cell recovery

Objectives

Immune recovery following haematopoietic cell transplantation (HCT) functions as a dynamical system. Reducing the duration of intense immune suppression and augmenting antigen presentation has the potential to optimise T-cell reconstitution, potentially influencing long-term outcomes.

Methods

Based on donor-derived T-cell recovery, 26 patients were adaptively randomised between mycophenolate mofetil (MMF) administered for 30-day post-transplant with filgrastim for cytokine support (MMF30 arm, N = 11), or MMF given for 15 days with sargramostim (MMF15 arm, N = 15). All patients underwent in vivo T-cell depletion with 5.1 mg kg^-1 antithymocyte globulin (administered over 3 days, Day -9 through to Day -7) and received reduced intensity 450 cGy total body irradiation (3 fractions on Day -1 and Day 0). Patients underwent HLA-matched related and unrelated donor haematopoietic cell transplantation (HCT).

Results

Clinical outcomes were equivalent between the two groups. The MMF15 arm demonstrated superior T-cell, as well as T-cell subset recovery and a trend towards superior T-cell receptor (TCR) diversity in the first month with this difference persisting through the first year. T-cell repertoire recovery was more rapid and sustained, as well as more diverse in the MMF15 arm.

Conclusion

The long-term superior immune recovery in the MMF15 arm, administered GMCSF, is consistent with a disproportionate impact of early interventions in HCT. Modifying the 'immune-milieu' following allogeneic HCT is feasible and may influence long-term T-cell recovery.

The Composite Immune Risk Score predicts overall survival after allogeneic hematopoietic stem cell transplantation: A retrospective analysis of 1838 cases

There has been little consensus on how to quantitatively assess immune reconstitution after hematopoietic stem cell transplantation (HSCT) as part of the standard of care. We retrospectively analyzed 11 150 post-transplant immune profiles of 1945 patients who underwent HSCT between 2012 and 2020. 1838 (94.5%) of the cases were allogeneic HSCT. Using the training set of patients (n = 729), we identified a composite immune signature (integrating neutrophil, total lymphocyte, natural killer, total T, CD4⁺ T, and B cell counts in the peripheral blood) during days 91-180 after allogeneic HSCT that was predictive of early mortality and moreover simplified it into a formula for a Composite Immune Risk Score. When we verified the Composite Immune Risk Score in the validation (n = 284) and test (n = 391) sets of patients, a high score value was found to be associated with hazard ratios (HR) of 3.64 (95% C.I. 1.55-8.51; p = .0014) and 2.44 (95% C.I., 1.22-4.87; p = .0087), respectively, for early mortality. In multivariate analysis, a high Composite Immune Risk Score during days 91-180 remained an independent risk factor for early mortality after allogeneic HSCT (HR, 1.80; 95% C.I., 1.28-2.55; p = .00085). In conclusion, the Composite Immune Risk Score is easy to compute and could identify the high-risk patients of allogeneic HSCT who require targeted effort for prevention and control of infection.

Blasts in context: the impact of the immune environment on acute myeloid leukemia prognosis and treatment

Acute myeloid leukemia (AML) is a cancer that originates from the bone marrow (BM). Under physiological conditions, the bone marrow supports the homeostasis of immune cells and hosts memory lymphoid cells. In this review, we summarize our present understanding of the role of the immune microenvironment on healthy bone marrow and on the development of AML, with a focus on T cells and other lymphoid cells. The types and function of different immune cells involved in the AML microenvironment as well as their putative role in the onset of disease and response to treatment are presented. We also describe how the immune context predicts the response to immunotherapy in AML and how these therapies modulate the immune status of the bone marrow. Finally, we focus on allogeneic stem cell transplantation and summarize the current understanding of the immune environment in the post-transplant bone marrow, the factors associated with immune escape and relevant strategies to prevent and treat relapse.

High dose valacyclovir for cytomegalovirus prophylaxis following allogeneic hematopoietic cell transplantation

INTRODUCTION

Cytomegalovirus (CMV) is one of the most common and clinically significant viral infections following allogeneic hematopoietic cell transplantation (HCT). Currently available options for CMV prophylaxis and treatment present challenges related to side effects and cost.

METHODS

In this retrospective medical record review, the incidence of clinically significant CMV infection (CMV disease or reactivation requiring preemptive treatment) following allogeneic HCT was compared in patients receiving valacyclovir 1 g three times daily versus acyclovir 400 mg every 12 h for viral prophylaxis.

RESULTS

Forty-five patients who received valacyclovir were matched based on propensity scoring to 35 patients who received acyclovir. All patients received reduced-intensity conditioning regimens containing anti-thymocyte globulin. Clinically significant CMV infection by day + 180 was lower in the valacyclovir group compared to the acyclovir group (18% vs. 57%, p = 0.0004). Patients receiving valacyclovir prophylaxis also had less severe infection evidenced by a reduction in CMV disease, lower peak CMV titers, delayed CMV reactivation, and less secondary neutropenia.

CONCLUSION

Prospective evaluation of valacyclovir 1 g three times daily for viral prophylaxis following allogeneic HCT is warranted. Due to valacyclovir's favorable toxicity profile and affordable cost, it has the potential to benefit patients on a broad scale as an option for CMV prophylaxis.

Dynamical Systems Modeling of Early-Term Immune Reconstitution with Different Antithymocyte Globulin Administration Schedules in Allogeneic Stem Cell Transplantation

Alloreactivity forms the basis of allogeneic hematopoietic cell transplantation (HCT), with donor-derived T cell response to recipient antigens mediating clinical responses either in part or entirely. These encompass the different manifestations of graft-versus-host disease (GVHD), infection risk, and disease response. While the latter is contingent on disease biology and thus may be less predictable, the former 2 manifestations are more likely to be directly proportional to the magnitude of donor-derived T cell recovery. Herein we explore the quantitative aspects of immune cell recovery following allogeneic HCT and clinical outcomes in 2 cohorts of HLA-matched allograft recipients who received rabbit antithymocyte globulin (ATG) on different schedules (days -9 to -7 versus days -3 to -1). Monocyte as well as donor-derived T cell (ddCD3) recovery was superior in those given ATG early in the course of disease (days -9/-7). This difference was related to a more rapid rate of ddCD3 recovery, driven largely by CD3⁺/CD8⁺ cells in the first month post-transplantation. Early monocyte recovery was associated with later T cell recovery and improved survival. In contrast, rapid and early ddCD3 expansion out of proportion to monocyte recovery was associated with a high likelihood of acute GVHD and poor survival. This analytic methodology demonstrates that modeling "early-term immune reconstitution" following HCT yields insights that may be useful in the management of post-transplantation immunosuppression and adaptive cellular therapy to optimize clinical outcomes. © 2021 American Society for Transplantation and Cellular Therapy. Published by Elsevier Inc.

Comprehensive Prognostication in Critically Ill Pediatric Hematopoietic Cell Transplant Patients: Results from Merging the Center for International Blood and Marrow Transplant Research (CIBMTR) and Virtual Pediatric Systems (VPS) Registries

Critically ill pediatric allogeneic hematopoietic cell transplant (HCT) patients may benefit from early and aggressive interventions aimed at reversing the progression of multiorgan dysfunction. Therefore, we evaluated 25 early risk factors for pediatric intensive care unit (PICU) mortality to improve mortality prognostication. We merged the Virtual Pediatric Systems and Center for International Blood and Marrow Transplant Research databases and analyzed 936 critically ill patients ≤21 years of age who had undergone allogeneic HCT and subsequently required PICU admission between January 1, 2009, and December 31, 2014. Of 1532 PICU admissions, the overall PICU mortality rate was 17.4% (95% confidence interval [CI], 15.6% to 19.4%) but was significantly higher for patients requiring mechanical ventilation (44.0%), renal replacement therapy (56.1%), or extracorporeal life support (77.8%). Mortality estimates increased significantly the longer that patients remained in the PICU. Of 25 HCT- and PICU-specific characteristics available at or near the time of PICU admission, moderate/severe pre-HCT renal injury, pre-HCT recipient cytomegalovirus seropositivity, <100-day interval between HCT and PICU admission, HCT for underlying acute myeloid leukemia, and greater admission organ dysfunction as approximated by the Pediatric Risk of Mortality 3 score were each independently associated with PICU mortality. A multivariable model using these components identified that patients in the top quartile of risk had 3 times greater mortality than other patients (35.1% versus 11.5%, P < .001, classification accuracy 75.2%; 95% CI, 73.0% to 77.4%). These data improve our working knowledge of the factors influencing the progression of critical illness in pediatric allogeneic HCT patients. Future investigation aimed at mitigating the effect of these risk factors is warranted.

Reduced plasma ascorbic acid levels in recipients of myeloablative conditioning and hematopoietic cell transplantation

Hematopoietic cell transplantation (HCT) conditioned using myeloablative conditioning (MAC) is complicated by end organ injury due to endothelial dysfunction and graft versus host disease. Mucositis and oxidant injury results in micronutrient deficiency. Ascorbic acid (AA) levels were measured in 15 patients undergoing HCT conditioned with MAC (11 allogeneic and four autologous HCT). Ascorbate levels declined postconditioning to 27.3 μMol/L (±14.1) by day 0 (P = .03 compared with pretransplant baseline), reaching a nadir level of 21.5 (±13.8) on day 14 (P = .003) post-transplant. Patients undergoing allogeneic HCT continued to have low AA levels to day 60 post-transplant. The role of AA in maintaining endothelial function and hematopoietic as well as T-cell recovery is provided, developing the rationale for repletion of vitamin C following HCT.

T Cell Repertoire Evolution after Allogeneic Bone Marrow Transplantation: An Organizational Perspective

High-throughput sequencing (HTS) of human T cell receptors has revealed a high level of complexity in the T cell repertoire, which makes it difficult to correlate T cell reconstitution with clinical outcomes. The associations identified thus far are of a broadly statistical nature, precluding precise modeling of outcomes based on T cell repertoire development following bone marrow transplantation (BMT). Previous work has demonstrated an inherent, mathematically definable order observed in the T cells from a diverse group of donors, which is perturbed in recipients following BMT. In this study, T cell receptor (TCR)-β sequences from HLA-matched related donor and recipient pairs are analyzed to further develop this methodology. TCR-β sequencing from unsorted and sorted T cell subsets isolated from the peripheral blood samples of BMT donors and recipients show conservation and symmetry of VJ segment usage in the clonal frequencies, linked to the organization of the gene segments along the TCR locus. This TCR-β VJ segment translational symmetry is preserved post-transplantation and even in cases of acute graft-versus-host disease (aGVHD), suggesting that GVHD occurrence represents a polyclonal donor T cell response to recipient antigens. The complexity of the repertoire is significantly diminished after BMT, and the T cell clonal hierarchy is altered post-transplantation. Low-frequency donor clones tended to take on a higher rank in the recipients following BMT, especially in patients with aGVHD. Over time, the repertoire evolves to a more donor-like state in the recipients who did not develop GVHD as opposed to those who did. The results presented here support new methods of quantifying and characterizing post-transplantation T cell repertoire reconstitution.

A dynamical systems perspective on chimeric antigen receptor T-cell dosing

Chimeric antigen receptor T cells (CAR T cells) are dosed similarly to donor lymphocyte infusions following hematopoietic cell transplantation. However, the mechanism driving proliferation in CAR T cells is distinct from conventional T cells. As such there are quantitative differences in the antigen response of these engineered cells when compared with conventional T cells. In this perspective paper the logistic equation of growth is used to develop a mathematical basis for understanding the difference between CAR T cell and conventional T cell response to antigen burden.

A Review on the Nonlinear Dynamical System Analysis of Electrocardiogram Signal

Electrocardiogram (ECG) signal analysis has received special attention of the researchers in the recent past because of its ability to divulge crucial information about the electrophysiology of the heart and the autonomic nervous system activity in a noninvasive manner. Analysis of the ECG signals has been explored using both linear and nonlinear methods. However, the nonlinear methods of ECG signal analysis are gaining popularity because of their robustness in feature extraction and classification. The current study presents a review of the nonlinear signal analysis methods, namely, reconstructed phase space analysis, Lyapunov exponents, correlation dimension, detrended fluctuation analysis (DFA), recurrence plot, Poincaré plot, approximate entropy, and sample entropy along with their recent applications in the ECG signal analysis.

Reconstituting donor T cells increase their biomass following hematopoietic stem cell transplantation

In this study, we used a rapid, highly-sensitive, single-cell biomass measurement method, Live Cell Interferometry (LCI), to measure biomass in populations of CD3 + T cells isolated from hematopoietic stem cell transplant (SCT) patients at various times pre- and post-transplant (days 0–100).

Dynamical system modeling to simulate donor T cell response to whole exome sequencing-derived recipient peptides: Understanding randomness in alloreactivity incidence following stem cell transplantation

Quantitative relationship between the magnitude of variation in minor histocompatibility antigens (mHA) and graft versus host disease (GVHD) pathophysiology in stem cell transplant (SCT) donor-recipient pairs (DRP) is not established. In order to elucidate this relationship, whole exome sequencing (WES) was performed on 27 HLA matched related (MRD), & 50 unrelated donors (URD), to identify nonsynonymous single nucleotide polymorphisms (SNPs). An average 2,463 SNPs were identified in MRD, and 4,287 in URD DRP (p<0.01); resulting peptide antigens that may be presented on HLA class I molecules in each DRP were derived in silico (NetMHCpan ver2.0) and the tissue expression of proteins these were derived from determined (GTex). MRD DRP had an average 3,670 HLA-binding-alloreactive peptides, putative mHA (pmHA) with an IC50 of <500 nM, and URD, had 5,386 (p<0.01). To simulate an alloreactive donor cytotoxic T cell response, the array of pmHA in each patient was considered as an operator matrix modifying a hypothetical cytotoxic T cell clonal vector matrix; each responding T cell clone’s proliferation was determined by the logistic equation of growth, accounting for HLA binding affinity and tissue expression of each alloreactive peptide. The resulting simulated organ-specific alloreactive T cell clonal growth revealed marked variability, with the T cell count differences spanning orders of magnitude between different DRP. Despite an estimated, uniform set of constants used in the model for all DRP, and a heterogeneously treated group of patients, higher total and organ-specific T cell counts were associated with cumulative incidence of moderate to severe GVHD in recipients. In conclusion, exome wide sequence differences and the variable alloreactive peptide binding to HLA in each DRP yields a large range of possible alloreactive donor T cell responses. Our findings also help understand the apparent randomness observed in the development of alloimmune responses.

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.

Methodologic Considerations in the Application of Next-Generation Sequencing of Human TRB Repertoires for Clinical Use

Next-generation sequencing (NGS) of immune receptors has become a standard tool to assess minimal residual disease (MRD) in patients treated for lymphoid malignancy, and it is being used to study the T-cell repertoire in many clinical settings. To better understanding the potential clinical utility and limitations of this application outside of MRD, we developed a BIOMED-2 primer-based NGS method and characterized its performance in controls and patients with graft-versus-host disease (GVHD) after allogeneic hematopoietic transplant. For controls and patients with GVHD, replicate sequencing of the same T-cell receptor β (TRB) libraries was highly reproducible. Higher variability was observed in sequencing of different TRB libraries made from the same DNA stock. Variability was increased in patients with GVHD compared with controls; patients with GVHD also had lower diversity than controls. In the T-cell repertoire of a healthy person, approximately 99.6% of the CDR3 clones were in low abundance, with frequency <10^-3. A single library could identify >93% of the clones with frequency ≥10^-3 in the repertoire. Sequencing in duplicate increased the average detection rate to >97%. This work demonstrates that NGS reliably and robustly characterizes TRB populations in healthy individuals and patients with GVHD with frequency ≥10^-3 and provides a methodologic framework for applying NGS immune repertoire methods to clinical testing applications beyond MRD.

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.

Donor lymphocyte infusion after allogeneic hematopoietic stem cell transplantation

Donor lymphocyte infusion, a rescue therapy after hematopoietic stem cell transplantation, has been increasingly adopted, as modalities of stem cell transplantation have widened. First described as donor lymphocyte transfusion or cell therapy, it consists of infusion of donor lymphocytes, collected in steady state or after growth factor enhancement. As in literature the most used name is donor lymphocyte infusion, we'll adopt it here. Its most striking efficacy is observed in patients with chronic myelogenous leukemia, who relapsed after allogeneic stem cells transplantation. However, graft-versus-host disease, its main complication, may still hamper its feasibility.