Exome chip analyses identify genes affecting mortality after HLA-matched unrelated-donor blood and marrow transplantation

Associations of minor histocompatibility antigens with outcomes following allogeneic hematopoietic cell transplantation

The role of minor histocompatibility antigens (mHAs) in mediating graft versus leukemia and graft versus host disease (GvHD) following allogeneic hematopoietic cell transplantation (alloHCT) is recognized but not well-characterized. By implementing improved methods for mHA prediction in two large patient cohorts, this study aimed to comprehensively explore the role of mHAs in alloHCT by analyzing whether (1) the number of predicted mHAs, or (2) individual mHAs are associated with clinical outcomes. The study population consisted of 2249 donor-recipient pairs treated for acute myeloid leukemia and myelodysplastic syndrome with alloHCT. A Cox proportional hazard model showed that patients with a class I mHA count greater than the population median had an increased hazard of GvHD mortality (hazard ratio [HR] = 1.39, 95% confidence interval [CI] = 1.01, 1.77, p = .046). Competing risk analyses identified the class I mHAs DLRCKYISL (GSTP), WEHGPTSLL (CRISPLD2), and STSPTTNVL (SERPINF2) were associated with increased GVHD mortality (HR = 2.84, 95% CI = 1.52, 5.31, p = .01), decreased leukemia-free survival (LFS) (HR = 1.94, 95% CI = 1.27, 2.95, p = .044), and increased disease-related mortality (DRM) (HR = 2.32, 95% CI = 1.5, 3.6, p = .008), respectively. One class II mHA YQEIAAIPSAGRERQ (TACC2) was associated with increased risk of treatment-related mortality (TRM) (HR = 3.05, 95% CI = 1.75, 5.31, p = .02). WEHGPTSLL and STSPTTNVL were both present within HLA haplotype B*40:01-C*03:04 and showed a positive dose-response relationship with increased all-cause mortality and DRM and decreased LFS, indicating these two mHAs contribute to the risk of mortality in an additive manner. Our study reports the first large-scale investigation of the associations of predicted mHA peptides with clinical outcomes following alloHCT.

Shared graft-versus-leukemia minor histocompatibility antigens in DISCOVeRY-BMT

T-cell responses to minor histocompatibility antigens (mHAs) mediate graft-versus-leukemia (GVL) effects and graft-versus-host disease (GVHD) in allogeneic hematopoietic cell transplantation. Therapies that boost T-cell responses improve allogeneic hematopoietic cell transplant (alloHCT) efficacy but are limited by concurrent increases in the incidence and severity of GVHD. mHAs with expression restricted to hematopoietic tissue (GVL mHAs) are attractive targets for driving GVL without causing GVHD. Prior work to identify mHAs has focused on a small set of mHAs or population-level single-nucleotide polymorphism-association studies. We report the discovery of a large set of novel GVL mHAs based on predicted immunogenicity, tissue expression, and degree of sharing among donor-recipient pairs (DRPs) in the DISCOVeRY-BMT data set of 3231 alloHCT DRPs. The total number of predicted mHAs varied by HLA allele, and the total number and number of each class of mHA significantly differed by recipient genomic ancestry group. From the pool of predicted mHAs, we identified the smallest sets of GVL mHAs needed to cover 100% of DRPs with a given HLA allele. We used mass spectrometry to search for high-population frequency mHAs for 3 common HLA alleles. We validated 24 predicted novel GVL mHAs that are found cumulatively within 98.8%, 60.7%, and 78.9% of DRPs within DISCOVeRY-BMT that express HLA-A∗02:01, HLA-B∗35:01, and HLA-C∗07:02, respectively. We confirmed the immunogenicity of an example novel mHA via T-cell coculture with peptide-pulsed dendritic cells. This work demonstrates that the identification of shared mHAs is a feasible and promising technique for expanding mHA-targeting immunotherapeutics.

Liquid Biopsy, a Potential New Detection Method in Heart Allograft Rejection

Allografts rejection remains the most important reason causing allograft dysfunction in heart transplantation recipients. Currently, the golden standard for detecting graft rejection is endomyocardial biopsy (EMB). As a new noninvasive technique, liquid biopsy emerges along with the great developments of droplet-based digital PCR and the various optimizations of next-generation sequencing technologies, which is also cheaper than EMB. This review introduces several types of liquid biopsy and its application in heart transplantation.

