HLA-C expression levels define permissible mismatches in hematopoietic cell transplantation

Development of data-driven models for the flow cytometric crossmatch

HLA laboratories use virtual crossmatching (VXM) to predict recipient and donor compatibility using HLA antibody data and donor HLA type. Increasingly, transplant centers are utilizing VXM as the final compatibility determination prior to transplant. However, the VXM interpretation is based on HLA experience of individual transplant centers. This study developed data-driven algorithms that predicted flow cytometric crossmatch (FCXM) outcomes using HLA antibody mean fluorescent intensity (MFI) data and donor HLA typing without the need for human interpretation.Two algorithms were evaluated; an MFI Optimal-Threshold model and a Least-Squares-Fitting model. The Optimal-Threshold model correctly determined between 81.5% and 85.5% of T or B-cell responses. A class I antibody MFI threshold of 4670 was optimal for predicting T-cell response while an antibody MFI threshold of 6180 was optimal for predicting B-cell responses. HLA class I antibodies had a 1.47-fold greater influence on FCXM outcomes than class II antibodies. HLA-B antibodies influenced T and B-cell responses more than HLA-A or -C (-B > -A > -C). The Least-Squares-Fitting model increased accuracy to 94.1% and 88.8% for T and B-cell responses, respectively. The algorithms described here provide enhanced FCXM prediction and novel insights into the influence of specific HLA antibodies on the crossmatch outcome.

HLA Matching in Unrelated Stem Cell Transplantation up to Date

Unrelated hematopoietic stem cell transplantation (HSCT) has evolved from an experimental protocol to a potentially curative first-line treatment in certain disease instances. Factors enabling this transformation were the optimization of treatment protocols and supportive care as well as the availability of a large number of donors worldwide along with the higher quality and reliability of HLA typing. The main criterion for donor selection is HLA compatibility. In this review we discuss the current clinical evidence of HLA matching in unrelated HSCT. In this context, we address methodical aspects of transplantation immunobiology research and discuss the impact of locus and resolution of HLA differences. Furthermore, we address special constellations such as unidirectional mismatches or the presence of nonexpressed alleles as well as HLA alloimmunization and describe the perspective for HLA typing and matching strategies in the future, given the implementation of novel complete or near-complete gene typing approaches using next-generation sequencing short read technology, which are now entering the standard of clinical care.

Selection of unrelated donors and cord blood units for hematopoietic cell transplantation: guidelines from the NMDP/CIBMTR

Allogeneic hematopoietic cell transplantation involves consideration of both donor and recipient characteristics to guide the selection of a suitable graft. Sufficient high-resolution donor-recipient HLA match is of primary importance in transplantation with adult unrelated donors, using conventional graft-versus-host disease prophylaxis. In cord blood transplantation, optimal unit selection requires consideration of unit quality, cell dose and HLA-match. In this summary, the National Marrow Donor Program (NMDP) and the Center for International Blood and Marrow Transplant Research, jointly with the NMDP Histocompatibility Advisory Group, provide evidence-based guidelines for optimal selection of unrelated donors and cord blood units.

Factors affecting HLA expression: A review

The detection and semiquantitative measurement of circulating human leucocyte antigen (HLA)-specific antibodies is essential for the management of patients before and after transplantation. In addition, the pretransplant cross-match to assess the reactivity of recipient HLA antibody against donor lymphocytes has long been the gold standard to prevent hyperacute rejection. Whilst both of these tests assume that recipient HLA-specific antibody is the only variable in the assessment of transplant risk, this is not the case. Transplant immunologists recognize that some HLA antigens are expressed at levels a magnitude lower than others (e.g., HLA-C, HLA-DQ), but within loci, and between different cell types there are many factors that influence HLA expression in both resting and activated cells. HLA is not usually expressed without the specific promoter proteins NLRC5, for HLA class I, and CIITA, for class II. The quantity of HLA protein production is then affected by factors including promoter region polymorphisms, alternative exon splice sites, methylation and microRNA-directed degradation. Different loci are influenced by multiple combinations of these control mechanisms making prediction of HLA regulation difficult, but an ability to measure the cellular expression of each HLA antigen, in conjunction with knowledge of circulating HLA-specific antibody, would lead to a more informed algorithm to assess transplant risk.

Identification of the rs9277534 HLA-DP expression marker by next generation sequencing for the selection of unrelated donors for hematopoietic cell transplantation

Mismatching of an unrelated donor against a high-expression HLA-DPB1 recipient allele is associated with a high risk of graft-versus-host disease and mortality. The Seattle Cancer Care Alliance (SCCA) and Fred Hutchinson Cancer Research Center transplant program employs an algorithm to match for HLA-A, B, C, DRB1, DQB1 and DPB1 alleles (12/12) and to avoid, whenever possible, donor mismatching against a recipient high-expression HLA-DPB1 allele. HLA-DPB1 expression is associated with the rs9277534 A/G polymorphism located in the 3'UTR of the HLA-DPB1 gene. Next generation sequencing of HLA-DPB1 using the Illumina TruSight HLA V2 Sequencing Panel and Conexio Assign software analyses provides information on rs9277534 variants without the need for any additional SNP testing. Here we present the molecular location of rs9277534 in NGS data and discuss the challenges to resolve HLA-DPB1 ambiguities.

Integrate CRISPR/Cas9 for protein expression of HLA-B*38:68Q via precise gene editing

