• BLEAT facilitates the selection of stem cell sources for single HLA-B mismatches.

  • BLEAT assessment shows 99.99% of more than 2 million HLA-B leader-type alleles have the same leader type as their HLA-B allele family.

Sequence variation in the HLA-B gene is critically linked to differential immune responses. A dimorphism at -21 of HLA-B exon 1 gives rise to leader peptides that are markers for risk of acute graft-versus-host disease, relapse, and mortality after unrelated donor and cord blood transplantation. To optimize the selection of stem cell transplant sources based on the HLA-B leader, an HLA-BLeader Assessment Tool (BLEAT) was developed to automate the assignment of leader genotypes, define HLA-B leader match statuses, and rank order candidate stem cell sources according to clinical risk. The base cohort consisted of 9 417 614 registered donors from the Be The Match Registry with HLA-B typing. Among these donors, the performance of BLEAT was assessed in 1 098 358 donors with sequence data for HLA-B exon 1 (2 196 716 haplotypes). The accuracy of leader assignment was then assessed in a second cohort of 1259 patients and their unrelated transplant donors. We furthermore established the frequencies of HLA-B leader genotype (MM, MT, TT) representations in broad racial categories in the 9.42 million donors. BLEAT has direct applications for the selection of optimal stem cell sources for transplantation and broad utility in basic and clinical research in pharmacogenomics, vaccine development, and cancer and infectious disease studies of human populations.

HLA-B is the most polymorphic locus in the human genome with >7561 unique alleles recognized currently.1  The HLA-B protein plays an essential role in T-cell and natural killer (NK) cell immunity.2  Antigen presentation of HLA-B to T cells is the basis for host defense against infectious agents and graft acceptance in transplantation.3  Its multidimensional role in both adaptive and innate immunity is the foundation for HLA-B associations with many human diseases and responses to pharmacologic agents.4-9 

HLA class I genes comprise 7 exons. Exon 1 encodes a nonamer leader peptide cleaved from the N-terminus of the HLA-B molecule during protein maturation. The resultant peptide can be presented by HLA-E molecules to stabilize their expression on the cell surface and enhance binding to the cognate receptor CD94/NKG2-A on NK cells. Polymorphisms in the sequence of the HLA-E-bound peptide are predicted to alter HLA-E stability and CD94/NKG2-A binding.10-12  Exons 2, 3, and 4 give rise to α1, α2, and α3 domains, respectively, and exons 5 through 7 encode the transmembrane and cytoplasmic tail. Variation within α1 and α2 domains provides structural and functional diversity to the peptide-binding region conventionally assessed for the matching of candidate stem cell donors and cord blood units for transplantation.

Beyond the peptide-binding region, functional sequence variation exists in exon 1.1  The -21 position of exon 1 (rs1050458) is invariant in HLA-A and HLA-C but dimorphic in HLA-B, where a cytosine/thymine change gives rise to leader peptides with threonine (T) or methionine (M) respectively at the second position (P2) of the nonamer peptide.1  Two ancient clades of HLA-B alleles defined by full-length gene sequences demonstrate a strong linkage between the M/T-leader dimorphism and coding variants of a given allele family with very few exceptions.13 

Recent data implicate a sentinel role for the HLA-B leader in transplant survivorship, NK responses to leukemia, and HIV control.13-18  In unrelated donor and cord blood transplantation, the HLA-B leader identifies donor and cord blood unit HLA mismatch scenarios that are better tolerated than others and inform transplant survival.13,16,17  The mechanisms through which the dimorphic HLA-B leader affects immune responses in transplantation are not known but may include differential expression of the non-classical class I molecule HLA-E, a ligand for T- and NK-cell receptors.10,11,19 

The use of the HLA-B leader to guide the judicious selection of stem cell sources has a particular impact in underserved populations in which HLA-matched stem cell sources are most often unavailable.20  However, a major roadblock for clinical translation of the leader to donor selection is the lack of complete exon 1 sequence data to relate leader assignment to the HLA-B allele. Currently, exon 1 sequences are publicly cataloged for 4148 (55%) of the currently recognized 7561 HLA-B alleles.1  Beyond incomplete HLA-B exon 1 sequence information, clinical laboratories lack tools for translation of HLA-B tissue typing to the leader defined by the P2 dimorphism to inform clinical decision making. This deficiency is particularly acute when HLA typing approaches do not include a direct assessment of exon 1.

To address these needs, we developed HLA-B Leader Assessment Tool (BLEAT), a program to automate leader assignment to each HLA-B allele, define leader genotypes (MM, MT, TT) of individuals, and indicate preferential order for selection of stem cell sources for transplantation based on prior outcome studies. Provision of the patient’s leader genotype provides information on transplant risks conferred by the patient’s germline HLA-B. Pairwise assessment of a patient with each candidate donor or cord blood unit provides information on whether mismatched HLA-B alleles have the same (“leader-matched”) or different (“leader-mismatched”) leaders and whether the matched allele has an M or a T leader. We leveraged phased HLA-B gene sequences in 1 098 358 volunteer US registry donors and a large transplant cohort to develop and validate the leader tool and define leader genotypes in racially diverse US populations. Beyond enhancing the evaluation of unrelated donors and cord blood units for transplantation, BLEAT provides a foundation for anthropologic, pharmacogenomic, and disease association studies.21-24 