Multi-faceted role of LRP1 in the immune system

Graft versus host disease (GVHD) represents the major complication after allogeneic hematopoietic stem cell transplantation (Allo-SCT). GVHD-prone patients rely on GVHD prophylaxis (e.g. methotrexate) and generalized anti-GVHD medical regimen (glucocorticoids). New anti-GVHD therapy strategies are being constantly explored, however there is an urgent need to improve current treatment, since GVHD-related mortality reaches 22% within 5 years in patients with chronic GVHD. This review is an attempt to describe a very well-known receptor in lipoprotein studies - the low-density lipoprotein receptor related protein 1 (LRP1) - in a new light, as a potential therapeutic target for GVHD prevention and treatment. Our preliminary studies demonstrated that LRP1 deletion in donor murine T cells results in significantly lower GVHD-related mortality in recipient mice with MHC (major histocompatibility complex) -mismatched HSCT. Given the importance of T cells in the development of GVHD, there is a significant gap in scientific literature regarding LRP1's role in T cell biology. Furthermore, there is limited research interest and publications on this classical receptor molecule in other immune cell types. Herein, we endeavor to summarize existing knowledge about LRP1's role in various immune cells to demonstrate the possibility of this receptor to serve as a novel target for anti-GVHD treatment.

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.

Novel genetic variants associated with mortality after unrelated donor allogeneic hematopoietic cell transplantation

BACKGROUND

Identification of non-human leukocyte antigen (HLA) genetic risk factors could improve survival after allogeneic blood or marrow transplant (BMT) through matching at additional loci or individualizing risk prediction. We hypothesized that non-HLA loci contributed significantly to 1-year overall survival (OS), disease related mortality (DRM) or transplant related mortality (TRM) after unrelated donor (URD)BMT.

METHODS

We performed a genome-wide association study (GWAS) in 2,887 acute myeloid leukemia (AML), myelodysplastic syndrome (MDS) and acute lymphoblastic leukemia (ALL) patients and their ≥8/8 HLA-matched URDs comprising two independent cohorts treated from 2000-2011.

FINDINGS

Using meta-analyses of both cohorts, genome-wide significant associations (p < 5 × 10^-8) were identified in: recipient genomes with OS at MBNL1 (rs9990017, HR = 1.4, 95% CI 1.24-1.56, p = 3.3 × 10^-8) and donor-recipient genotype mismatch with OS at LINC02774 (rs10927108, HR = 1.34, 95% CI 1.21-1.48, p = 2.0 × 10^-8); donor genomes with DRM at PCNX4 (rs79076914, HR = 1.7, 95% CI 1.41-2.05, p = 3.15 × 10^-8), LINC01194 (rs79498125, HR = 1.86, 95% CI 1.49-2.31, p = 2.84 × 10^-8), ARID5B (rs2167710, HR = 1.5, 95% CI 1.31-1.73, p = 6.9 × 10^-9) and CT49 (rs32250, HR = 1.44, 95% CI1.26-1.64, p = 2.6 × 10^-8); recipient genomes at PILRB with TRM (rs141591562, HR = 2.33, 95% CI 1.74-3.12, p = 1.26 × 10^-8) and donor-recipient genotype mismatch between EPGN and MTHF2DL with TRM (rs75868097, HR = 2.66, 95% CI 1.92-3.58, p = 4.6 × 10^-9). Results publicly available at https://fuma.ctglab.nl/browse.

INTERPRETATION

These data provide the first evidence that non-HLA common genetic variation at novel loci with biochemical function significantly impacts 1-year URD-BMT survival. Our findings have implications for donor selection, could guide treatment strategies and provide individualized risk prediction after future validation and functional studies.

FUNDING

This project was funded by grants from the National Institutes of Health, USA.

Genetic susceptibility to acute graft versus host disease in pediatric patients undergoing HSCT