The determination of null- or low-expressed HLA alleles is clinically relevant in both hematopoietic stem cell transplantation and solid organ transplantation. We studied the expression level of a questionable (Q) HLA-B*38:68Q allele, which carries a 9-nucleotide (nt) deletion at codon 230–232 in exon 4 of HLA-B*38:01:01:01 using CRISPR/Cas9 gene editing technology. CRISPR/Cas9 gene editing of HLA-B*38:01:01:01 homozygous EBV B cell line resulted in one HLA-B*38:68Q/B*38:01:01:01 heterozygous and one HLA-B*38:68Q homozygous clone. Flow cytometric analysis of monoclonal anti-Bw4 antibody showed the protein expression of HLA-B*38:01:01:01 in homozygous cells was 2.2 fold higher than HLA-B*38:68Q/B*38:01:01:01 heterozygous cells, and the expression of HLA-B*38:68Q/B*38:01:01:01 heterozygous cells was over 2.0 fold higher than HLA-B*38:68Q homozygous cells. The HLA-B*38:68Q expression was further confirmed using anti-B38 polyclonal antibody. Similarly, the expression of the HLA-B*38:01:01:01 homozygous cells was 1.5 fold higher than that of HLA-B*38:68Q/B*38:01:01:01 heterozygous cells, and the HLA-B*38:68Q/B*38:01:01:01 heterozygous cells was over 1.6 fold higher than that of HLA-B*38:68Q homozygous cells. The treatment of HLA-B*38:68Q homozygous cells with IFN-γ significantly increased its expression. In conclusion, we demonstrate that HLA-B*38:68Q is a low-expressing HLA allele. The CRISPR/Cas9 technology is a useful tool to induce precise gene editing in HLA genes to enable the characterization of HLA gene variants on expression and function.

Expression estimation and eQTL mapping for HLA genes with a personalized pipeline

The HLA (Human Leukocyte Antigens) genes are well-documented targets of balancing selection, and variation at these loci is associated with many disease phenotypes. Variation in expression levels also influences disease susceptibility and resistance, but little information exists about the regulation and population-level patterns of expression. This results from the difficulty in mapping short reads originated from these highly polymorphic loci, and in accounting for the existence of several paralogues. We developed a computational pipeline to accurately estimate expression for HLA genes based on RNA-seq, improving both locus-level and allele-level estimates. First, reads are aligned to all known HLA sequences in order to infer HLA genotypes, then quantification of expression is carried out using a personalized index. We use simulations to show that expression estimates obtained in this way are not biased due to divergence from the reference genome. We applied our pipeline to the GEUVADIS dataset, and compared the quantifications to those obtained with reference transcriptome. Although the personalized pipeline recovers more reads, we found that using the reference transcriptome produces estimates similar to the personalized pipeline (r ≥ 0.87) with the exception of HLA-DQA1. We describe the impact of the HLA-personalized approach on downstream analyses for nine classical HLA loci (HLA-A, HLA-C, HLA-B, HLA-DRA, HLA-DRB1, HLA-DQA1, HLA-DQB1, HLA-DPA1, HLA-DPB1). Although the influence of the HLA-personalized approach is modest for eQTL mapping, the p-values and the causality of the eQTLs obtained are better than when the reference transcriptome is used. We investigate how the eQTLs we identified explain variation in expression among lineages of HLA alleles. Finally, we discuss possible causes underlying differences between expression estimates obtained using RNA-seq, antibody-based approaches and qPCR.

The level at which a gene is expressed can have important influence on the phenotype of an organism, including its predisposition to develop diseases. One way to estimate gene expression is by quantifying the abundance of RNA. RNA-seq has become the method of choice to provide such estimates at the genomewide scale. However, the application of RNA-seq to HLA genes —key players in the immune adaptive response— has remained a rarely explored approach. This is due to the problem of mapping bias, which causes deficient read alignment at genes which are very polymorphic and different from the reference genome. This has motivated approaches that replace the single reference genome with personalized sequences, comprised of the individual’s specific HLA genotype. Here we explore the use of computational frameworks to obtain reliable expression levels for HLA genes from RNA-seq datasets. We present a pipeline in which the quantification of HLA expression is carried out using methods which account for HLA diversity, avoiding the biases of standard approaches. We then evaluate the impact of this form of quantifying HLA expression on downstream analyses. The pipeline also allows us to integrate information on eQTLs with expression levels at the HLA allele-level, which can help disentangle different contributions to disease phenotypes and help understand the regulatory architecture at the HLA region.

The MHC in the era of next-generation sequencing: Implications for bridging structure with function

The MHC continues to have the most disease-associations compared to other regions of the human genome, even in the genome-wide association study (GWAS) and single nucleotide polymorphism (SNP) era. Analysis of non-coding variation and their impact on the level of expression of HLA allotypes has shed new light on the potential mechanisms underlying HLA disease associations and alloreactivity in transplantation. Next-generation sequencing (NGS) technology has the capability of delineating the phase of variants in the HLA antigen-recognition site (ARS) with non-coding regulatory polymorphisms. These relationships are critical for understanding the qualitative and quantitative implications of HLA gene diversity. This article summarizes current understanding of non-coding region variation of HLA loci, the consequences of regulatory variation on HLA expression, the role for evolution in shaping lineage-specific expression, and the impact of HLA expression on disease susceptibility and transplantation outcomes. A role for phased sequencing methods for the MHC, and perspectives for future directions in basic and applied immunogenetic studies of the MHC are presented.

Prevention of relapse after allogeneic hematopoietic cell transplantation by donor and cell source selection

Allogeneic hematopoietic cell transplantation (HCT) is the most established form of cancer immunotherapy and has been successfully applied for the treatment and cure of otherwise lethal neoplastic blood disorders. Cancer immune surveillance is mediated to a large extent by alloreactive T and natural killer (NK) cells recognizing genetic differences between patient and donor. Profound insights into the biology of these effector cells has been obtained over recent years and used for the development of innovative strategies for intelligent donor selection, aiming for improved graft-versus-leukemia effect without unmanageable graft-versus-host disease. The cellular composition of the stem cell source plays a major role in modulating these effects. This review summarizes the current state-of the-art of donor selection according to HLA, NK alloreactivity and stem cell source.

Dissecting Genetic Control of HLA-DPB1 Expression and Its Relation to Structural Mismatch Models in Hematopoietic Stem Cell Transplantation

