Impact of the HLA Immunopeptidome on Survival of Leukemia Patients After Unrelated Donor Transplantation

PURPOSE

Immunopeptidome divergence between mismatched HLA-DP is a determinant of T-cell alloreactivity and clinical tolerability after fully HLA-A, -B, -C, -DRB1, -DQB1 matched unrelated donor hematopoietic cell transplantation (UD-HCT). Here, we tested this concept in HLA-A, -B, and -C disparities after single class I HLA-mismatched UD-HCT.

PATIENTS AND METHODS

We studied 2,391 single class I HLA-mismatched and 14,426 fully HLA-matched UD-HCT performed between 2008 and 2018 for acute leukemia or myelodysplastic syndromes. Hierarchical clustering of experimentally determined peptide-binding motifs (PBM) was used as a proxy for immunopeptidome divergence of HLA-A, -B, or -C disparities, allowing us to classify 1,629/2,391 (68.1%) of the HLA-mismatched UD-HCT as PBM-matched or PBM-mismatched. Risks associated with PBM-matching status were assessed by Cox proportional hazards models, with overall survival (OS) as the primary end point.

RESULTS

Relative to full matches, bidirectional or unidirectional PBM mismatches in graft-versus-host (GVH) direction (PBM-GVH mismatches, 60.7%) were associated with significantly lower OS (hazard ratio [HR], 1.48; P < .0001), while unidirectional PBM mismatches in host-versus-graft direction or PBM matches (PBM-GVH matches, 39.3%) were not (HR, 1.13; P = .1017). PBM-GVH mismatches also had significantly lower OS than PBM-GVH matches in direct comparison (HR, 1.32; P = .0036). The hazards for transplant-related mortality and acute and chronic graft-versus-host disease but not relapse increased stepwise from full HLA matches to single PBM-GVH matches, and single PBM-GVH mismatches. A webtool for PBM-matching of single class I HLA-mismatched donor-recipient pairs was developed.

CONCLUSION

PBM-GVH mismatches inform mortality risks after single class I HLA-mismatched UD-HCT, suggesting that prospective consideration of directional PBM-matching status might improve outcome. These findings highlight immunopeptidome divergence between mismatched HLA as a driver of clinical tolerability in UD-HCT.

Time series clustering of T cell subsets dissects heterogeneity in immune reconstitution and clinical outcomes among MUD-HCT patients receiving ATG or PTCy

Introduction

Anti-T-lymphocyte globulin (ATG) or post-transplant cyclophosphamide (PTCy) prevent graft-versus-host disease (GVHD) after hematopoietic cell transplantation (HCT), yet individual patients benefit differentially.

Methods

Given the sparse comparative data on the impact of cellular immune reconstitution in this setting, we studied flow cytometry and clinical outcomes in 339 recipients of 10/10 matched-unrelated donor (MUD) HCT using either ATG (n=304) or PTCy (n=35) for in vivo T cell manipulation along with a haploidentical PTCy control cohort (n=45). Longitudinal cellular immune reconstitution data were analyzed conventionally and with a data science approach using clustering with dynamic time warping to determine the similarity between time-series of T cell subsets.

Results

Consistent with published studies, no significant differences in clinical outcomes were observed at the cohort level between MUD-ATG and MUD-PTCy. However, cellular reconstitution revealed preferences for distinct T cell subpopulations associating with GVHD protection in each setting. Starting early after HCT, MUD-PTCy patients had higher regulatory T cell levels after HCT (p &lt;0.0001), while MUD-ATG patients presented with higher levels of γδ T- or NKT cells (both p &lt;0.0001). Time-series clustering further dissected the patient population’s heterogeneity revealing distinct immune reconstitution clusters. Importantly, it identified phenotypes that reproducibly associated with impaired clinical outcomes within the same in vivo T cell manipulation platform. Exemplarily, patients with lower activated- and αβ T cell counts had significantly higher NRM (p=0.032) and relapse rates (p =0.01).

Discussion

The improved understanding of the heterogeneity of cellular reconstitution in MUD patients with T cell manipulation both at the cohort and individual level may support clinicians in managing HCT complications.

An HLA map of the world: A comparison of HLA frequencies in 200 worldwide populations reveals diverse patterns for class I and class II