Base cohort

As determined by the National Marrow Donor Program (NMDP) institutional review board staff, this study is not human subjects research as defined by regulations. The base cohort includes 9.42 million US-based donors with HLA-B typing from the Be The Match Registry as of May 29, 2020, who self-identified using NMDP race categorizations through either miscellaneous categories (ie, unknown, other, declined, multirace) or broad race categories (ie, African American [AFA], Asian or Pacific Islander [API], White [CAU], Hispanic or Latin American [HIS], American Indian or Alaskan Native [NAM]).25  NMDP broad race categories were mapped to US Census categories as Black or AFA, API, CAU, HIS, and NAM, respectively.26  Specific HLA-B alleles or availability of the leader sequence in International ImMunogGeneTics (IMGT) were not used to define the base cohort.

Cohort 1

Within the base cohort of 9.42 million donors, 1 098 358 donors (cohort 1) were sequenced for HLA-B exons 1 through 7 using single-molecule real-time PacBio sequencing.27  HLA sequencing laboratories prioritize coding sequences but not necessarily noncoding regions. Resulting ambiguity may be represented using any number of typing fields from 1 to 4, which are delimited via colons in an allele name (mock allele: HLA-B*01:02:03:04). Each donor in cohort 1 contributes 2 alleles, resulting in 2 196 716 total alleles. There were 2259 unique HLA-B alleles including 44 G groups (eg, HLA-B*07:02:01G) and 1 multiple allele code (MAC) allele (ie, HLA-B*39:BMFM). G groups represent ambiguous strings of alleles with the same nucleotide sequences across peptide-binding domains (exons 2 and 3 for HLA-B). MAC alleles are ambiguous strings encoded and decoded using a MAC service (HLA-B*39:BMFM encodes 46 alleles).28 

Long-read sequences were reported to NMDP from March 2016 to July 2020. The 1 098 358 donors in cohort 1 served as the foundation for analysis of 2 196 716 HLA-B alleles (2259 total unique HLA-B alleles observed in cohort 1). Frequencies of the P2 M- and T-leader alleles and TT, MT, and MM leader genotype for each subject were identified for each broad race group.

Cohort 2

Cohort 2 consisted of 1269 patients with HLA typing in the Center for International Blood and Marrow Transplant Research database who received an unrelated donor transplant. An input supplemental File 1 was imported within BLEAT, leader-assigned, and exported. The resulting output dataset is supplemental File 2. This dataset served as a blinded validation of leader genotype assignment, shared leader P2, and match status of the leader alleles in patients and donors. High-resolution HLA typing was captured for patients and donors as previously described.13  Informed consent was obtained from transplant patients and donors. Protocols were approved by the institutional review boards of the National Institutes of Health Office for Human Research Protections, Be The Match/NMDP, and the Fred Hutchinson Cancer Research Center.

HLA-B exon 1 sequencing analysis

The 2 196 716 exon 1 sequences from 932 692 donors in cohort 1 were analyzed using a Jupyter notebook. The analysis leverages ImmunoPolymorphism Database (IPD)-IMGT/HLA database (version 3.42.0) exon 1 reference sequences1  to verify the exon 1 sequences of alleles with various resolutions: G groups, MAC, 2-to-4-field typing. Analysis scripts are available at https://github.com/nmdp-bioinformatics/b-leader/tree/master/analysis.

Statistical analysis

The genotype frequencies of HLA-B leader genotypes presented in Table 3 were tested for significant differences across 5 distinct broad race categories using Pearson’s χ2 test of independence with 8 degrees of freedom using R (version 4.0.2).29  Residual ambiguity at the level of leader allele polymorphism was resolved by Monte Carlo sampling of the likelihood distribution from statistical imputation based on haplotype frequencies derived from the corresponding US population.

HLA-B leader sequence variation

The rs1050458C/T exon 1 dimorphism of the HLA-B leader is present in all HLA-B alleles recognized to date with very rare exceptions.1  Of the 9.42 million HLA-typed donors (base cohort) evaluated in the current study, a subset of 1 098 358 (cohort 1) were sequenced across exons 1 through 7 and introns 1 through 6, inclusively, with a median of 2724 base pairs and an interquartile range of 13 (2712, 2725) base pairs. HLA-B alleles were assigned based on long-read, phased sequencing of exons 1 through 7,27  with most sequences initiating 33 base pairs downstream of exon 1 in the 5' untranslated region and terminating at the end of exon 7.

Sequence variation at 12 of the 27 base pairs that encode the nonamer peptide (chromosome 6: 31 357,125-31 357,152 [GRCh38/hg38])30  was identified in cohort 1. These variant positions gave rise to 35 unique nucleotide sequences expected to encode 26 unique nonamer leader peptides (Table 1). The 5 most frequently observed sequences were found in 99.99% of the 2 196 716 haplotypes and encoded either T (68.9%) or M (31.1%) at P2. Rare variants were found, including 3 observations of a leader with arginine (R) and 1 observation with valine (V) at P2. Of the 3413 HLA-B alleles that currently lack exon 1 data in IMGT, 510 were present in cohort 1. Most of the 510 HLA-B alleles (99.99%) encoded a leader allele whose P2 aligned with the first-field HLA-B allele family, as demonstrated in Table 2.