HLA expression levels have been suggested to be genetically controlled by single nucleotide polymorphisms (SNP) in the untranslated regions (UTR), and expression variants have been associated with the outcome of chronic viral infection and hematopoietic stem cell transplantation (HSCT). In particular, the 3′UTR rs9277534-G/A SNP in HLA-DPB1 has been associated with graft-versus-host-disease after HSCT (Expression model); however its relevance in different immune cells and its mode of action have not been systematically addressed. In addition, there is a strong though not complete overlap between the rs9277534-G/A SNP and structural HLA-DPB1 T cell epitope (TCE) groups which have also been associated with HSCT outcome (TCE Structural model). Here we confirm and extend previous findings of significantly higher HLA-DPB1 expression in B cell lines, unstimulated primary B cells, and monocytes homozygous for rs9277534-G compared to those homozygous for rs9277534-A. However, these differences were abrogated by interferon-γ stimulation or differentiation into dendritic cells. We identify at least seven 3′UTR rs9277534-G/A haplotypes differing by a total of 37 SNP, also characterized by linkage to length variants of a short tandem repeat (STR) in intron 2 and TCE group assignment. 3′UTR mapping did not show any significant differences in post-transcriptional regulation assessed by luciferase assays between two representative rs9277534-G/A haplotypes for any of eight overlapping fragments. Moreover, no evidence for alternative splicing associated with the intron 2 STR was obtained by RT-PCR. In an exemplary cohort of 379 HLA-DPB1 mismatched donor-recipient pairs, risk prediction by the Expression model and the Structural TCE model was 36.7% concordant, with the majority of discordances due to non-applicability of the Expression model. HLA-DPB1 from different TCE groups expressed in the absence of the 3′UTR at similar levels by transfected HeLa cells elicited significantly different mean alloreactive CD4+ T-cell responses, as assessed by CD137 upregulation assays in 178 independent cultures. Taken together, our data provide new insights into the cell type-specific and mechanistic basis of the association between the rs9277534-G/A SNP and HLA-DPB1 expression, and show that, despite partial overlap between both models in HSCT risk-prediction, differential alloreactivity determined by the TCE structural model occurs independently from HLA-DPB1 differential expression.

Genetics of antigen processing and presentation

Immune response to disease requires coordinated expression of an army of molecules. The highly polymorphic MHC class I and class II molecules are key to control of specificity of antigen presentation. Processing of the antigen, to peptides or other moieties, requires other sets of molecules. For classical class I, this includes TAP peptide transporters, proteasome components and Tapasin, genes which are encoded within the MHC. Similarly, HLA-DO and -DM, which influence presentation by HLA class II molecules, are encoded in the MHC region. Analysis of MHC mutants, including point mutations and large deletions, has been central to understanding the roles of these genes. Mouse genetics has also played a major role. Many other genes have been identified including those controlling expression of HLA class I and class II at the transcriptional level. Another genetic approach that has provided insight has been the analysis of microorganisms, including viruses and bacteria that escape immune recognition by blocking these antigen processing and presentation pathways. Here, we provide a brief history of the genetic approaches, both traditional and modern, that have been used in the quest to understand antigen processing and presentation.

The prevalence of antibodies against the HLA-DRB3 protein in kidney transplantation and the correlation with HLA expression

Human leukocyte antigen (HLA)-DRB3 is a functional HLA class II gene, which has a limited allele diversity in the human population. Furthermore, the HLA-DRB3 gene is only present in a subset of individuals. Therefore, in organ transplantation, this HLA molecule is frequently mismatched between patient and graft donor and thus antibodies against this mismatched HLA molecule can develop. In this study, we aimed to evaluate the prevalence and reactivity of these antibodies and aimed to identify factors that underlie antibody formation against HLA-DRB3. We showed in our patient cohort that HLA-DRB3 antibodies are identified in about 7% of all patients that were screened with solid phase assays. In these assays, we observed multiple antibody reactivity patterns indicating that HLA-DRB3 harbours multiple epitopes. In those cases, where we succeeded at tracing back the induction of these antibodies to the molecular HLA typing of the immunogenic event, we noticed a different frequency of HLA-DRB1 allele groups in the donors as compared to a control group. To a certain extent this distribution (e.g. HLA-DRB1*11 individuals) could be linked to an altered expression level. However, it also appears that different HLA-DRB3 alleles (e.g. HLA-DRB3*01 group) vary in their immunogenicity without having an expression difference. In conclusion, our study provides information on the immunogenicity and reactivity patterns of antibodies against HLA-DRB3 in kidney transplantation, and it points towards the possibility of HLA expression as a factor underlying antibody formation.

In celebration of Ruggero Ceppellini: HLA in transplantation

The availability of hematopoietic cell transplantation as curative therapy for blood disorders has been dramatically improved through a better understanding of the human leukocyte antigen (HLA) barrier. Although a fully compatible unrelated donor is preferable, transplantation from donors with a limited degree of HLA mismatching is associated with acceptable outcomes in many cases. Research on the limits of HLA mismatching, and the features that define permissible HLA mismatches will continue to enable transplantation to be more broadly available to patients in need.

Clinical implications of HLA locus mismatching in unrelated donor hematopoietic cell transplantation: a meta-analysis

It remains controversial that the impacts of individual HLA locus mismatches on clinical outcomes of patients receiving unrelated-donor hematopoietic cell transplantation (HCT), as compared to HLA allele matched controls. We conducted a meta-analysis to address these issues. Four databases (PubMed, Embase, Web of Science and the Cochrane Library) were searched to select eligible studies. All donor-recipient pairs were high-resolution typing for HLA-A, -B, -C, -DRB1, DQB1 and DPB1 loci. Multivariate-adjusted hazard ratios (HRs) were extracted and pooled using a random-effects model. A total of 36 studies were included, with 100,072 patients receiving HCT. Surprisingly, we found that HLA-DQB1 locus mismatches had no significantly increased risk of multiple outcomes including acute and chronic graft-versus-host disease (GVHD), overall mortality and disease relapse (HR, 1.07; P = .153; HR, 1.07; P = .271; HR, 1.09; P = .230; HR, 1.07; P = .142 and HR, 1.02; P = .806, respectively). Mismatched HLA-DPB1 was significantly associated with a reduced risk of disease relapse (HR, 0.74; P < .001) but not with increased risks of transplant-related mortality (TRM) and overall mortality (HR, 1.09; P = .591; I² = 74.2% and HR, 1.03; P = .460, respectively). In conclusion, HLA-DQB1 locus mismatches is a permissive mismatching. HLA-DPB1 locus mismatches significantly protect against leukemia relapse. Refining effects of individual HLA locus mismatches contributes to predicting prognosis of patients receiving unrelated donor HCT.

