Significant variation between SNP-based HLA imputations in diverse populations: the last mile is the hardest.

Four single nucleotide polymorphism (SNP)-based human leukocyte antigen (HLA) imputation methods (e-HLA, HIBAG, HLA*IMP:02 and MAGPrediction) were trained using 1000 Genomes SNP and HLA genotypes and assessed for their ability to accurately impute molecular HLA-A, -B, -C and -DRB1 genotypes in the Human Genome Diversity Project cell panel. Imputation concordance was high (>89%) across all methods for both HLA-A and HLA-C, but HLA-B and HLA-DRB1 proved generally difficult to impute. Overall, <27.8% of subjects were correctly imputed for all HLA loci by any method. Concordance across all loci was not enhanced via the application of confidence thresholds; reliance on confidence scores across methods only led to noticeable improvement (+3.2%) for HLA-DRB1. As the HLA complex is highly relevant to the study of human health and disease, a standardized assessment of SNP-based HLA imputation methods is crucial for advancing genomic research. Considerable room remains for the improvement of HLA-B and especially HLA-DRB1 imputation methods, and no imputation method is as accurate as molecular genotyping. The application of large, ancestrally diverse HLA and SNP reference data sets and multiple imputation methods has the potential to make SNP-based HLA imputation methods a tractable option for determining HLA genotypes.

Revealing complete complex KIR haplotypes phased by long-read sequencing technology.

The killer cell immunoglobulin-like receptor (KIR) region of human chromosome 19 contains up to 16 genes for natural killer (NK) cell receptors that recognize human leukocyte antigen (HLA)/peptide complexes and other ligands. The KIR proteins fulfill functional roles in infections, pregnancy, autoimmune diseases and transplantation. However, their characterization remains a constant challenge. Not only are the genes highly homologous due to their recent evolution by tandem duplications, but the region is structurally dynamic due to frequent transposon-mediated recombination. A sequencing approach that precisely captures the complexity of KIR haplotypes for functional annotation is desirable. We present a unique approach to haplotype the KIR loci using single-molecule, real-time (SMRT) sequencing. Using this method, we have-for the first time-comprehensively sequenced and phased sixteen KIR haplotypes from eight individuals without imputation. The information revealed four novel haplotype structures, a novel gene-fusion allele, novel and confirmed insertion/deletion events, a homozygous individual, and overall diversity for the structural haplotypes and their alleles. These KIR haplotypes augment our existing knowledge by providing high-quality references, evolutionary informers, and source material for imputation. The haplotype sequences and gene annotations provide alternative loci for the KIR region in the human genome reference GrCh38.p8.

Validation of statistical imputation of allele-level multilocus phased genotypes from ambiguous HLA assignments.

Genetic matching for loci in the human leukocyte antigen (HLA) region between a donor and a patient in hematopoietic stem cell transplantation (HSCT) is critical to outcome; however, methods for HLA genotyping of donors in unrelated stem cell registries often yield results with allelic and phase ambiguity and/or do not query all clinically relevant loci. We present and evaluate a statistical method for in silico imputation of HLA alleles and haplotypes in large ambiguous population data from the Be The Match(®) Registry. Our method builds on haplotype frequencies estimated from registry populations and exploits patterns of linkage disequilibrium (LD) across HLA haplotypes to infer high resolution HLA assignments. We performed validation on simulated and real population data from the Registry with non-trivial ambiguity content. While real population datasets caused some predictions to deviate from expectation, validations still showed high percent recall for imputed results with average recall >76% when imputing HLA alleles from registry data. We simulated ambiguity generated by several HLA genotyping methods to evaluate the imputation performance on several levels of typing resolution. On average, imputation percent recall of allele-level HLA haplotypes was >95% for allele-level typing, >92% for intermediate resolution typing and >58% for serology (low-resolution) typing. Thus, allele-level HLA assignments can be imputed through the application of a set of statistical and population genetics inferences and with knowledge of haplotype frequencies and self-identified race and ethnicities.

Extensive haplotype diversity in African American mothers and their cord blood units.

HLA-A, -B, -C, -DRB1, -DQB1 assignments were obtained for 374 pairs of African American mothers and their umbilical cord blood units (CBU) by DNA sequencing. An algorithm developed by the National Marrow Donor Program was used to assign 1122 haplotypes by segregation. Seventy percent of the haplotypes carried assignments at all five loci. In the remainder, alleles at various loci, most often DQB1 in 48% of the haplotypes with a missing assignment, could not be assigned due to sharing of both alleles by mother and CBU. There were 652 haplotypes carrying a unique combination of alleles at the five loci; the majority (74%) were singletons. Novel B∼C and DRB1~DQB1 associations were observed. The results show the genetic diversity in this population and provide validation for a publically available tool for pedigree analysis. Our observations underscore the need for procurement of increased numbers of units in the national cord blood inventory in order to identify matching donors for all patients requiring hematopoietic stem cell transplantation.