HLA frequencies show widespread variation across human populations. Demographic factors as well as selection are thought to have shaped HLA variation across continents. In this study, a worldwide comparison of HLA class I and class II diversity was carried out. Multidimensional scaling techniques were applied to 50 HLA-A and HLA-B (class I) as well as 13 HLA-DRB1 (class II) first-field frequencies in 200 populations from all continents. Our results confirm a strong effect of geography on the distribution of HLA class I allele groups, with principal coordinates analysis closely resembling geographical location of populations, especially those of Africa-Eurasia. Conversely, class II frequencies stratify populations along a continuum of differentiation less clearly correlated to actual geographic location. Double clustering analysis revealed finer intra-continental sub-clusters (e.g., Northern and Western Europe vs. South East Europe, North Africa and Southwest Asia; South and East Africa vs. West Africa), and HLA allele group patterns characteristic of these clusters. Ancient (Austronesian expansion) and more recent (Romani people in Europe) migrations, as well as extreme differentiation (Taiwan indigenous peoples, Native Americans), and interregional gene flow (Sámi, Egyptians) are also reflected by the results. Barrier analysis comparing D_ST and geographic location identified genetic discontinuities caused by natural barriers or human behavior explaining inter and intra-continental HLA borders for class I and class II. Overall, a progressive reduction in HLA diversity from African to Oceanian and Native American populations is noted. This analysis of HLA frequencies in a unique set of worldwide populations confirms previous findings on the remarkable similarity of class I frequencies to geography, but also shows a more complex development for class II, with implications for both human evolutionary studies and biomedical research.

The immunogenetic diversity of the HLA system in Mexico correlates with underlying population genetic structure