Major Histocompatibility Complex and Hematopoietic Stem Cell Transplantation: Beyond the Classical HLA Polymorphism

Allogeneic hematopoietic stem cell transplantation (HSCT) represents a curative treatment for many patients with hematological malignant or non-malignant disorders. Evaluation of potential donors for HSCT includes a rigorous assessment of the human leukocyte antigens (HLA) match status of family members, and the identification of suitable unrelated donors. Genes encoding transplantation antigens are placed both within and outside the major histocompatibility complex (MHC). The human MHC is located on the short arm of chromosome 6 and contains a series of genes encoding two distinct types of highly polymorphic cell surface glycoproteins. Donors for HSCT are routinely selected based on the level of matching for HLA-A, -B, -C, -DRB1, and -DQB1 loci. However, disease relapse, graft-versus-host-disease, and infection remain significant risk factors of morbidity and mortality. In the same breath, in high-risk patients, graft-versus-leukemia effects inherent in HLA mismatching play a substantial immunological role to limit the recurrence of post-transplant disease. The definition of a suitable donor is ever changing, shaped not only by current typing technology, but also by the specific transplant procedure. Indeed, a more complete understanding of permissible HLA mismatches and the role of Killer Immunoglobulin-like receptors' genes increases the availability of HLA-haploidentical and unrelated donors.

Matching donor and recipient based on predicted indirectly recognizable human leucocyte antigen epitopes

The predicted indirectly recognizable human leucocyte antigen (HLA) epitopes (PIRCHE) algorithm is a novel in silico algorithm to determine donor-recipient compatibility. The PIRCHE algorithm determines donor-recipient compatibility by counting the number of mismatched HLA-derived epitopes that are involved in indirect T-cell alloimmune responses; these epitopes are designated as PIRCHE. Over the last few years, the PIRCHE algorithm has been investigated in both hematopoietic stem cell transplantation and solid organ transplantation. This review describes the theory of the algorithm, its application in transplantation, and highlights the future perspectives on the clinical application of the PIRCHE algorithm.

Development of an Unrelated Donor Selection Score Predictive of Survival after HCT: Donor Age Matters Most

Donor factors, in addition to HLA matching status, have been associated with recipient survival in unrelated donor (URD) hematopoietic cell transplantation (HCT); however, there is no hierarchical algorithm that weights the characteristics of individual donors against each other in a quantitative manner to facilitate donor selection. The goal of this study was to develop and validate a donor selection score that prioritizes donor characteristics associated with better survival in 8/8 HLA-matched URDs. Two separate patient/donor cohorts, the first receiving HCT between 1999 and 2011 (n = 5952, c1), and the second between 2012 and 2014 (n = 4510, c2) were included in the analysis. Both cohorts were randomly spilt, 2:1, into training and testing sets. Despite studying over 10,000 URD transplants, we were unable to validate a donor selection score. The only donor characteristic associated with better survival was younger age, with 2-year survival being 3% better when a donor 10 years younger is selected. These results support previous studies suggesting prioritization of a younger 8/8 HLA-matched donor. This large dataset also shows that none of the other donor clinical factors tested were reproducibly associated with survival, and hence flexibility in selecting URDs based on other characteristics is justified. These data support a simplified URD selection process and have significant implications for URD registries.

Elevated HLA-A expression impairs HIV control through inhibition of NKG2A-expressing cells

The highly polymorphic human leukocyte antigen (HLA) locus encodes cell surface proteins that are critical for immunity. HLA-A expression levels vary in an allele-dependent manner, diversifying allele-specific effects beyond peptide-binding preference. Analysis of 9763 HIV-infected individuals from 21 cohorts shows that higher HLA-A levels confer poorer control of HIV. Elevated HLA-A expression provides enhanced levels of an HLA-A–derived signal peptide that specifically binds and determines expression levels of HLA-E, the ligand for the inhibitory NKG2A natural killer (NK) cell receptor. HLA-B haplotypes that favor NKG2A-mediated NK cell licensing (i.e., education) exacerbate the deleterious effect of high HLA-A on HIV control, consistent with NKG2A-mediated inhibition impairing NK cell clearance of HIV-infected targets. Therapeutic blockade of HLA-E:NKG2A interaction may yield benefit in HIV disease.

HLA variation and disease

Fifty years since the first description of an association between HLA and human disease, HLA molecules have proven to be central to physiology, protective immunity and deleterious, disease-causing autoimmune reactivity. Technological advances have enabled pivotal progress in the determination of the molecular mechanisms that underpin the association between HLA genetics and functional outcome. Here, we review our current understanding of HLA molecules as the fundamental platform for immune surveillance and responsiveness in health and disease. We evaluate the scope for personalized antigen-specific disease prevention, whereby harnessing HLA-ligand interactions for clinical benefit is becoming a realistic prospect.

Quantification of classical HLA class I mRNA by allele-specific, real-time polymerase chain reaction for most Han individuals

