Novel HLA-DP region susceptibility loci associated with severe acute GvHD

A genome-wide association study on hematopoietic stem cell transplantation reveals novel genomic loci associated with transplant outcomes

Introduction

Data on genomic susceptibility for adverse outcomes after hematopoietic stem cell transplantation (HSCT) for recipients are scarce.

Methods

We performed a genome wide association study (GWAS) to identify genes associated with survival/mortality, relapse, and severe graft-versus-host disease (sGvHD), fitting proportional hazard and subdistributional models to data of n=1,392 recipients of European ancestry from three centres.

Results

The single nucleotide polymorphism (SNP) rs17154454, intronic to the neuronal growth guidant semaphorin 3C gene (SEMA3C), was genome-wide significantly associated with event-free survival (p=7.0x10-8) and sGvHD (p=7.5x10-8). Further associations were detected for SNPs in the Paxillin gene (PXN) with death without prior relapse or sGvHD, as well as for SNPs of the Plasmacytoma Variant Translocation 1 gene (PVT1, a long non-coding RNA gene), the Melanocortin 5 Receptor (MC5R) gene and the WW Domain Containing Oxidoreductase gene (WWOX), all associated with the occurrence of sGvHD. Functional considerations support the observed associations.

Discussion

Thus, new genes were identified, potentially influencing the outcome of HSCT.

Fifty novel HLA-DPB1 alleles identified in a UK cohort of unrelated hematopoietic cell donors and recipients

HLA-DPB1 is the classical HLA class II genes with the least recorded variation on the IPD-IMGT/HLA Database, suggesting the full extent of its diversity is perhaps yet to be characterized. Here, a full-gene typing strategy was employed to genotype a UK cohort of 1470 HCT recipients (n = 744) and donors (n = 726). In total, 2940 full-length HLA-DPB1 sequences were generated, comprising 193 distinct alleles. Of these, 107 sequences contained novel variation, totaling 49 unique intronic HLA-DPB1 alleles, and one coding variant (HLA-DPB1*1188:01). Full-gene sequencing resulted in zygosity changes for 129 individuals by identifying two distinct intronic variants of the same coding allele. We verified the existence of nine unconfirmed alleles and extended the sequence of two existing alleles on the IPD-IMGT/HLA Database.

Sequence complementarity between human noncoding RNAs and SARS-CoV-2 genes: What are the implications for human health?

Objectives

To investigate in silico the presence of nucleotide sequence complementarity between the RNA genome of Severe Acute Respiratory Syndrome CoronaVirus-2 (SARS-CoV-2) and human non-coding (nc)RNA genes.

Methods

The FASTA sequence (NC_045512.2) of each of the 11 SARS-CoV-2 isolate Wuhan-Hu-1 genes was retrieved from NCBI.nlm.nih.gov/gene and the Ensembl.org library interrogated for any base-pair match with human ncRNA genes. SARS-CoV-2 gene-matched human ncRNAs were screened for functional activity using bioinformatic analysis. Finally, associations between identified ncRNAs and human diseases were searched in GWAS databases.

Results

A total of 252 matches were found between the nucleotide sequence of SARS-CoV-2 genes and human ncRNAs. With the exception of two small nuclear RNAs, all of them were long non-coding (lnc)RNAs expressed mainly in testis and central nervous system under physiological conditions. The percentage of alignment ranged from 91.30% to 100% with a mean nucleotide alignment length of 17.5 ± 2.4. Thirty-three (13.09%) of them contained predicted R-loop forming sequences, but none of these intersected the complementary sequences of SARS-CoV-2. However, in 31 cases matches fell on ncRNA regulatory sites, whose adjacent coding genes are mostly involved in cancer, immunological and neurological pathways. Similarly, several polymorphic variants of detected non-coding genes have been associated with neuropsychiatric and proliferative disorders.

Conclusion

This pivotal in silico study shows that SARS-CoV-2 genes have Watson-Crick nucleotide complementarity to human ncRNA sequences, potentially disrupting ncRNA epigenetic control of target genes. It remains to be elucidated whether this could result in the development of human disease in the long term.

Genetic Associations with Immune-mediated Outcomes after Allogeneic Hematopoietic Cell Transplantation

A recipient variant correlated with DPB1 expression is associated with the risks of GVHD and relapse after unrelated HCT.