Comparative validation of computer programs for haplotype frequency estimation from donor registry data.

Estimation of human leukocyte antigen (HLA) haplotype frequencies from unrelated stem cell donor registries presents a challenge because of large sample sizes and heterogeneity of HLA typing data. For the 14th International HLA and Immunogenetics Workshop, five bioinformatics groups initiated the 'Registry Diversity Component' aiming to cross-validate and improve current haplotype estimation tools. Five datasets were derived from different donor registries and then used as input for five different computer programs for haplotype frequency estimation. Because of issues related to heterogeneity and complexity of HLA typing data identified in the initial phase, the same five implementations, and two new ones, were used on simulated datasets in a controlled experiment where the correct results were known a priori. These datasets contained various fractions of missing HLA-DR modeled after European haplotype frequencies. We measured the contribution of sampling fluctuation and estimation error to the deviation of the frequencies from their true values, finding equivalent contributions of each for the chosen samples. Because of patient-directed activities, selective prospective typing strategies and the variety and evolution of typing technology, some donors have more complete and better HLA data. In this setting, we show that restricting estimation to fully typed individuals introduces biases that could be overcome by including all donors in frequency estimation. Our study underlines the importance of critical review and validation of tools in registry-related activity and provides a sustainable framework for validating the computational tools used. Accurate frequencies are essential for match prediction to improve registry operations and to help more patients identify suitably matched donors.

Genotype List String: a grammar for describing HLA and KIR genotyping results in a text string.

Knowledge of an individual's human leukocyte antigen (HLA) genotype is essential for modern medical genetics, and is crucial for hematopoietic stem cell and solid-organ transplantation. However, the high levels of polymorphism known for the HLA genes make it difficult to generate an HLA genotype that unambiguously identifies the alleles that are present at a given HLA locus in an individual. For the last 20 years, the histocompatibility and immunogenetics community has recorded this HLA genotyping ambiguity using allele codes developed by the National Marrow Donor Program (NMDP). While these allele codes may have been effective for recording an HLA genotyping result when initially developed, their use today results in increased ambiguity in an HLA genotype, and they are no longer suitable in the era of rapid allele discovery and ultra-high allele polymorphism. Here, we present a text string format capable of fully representing HLA genotyping results. This Genotype List (GL) String format is an extension of a proposed standard for reporting killer-cell immunoglobulin-like receptor (KIR) genotype data that can be applied to any genetic data that use a standard nomenclature for identifying variants. The GL String format uses a hierarchical set of operators to describe the relationships between alleles, lists of possible alleles, phased alleles, genotypes, lists of possible genotypes, and multilocus unphased genotypes, without losing typing information or increasing typing ambiguity. When used in concert with appropriate tools to create, exchange, and parse these strings, we anticipate that GL Strings will replace NMDP allele codes for reporting HLA genotypes.

Common and well-documented HLA alleles: 2012 update to the CWD catalogue.

We have updated the catalogue of common and well-documented (CWD) human leukocyte antigen (HLA) alleles to reflect current understanding of the prevalence of specific allele sequences. The original CWD catalogue designated 721 alleles at the HLA-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, and -DPB1 loci in IMGT (IMmunoGeneTics)/HLA Database release 2.15.0 as being CWD. The updated CWD catalogue designates 1122 alleles at the HLA-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, -DPA1 and -DPB1 loci as being CWD, and represents 14.3% of the HLA alleles in IMGT/HLA Database release 3.9.0. In particular, we identified 415 of these alleles as being 'common' (having known frequencies) and 707 as being 'well-documented' on the basis of ~140,000 sequence-based typing observations and available HLA haplotype data. Using these allele prevalence data, we have also assigned CWD status to specific G and P designations. We identified 147/151 G groups and 290/415 P groups as being CWD. The CWD catalogue will be updated on a regular basis moving forward, and will incorporate changes to the IMGT/HLA Database as well as empirical data from the histocompatibility and immunogenetics community. This version 2.0.0 of the CWD catalogue is available online at cwd.immunogenomics.org, and will be integrated into the Allele Frequencies Net Database, the IMGT/HLA Database and National Marrow Donor Program's bioinformatics web pages.

16(th) IHIW: immunogenomic data-management methods. report from the immunogenomic data analysis working group (IDAWG).