Recent studies have shown that expression levels of different alleles at the same HLA class I locus can vary dramatically, which might have a broad influence on human disease. However, precise quantification of the relative expression level of each HLA allele is challenging, because distinguishing different alleles on the same locus is difficult. Here, we developed a series of allele-specific, real-time polymerase chain reaction assays for quantifying HLA class I allele mRNA in most Han individuals. The alleles of almost all heterozygous genotypes with a frequency higher than 0.5% in our population (78 alleles on HLA-A locus, 124 alleles on HLA-B locus, and 74 alleles on HLA-C locus) were specifically amplified. The specificity of the amplification was strictly validated by setting the corresponding negative control for each allele of each genotype. The amplification efficiency of each reaction was determined, and the slopes of the reactions were compared. This study provides a tool for detecting the comprehensive expression profile of HLA class I alleles and will be useful not only for the investigation of the molecular mechanism underlying HLA allele expression regulation but also for exploration of immunological mechanisms involving HLA expression in the fields of tumour immune evasion, viral infection, auto-immune disorders, and graft vs host disease after haematopoietic stem cell transplantation.

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.

Which factors influence the development of GVHD in HLA-matched or mismatched transplants?

The sheer diversity of HLA alleles makes the probability of finding matched unrelated donors for patients requiring hematopoietic cell transplantation (HCT) a complex situation. New evidence suggests that mismatching at certain HLA loci may provide a greater benefit in terms of graft-versus-leukemia effect than other mismatches when HLA-matched donors are not available. This review summarizes the current understanding of HLA matching requirements for unrelated donor HCT.

Posttranscriptional Regulation of HLA-A Protein Expression by Alternative Polyadenylation Signals Involving the RNA-Binding Protein Syncrip

Genomic variation in the untranslated region (UTR) has been shown to influence HLA class I expression level and associate with disease outcomes. Sequencing of the 3'UTR of common HLA-A alleles indicated the presence of two polyadenylation signals (PAS). The proximal PAS is conserved, whereas the distal PAS is disrupted within certain alleles by sequence variants. Using 3'RACE, we confirmed expression of two distinct forms of the HLA-A 3'UTR based on use of either the proximal or the distal PAS, which differ in length by 100 bp. Specific HLA-A alleles varied in the usage of the proximal versus distal PAS, with some alleles using only the proximal PAS, and others using both the proximal and distal PAS to differing degrees. We show that the short and the long 3'UTR produced similar mRNA expression levels. However, the long 3'UTR conferred lower luciferase activity as compared with the short form, indicating translation inhibition of the long 3'UTR. RNA affinity pull-down followed by mass spectrometry analysis as well as RNA coimmunoprecipitation indicated differential binding of Syncrip to the long versus short 3'UTR. Depletion of Syncrip by small interfering RNA increased surface expression of an HLA-A allotype that uses primarily the long 3'UTR, whereas an allotype expressing only the short form was unaffected. Furthermore, specific blocking of the proximal 3'UTR reduced surface expression without decreasing mRNA expression. These data demonstrate HLA-A allele-specific variation in PAS usage, which modulates their cell surface expression posttranscriptionally.

Shared peptide binding of HLA Class I and II alleles associate with cutaneous nevirapine hypersensitivity and identify novel risk alleles

Genes of the human leukocyte antigen (HLA) system encode cell-surface proteins involved in regulation of immune responses, and the way drugs interact with the HLA peptide binding groove is important in the immunopathogenesis of T-cell mediated drug hypersensitivity syndromes. Nevirapine (NVP), is an HIV-1 antiretroviral with treatment-limiting hypersensitivity reactions (HSRs) associated with multiple class I and II HLA alleles. Here we utilize a novel analytical approach to explore these multi-allelic associations by systematically examining HLA molecules for similarities in peptide binding specificities and binding pocket structure. We demonstrate that primary predisposition to cutaneous NVP HSR, seen across ancestral groups, can be attributed to a cluster of HLA-C alleles sharing a common binding groove F pocket with HLA-C*04:01. An independent association with a group of class II alleles which share the HLA-DRB1-P4 pocket is also observed. In contrast, NVP HSR protection is afforded by a cluster of HLA-B alleles defined by a characteristic peptide binding groove B pocket. The results suggest drug-specific interactions within the antigen binding cleft can be shared across HLA molecules with similar binding pockets. We thereby provide an explanation for multiple HLA associations with cutaneous NVP HSR and advance insight into its pathogenic mechanisms.

What compatibility in 2017 for the haematopoietic stem cell transplantation?

The diversification of potential donors to perform stem cell allografts now enables to propose a compatible graft cell source adapted to the different clinical situations. Transplants with a geno-identical sibling donor, otherwise with the most HLA-compatible unrelated donor, remain the first-line solutions. Alternative transplants allow to graft patients having no donors in international registries, owing to the rarity of their HLA typing. They are carried out with fairly incompatible grafts and are therefore limited by the existence in the recipient of preformed anti-HLA antibodies which predispose to their rejection. The simple prevention of acute Graft-versus-host disease in haplo-identical transplants, as well as the availability of donors, explain why they have very often replaced placental stem cell transplants. These latter remain useful for pediatric patients or in the absence of family donors.

Structural and regulatory diversity shape HLA-C protein expression levels

Expression of HLA-C varies widely across individuals in an allele-specific manner. This variation in expression can influence efficacy of the immune response, as shown for infectious and autoimmune diseases. MicroRNA binding partially influences differential HLA-C expression, but the additional contributing factors have remained undetermined. Here we use functional and structural analyses to demonstrate that HLA-C expression is modulated not just at the RNA level, but also at the protein level. Specifically, we show that variation in exons 2 and 3, which encode the α1/α2 domains, drives differential expression of HLA-C allomorphs at the cell surface by influencing the structure of the peptide-binding cleft and the diversity of peptides bound by the HLA-C molecules. Together with a phylogenetic analysis, these results highlight the diversity and long-term balancing selection of regulatory factors that modulate HLA-C expression.

Killer Immunoglobulin-Like Receptor Allele Determination Using Next-Generation Sequencing Technology