The most frequent complication of allogeneic hematopoietic stem cell transplantation is acute Graft versus Host Disease (aGVHD). Proliferation and differentiation of donor T cells initiate inflammatory response affecting the skin, liver, and gastrointestinal tract. Besides recipient-donor HLA disparities, disease type, and the conditioning regimen, variability in the non-HLA genotype have an impact on aGVHD onset, and genetic variability of key cytokines and chemokines was associated with increased risk of aGVHD. To get further insight into the recipient genetic component of aGVHD grades 2-4 in pediatric patients, we performed an exome-wide association study in a discovery cohort (n = 87). Nine loci sustained correction for multiple testing and were analyzed in a validation group (n = 168). Significant associations were replicated for ERC1 rs1046473, PLEK rs3816281, NOP9 rs2332320 and SPRED1 rs11634702 variants through the interaction with non-genetic factors. The ERC1 variant was significant among patients that received the transplant from HLA-matched related individuals (p = 0.03), bone marrow stem cells recipients (p = 0.007), and serotherapy-negative patients (p = 0.004). NOP9, PLEK, and SPRED1 effects were modulated by stem cell source, and serotherapy (p < 0.05). Furthermore, ERC1 and PLEK SNPs correlated with aGVHD 3-4 independently of non-genetic covariates (p = 0.02 and p = 0.003). This study provides additional insight into the genetic component of moderate to severe aGVHD.

Neither Donor nor Recipient Mitochondrial Haplotypes Are Associated with Unrelated Donor Transplant Outcomes: A Validation Study from the CIBMTR

Graft-versus-host-disease (GVHD) is a multistep process that involves T-cell recognition and priming toward alloantigen, expansion, acquisition of effector function, and repeated tissue injury, resulting in clinical manifestations of the disease. All of these processes have considerable metabolic demands and understanding the key role of mitochondria in cellular metabolism as it relates to GVHD has increased significantly. Mitochondrial DNA (mtDNA) haplotypes have been linked to functional differences in vitro, suggesting they have functional differences at an organismal level. We previously used mtDNA typing to assess the impact of mtDNA haplotypes on outcomes of ~400 allo-HCT patients. This pilot study identified uncommon mtDNA haplotypes potentially associated with inferior outcomes. We sought to validate pilot findings of associations between donor and recipient mitochondrial haplotypes and transplant outcome. We examined a cohort of 4143 donor-recipient pairs obtained from the Center for International Blood and Marrow Transplant Research. MtDNA was extracted from whole blood or peripheral blood mononuclear cells from donors and recipients and sequenced to discern haplotype. We used multiple regression analysis to examine the independent association of mtDNA haplotype with overall survival and grade III-IV acute GVHD (aGVHD) adjusting for known risk factors for poor transplant outcome. Neither recipient nor donor mtDNA haplotype reached groupwise significance for overall survival (P =.26 and .39, respectively) or grade III-IV aGVHD (P = .68 and.57, respectively). Adjustment for genomically determined ancestry in the subset of donor-recipient pairs for which this was available did not materially change results. We conclude that our original finding was due to chance in a small sample size and that there is essentially no evidence that mtDNA haplotype or haplotype mismatch contributes to risk of serious outcomes after allogeneic transplantation.

Genome-Wide Association Analyses Identify Variants in IRF4 Associated With Acute Myeloid Leukemia and Myelodysplastic Syndrome Susceptibility

The role of common genetic variation in susceptibility to acute myeloid leukemia (AML), and myelodysplastic syndrome (MDS), a group of rare clonal hematologic disorders characterized by dysplastic hematopoiesis and high mortality, remains unclear. We performed AML and MDS genome-wide association studies (GWAS) in the DISCOVeRY-BMT cohorts (2,309 cases and 2,814 controls). Association analysis based on subsets (ASSET) was used to conduct a summary statistics SNP-based analysis of MDS and AML subtypes. For each AML and MDS case and control we used PrediXcan to estimate the component of gene expression determined by their genetic profile and correlate this imputed gene expression level with risk of developing disease in a transcriptome-wide association study (TWAS). ASSET identified an increased risk for de novo AML and MDS (OR = 1.38, 95% CI, 1.26-1.51, Pmeta = 2.8 × 10-12) in patients carrying the T allele at s12203592 in Interferon Regulatory Factor 4 (IRF4), a transcription factor which regulates myeloid and lymphoid hematopoietic differentiation. Our TWAS analyses showed increased IRF4 gene expression is associated with increased risk of de novo AML and MDS (OR = 3.90, 95% CI, 2.36-6.44, Pmeta = 1.0 × 10-7). The identification of IRF4 by both GWAS and TWAS contributes valuable insight on the role of genetic variation in AML and MDS susceptibility.