Larger cohorts would be needed to detect other genotypic associations with aGVHD, cGVHD, and relapse after allogeneic HCT.

Previous studies have identified more than 200 genetic variants associated with acute or chronic graft-versus-host disease (aGVHD; cGVHD) or recurrent malignancy after allogeneic hematopoietic cell transplantation (HCT). We tested these candidate donor and recipient variants in a cohort of 4270 HCT recipients of European ancestry and in subcohorts of 1827 sibling and 1447 unrelated recipients who had 10/10 HLA-A, B, C, DRB1, and DQB1-matched donors. We also carried out a genome-wide association study (GWAS) for these same outcomes. The discovery and replication analysis of candidate variants identified a group of closely linked recipient HLA-DPB1 single-nucleotide polymorphisms (SNPs) associated with an increased risk of aGVHD and a corresponding decreased risk of recurrent malignancy after unrelated HCT. These results reflect a correlation with the level of HLA-DPB1 expression previously shown to affect the risks of aGVHD and relapse in unrelated recipients. Our GWAS identified an association of cGVHD with a locus of X-linked recipient intron variants in NHS, a gene that regulates actin remodeling and cell morphology. Evaluation of this association in a second replication cohort did not confirm the original replication results, and we did not reach any definitive conclusion regarding the validity of this discovery. The cohort used for our study is larger than those used in most previous HCT studies but is smaller than those typically used for other genotype-phenotype association studies. Genomic and disease data from our study are available for further analysis in combination with data from other cohorts.

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.

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.

Chromosome Y-encoded antigens associate with acute graft-versus-host disease in sex-mismatched stem cell transplant

Allogeneic hematopoietic stem cell transplantation (allo-HCT) is a curative option for blood cancers, but the coupled effects of graft-versus-tumor and graft-versus-host disease (GVHD) limit its broader application. Outcomes improve with matching at HLAs, but other factors are required to explain residual risk of GVHD. In an effort to identify genetic associations outside the major histocompatibility complex, we conducted a genome-wide clinical outcomes study on 205 acute myeloid leukemia patients and their fully HLA-A-, HLA-B-, HLA-C-, HLA-DRB1-, and HLA-DQB1-matched (10/10) unrelated donors. HLA-DPB1 T-cell epitope permissibility mismatches were observed in less than half (45%) of acute GVHD cases, motivating a broader search for genetic factors affecting clinical outcomes. A novel bioinformatics workflow adapted from neoantigen discovery found no associations between acute GVHD and known, HLA-restricted minor histocompatibility antigens (MiHAs). These results were confirmed with microarray data from an additional 988 samples. On the other hand, Y-chromosome-encoded single-nucleotide polymorphisms in 4 genes (PCDH11Y, USP9Y, UTY, and NLGN4Y) did associate with acute GVHD in male patients with female donors. Males in this category with acute GVHD had more Y-encoded variant peptides per patient with higher predicted HLA-binding affinity than males without GVHD who matched X-paralogous alleles in their female donors. Methods and results described here have an immediate impact for allo-HCT, warranting further development and larger genomic studies where MiHAs are clinically relevant, including cancer immunotherapy, solid organ transplant, and pregnancy.

Hidden genomic MHC disparity between HLA-matched sibling pairs in hematopoietic stem cell transplantation

Matching classical HLA alleles between donor and recipient is an important factor in avoiding adverse immunological effects in HSCT. Siblings with no differences in HLA alleles, either due to identical-by-state or identical-by-descent status, are considered to be optimal donors. We carried out a retrospective genomic sequence and SNP analysis of 336 fully HLA-A, -B, -DRB1 matched and 14 partially HLA-matched sibling HSCT pairs to determine the level of undetected mismatching within the MHC segment as well as to map their recombination sites. The genomic sequence of 34 genes locating in the MHC region revealed allelic mismatching at 1 to 8 additional genes in partially HLA-matched pairs. Also, fully matched pairs were found to have mismatching either at HLA-DPB1 or at non-HLA region within the MHC segment. Altogether, 3.9% of fully HLA-matched HSCT pairs had large genomic mismatching in the MHC segment. Recombination sites mapped to certain restricted locations. The number of mismatched nucleotides correlated with the risk of GvHD supporting the central role of full HLA matching in HSCT. High-density genome analysis revealed that fully HLA-matched siblings may not have identical MHC segments and even single allelic mismatching at any classical HLA gene often implies larger genomic differences along MHC.