The impact of natural killer (NK) cell alloreactivity on hematopoietic stem cell transplantation (HSCT) outcome is still debated due to the complexity of graft parameters, HLA class I environment, the nature of killer cell immunoglobulin-like receptor (KIR)/KIR ligand genetic combinations studied, and KIR⁺ NK cell repertoire size. KIR genes are known to be polymorphic in terms of gene content, copy number variation, and number of alleles. These allelic polymorphisms may impact both the phenotype and function of KIR⁺ NK cells. We, therefore, speculate that polymorphisms may alter donor KIR⁺ NK cell phenotype/function thus modulating post-HSCT KIR⁺ NK cell alloreactivity. To investigate KIR allele polymorphisms of all KIR genes, we developed a next-generation sequencing (NGS) technology on a MiSeq platform. To ensure the reliability and specificity of our method, genomic DNA from well-characterized cell lines were used; high-resolution KIR typing results obtained were then compared to those previously reported. Two different bioinformatic pipelines were used allowing the attribution of sequencing reads to specific KIR genes and the assignment of KIR alleles for each KIR gene. Our results demonstrated successful long-range KIR gene amplifications of all reference samples using intergenic KIR primers. The alignment of reads to the human genome reference (hg19) using BiRD pipeline or visualization of data using Profiler software demonstrated that all KIR genes were completely sequenced with a sufficient read depth (mean 317× for all loci) and a high percentage of mapping (mean 93% for all loci). Comparison of high-resolution KIR typing obtained to those published data using exome capture resulted in a reported concordance rate of 95% for centromeric and telomeric KIR genes. Overall, our results suggest that NGS can be used to investigate the broad KIR allelic polymorphism. Hence, these data improve our knowledge, not only on KIR⁺ NK cell alloreactivity in HSCT but also on the role of KIR⁺ NK cell populations in control of viral infections and diseases.

Role of major histocompatibility complex variation in graft-versus-host disease after hematopoietic cell transplantation

Graft-versus-host disease (GVHD) remains a significant potentially life-threatening complication of allogeneic hematopoietic cell transplantation (HCT). Since the discovery of the human leukocyte antigen (HLA) system over 50 years ago, significant advances have clarified the nature of HLA variation between transplant recipients and donors as a chief etiology of GVHD. New information on coding and non-coding gene variation and GVHD risk provides clinicians with options to consider selected mismatched donors when matched donors are not available. These advances have increased the availability of unrelated donors for patients in need of a transplant and have lowered the overall morbidity and mortality of HCT.

Sequence and Phylogenetic Analysis of the Untranslated Promoter Regions for HLA Class I Genes

Polymorphisms located within the MHC have been linked to many disease outcomes by mechanisms not yet fully understood in most cases. Variants located within untranslated regions of HLA genes are involved in allele-specific expression and may therefore underlie some of these disease associations. We determined sequences extending nearly 2 kb upstream of the transcription start site for 68 alleles from 57 major lineages of classical HLA class I genes. The nucleotide diversity within this promoter segment roughly follows that seen within the coding regions, with HLA-B showing the highest (∼1.9%), followed by HLA-A (∼1.8%), and HLA-C showing the lowest diversity (∼0.9%). Despite its greater diversity, HLA-B mRNA expression levels determined in 178 European Americans do not vary in an allele- or lineage-specific manner, unlike the differential expression levels of HLA-A or HLA-C reported previously. Close proximity of promoter sequences in phylogenetic trees is roughly reflected by similarity of expression pattern for most HLA-A and -C loci. Although promoter sequence divergence might impact promoter activity, we observed no clear link between the phylogenetic structures as represented by pairwise nucleotide differences in the promoter regions with estimated differences in mRNA expression levels for the classical class I loci. Further, no pair of class I loci showed coordinated expression levels, suggesting that distinct mechanisms across loci determine their expression level under nonstimulated conditions. These data serve as a foundation for more in-depth analysis of the functional consequences of promoter region variation within the classical HLA class I loci.

Pharmacokinetics, Pharmacodynamics and Pharmacogenomics of Immunosuppressants in Allogeneic Haematopoietic Cell Transplantation: Part I

Although immunosuppressive treatments and target concentration intervention (TCI) have significantly contributed to the success of allogeneic haematopoietic cell transplantation (alloHCT), there is currently no consensus on the best immunosuppressive strategies. Compared with solid organ transplantation, alloHCT is unique because of the potential for bidirectional reactions (i.e. host-versus-graft and graft-versus-host). Postgraft immunosuppression typically includes a calcineurin inhibitor (cyclosporine or tacrolimus) and a short course of methotrexate after high-dose myeloablative conditioning, or a calcineurin inhibitor and mycophenolate mofetil after reduced-intensity conditioning. There are evolving roles for the antithymyocyte globulins (ATGs) and sirolimus as postgraft immunosuppression. A review of the pharmacokinetics and TCI of the main postgraft immunosuppressants is presented in this two-part review. All immunosuppressants are characterized by large intra- and interindividual pharmacokinetic variability and by narrow therapeutic indices. It is essential to understand immunosuppressants' pharmacokinetic properties and how to use them for individualized treatment incorporating TCI to improve outcomes. TCI, which is mandatory for the calcineurin inhibitors and sirolimus, has become an integral part of postgraft immunosuppression. TCI is usually based on trough concentration monitoring, but other approaches include measurement of the area under the concentration-time curve (AUC) over the dosing interval or limited sampling schedules with maximum a posteriori Bayesian personalization approaches. Interpretation of pharmacodynamic results is hindered by the prevalence of studies enrolling only a small number of patients, variability in the allogeneic graft source and variability in postgraft immunosuppression. Given the curative potential of alloHCT, the pharmacodynamics of these immunosuppressants deserves to be explored in depth. Development of sophisticated systems pharmacology models and improved TCI tools are needed to accurately evaluate patients' exposure to drugs in general and to immunosuppressants in particular. Sequential studies, first without and then with TCI, should be conducted to validate the clinical benefit of TCI in homogenous populations; randomized trials are not feasible, because there are higher-priority research questions in alloHCT. In Part I of this article, we review the alloHCT process to facilitate optimal design of pharmacokinetic and pharmacodynamics studies. We also review the pharmacokinetics and TCI of calcineurin inhibitors and methotrexate.