SUMMARY: The goal of the immunogenomic data analysis working group (IDAWG) is to facilitate the consistent analysis of HLA and KIR data, and the sharing of those data among the immunogenomic and larger genomic communities. However, the data management approaches currently applied by immunogenomic researchers are not widely discussed or reported in the literature, and the effect of different approaches on data analyses is not known. With ASHI's support, the IDAWG developed a 45 question survey on HLA and KIR data generation, data management and data analysis practices. Survey questions detailed the loci genotyped, typing systems used, nomenclature versions reported, computer operating systems and software used to manage and transmit data, the approaches applied to resolve HLA ambiguity and the methods used for basic population-level analyses. Respondents were invited to demonstrate their HLA ambiguity resolution approaches in simulated data sets. By May 2012, 156 respondents from 35 nations had completed the survey. These survey respondents represent a broad sampling of the Immunogenomic community; 52% were European, 30% North American, 10% Asian, 4% South American and 4% from the Pacific. The project will continue in conjunction with the 17th Workshop, with the aim of developing community data sharing standards, ambiguity resolution documentation formats, single-task data Management tools and novel data analysis methods and applications. While additional project details and plans for the 17th IHIW will be forthcoming, we welcome the input and participation in these projects from the histocompatibility and immunogenetics community.

16(th) IHIW: global analysis of registry HLA haplotypes from 20 million individuals: report from the IHIW Registry Diversity Group.

This project has the goal to validate bioinformatics methods and tools for HLA haplotype frequency analysis specifically addressing unique issues of haematopoietic stem cell registry data sets. In addition to generating new methods and tools for the analysis of registry data sets, the intent is to produce a comprehensive analysis of HLA data from 20 million donors from the Bone Marrow Donors Worldwide (BMDW) database. This report summarizes the activity on this project as of the 16IHIW meeting in Liverpool.

A community standard for immunogenomic data reporting and analysis: proposal for a STrengthening the REporting of Immunogenomic Studies statement.

Modern high-throughput HLA and KIR typing technologies are generating a wealth of immunogenomic data with the potential to revolutionize the fields of histocompatibility and immune-related disease association and population genetic research, much as SNP-based approaches have revolutionized association research. The STrengthening the REporting of Genetic Association studies (STREGA) statement provides community-based data reporting and analysis standards for genomic disease-association studies, identifying specific areas in which adoption of reporting guidelines can improve the consistent interpretation of genetic studies. While aspects of STREGA can be applied to immunogenomic studies, HLA and KIR research requires additional consideration, as the high levels of polymorphism associated with immunogenomic data pose unique methodological and computational challenges to the synthesis of information across datasets. Here, we outline the principle challenges to consistency in immunogenomic studies, and propose that an immunogenomic-specific analog to the STREGA statement, a STrengthening the REporting of Immunogenomic Studies (STREIS) statement, be developed as part of the 16th International HLA and Immunogenetics Workshop. We propose that STREIS extends at least four of the 22 elements of the STREGA statement to specifically address issues pertinent to immunogenomic data: HLA and KIR nomenclature, data-validation, ambiguity resolution, and the analysis of highly polymorphic genetic systems. As with the STREGA guidelines, the intent behind STREIS is not to dictate the design of immunogenomic studies, but to ensure consistent and transparent reporting of research, facilitating the synthesis of HLA and KIR data across studies.

Genetic differentiation of Jewish populations.

The Jewish diaspora can be viewed as a natural process in population dispersion and differentiation. We extend genetic studies on the Jewish diaspora to an analysis of human leukocyte antigen (HLA) haplotype distributions in the Jewish peoples, and show the value of this information for the design of Jewish marrow donor registries. HLA data from the Hadassah Bone Marrow Registry having parental country-of-origin information comprise samples of geographically discrete regions. We analyzed the HLA allele and haplotype frequencies for each national sample using population genetic and clustering methods. Population differentiation among diaspora populations was shown on the basis of HLA haplotype frequencies, including differences within the more recently diverged European groups. A method of haplotype and population clustering showed patterns of unique haplotype affinities associated with specific Jewish populations. The evidence showed that diaspora Jewish populations can be sorted into distinct clades of which the Ashkenazi are but one. Relationships among Jewish populations are interpretable in light of the historical record. We suggest that a major contributing factor to the genetic divergence between Jewish groups may have been admixture with local host populations, while, at the same time, threads of Eastern Mediterranean ancestry remain evident.

The HLA dictionary 2008: a summary of HLA-A, -B, -C, -DRB1/3/4/5, and -DQB1 alleles and their association with serologically defined HLA-A, -B, -C, -DR, and -DQ antigens.