In summary, phased sequencing of exon 1 of HLA-B confirmed the general pattern of 1-to-1 linkage for the HLA-B coding region with the M and T dimorphism at P2 of the leader allele. We leveraged the complete sequence data on leader-allele haplotypes in cohort 1 and designed BLEAT to be informative for all currently known exon 1 through 4 HLA-B haplotypes.

Overview of BLEAT

More than 99.9% of HLA-B allele leader haplotypes encode either an M or T leader. Three steps are performed by BLEAT: (1) mapping of subject HLA-B allele typing to M- and T-leader designations (Figure 1); (2) assignment of leader genotype (MM, MT, TT) and match status for transplant pairs with 1 HLA-B mismatch (Figure 2); and (3) rank ordering of candidate donors based on previously defined risks associated with leader matching (Figure 2).

Figure 1.

BLEAT leader assignment. Two scenarios are depicted using the 2 HLA-B alleles from 1 subject to illustrate a simple use case (HLA-B*07:06:01) and a more complex use case (HLA-B*56:01:01G). For clarity, only the last 5 alleles (by number) are displayed for HLA-B*56:01:01G.36  (A) The leader peptide is deduced from the exon 1 sequence. (B) The deduced leader peptide is available from IMGT in the Anthony Nolan HLA Informatics Group GitHub repository.35  (C) This information is processed and reflected on BLEAT’s user interface to display P2 information. The user may hover over the P2 “Leader” icon to reveal a tooltip (D) that organizes all potential HLA-B alleles from the provided allele into known and unknown major (and any minor, if applicable) leader alleles. Rare minor leader alleles are highlighted with red on the user interface.

Figure 1.

BLEAT leader assignment. Two scenarios are depicted using the 2 HLA-B alleles from 1 subject to illustrate a simple use case (HLA-B*07:06:01) and a more complex use case (HLA-B*56:01:01G). For clarity, only the last 5 alleles (by number) are displayed for HLA-B*56:01:01G.36  (A) The leader peptide is deduced from the exon 1 sequence. (B) The deduced leader peptide is available from IMGT in the Anthony Nolan HLA Informatics Group GitHub repository.35  (C) This information is processed and reflected on BLEAT’s user interface to display P2 information. The user may hover over the P2 “Leader” icon to reveal a tooltip (D) that organizes all potential HLA-B alleles from the provided allele into known and unknown major (and any minor, if applicable) leader alleles. Rare minor leader alleles are highlighted with red on the user interface.

Close modal
Figure 2.

BLEAT provides a rank order for candidate stem cell sources (“donors”) to optimize therapy selection and patient outcomes. The BLEAT user interface (left) has an infographic (right) that is toggled via the Help button (A). On the interface, users enter HLA-B typing information for each subject. (A) The patient’s leader genotype (MM, MT, or TT) is automatically classified upon entering the patient’s HLA-B typing. Ambiguous alleles can sometimes contain rare leader variants; the tool highlights this possibility in red. (B) HLA-B allele typing for donors can be manually entered or automatically imported. Before or during sorting, (C) the leader match status for each donor (matched with the patient) is calculated. Once calculated, a donor selection guide (D) highlights each selected donor’s relative level of risk for acute GVHD based on the patient’s leader genotype and published outcomes.13  The match status for each donor (E) is displayed through HapLogic match grades (A, allele match; P, potential match; L, allele mismatch; M, antigen mismatch)37  and leader match status (eg, MTT, a 3-letter nomenclature used to designate the leader of the patient’s mismatched HLA-B [first letter], the leader of the donor’s mismatched HLA-B [second letter], and the leader of the shared/matched HLA-B [third letter]).13 

Figure 2.

BLEAT provides a rank order for candidate stem cell sources (“donors”) to optimize therapy selection and patient outcomes. The BLEAT user interface (left) has an infographic (right) that is toggled via the Help button (A). On the interface, users enter HLA-B typing information for each subject. (A) The patient’s leader genotype (MM, MT, or TT) is automatically classified upon entering the patient’s HLA-B typing. Ambiguous alleles can sometimes contain rare leader variants; the tool highlights this possibility in red. (B) HLA-B allele typing for donors can be manually entered or automatically imported. Before or during sorting, (C) the leader match status for each donor (matched with the patient) is calculated. Once calculated, a donor selection guide (D) highlights each selected donor’s relative level of risk for acute GVHD based on the patient’s leader genotype and published outcomes.13  The match status for each donor (E) is displayed through HapLogic match grades (A, allele match; P, potential match; L, allele mismatch; M, antigen mismatch)37  and leader match status (eg, MTT, a 3-letter nomenclature used to designate the leader of the patient’s mismatched HLA-B [first letter], the leader of the donor’s mismatched HLA-B [second letter], and the leader of the shared/matched HLA-B [third letter]).13 

Close modal

Mapping subject alleles to leader designations

