Computational Analysis of HLA-presentation of Non-synonymous Recipient Mismatches Indicates Effect on the Risk of Chronic Graft-vs.-Host Disease After Allogeneic HSCT

Comparable relapse incidence after unrelated allogeneic stem cell transplantation with post-transplant cyclophosphamide versus conventional anti-graft versus host disease prophylaxis in patients with acute myeloid leukemia: A study on behalf of the Acute Leukemia Working Party of the European Society for Blood and Marrow Transplantation

We compared relapse incidence (RI) post-unrelated transplantation with post-transplant cyclophosphamide (PTCy) versus no PTCy graft-versus-host disease (GVHD) prophylaxis, in 7049 acute myeloid leukemia (AML) patients in remission, 707 with PTCy, and 6342 without (No PTCy). The patients in the PTCy group were younger, 52.7 versus 56.6 years (p < .001). There were more 9/10 donors in the PTCy group, 33.8% versus 16.4% (p < .001), and more received myeloablative conditioning, 61.7% versus 50.2% (p < .001). In the No PTCy group, 87.7% of patients received in vivo T-cell depletion. Neutrophil and platelet engraftment were lower in the PTCy versus No PTCy group, 93.8% and 80.9% versus 97.6% and 92.6% (p < .001). RI was not significantly different in the PTCy versus the No PTCy group, hazard ratio (HR) of 1.11 (95% confidence interval [CI] 0.9-1.37) (p = .31). Acute GVHD grades II-IV and III-IV, were significantly lower in the PTCy versus the No PTCy group, HR of 0.74 (95% CI 0.59-0.92, p = .007) and HR = 0.56 (95% CI 0.38-0.83, p = .004), as were total and extensive chronic GVHD, HRs of 0.5 (95% CI 0.41-0.62, p < .001) and HR = 0.31 (95% CI 0.22-0.42, p < .001). Non-relapse mortality (NRM) was significantly lower with PTCy versus the No PTCy group, HR of 0.67 (95% CI 0.5-0.91, p = .007). GVHD-free, relapse-free survival (GRFS) was higher in the PTCy versus the No PTCy group, HR of 0.69 (95% CI 0.59-0.81, p = .001). Leukemia-free survival (LFS) and overall survival (OS) did not differ between the groups. In summary, we observed comparable RI, OS, and LFS, significantly lower incidences of GVHD and NRM, and significantly higher GRFS in AML patients undergoing unrelated donor-hematopoietic stem cell transplantation with PTCy versus No PTCy GVHD prophylaxis.

Identification of single nucleotide polymorphisms (SNPs) associated with chronic graft-versus-host disease in patients undergoing allogeneic hematopoietic cell transplantation

INTRODUCTION

Chronic graft-versus-host disease (cGVHD) is a debilitating side effect of allogeneic hematopoietic cell transplantation (HCT), affecting the quality of life of patients. We used whole exome sequencing to identify candidate SNPs and complete a multi-marker gene-level analysis using a cohort of cGVHD( +) (N = 16) and cGVHD( -) (N = 66) HCT patients.

METHODS

Saliva samples were collected from HCT patients (N = 82) pre-conditioning in a multi-center study from March 2011 to May 2018. Exome sequencing was performed and FASTQ files were processed for sequence alignments. Significant SNPs were identified by logistic regression using PLINK2v3.7 and Fisher's exact test. One cGVHD( -) patient sample was excluded from further analysis since no SNP was present in at least 10% of the sample population. The FUMA platform's SNP2GENE was utilized to annotate SNPs and generate a MAGMA output. Chromatin state visualization of lead SNPs was completed using Epilogos tool. FUMA's GENE2FUNC was used to obtain gene function and tissue expression from lead genomic loci.

RESULTS

Logistic regression classified 986 SNPs associated with cGVHD( +). SNP2GENE returned three genomic risk loci, four lead SNPs, 48 candidate SNPs, seven candidate GWAS tagged SNPs, and four mapped genes. Fisher's exact test identified significant homozygous genotypes of four lead SNPs (p < 0.05). GENE2FUNC analysis of multi-marker SNP sets identified one positional gene set including lead SNPs for KANK1 and KDM4C and two curated gene sets including lead SNPs for PTPRD, KDM4C, and/or KANK1.