The 2008 report of the human leukocyte antigen (HLA) data dictionary presents serologic equivalents of HLA-A, -B, -C, -DRB1, -DRB3, -DRB4, -DRB5, and -DQB1 alleles. The dictionary is an update of the one published in 2004. The data summarize equivalents obtained by the World Health Organization Nomenclature Committee for Factors of the HLA System, the International Cell Exchange, UCLA, the National Marrow Donor Program, recent publications, and individual laboratories. The 2008 edition includes information on 832 new alleles (685 class I and 147 class II) and updated information on 766 previously listed alleles (577 class I and 189 class II). The tables list the alleles with remarks on the serologic patterns and the equivalents. The serological equivalents are listed as expert assigned types, and the data are useful for identifying potential stem cell donors who were typed by either serology or DNA-based methods. The tables with HLA equivalents are available as a searchable form on the IMGT/HLA database Web site (http://www.ebi.ac.uk/imgt/hla/dictionary.html).

Four-locus high-resolution HLA typing in a sample of Mexican Americans.

Mexicans are the most common minority population of the United States. From a sample of 553 bone marrow donor registrants of self-described Mexican ancestry, human leukocyte antigen (HLA) loci A, C, B and DRB1 were typed by high resolution sequence based typing (SBT) methods. A total of 47, 34, 76 and 46 distinct alleles at A, C, B and DRB1 respectively were identified, including 3 new alleles. The four-locus haplotype frequency distribution was extremely skewed with only 53.9% of 1106 chromosomes present with more than one estimated copy. Haplotypes of Native American origin were identified. These data form an initial basis for determining the requirements for an adequate donor pool for stem cell transplantation in this population.

Next generation sequencing characterizes the extent of HLA diversity in an Argentinian registry population.

Next generation DNA sequencing is used to determine the HLA-A, -B, -C, -DRB1, and -DQB1 assignments of 1472 unrelated volunteers for the unrelated donor registry in Argentina. The analysis characterized all HLA exons and introns for class I alleles; at least exons 2, 3 for HLA-DRB1; and exons 2 to 6 for HLA-DQB1. Of the distinct alleles present, there are 330 class I and 98 class II. The majority (~98%) of the cumulative allele frequency at each locus is contributed by alleles that appear at a frequency of at least 1 in 1000. Fourteen (18.2%) of the 77 novel class I and II alleles carry nonsynonymous variation within their exons; 52 (75.4%) class I novel alleles carry only single, apparently random, nucleotide variation within their introns/untranslated regions. Alleles encoding protein variation not usually detected by typing focused only on the exons encoding the antigen recognition domain are 1.0% of the class I assignments and 7.3% of the class II assignments (predominantly DQB1*02:02:01, DQB1*03:19:01, and DRB1*14:54:01). Updates to the common and well documented list of alleles include 10 alleles previously thought to be uncommon but that are found at least 30 times. Five locus haplotypes estimated using the expectation-maximization algorithm as present 3 or more times total 187. While the known HLA diversity continues to increase, the conservation of known allele sequences is remarkable. Overall, the HLA diversity observed in the Argentinian population reflects its European and Native American ancestry.

World Marrow Donor Association guidelines for use of HLA nomenclature and its validation in the data exchange among hematopoietic stem cell donor registries and cord blood banks.

Since the advent of the European Marrow Donor Information System in the first half of the last decade, fully automated data exchange between registry computer systems has been playing an ever-increasing role in the international search for unrelated donors of blood progenitor cells. This exchange, however, was hampered by different local conventions used to present HLA data and complicated by the need to extend the official WHO nomenclature to accommodate the registries' information systems and to cross-validate HLA data obtained with different methods and/or at different loci. The guidelines presented here have been developed by the World Marrow Donor Association to standardize the nomenclature to be used and the validation checks to be applied in the international electronic exchange of HLA-typing data among unrelated volunteer hematopoietic stem cell donor registries and umbilical cord blood banks. Two reference web sites have been designated to maintain and update the approved HLA nomenclature and all the ancillary information needed by the conventions described here.

Unrelated donor search prognostic score to support early HLA consultation and clinical decisions.

A simple scoring system that can provide a quick search prognosis at the onset of an adult unrelated donor (URD) search could be a useful tool for transplant physicians. We aimed to determine whether patient human leukocyte Ag genotype frequency (GF) could be used as a surrogate measure of whether or not a potential 10/10 and/or 9/10 URD in the Be The Match Registry (BTMR) can be identified for the patient. GF was assigned on a training data set of 2410 patients that searched the BTMR using the reported ethnic group. A proportional odds model was used to correlate GF with defined search productivity categories as follows: 'Good' (>2 10/10), 'Fair' (1-2 10/10 or No 10/10 and >2 9/10) or 'Poor' (No 10/10 and <3 9/10). A second cohort (n=2411) was used to calculate the concordance by the ethnic group in all three categories. In addition, we validated against an independent cohort (n=1344) resolved as having a 10/10 or 9/10 matched URD. Across the ethnic groups, >90% of cases with 'Good' GF prognosis, 20-26% 'Fair' and <10% 'Poor' had a 10/10 URD. Although not a replacement for an actual URD search, GF offers a quick way for transplant physicians to get an indication of the likely search outcome.

Identification of high-risk amino-acid substitutions in hematopoietic cell transplantation: a challenging task.

Allogeneic hematopoietic cell transplantation (HCT) offers the potential to cure hematologic malignancies. In the absence of an HLA-matched donor, HLA mismatched unrelated donors may be used, although risks of GvHD and treatment-related mortality (TRM) are higher. Identification and avoidance of amino-acid substitution and position types (AASPT) conferring higher risks of TRM and GvHD would potentially improve the success of transplantation from single HLA mismatched unrelated donors. Using random forest and logistic regression analyses, we identified 19 AASPT associated with greater risks for at least one adverse transplant outcome: grade III-IV acute GvHD, TRM, lower disease-free survival or worse overall survival relative to HLA-matched unrelated donors and to other AASPT. When tested in an independent validation cohort of 3530 patients, none of the AASPT from the training set were validated as high risk, however. Review of the literature shows that failure to validate original observations is the rule and not the exception in immunobiology and emphasizes the importance of independent validation before clinical application. Our current data do not support avoiding any specific class I AASPT for unrelated donors. Additional studies should be performed to fully understand the role of AASPT in HCT outcomes.

A comparative reference study for the validation of HLA-matching algorithms in the search for allogeneic hematopoietic stem cell donors and cord blood units.

The accuracy of human leukocyte antigen (HLA)-matching algorithms is a prerequisite for the correct and efficient identification of optimal unrelated donors for patients requiring hematopoietic stem cell transplantation. The goal of this World Marrow Donor Association study was to validate established matching algorithms from different international donor registries by challenging them with simulated input data and subsequently comparing the output. This experiment addressed three specific aspects of HLA matching using different data sets for tasks of increasing complexity. The first two tasks targeted the traditional matching approach identifying discrepancies between patient and donor HLA genotypes by counting antigen and allele differences. Contemporary matching procedures predicting the probability for HLA identity using haplotype frequencies were addressed by the third task. In each task, the identified disparities between the results of the participating computer programs were analyzed, classified and quantified. This study led to a deep understanding of the algorithms participating and finally produced virtually identical results. The unresolved discrepancies total to less than 1%, 4% and 2% for the three tasks and are mostly because of individual decisions in the design of the programs. Based on these findings, reference results for the three input data sets were compiled that can be used to validate future matching algorithms and thus improve the quality of the global donor search process.

The HLA dictionary 1999: a summary of HLA-A, -B, -C, -DRB1/3/4/5, -DQB1 alleles and their association with serologically defined HLA-A, -B, -C, -DR, and -DQ antigens.

This report presents serologic equivalents of 90 HLA-A, 190 HLA-B, and 145 HLA-DRB1 alleles. The equivalents cover over 70% of the presently identified HLA-A, -B, and -DRB1 alleles. The dictionary is an update of the one published in 1997 and now also includes equivalents for HLA-C, DRB3, DRB4, DRB5, and DQB1 alleles. The data summarize information obtained by the WHO HLA Nomenclature Committee, the International Cell Exchange (UCLA), the National Marrow Donor Program (NMDP), and by individual laboratories. In addition, a listing is provided of alleles that are expressed as antigens with serologic reaction patterns that differ from the well-established HLA specificities and that often lack official WHO nomenclature. The provided equivalents will be useful in guiding searches for unrelated donors in which patients and/or potential donors are typed by either serology or DNA-based methods. These equivalents will also serve typing and matching procedures for organ transplant programs where HLA typings from donors and from recipients on waiting lists represent mixtures of serologic and molecular typings. Some guidelines are provided for the use of appropriate WHO HLA nomenclature for serologic typings and for generic and allele specific typings obtained with molecular methods. The tables with HLA equivalents and the questionnaire for submission of serology on poorly identified alleles will also be available at the WMDA web page: www.bmdw.org/wmda.