The user may input HLA-B typing for a patient and any number of unrelated donor(s), cord blood unit(s), or other cellular therapy sources (“subjects”). BLEAT assigns a leader designation of M (methionine) or T (threonine) to each input allele based on alignment of reference HLA-B sequences from IMGT1  and primary exon 1 sequence data from this study.

For HLA-B alleles that lack exon 1 information, BLEAT maps the allele to the first-field allele family. A second well-known phenomenon in HLA typing is the presence of ambiguous combinations of HLA-B alleles depending on typing methodology and platform output resolution. For example, HLA-B*07:37:01G is a name that indicates a group of alleles consisting of HLA-B*07:37:01 and HLA-B*07:37:02, both are alleles that encode the M-leader, typical of the HLA-B*07 family. However, a few exceptions require careful consideration. For example, in Figure 1, HLA-B*56:01:01G contains HLA-B*56:68, a rare B*56 M-leader variant.

Leader genotype and leader match status assignment of patient and stem cell sources

As described, BLEAT categorizes >7561 HLA-B alleles by the -21 leader position of each HLA-B allele for an unlimited number of subjects. In transplantation where a patient and any number of candidate donors or cord blood units are under consideration, BLEAT provides the user with the leader genotype (MM, MT, TT) of the patient and each stem cell source and compares the leader genotype of each candidate stem cell source with that of the patient. Pairwise comparisons are then translated into leader match statuses. When the patient and the stem cell source have the same HLA-B alleles (HLA-B-matched), their leader sequences are “leader-matched.” For patient/stem cell sources with 1 HLA-B mismatch, a 3-letter leader nomenclature is used to describe the leader of the patient’s mismatched HLA-B allele, the leader of the stem cell source’s mismatched allele, and the leader of the shared (matched) allele. For single HLA-B mismatched pairs, a total of 8 unique leader match status combinations are possible: TTT, TTM, TMT, MTT, MMM, MMT, MTM, TMM.13  BLEAT uses tooltips and interactive color-coding to help interpret this nomenclature.

Rank ordering of candidate unrelated donor stem cell sources

Finally, BLEAT provides an ordered list of the preferred unrelated donor stem cell sources based on published patient outcomes.13,16,17  Selection order is defined by the patient’s leader genotype: MM and TT patients benefit from leader-matched donors. MT patients have a choice in donors based on leader match status and the leader of the shared allele; graft-versus-host-disease (GVHD) risk is lowest with leader-matched donors who share a T-leader allele (MT donors). Donors who are leader mismatched and/or share an M-leader allele are associated with higher clinical risk. The prioritization of suitable donors for MT patients is illustrated (Figure 2D). For instance, leader-matched MT donors sharing an HLA-B*58 (T-leader) allele (MMT leader match status) are sorted to the top of the donor list, followed by leader-matched MT donors sharing an HLA-B*07 (M-leader) allele (TTM), then leader-mismatched donors (MTT, TTM). Additionally, the possibility of rare, minor P2 alleles within the HLA-B*07:02:01G (HLA-B*07:65, HLA-B*07:271, HLA-B*07:371, HLA-B*07:390) and HLA-B*58:01:01G (HLA-B*58:112) allele groups are highlighted in Figure 2.

Evaluation of BLEAT

We evaluated the performance of BLEAT for the assignment of leader genotype and patient/donor leader match status by testing a cohort of 1269 patient/unrelated donor transplant pairs (cohort 2) whose HLA typing, leader genotypes, and match statuses were previously validated using a T-SQL script on a SQL server.13  In the current study, the user was blinded to the previously assigned leader genotype and match status and applied BLEAT to independently assign the leader genotype to each patient and donor. Supplemental File 1  provides the input of cohort 2 for BLEAT. The leader match status of each transplant pair was determined via BLEAT and exported as supplemental File 2. BLEAT-assigned leader genotypes were concordant as was leader match status (leader-matched vs leader-mismatched; shared M vs shared T leader).

Broad race frequencies of HLA-B leader peptides

Study subjects for clinical and basic research may possess HLA-B typing at various degrees of resolution (fields) for alleles with known or unknown exon 1 sequences. We assessed the capability of BLEAT to assign leader genotypes for large-scale analysis of racially diverse subjects in the base cohort of 9.42 million registered donors selected solely on the basis of availability of HLA-B typing. First, we noted that the major leader allele T is observed in >68% of US registry donors in all race categories. We then examined the distribution of leader alleles and genotypes according to self-described race categories. The distribution of HLA-B leader genotypes differed significantly by donor race (P < 2.2 × 10-16) (Table 3). The frequency of the TT genotype was highest among Asians.

Table 4 is a race-based extension of Table 1 from cohort 1. Leader allele frequencies are similar between the base cohort (Table 3) and its subset cohort 1 (Table 4). M-leader allele frequency percentages are similar for White (32.2% [base cohort] vs 34.2% [cohort 1]), Hispanic or Latino (27.3% vs 30.3%), Asian or Pacific Islander (14.4% vs 14.8%), Black or African American (24.6% vs 25.5%), and American Indian or Alaskan Native (28.9% vs 27.7%) groups.

Access to BLEAT