CONCLUSIONS

Our data suggest that SNPs in three genes located on chromosome 9 confer genetic susceptibility to cGVHD in HCT patients. These genes modulate STAT3 expression and phosphorylation in cancer pathogenesis. The findings may have implications in the modulation of pathways currently targeted by JAK inhibitors in cGVHD clinical trials.

A Model of Minor Histocompatibility Antigens in Allogeneic Hematopoietic Cell Transplantation

Minor histocompatibility antigens (mHAg) composed of peptides presented by HLA molecules can cause immune responses involved in graft-versus-host disease (GVHD) and graft-versus-leukemia effects after allogeneic hematopoietic cell transplantation (HCT). The current study was designed to identify individual graft-versus-host genomic mismatches associated with altered risks of acute or chronic GVHD or relapse after HCT between HLA-genotypically identical siblings. Our results demonstrate that in allogeneic HCT between a pair of HLA-identical siblings, a mHAg manifests as a set of peptides originating from annotated proteins and non-annotated open reading frames, which i) are encoded by a group of highly associated recipient genomic mismatches, ii) bind to HLA allotypes in the recipient, and iii) evoke a donor immune response. Attribution of the immune response and consequent clinical outcomes to individual peptide components within this set will likely differ from patient to patient according to their HLA types.

LMP7 polymorphism may modify the presentation and clinical impact of minor histocompatibility antigens in matched related hematopoietic stem cell transplantation

Differential expression of minor histocompatibility antigens between the recipient and donor determines their disparity and can be modified by immunoproteasomes that regulate their processing and presentation. We examined the impact of HA-1 and HA-8 disparity, and immunoproteasome LMP7 polymorphism in 130 pairs. In multivariate analysis, HA-1 disparity showed a statistically significant association with an increased incidence of acute graft-versus-host disease (aGVHD) II-IV (p = 0.043, HR: 3.71, 95%CI = 1.04-13.26), while LMP7-Q/Q showed a trend toward increased incidence of aGVHD compared to LMP7-Q/K and K/K genotypes (p = 0.087, HR: 2.36, 95%CI = 0.88-6.31). All HA-1 and HA-8 disparate patients who developed aGVHD had the LMP7-Q/Q genotype. No significant association could be detected between HA-1, HA-8, or LMP7 and chronic GVHD, relapse-free survival (RFS), overall survival (OS), or transplant-related mortality (TRM). In conclusion, we suggested an association between the HA-1 disparity and the risk of developing aGVHD with a possible modifying effect of LMP7.

Review of Genetic Variation as a Predictive Biomarker for Chronic Graft-Versus-Host-Disease After Allogeneic Stem Cell Transplantation

Chronic graft-versus-host disease (cGvHD) is one of the major complications of allogeneic stem cell transplantation (HSCT). cGvHD is an autoimmune-like disorder affecting multiple organs and involves a dermatological rash, tissue inflammation and fibrosis. The incidence of cGvHD has been reported to be as high as 30% to 60% and there are currently no reliable tools for predicting the occurrence of cGvHD. There is therefore an important unmet clinical need for predictive biomarkers. The present review summarizes the state of the art for genetic variation as a predictive biomarker for cGvHD. We discuss three different modes of action for genetic variation in transplantation: genetic associations, genetic matching, and pharmacogenetics. The results indicate that currently, there are no genetic polymorphisms or genetic tools that can be reliably used as validated biomarkers for predicting cGvHD. A number of recommendations for future studies can be drawn. The majority of studies to date have been under-powered and included too few patients and genetic markers. Like in all complex multifactorial diseases, large collaborative genome-level studies are now needed to achieve reliable and unbiased results. Some of the candidate genes, in particular, CTLA4, HSPE, IL1R1, CCR6, FGFR1OP, and IL10, and some non-HLA variants in the HLA gene region have been replicated to be associated with cGvHD risk in independent studies. These associations should now be confirmed in large well-characterized cohorts with fine mapping. Some patients develop cGvHD despite very extensive immunosuppression and other treatments, indicating that the current therapeutic regimens may not always be effective enough. Hence, more studies on pharmacogenetics are also required. Moreover, all of these studies should be adjusted for diagnostic and clinical features of cGvHD. We conclude that future studies should focus on modern genome-level tools, such as machine learning, polygenic risk scores and genome-wide association study-transcription meta-analyses, instead of focusing on just single variants. The risk of cGvHD may be related to the summary level of immunogenetic differences, or whole genome histocompatibility between each donor-recipient pair. As the number of genome-wide analyses in HSCT is increasing, we are approaching an era where there will be sufficient data to incorporate these approaches in the near future.