Allorecognition of HLA-C Mismatches by CD8+ T Cells in Hematopoietic Stem Cell Transplantation Is a Complex Interplay between Mismatched Peptide-Binding Region Residues, HLA-C Expression, and HLA-DPB1 Disparities

HLA-C locus mismatches (MMs) are the most frequent class I disparities in unrelated hematopoietic stem cell transplantation (HSCT) and have a detrimental impact on clinical outcome. Recently, a few retrospective clinical studies have reported some variability in the immunogenicity of HLA-C incompatibilities. To get better insight into presumably permissive HLA-C MMs, we have developed a one-way in vitro mixed lymphocyte reaction (MLR) assay allowing to quantify activated CD56^-CD137⁺CD8⁺ lymphocytes in HLA-C incompatible combinations. T cell-mediated alloresponses were correlated with genetic markers such as HLA-C mRNA expression and the number of amino acid (aa) MMs in the α1/α2 domains (peptide-binding region). Because of the high rate of HLA-DPB1 incompatibilities in HLA-A-, B-, C-, DRB1-, and DQB1-matched unrelated HSCT patient/donor pairs, the impact of HLA-DPB1 mismatching, a potential bystander of CD4⁺ T cell activation, was also considered. Heterogeneous alloresponses were measured in 63 HLA-C-mismatched pairs with a positive assay in 52% of the combinations (2.3-18.6% activated CTLs), representing 24 different HLA-A~B~DRB1~DQB1 haplotypes. There was no correlation between measured alloresponses and mRNA expression of the mismatched HLA-C alleles. The HLA-C*03:03/03:04 MM did not induce any positive alloresponse in five MLRs. We also identified HLA-C*02:02 and HLA-C*06:02 as mismatched alleles with lower immunogenicity, and HLA-C*14:02 as a more immunogenic MM. A difference of at least 10 aa residues known to impact peptide/T cell receptor (TCR) binding and a bystander HLA-DPB1 incompatibility had a significant impact on CTL alloreactivity (p = 0.021). The same HLA-C MM, when recognized by two different responders with the same HLA haplotypes, was recognized differently, emphasizing the role of the T-cell repertoire of responding cells. In conclusion, mismatched HLA-C alleles differing by 10 or more aas in the peptide/TCR-binding region, when occurring together with HLA-DPB1 incompatibilities, should be considered as high-risk MMs in unrelated HSCT.

Low Constitutive Cell Surface Expression of HLA-B Is Caused by a Posttranslational Mechanism Involving Glu180 and Arg239

HLA class I cell surface expression is crucial for normal immune responses, and variability in HLA expression may influence the course of infections. We have previously shown that classical HLA class I expression on many human cell types is biased with greatly reduced expression of HLA-B compared with HLA-A in the absence of inflammatory signals. In the search for the mechanisms responsible for this discrepancy, we have recently reported that the regulation is mainly posttranslational and that the C-terminal part of the α2 domain and the α3 domain contain the molecular determinants that explain most of the variability of expression between common HLA-A and -B allomorphs. In this study, we present a fine mapping of the structural determinants that allow such variability by exchanging key amino acids located within the C-terminal part of the α2 domain and the α3 domain of HLA-A2 and -B8, including Glu/Asp at position 177, Gln/Glu at position 180, Gly/Arg at position 239, and Pro/Ser at position 280. We found that the HLA-A2 and -B8 expression profiles could be interconverted to a large extent by mutual exchange of Gln/Glu at position 180 or by Gly/Arg at position 239. The presence of Gln¹⁸⁰ and Gly²³⁹, as in HLA-A2, led to higher cell surface expression levels when compared with the presence of Glu¹⁸⁰ and Arg²³⁹, as in HLA-B8. This indicates that the amino acids at positions 180 and 239 determine the level of cell surface expression of common HLA-A and -B allomorphs, probably by affecting HLA processing in the Ag presentation pathway.

Application of High-Throughput Next-Generation Sequencing for HLA Typing on Buccal Extracted DNA: Results from over 10,000 Donor Recruitment Samples

Background

Unambiguous HLA typing is important in hematopoietic stem cell transplantation (HSCT), HLA disease association studies, and solid organ transplantation. However, current molecular typing methods only interrogate the antigen recognition site (ARS) of HLA genes, resulting in many cis-trans ambiguities that require additional typing methods to resolve. Here we report high-resolution HLA typing of 10,063 National Marrow Donor Program (NMDP) registry donors using long-range PCR by next generation sequencing (NGS) approach on buccal swab DNA.

Methods

Multiplex long-range PCR primers amplified the full-length of HLA class I genes (A, B, C) from promotor to 3’ UTR. Class II genes (DRB1, DQB1) were amplified from exon 2 through part of exon 4. PCR amplicons were pooled and sheared using Covaris fragmentation. Library preparation was performed using the Illumina TruSeq Nano kit on the Beckman FX automated platform. Each sample was tagged with a unique barcode, followed by 2×250 bp paired-end sequencing on the Illumina MiSeq. HLA typing was assigned using Omixon Twin software that combines two independent computational algorithms to ensure high confidence in allele calling. Consensus sequence and typing results were reported in Histoimmunogenetics Markup Language (HML) format. All homozygous alleles were confirmed by Luminex SSO typing and exon novelties were confirmed by Sanger sequencing.

Results

Using this automated workflow, over 10,063 NMDP registry donors were successfully typed under high-resolution by NGS. Despite known challenges of nucleic acid degradation and low DNA concentration commonly associated with buccal-based specimens, 97.8% of samples were successfully amplified using long-range PCR. Among these, 98.2% were successfully reported by NGS, with an accuracy rate of 99.84% in an independent blind Quality Control audit performed by the NDMP. In this study, NGS-HLA typing identified 23 null alleles (0.023%), 92 rare alleles (0.091%) and 42 exon novelties (0.042%).

Conclusion