Recipient and donor genetic variants associated with mortality after allogeneic hematopoietic cell transplantation

Many studies have suggested that genetic variants in donors and recipients are associated with survival-related outcomes after allogeneic hematopoietic cell transplantation (HCT), but these results have not been confirmed. Therefore, the utility of testing genetic variants in donors and recipients for risk stratification or understanding mechanisms leading to mortality after HCT has not been established. We tested 122 recipient and donor candidate variants for association with nonrelapse mortality (NRM) and relapse mortality (RM) in a cohort of 2560 HCT recipients of European ancestry with related or unrelated donors. Associations discovered in this cohort were tested for replication in a separate cohort of 1710 HCT recipients. We found that the donor rs1051792 A allele in MICA was associated with a lower risk of NRM. Donor and recipient rs1051792 genotypes were highly correlated, making it statistically impossible to determine whether the donor or recipient genotype accounted for the association. Risks of grade 3 to 4 graft-versus-host disease (GVHD) and NRM in patients with grades 3 to 4 GVHD were lower with donor MICA-129Met but not with MICA-129Val, implicating MICA-129Met in the donor as an explanation for the decreased risk of NRM after HCT. Our analysis of candidate variants did not show any other association with NRM or RM. A genome-wide association study did not identify any other variants associated with NRM or RM.

Pre-HCT mosaicism increases relapse risk and lowers survival in acute lymphoblastic leukemia patients post-unrelated HCT

Pre-HCT mosaicism is related to increased relapse risk and lower survival after unrelated HCT, independent of cytogenetics at diagnosis. Pre-HCT mosaicism could be a useful clinical tool to guide risk stratification in acute lymphoblastic leukemia patients.

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.

Multiple functional variants in the IL1RL1 region are pretransplant markers for risk of GVHD and infection deaths

Graft-versus-host disease (GVHD) and infections are the 2 main causes of death without relapse after allogeneic hematopoietic cell transplantation (HCT). Elevated soluble serum simulation-2 (sST2), the product of IL1RL1 in plasma/serum post-HCT, is a validated GVHD biomarker. Hundreds of SNPs at 2q12.1 have been shown to be strongly associated with sST2 concentrations in healthy populations. We therefore hypothesized that the donor genetic variants in IL1RL1 correlate with sST2 protein levels associated with patient survival outcomes after HCT. We used DISCOVeRY-BMT (Determining the Influence of Susceptibility Conveying Variants Related to 1-Year Mortality after Blood and Marrow Transplantation), a genomic study of >3000 donor-recipient pairs, to inform our hypothesis. We first measured pre-HCT plasma/serum sST2 levels in a subset of DISCOVeRY-BMT donors (n = 757) and tested the association of donor sST2 levels with donor single nucleotide polymorphisms (SNPs) in the 2q12.1 region. Donor SNPs associated with sST2 levels were then tested for association with recipient death caused by acute GVHD (aGVHD)-, infection-, and transplant-related mortality in cohorts 1 and 2. Meta-analyses of cohorts 1 and 2 were performed using fixed-effects inverse variance weighting, and P values were corrected for multiple comparisons. Donor risk alleles in rs22441131 (P _meta = .00026) and rs2310241 (P _meta = .00033) increased the cumulative incidence of aGVHD death up to fourfold and were associated with high sST2 levels. Donor risk alleles at rs4851601 (P _meta = 9.7 × 10^-7), rs13019803 (P _meta = 8.9 × 10^-6), and rs13015714 (P _meta = 5.3 × 10^-4) increased cumulative incidence of infection death to almost sevenfold and were associated with low sST2 levels. These functional variants are biomarkers of infection or aGVHD death and could facilitate donor selection, prophylaxis, and a conditioning regimen to reduce post-HCT mortality.