The BLEAT user interface is accessible through https://bleader.nmdp.org. The codebase and application programming interface are provided at https://github.com/nmdp-bioinformatics/b-leader for flexible integration into any computational workflow and the ability to incorporate future added functionality. This tool has been tested by independent users since early 2020 and is supported as a service to the global transplant center and research community.

Elucidation of the HLA-B leader as a functional determinant in HIV, leukemia relapse, and stem cell transplantation has major implications for advancing knowledge of the immunobiology of HLA-B.13-18  For unrelated donor searches that yield HLA-B-mismatched donors, BLEAT is designed to assess whether the HLA-B mismatch is associated with lower (leader matched) or higher (leader mismatched) risk of clinically significant acute GVHD.10  When HLA-B is matched (patient and donor leader genotype are identical) and the donor mismatch is located at either HLA-A, -C, -DRB1, or -DQB1, the clinical significance of an HLA-A, -C, -DRB1, or -DQB1 mismatch depends on the patient/donor leader genotype.14 HLA-B-matched patients and donors with the MM leader genotype have increased risks of mortality, nonrelapse mortality, and acute GVHD relative to HLA-B-matched patients and donors with the TT leader genotype; single HLA-DRB1 and HLA-DQB1 mismatches are better tolerated than other mismatched loci when the patient and donor have the TT leader genotype.14  BLEAT may be used to determine HLA-B leader genotypes to compare risks and enhance donor selection.

Although >7561 unique HLA-B alleles are currently recognized,1  only 55% of these alleles have sequence information on the exon 1 rs1050458 dimorphism. Furthermore, exon 1 sequence data are not readily captured to facilitate clinical decision-making. We approached this unmet need via the validation of HLA-B leader sequences to known HLA-B alleles and the addition of new information on leader sequences for novel or rare HLA-B alleles. With a validated sequence library, we developed a tool to automate leader assignment to an HLA-B allele, derive leader genotypes for patients and donors, and compare leader match status in the patient to candidate transplant donors.

The current study was designed to address unmet needs in transplantation, namely, a user-friendly tool to accurately and efficiently determine the HLA-B leader to facilitate the selection of stem cell sources for patients. To this end, we mandated rigorous quality control for leader assignment and ascertainment of leader match status by the blinded study of a previously characterized cohort of patients and their transplant donors. We also leveraged a unique pool of 1 098 358 registered US donors with phased sequence data for HLA-B to verify the physical linkage of the -21 variant to the coding region and assigned HLA-B allele type. This cohort provided a unique opportunity to explore and extend known sequence information within exon 1 of HLA-B and better understand genotype frequencies of the leader across racially diverse US populations. We confirmed that >99.9% of 1 098 358 donors have leader alleles with either M or T at P2. We found striking differences in major and minor M- and T-leader frequencies across 1 098 358 donors representative of US continental populations. Although we present data on a large US population for HLA-B, the frequencies of HLA-B M and T leaders remain to be validated independently. The extension of these studies to other diverse populations remains an important research objective.

Although complete donor HLA matching at HLA-A, -B, -C, -DRB1 (HLA-8/8), and -DQB1 (HLA-10/10) lowers acute and chronic GVHD and improves survival after unrelated donor hematopoietic cell transplantation,10  many patients only have mismatched donor options.19  Current donor criteria consider the total number of HLA mismatches with preference given to donors mismatched at only 1 HLA gene (HLA-7/8 or HLA-9/10) over donors mismatched at 2 HLA loci (HLA-6/8 or HLA-8/10) to lower risks. When a search yields donors mismatched at HLA-A, -C, -DRB1, or -DQB1, the patient leader genotype defines the risks of mismatching; for all leader genotypes, a single HLA-DQB1 or -DRB1 mismatch is better tolerated compared with an HLA-A or -C mismatch, but risks depend on the combination of the leader genotype as well as the mismatched locus.14  Use of BLEAT enhances donor selection to lower overall risks to patients with HLA-mismatched donors. Although global leader genotype frequencies may differ across diverse populations, the HLA-B leader has clinical significance in all patients and donors; hence, the integration of BLEAT into clinical practice is expected to improve survival and lower transplant-associated risks for patients of all race/ethnic backgrounds.10,14  The specific mechanisms by which the leader influences GVHD and survival are the subjects of current research efforts; however, mechanisms include the potential role for NK and T-cell recognition of HLA-E. Class I leader peptides promote cell surface expression of HLA-E, where it may interact with T-cell and NK-cell receptors.11,12,16-18 

Beyond applications in transplantation, BLEAT has utility in population-based studies and research on diseases with strong associations to HLA-B (Figure 3). BLEAT may be used by clinical and basic researchers and stem cell registries to automate leader assignment in small- or large-scale studies to define the population genetics of HLA-B and help inform recruitment of future stem cell donors to support transplantation. Clinical and basic researchers may also use BLEAT to understand the functional significance of the leader particularly when information on HLA-B exon 1-linked variation is desirable.15,18  Because BLEAT requires HLA-B typing, it is subject to the limitations of both data captured in reference databases and data supplied by the user. High-resolution (allele) typing data are ideal; however, BLEAT accommodates typing at lower levels of resolution and when exon 1 is not directly sequenced. We demonstrate that general concordance of HLA-B leader assignments using first-field typing holds true for 2259 unique alleles encoded by 2.2 million haplotypes, but exceptions do occur. These exceptions give rise to both major and minor leader alleles in some allele families more than others. When HLA-B typing is acquired at low-resolution (ie, first-field), BLEAT signals the potential existence of an alternative leader based on the frequencies of major and minor leader alleles in a given allele family. As new HLA-B alleles are discovered, delineation of their exon 1 sequences will continue to be an important element of annotating the full spectrum of HLA-B variation.