Long-range, unambiguous HLA genotyping is achievable on clinical buccal swab-extracted DNA. Importantly, full-length gene sequencing and the ability to curate full sequence data will permit future interrogation of the impact of introns, expanded exons, and other gene regulatory sequences on clinical outcomes in transplantation.

The effect of HLA allele and haplotype polymorphisms on donor matching in hematopoietic stem cell transplantation - Croatian experience

The knowledge of HLA characteristics of a patient's population helps to predict the probability of finding a MUD. The study included 170 transplanted patients for whom a search for a MUD in BMDW was performed and a sample of 4000 volunteer unrelated donors from the Croatian Bone Marrow Donor Registry (CBMDR). Patients and their MUDs were typed for HLA-A, -B, -C, -DRB1, and -DQB1 loci using PCR-SSO and PCR-SSP methods while donors were typed for HLA-A, -B, -C, and -DRB1 loci using the PCR-SSO method. A comparison of allele frequencies at tested HLA loci between patients and donors from CBMDR did not reveal significant differences. The majority of patients (117, 68.8%) had a 10/10 MUD, 45 (26.5%) patients had a 9/10 MUD and eight (4.7%) patients had an 8/10 MUD. The highest number of mismatches (MM) was present at HLA-DRB1 (19; 31.1%). The presence of DRB1*11 and DRB1*04 allelic groups among patients caused allelic MMs at HLA-DRB1 in most cases. The presence of an infrequent HLA-B∼C haplotype resulted in the HLA-C MM at antigen level in the majority of cases. The present study clarified HLA factors that cause difficulties in searching for a 10/10 MUD for Croatian patients.

Mismatched unrelated donor transplantation

There are now more than 25 million volunteer donors registered worldwide for patients in need of a life-saving hematopoietic cell transplant to cure blood disorders. Although a human leukocyte antigen (HLA)-matched donor remains the preferred stem cell source for transplantation, the use of a donor with limited HLA mismatching may be considered. Significant advances in clinical and basic research have been instrumental in furthering the understanding of donor-recipient HLA mismatches that are better tolerated than other mismatches. An increased appreciation of the role of regulatory region variation that affects the level of HLA expression provides new approaches for the selection of HLA-mismatched donors.

How to select the best available related or unrelated donor of hematopoietic stem cells?

Recognition of HLA incompatibilities by the immune system represents a major barrier to allogeneic hematopoietic stem cell transplantation. HLA genotypically identical sibling donors are, therefore, the gold standard for transplantation purposes, but only 30% patients have such a donor. For the remaining 70% patients alternative sources of stem cells are a matched unrelated adult volunteer donor, a haploidentical donor or a cord blood unit. The definition of 'HLA matching' depends on the level of resolution and on which loci are tested. The development of HLA molecular typing technologies and the availability of more than 27 million donors in the international database has greatly facilitated unrelated donor searches. The gold standard is high resolution typing at the HLA-A, -B, -C, -DRB1, and -DQB1 loci (10/10 match). Single disparities for HLA-A, -B, - C, or -DRB1 are associated with increased risk of post-transplant complications, but less so in patients with advanced disease, and in those undergoing T-cell-depleted allografting. HLA-DQB1 mismatches seem to be better tolerated and some HLA-C, -DRB1 and -DPB1 disparities are potentially less immunogenic. HLA typing by next-generation sequencing methods is likely to change matching algorithms by providing full sequence information on all HLA loci in a single step. In most European populations a 10/10 matched donor can be found for at least 50% of patients and an additional 20-30% patients may have a 9/10 matched donor. Genetic factors that help in identifying donors with less immunogenic mismatches are discussed. Haploidentical donors are increasingly used as an alternative source of stem cells for those patients lacking a matched unrelated donor.

Identification of a novel HLA-C allele, HLA-C*07:477 in a Drung Chinese individual

HLA-C*07:477 has one nonsynonymous nucleotide change from HLA-C*07:02:01:01 in codon 176, exon 3 (AAG>GAG).

HIV-1 Vpu Mediates HLA-C Downregulation

Many pathogens evade cytotoxic T lymphocytes (CTLs) by downregulating HLA molecules on infected cells, but the loss of HLA can trigger NK cell-mediated lysis. HIV-1 is thought to subvert CTLs while preserving NK cell inhibition by Nef-mediated downregulation of HLA-A and -B but not HLA-C molecules. We find that HLA-C is downregulated by most primary HIV-1 clones, including transmitted founder viruses, in contrast to the laboratory-adapted NL4-3 virus. HLA-C reduction is mediated by viral Vpu and reduces the ability of HLA-C restricted CTLs to suppress viral replication in CD4+ cells in vitro. HLA-A/B are unaffected by Vpu, and primary HIV-1 clones vary in their ability to downregulate HLA-C, possibly in response to whether CTLs or NK cells dominate immune pressure through HLA-C. HIV-2 also suppresses HLA-C expression through distinct mechanisms, underscoring the immune pressure HLA-C exerts on HIV. This viral immune evasion casts new light on the roles of CTLs and NK cells in immune responses against HIV.

Isolation of full-length HLA-C*18:02 allele in an individual from Sub-Saharan Africa

The full-length sequence of HLA-C*18:02 differs from that of HLA-C*18:01 by a single nucleotide polymorphism.

Expression of classical HLA class I molecules: regulation and clinical impacts: Julia Bodmer Award Review 2015

Human leukocyte antigen (HLA) class I genes are ubiquitously expressed, but in a tissue specific-manner. Their expression is primarily regulated at the transcriptional level and can be modulated both positively and negatively by different stimuli. Advances in sequencing technologies led to the identification of new regulatory variants located in the untranslated regions (UTRs), which could influence the expression. After a brief description of the mechanisms underlying the transcriptional regulation of HLA class I genes expression, we will review how the expression levels of HLA class I genes could affect biological and pathological processes. Then, we will discuss on the differential expression of HLA class I genes according to the locus, allele and UTR polymorphisms and its clinical impact. This interesting field of study led to a new dimension of HLA typing, going beyond a qualitative aspect.