We studied HLA class I (HLA-A, -B) and class II (HLA-DRB1, -DQB1) allele groups and alleles by PCR-SSP based typing in a total of 15,318 mixed ancestry Mexicans from all the states of the country divided into 78 sample sets, providing information regarding allelic and haplotypic frequencies and their linkage disequilibrium, as well as admixture estimates and genetic substructure. We identified the presence of 4268 unique HLA extended haplotypes across Mexico and find that the ten most frequent (HF > 1%) HLA haplotypes with significant linkage disequilibrium (Δ'≥0.1) in Mexico (accounting for 20% of the haplotypic diversity of the country) are of primarily Native American ancestry (A*02~B*39~DRB1*04~DQB1*03:02, A*02~B*35~DRB1*08~DQB1*04, A*68~B*39~DRB1*04~DQB1*03:02, A*02~B*35~DRB1*04~DQB1*03:02, A*24~B*39~DRB1*14~DQB1*03:01, A*24~B*35~DRB1*04~DQB1*03:02, A*24~B*39~DRB1*04~DQB1*03:02, A*02~B*40:02~DRB1*04~DQB1*03:02, A*68~B*35~DRB1*04~DQB1*03:02, A*02~B*15:01~DRB1*04~DQB1*03:02). Admixture estimates obtained by a maximum likelihood method using HLA-A/-B/-DRB1 as genetic estimators revealed that the main genetic components in Mexico as a whole are Native American (ranging from 37.8% in the northern part of the country to 81.5% in the southeastern region) and European (ranging from 11.5% in the southeast to 62.6% in northern Mexico). African admixture ranged from 0.0 to 12.7% not following any specific pattern. We were able to detect three major immunogenetic clusters correlating with genetic diversity and differential admixture within Mexico: North, Central and Southeast, which is in accordance with previous reports using genome-wide data. Our findings provide insights into the population immunogenetic substructure of the whole country and add to the knowledge of mixed ancestry Latin American population genetics, important for disease association studies, detection of demographic signatures on population variation and improved allocation of public health resources.

High-resolution HLA allele and haplotype frequencies in majority and minority populations of Costa Rica and Nicaragua: Differential admixture proportions in neighboring countries

The HLA system shows the most extensive polymorphism in the human genome. Allelic and haplotypic frequencies of HLA genes vary dramatically across human populations. Due to a complex history of migration, populations in Latin America show a broad variety of admixture proportions, usually varying not only between countries, but also within countries. Knowledge of HLA allele and haplotype frequencies is essential for medical fields such as transplantation, but also serves as a means to assess genetic diversity and ancestry in human populations. Here, we have determined high-resolution HLA-A, -B, -C, and -DRB1 allele and haplotype frequencies in a sample of 713 healthy subjects from three Mestizo populations, one population of African descent, and Amerindians of five different groups from Costa Rica and Nicaragua and compared their profiles to a large set of indigenous populations from Iberia, Sub-Saharan Africa, and the Americas. Our results show a great degree of allelic and haplotypic diversity within and across these populations, with most extended haplotypes being private. Mestizo populations show alleles and haplotypes of putative European, Amerindian, and Sub-Saharan African origin, albeit with differential proportions. Despite some degree of gene flow, Amerindians and Afro-descendants show great similarity to other Amerindian and West African populations, respectively. This is the first comprehensive study reporting high-resolution HLA diversity in Central America, and its results will shed light into the genetic history of this region while also supporting the development of medical programs for organ and stem cell transplantation.

4-Locus high-resolution HLA allele and haplotype frequencies in Amerindians from Costa Rica

A total of 125 Costa Ricans of Amerindian descent were genotyped at high-resolution for the human leukocyte antigen loci HLA-A, -B, -C, and -DRB1 using sequence-based typing methods. The respective allele and extended haplotype frequencies, as well as Hardy-Weinberg proportions were calculated. The most frequent extended haplotype identified was A*24:02:01-B*40:02:01-C*03:05-DRB1*04:07:01G, with an estimated frequency of 8.26%. A deviation from Hardy-Weinberg Equilibrium was detected at the DRB1 locus (p = 0.099). The HLA genotypic data of the population sample reported here are available publicly in the Allele Frequencies Net Database under the population name "Costa Rica Amerindians" and the identifier AFN3608.

5-Locus high-resolution HLA allele and haplotype frequencies in Costa Ricans from the Central Valley

A total of 221 Costa Rican Mestizos from the Central Valley were genotyped at high-resolution for the human leukocyte antigen loci HLA-A, -B, -C, -DRB1, and -DQB1 using sequence-based typing methods. The respective allele and extended haplotype frequencies, as well as Hardy-Weinberg proportions were calculated. The most frequent extended haplotype identified was A*24:02:01-B*40:02:01-C*03:05-DRB1*08:02:01-DQB1*04:02:01, with an estimated frequency of 2.04%. No deviation from Hardy-Weinberg Equilibrium was detected at any of the loci studied. The HLA genotypic data of the population sample reported here are available publicly in the Allele Frequencies Net Database under the population name "Costa Rica Central Valley Mestizo" and the identifier AFN3606.

4-Locus high-resolution HLA allele and haplotype frequencies in Costa Ricans from Guanacaste

A total of 110 Costa Rican Mestizos from the province of Guanacaste were genotyped at high-resolution for the human leukocyte antigen loci HLA-A, -B, -C, and -DRB1 using sequence-based typing methods. The respective allele and extended haplotype frequencies, as well as Hardy-Weinberg proportions were calculated. The most frequent extended haplotype identified was A*24:02:01-B*35:12:01-C*04:01:01-DRB1*04:07:01G, with an estimated frequency of 2.73%. No deviation from Hardy-Weinberg Equilibrium was detected at any of the loci studied. The HLA genotypic data of the population sample reported here are available publicly in the Allele Frequencies Net Database under the population name "Costa Rica Guanacaste Mestizo" and the identifier AFN3609.

4-Locus high-resolution HLA allele and haplotype frequencies in admixed population from Nicaragua

A total of 155 Nicaraguan Mestizos from across the country were genotyped at high-resolution for the human leukocyte antigen loci HLA-A, -B, -C, and -DRB1 using sequence-based typing methods. The respective allele and extended haplotype frequencies, as well as Hardy-Weinberg proportions were calculated. The most frequent extended haplotype identified was A*24:02:01-B*40:02:01-C*03:05-DRB1*04:07:01G, with an estimated frequency of 2.26%. No deviation from Hardy-Weinberg Equilibrium was detected at any of the loci studied. The HLA genotypic data of the population sample reported here are available publicly in the Allele Frequencies Net Database under the population name "Nicaragua Mestizo" and the identifier AFN3610.

Multiple Knockout of Classical HLA Class II β-Chains by CRISPR/Cas9 Genome Editing Driven by a Single Guide RNA

Comprehensive knockout of HLA class II (HLA-II) β-chain genes is complicated by their high polymorphism. In this study, we developed CRISPR/Cas9 genome editing to simultaneously target HLA-DRB, -DQB1, and -DPB1 through a single guide RNA recognizing a conserved region in exon 2. Abrogation of HLA-II surface expression was achieved in five different HLA-typed, human EBV-transformed B lymphoblastoid cell lines (BLCLs). Next-generation sequencing-based detection confirmed specific genomic insertion/deletion mutations with 99.5% penetrance in sorted cells for all three loci. No alterations were observed in HLA-I genes, the HLA-II peptide editor HLA-DMB, or its antagonist HLA-DOB, showing high on-target specificity. Transfection of full-length HLA-DPB1 mRNA into knockout BLCLs fully restored HLA-DP surface expression and recognition by alloreactive human CD4 T cells. The possibility to generate single HLA-II-expressing BLCLs by one-shot genome editing opens unprecedented opportunities for mechanistically dissecting the interaction of individual HLA variants with the immune system.

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.

Alloreactive T Cell Receptor Diversity against Structurally Similar or Dissimilar HLA-DP Antigens Assessed by Deep Sequencing

T cell alloreactivity is mediated by a self-human leukocyte antigen (HLA)-restricted T cell receptor (TCR) repertoire able to recognize both structurally similar and dissimilar allogeneic HLA molecules (i.e., differing by a single or several amino acids in their peptide-binding groove). We hypothesized that thymic selection on self-HLA molecules could have an indirect impact on the size and diversity of the alloreactive response. To test this possibility, we used TCR Vβ immunophenotyping and immunosequencing technology in a model of alloreactivity between self-HLA selected T cells and allogeneic HLA-DPB1 (DPB1) differing from self-DPB1*04:02 by a single (DPB1*02:01) or several (DPB1*09:01) amino acids in the peptide-binding groove. CD4+ T cells from three different self-DPB1*04:01,*04:02 individuals were stimulated with HeLa cells stably transduced with the relevant peptide processing machinery, co-stimulatory molecules, and HLA-DP. Flow cytometric quantification of the DPB1-specific T cell response measured as upregulation of the activation marker CD137 revealed significantly lower levels of alloreactivity against DPB1*02:01 compared with DPB1*09:01 (mean CD4+CD137+ frequency 35.2 ± 9.9 vs. 61.5 ± 7.7%, respectively, p < 0.0001). These quantitative differences were, however, not reflected by differences in the breadth of the alloreactive response at the Vβ level, with both alloantigens eliciting specific responses from all TCR-Vβ specificities tested by flow cytometry, albeit with higher levels of reactivity from most Vβ specificities against DPB1*09:01. In line with these observations, TCRB-CDR3 immunosequencing showed no significant differences in mean clonality of sorted CD137+CD4+ cells alloreactive against DPB1*02:01 or DPB1*09:01 [0.39 (0.36–0.45) and 0.39 (0.30–0.46), respectively], or in the cumulative frequencies of the 10 most frequent responding clones (55–67 and 58–62%, respectively). Most of the clones alloreactive against DPB1*02:01 (68.3%) or DPB1*09:01 (75.3%) were characterized by low-abundance (i.e., they were not appreciable among the pre-culture T cells). Interestingly, however, their cumulative frequency was lower against DPB1*02:01 compared with DPB1*09:01 (mean cumulative frequency 35.3 vs. 50.6%, respectively). Our data show that, despite lower levels of alloreactivity, a similar clonal diversity can be elicited by structurally similar compared with structurally dissimilar HLA-DPB1 alloantigens and demonstrate the power of TCRB immunosequencing in unraveling subtle qualitative changes not appreciable by conventional methods.

Prevalence and geographic distribution of haemophilia in Costa Rica

Haemophilia is the most frequent hereditary haemorrhagic illness and it is due to the deficiency of coagulation factors VIII (haemophilia A, HA) or IX (haemophilia B, HB).The prevalence of this disease varies according to the country, those having better survival rates having also higher prevalences. Specifically in Costa Rica, there are around 130 HA and 30 HB families. This study reports the prevalence and a spatial distribution analysis of both types of the disease in this country. The prevalence of haemophilia in this country is 7 cases per 100 000 men, for HA it is 6 cases per 100 000 and for HB it is 1 case per 100 000 male inhabitants. The prevalence of this disease is low when compared with other populations. This low prevalence could be due to the many patients that have died because of infection with human immunodeficiency virus during the 1980s. The prevalence of haemophilia in Costa Rica is almost one half of that present in developed countries. Nevertheless, the ratio between HA and HB follows world tendency: 5 : 1. In this study, nationwide geographical distribution maps were drawn in order to visualize the origin of severe cases and how this influences the pattern of distribution for both types of haemophilia. By means of these maps, it was possible to state that there is no association between the sites of maximum prevalence of mutated alleles and ethnicity. With this study, haemophilia prevalence distribution maps can be used to improve efforts for the establishment of hemophilia clinics or specialized health centers in those areas which hold the highest prevalences in this country. Also, this knowledge can be applied to improve treatment skills and offer the possibility of developing focused genetic counseling for these populations.

Polymorphism in TGFB1 is associated with worse non-relapse mortality and overall survival after stem cell transplantation with unrelated donors

Transforming growth factor β-1, encoded by the TGFB1 gene, is a cytokine that plays a central role in many physiological and pathogenic processes. We have sequenced TGFB1 regulatory region and assigned allelic genotypes in a large cohort of hematopoietic stem cell transplantation patients and donors. In this study, we analyzed 522 unrelated donor-patient pairs and examined the combined effect of all the common polymorphisms in this genomic region. In univariate analysis, we found that patients carrying a specific allele, 'p001', showed significantly reduced overall survival (5-year overall survival 30.7% for p001/p001 patients vs. 41.6% others; P=0.032) and increased non-relapse mortality (1-year non-relapse mortality: 39.0% vs. 25.4%; P=0.039) after transplantation. In multivariate analysis, the presence of a p001/p001 genotype in patients was confirmed as an independent factor for reduced overall survival [hazard ratio=1.53 (1.04-2.24); P=0.031], and increased non-relapse mortality [hazard ratio=1.73 (1.06-2.83); P=0.030]. In functional experiments we found a trend towards a higher percentage of surface transforming growth factor β-1-positive regulatory T cells after activation when the cells had a p001 allele (P=0.07). Higher or lower production of transforming growth factor β-1 in the inflammatory context of hematopoietic stem cell transplantation may influence the development of complications in these patients. Findings indicate that TGFB1 genotype could potentially be of use as a prognostic factor in hematopoietic stem cell transplantation risk assessment algorithms.

Novel alleles of the transforming growth factor β-1 regulatory region and exon 1

Transforming growth factor β-1, encoded by the TGFB1 gene, is a cytokine that plays a central role in many physiologic and pathogenic processes with pleiotropic effects. Regulatory activity for this gene has been shown for 3.0 kb between positions -2665 and +423 from its translational start site. At least 17 TGFB1 regulatory region and exon 1 alleles have been defined on the basis of 18 polymorphic sites. Polymorphisms in TGFB1's regulatory region have been associated with differential levels of expression of this cytokine and to genetic risk in cancer and transplantation. In this report, we present 19 novel TGFB1 regulatory region and exon 1 alleles: p018-p036. Amplification of TGFB1's regulatory region was performed with an in-house protocol, and novel alleles were defined by either allele-specific amplification and/or molecular cloning of the amplicons, followed by sequencing in isolation. Three of these novel alleles (p018, p019, and p020) are shown to be formed by novel combinations of the aforementioned known polymorphic positions. Another 16 novel alleles are shown to carry additional known and unknown single-nucleotide polymorphisms. Polymorphism in TGFB1's regulatory region could have an impact on important processes, including embryogenesis, hematopoiesis, carcinogenesis, angiogenesis, fibrosis, immune responses, and transplantation, making its characterization necessary.

Prevalence and geographic distribution of haemophilia in Costa Rica

Haemophilia is the most frequent hereditary haemorrhagic illness and it is due to the deficiency of coagulation factors VIII (haemophilia A, HA) or IX (haemophilia B, HB). The prevalence of this disease varies according to the country, those having better survival rates having also higher prevalences. Specifically in Costa Rica, there are around 130 HA and 30 HB families. This study reports the prevalence and a spatial distribution analysis of both types of the disease in this country. The prevalence of haemophilia in this country is 7 cases per 100000 men, for HA it is 6 cases per 100000 and for HB it is 1 case per 100000 male inhabitants. The prevalence of this disease is low when compared with other populations. This low prevalence could be due to the many patients that have died because of infection with human immunodeficiency virus during the 1980s. The prevalence of haemophilia in Costa Rica is almost one half of that present in developed countries. Nevertheless, the ratio between HA and HB follows world tendency: 5:1. In this study, nationwide geographical distribution maps were drawn in order to visualize the origin of severe cases and how this influences the pattern of distribution for both types of haemophilia. By means of these maps, it was possible to state that there is no association between the sites of maximum prevalence of mutated alleles and ethnicity. With this study, haemophilia prevalence distribution maps can be used to improve efforts for the establishment of hemophilia clinics or specialized health centers in those areas which hold the highest prevalences in this country. Also, this knowledge can be applied to improve treatment skills and offer the possibility of developing focused genetic counseling for these populations.