Figure 3.

Utility of BLEAT for diverse applications. The assignment of the leader genotype to the HLA-B allele may be useful for downstream applications by HLA laboratories, clinical researchers, donor registries, cell source banks, and health care workers. For HLA laboratories involved in basic research, these include exploration of allele diversity and haplotype linkage involving HLA-B to examine HLA structure and function. Laboratories may also supply this information for clinical research pursuits in understanding historical patient risk or modeling future risk to improve outcomes. Health care workers can apply BLEAT toward the selection of optimal donors based on transplant risk that can also be facilitated by donor registries or cell source banks. The assessment of HLA-B allele diversity can also aid in ensuring that therapies are available for diverse populations and enable individualized treatment.

Figure 3.

Utility of BLEAT for diverse applications. The assignment of the leader genotype to the HLA-B allele may be useful for downstream applications by HLA laboratories, clinical researchers, donor registries, cell source banks, and health care workers. For HLA laboratories involved in basic research, these include exploration of allele diversity and haplotype linkage involving HLA-B to examine HLA structure and function. Laboratories may also supply this information for clinical research pursuits in understanding historical patient risk or modeling future risk to improve outcomes. Health care workers can apply BLEAT toward the selection of optimal donors based on transplant risk that can also be facilitated by donor registries or cell source banks. The assessment of HLA-B allele diversity can also aid in ensuring that therapies are available for diverse populations and enable individualized treatment.

Close modal

At the population level, the assignment of the HLA-B leader genotypes for registered unrelated donors and cord blood units can be automated upon entry of validated HLA data into the registry systems. Such application is already underway for the automatic interpretation of HLA-B typing entries in the Stem Cell Therapeutic Outcomes Database at the Center for International Blood and Marrow Transplant Research and for the Be The Match Registry operated by the NMDP.

The intended application of BLEAT in transplantation is to aid the user in handling HLA genetic information rather than replace clinical decision-making. Selection of stem cell sources for hematopoietic cell transplantation is complex. Individual HLA features (eg, the HLA-B leader and HLA-DPB1 mismatching)31  and non-HLA characteristics (eg, donor age)32  affect the success of transplantation; however, the impact of HLA and non-HLA factors in combination remain important research questions. As potential synergistic effects of HLA and non-HLA features are elucidated, characteristics that define an “ideal” donor for a given patient provide an important basis for building the next generation of tools to facilitate donor selection. To this end, BLEAT components are provided for connection to existing pipelines and interfaces.

We envision a flexible user interface in which therapy selection criteria and sorting are provided based on evidence-based research guidelines in the HLA-mismatch setting. Exact criteria and sorting rules can be designated by the user as desired for different disease categories, therapy options, and other needs in the application. Future refinements to BLEAT to accommodate prioritization of preferable cord blood units based on the leader are forthcoming. In addition, a display of the results from the HLA-DPB1 TCE Prediction Tool for assessment of HLA-DPB1 permissive mismatches31  along with BLEAT may connect to a consolidated interface or other platform and enable users to view both HLA-B and HLA-DPB1 risk calculations among available stem cell source options. The utility of a comprehensive HLA genetic tool has particular relevance when there is clinical urgency to proceed to transplantation and when the only available stem cell sources are HLA-mismatched. With the continued development of new transplant regimens, the impact of HLA and non-HLA factors in different clinical settings will also be an important area for future investigation.

Each year, approximately 18 000 patients need an unrelated or cord blood hematopoietic cell transplant to cure a life-threatening blood disorder.33,34  Many of these patients will have only HLA-mismatched stem cell sources.20  New information on the role of the HLA-B leader in HLA-mismatched transplantation showcases the clinical importance of variation located outside of the peptide-binding region. The availability of BLEAT for the direct translation of HLA-B sequence variation to donor selection and pretransplant risk assessment represents a patient-directed approach for improving clinical outcomes.

The authors thank Caroline McKallor for data support and Pradeep Bashyal for software infrastructure.

This work was supported by grants from the National Institutes of Health National Institute of Allergy and Infectious Diseases (AI069197) and National Cancer Institute (CA231838 and CA100019) to E.W.P. and C.M. Bioinformatics methods development was supported in part by the Office of Naval Research Grant N00014-19-1-2705 to National Marrow Donor Program/Be The Match. Funding agencies had no role in study design, data collection and analysis, the decision to submit the manuscript for publication, or preparation of the manuscript.

Contribution: E.W.P. and Y.-T.B. conceived and supervised the study; R.S. developed BLEAT; R.S. conducted sequence analysis on cohort 1 and leader annotation on cohort 2; M.J.M. conducted frequency analysis on the base cohort; and all authors contributed to the manuscript figures, tables, and text.