Optimized Whole Genome Association Scanning for Discovery of HLA Class I-Restricted Minor Histocompatibility Antigens

Patients undergoing allogeneic stem cell transplantation as treatment for hematological diseases face the risk of Graft-versus-Host Disease as well as relapse. Graft-versus-Host Disease and the favorable Graft-versus-Leukemia effect are mediated by donor T cells recognizing polymorphic peptides, which are presented on the cell surface by HLA molecules and result from single nucleotide polymorphism alleles that are disparate between patient and donor. Identification of polymorphic HLA-binding peptides, designated minor histocompatibility antigens, has been a laborious procedure, and the number and scope for broad clinical use of these antigens therefore remain limited. Here, we present an optimized whole genome association approach for discovery of HLA class I minor histocompatibility antigens. T cell clones isolated from patients who responded to donor lymphocyte infusions after HLA-matched allogeneic stem cell transplantation were tested against a panel of 191 EBV-transformed B cells, which have been sequenced by the 1000 Genomes Project and selected for expression of seven common HLA class I alleles (HLA-A^∗01:01, A^∗02:01, A^∗03:01, B^∗07:02, B^∗08:01, C^∗07:01, and C^∗07:02). By including all polymorphisms with minor allele frequencies above 0.01, we demonstrated that the new approach allows direct discovery of minor histocompatibility antigens as exemplified by seven new antigens in eight different HLA class I alleles including one antigen in HLA-A^∗24:02 and HLA-A^∗23:01, for which the method has not been originally designed. Our new whole genome association strategy is expected to rapidly augment the repertoire of HLA class I-restricted minor histocompatibility antigens that will become available for donor selection and clinical use to predict, follow or manipulate Graft-versus-Leukemia effect and Graft-versus-Host Disease after allogeneic stem cell transplantation.

Meta-Analysis of Genome-Wide Association and Gene Expression Studies Implicates Donor T Cell Function and Cytokine Pathways in Acute GvHD

Graft-vs.-host disease (GvHD) is a major complication after allogeneic hematopoietic stem cell transplantation that causes mortality and severe morbidity. Genetic disparities in human leukocyte antigens between the recipient and donor are known contributors to the risk of the disease. However, the overall impact of genetic component is complex, and consistent findings across different populations and studies remain sparse. To gain a comprehensive understanding of the genes responsible for GvHD, we combined genome-wide association studies (GWAS) from two distinct populations with previously published gene expression studies on GvHD in a single gene-level meta-analysis. We hypothesized that genes driving GvHD should be associated in both data modalities and therefore could be detected more readily through their combined effects in the integrated analysis rather than in separate analyses. The meta-analysis yielded a total of 51 acute GvHD-associated genes (false detection rate [FDR] <0.1). In support of our hypothesis, this number was significantly higher than that in a permutation meta-analysis involving the whole data set, as well as in separate meta-analyses on the GWAS and gene expression data sets. The genes indicated by the meta-analysis were significantly enriched in 277 Gene Ontology terms (FDR < 0.05), such as T cell function and cytokine-mediated signaling pathways, and the results highlighted several established immune mediators, such as interleukins and JAK-STAT signaling, and presented TRAF6 and TERT as potential effector candidates. Altogether, the results support the chosen methodological approach, implicate a role of gene-level variation in donors' key immunological regulators predisposing patients to acute GVHD, and present potential targets for therapeutic intervention.