Conflict-of-interest disclosure: The authors declare no competing financial interests.

Data sharing statement: The user interface and code for BLEAT and leader analysis are available at https://bleader.nmdp.org and https://github.com/nmdp-bioinformatics/b-leader, and supporting data and details are provided in supplemental Materials and Methods. Genbank accession numbers are MH173353, MG756798, MH973951, and MG769755 for novel sequences. For more information, please contact bioinformatics-web@nmdp.org.

Correspondence: Yung-Tsi Bolon, Immunobiology and Bioinformatics Research, National Marrow Donor Program/Be The Match, Center for International Blood and Marrow Transplant Research, 500 N 5th St, Minneapolis, MN 55455; e-mail: ybolon@nmdp.org.

1.
Robinson
J
,
Halliwell
JA
,
Hayhurst
JD
,
Flicek
P
,
Parham
P
,
Marsh
SGE
.
The IPD and IMGT/HLA database: allele variant databases
.
Nucleic Acids Res.
2015
;
43
(
database issue
):
D423
-
D431
.
2.
Parham
P
,
Adams
EJ
,
Arnett
KL
.
The origins of HLA-A,B,C polymorphism
.
Immunol Rev.
1995
;
143
(
1
):
141
-
180
.
3.
Parham
P
,
Guethlein
LA
.
Genetics of natural killer cells in human health, disease, and survival
.
Annu Rev Immunol.
2018
;
36
(
1
):
519
-
548
.
4.
Greco
A
,
De Virgilio
A
,
Ralli
M
, et al
.
Behçet’s disease: new insights into pathophysiology, clinical features and treatment options
.
Autoimmun Rev.
2018
;
17
(
6
):
567
-
575
.
5.
Hanscombe
KB
,
Morris
DL
,
Noble
JA
, et al
.
Genetic fine mapping of systemic lupus erythematosus MHC associations in Europeans and African Americans
.
Hum Mol Genet.
2018
;
27
(
21
):
3813
-
3824
.
6.
Hanson
AL
,
Vukcevic
D
,
Leslie
S
, et al;
International Genetics of Ankylosing Spondylitis Consortium
.
Epistatic interactions between killer immunoglobulin-like receptors and human leukocyte antigen ligands are associated with ankylosing spondylitis
.
PLoS Genet.
2020
;
16
(
8
):
e1008906
.
7.
Illing
PT
,
Purcell
AW
,
McCluskey
J
.
The role of HLA genes in pharmacogenomics: unravelling HLA associated adverse drug reactions
.
Immunogenetics.
2017
;
69
(
8-9
):
617
-
630
.
8.
Leccese
P
,
Alpsoy
E
.
Behçet’s disease: an overview of etiopathogenesis
.
Front Immunol.
2019
;
10
:
1067
.
9.
Terao
C
,
Brynedal
B
,
Chen
Z
, et al
.
Distinct HLA associations with rheumatoid arthritis subsets defined by serological subphenotype [published correction appears in Am J Hum Genet. 2019;105(4):880]
.
Am J Hum Genet.
2019
;
105
(
3
):
616
-
624
.
10.
Braud
V
,
Jones
EY
,
McMichael
A
.
The human major histocompatibility complex class Ib molecule HLA-E binds signal sequence-derived peptides with primary anchor residues at positions 2 and 9
.
Eur J Immunol.
1997
;
27
(
5
):
1164
-
1169
.
11.
Horowitz
A
,
Djaoud
Z
,
Nemat-Gorgani
N
, et al
.
Class I HLA haplotypes form two schools that educate NK cells in different ways
.
Sci Immunol.
2016
;
1
(
3
):
eaag1672
.
12.
Braud
VM
,
Allan
DS
,
O’Callaghan
CA
, et al
.
HLA-E binds to natural killer cell receptors CD94/NKG2A, B and C
.
Nature.
1998
;
391
(
6669
):
795
-
799
.
13.
Petersdorf
EW
,
Carrington
M
,
O’hUigin
C
, et al;
International Histocompatibility Working Group in Hematopoietic Cell Transplantation
.
Role of HLA-B exon 1 in graft-versus-host disease after unrelated haemopoietic cell transplantation: a retrospective cohort study
.
Lancet Haematol.
2020
;
7
(
1
):
e50
-
e60
.
14.
Hallner
A
,
Bernson
E
,
Hussein
BA
, et al
.
The HLA-B -21 dimorphism impacts on NK cell education and clinical outcome of immunotherapy in acute myeloid leukemia
.
Blood.
2019
;
133
(
13
):
1479
-
1488
.
15.
Merino
AM
,
Song
W
,
He
D
, et al
.
HLA-B signal peptide polymorphism influences the rate of HIV-1 acquisition but not viral load
.
J Infect Dis.
2012
;
205
(
12
):
1797
-
1805
.
16.
Petersdorf
EW
,
Gooley
T
,
Volt
F
, et al
.
Use of the HLA-B leader to optimize cord-blood transplantation [published online ahead of print 29 October 2020]
.
Haematologica.
2020
.
17.
Petersdorf
EW
,
Stevenson
P
,
Bengtsson
M
, et al
.
HLA-B leader and survivorship after HLA-mismatched unrelated donor transplantation
.
Blood.
2020
;
136
(
3
):
362
-
369
.
18.
Ramsuran
V
,
Naranbhai
V
,
Horowitz
A
, et al
.
Elevated HLA-A expression impairs HIV control through inhibition of NKG2A-expressing cells
.
Science.
2018
;
359
(
6371
):
86
-
90
.
19.
Borrego
F
,
Ulbrecht
M
,
Weiss
EH
,
Coligan
JE
,
Brooks
AG
.
Recognition of human histocompatibility leukocyte antigen (HLA)-E complexed with HLA class I signal sequence-derived peptides by CD94/NKG2 confers protection from natural killer cell-mediated lysis
.
J Exp Med.
1998
;
187
(
5
):
813
-
818
.
20.
Gragert
L
,
Eapen
M
,
Williams
E
, et al
.
HLA match likelihoods for hematopoietic stem-cell grafts in the U.S. registry
.
N Engl J Med.
2014
;
371
(
4
):
339
-
348
.
21.
Dos Santos Francisco
R
,
Buhler
S
,
Nunes
JM
, et al
.
HLA supertype variation across populations: new insights into the role of natural selection in the evolution of HLA-A and HLA-B polymorphisms
.
Immunogenetics.
2015
;
67
(
11-12
):
651
-
663
.
22.
Fernandez Vina
MA
,
Hollenbach
JA
,
Lyke
KE
, et al
.
Tracking human migrations by the analysis of the distribution of HLA alleles, lineages and haplotypes in closed and open populations
.
Philos Trans R Soc Lond B Biol Sci.
2012
;
367
(
1590
):
820
-
829
.
23.
Lienert
K
,
Parham
P
.
Evolution of MHC class I genes in higher primates
.
Immunol Cell Biol.
1996
;
74
(
4
):
349
-
356
.
24.
Sanchez-Mazas
A
,
Černý
V
,
Di
D
, et al
.
The HLA-B landscape of Africa: signatures of pathogen-driven selection and molecular identification of candidate alleles to malaria protection
.
Mol Ecol.
2017
;
26
(
22
):
6238
-
6252
.
25.
Gragert
L
,
Madbouly
A
,
Freeman
J
,
Maiers
M
.
Six-locus high resolution HLA haplotype frequencies derived from mixed-resolution DNA typing for the entire US donor registry
.
Hum Immunol.
2013
;
74
(
10
):
1313
-
1320
.
26.
ProQuest
.
United States. Bureau of the Census., United States. Economics and Statistics Administration. ProQuest statistical abstract of the United States: 2015.
Lanham, MD
:
Bernan Press
;
2015
.
27.
Rhoads
A
,
Au
KF
.
PacBio sequencing and its applications
.
Genomics Proteomics Bioinformatics.
2015
;
13
(
5
):
278
-
289
.
28.
National Marrow Donor Program
. Multiple allele code service API documentation. https://hml.nmdp.org/mac/raml/allele-code.html. Published
2020
. Accessed 2 August 2021.
29.
R Core Team
. R: a language and environment for statistical computing.
2019
.
30.
Kent
WJ
,
Sugnet
CW
,
Furey
TS
, et al
.
The human genome browser at UCSC
.
Genome Res.
2002
;
12
(
6
):
996
-
1006
.
31.
Pidala
J
,
Lee
SJ
,
Ahn
KW
, et al
.
Nonpermissive HLA-DPB1 mismatch increases mortality after myeloablative unrelated allogeneic hematopoietic cell transplantation
.
Blood.
2014
;
124
(
16
):
2596
-
2606
.
32.
Shaw
BE
,
Logan
BR
,
Spellman
SR
, et al
.
Development of an unrelated donor selection score predictive of survival after HCT: donor age matters most
.
Biol Blood Marrow Transplant.
2018
;
24
(
5
):
1049
-
1056
.
33.
Health Resources & Services Administration
. Frequently asked questions (FAQ). https://bloodstemcell.hrsa.gov/about/faqs. Published
2020
. Accessed 2 August 2021.
34.
Nietfeld
JJ
,
Pasquini
MC
,
Logan
BR
,
Verter
F
,
Horowitz
MM
.
Lifetime probabilities of hematopoietic stem cell transplantation in the U.S
.
Biol Blood Marrow Transplant.
2008
;
14
(
3
):
316
-
322
.
35.
Marsh
SGEB.
_prot. https://raw.githubusercontent.com/ANHIG/IMGTHLA/Latest/alignments/B_prot.txt. Updated 15 October 15,
2020
. Accessed 2 August 2021.
36.
National Marrow Donor Program
. Multiple allele code API. https://hml.nmdp.org/mac/api/expand?typing=B*56:01:01G. Accessed 2 August 2021.
37.
Dehn
J
,
Setterholm
M
,
Buck
K
, et al
.
HapLogic: a predictive human leukocyte antigen-matching algorithm to enhance rapid identification of the optimal unrelated hematopoietic stem cell sources for transplantation
.
Biol Blood Marrow Transplant.
2016
;
22
(
11
):
2038
-
2046
.

Author notes

The full-text version of this article contains a data supplement.

Supplemental data