Standard reference sequences for submission of HLA genotyping for the 18th International HLA and Immunogenetics Workshop

Repeated local ellipsoid protrusion supplements HLA surface characterization

Allorecognition of donor HLA is a major risk factor for long-term kidney graft survival. Although several molecular matching algorithms have been proposed that compare physiochemical and structural features of the donors' and recipients' HLA proteins in order to predict their compatibility, the exact underlying mechanisms are still not fully understood. We hypothesized that the ElliPro approach of single ellipsoid fitting and protrusion ranking lacks sensitivity for the characteristic shape of HLA molecules and developed a prediction pipeline named Snowball that is fitting smaller ellipsoids iteratively to substructures. Aggregated protrusion ranks of locally fitted ellipsoids were calculated for 712 publicly available HLA structures and 78 predicted structures using AlphaFold 2. Amino-acid sequence and protrusion ranks were used to train deep neural network predictors to infer protrusion ranks for all known HLA sequences. Snowball protrusion ranks appear to be more sensitive than ElliPro scores in fine parts of the HLA such as the helix structures forming the HLA binding groove in particular when the ellipsoids are fitted to substructures considering atoms within a 15 Å radius. A cloud-based web service was implemented based on amino-acid matching considering both protein- and position-specific surface area and protrusion ranks extending the previously presented Snowflake prediction pipeline.

The genotype list string code syntax for exchanging nomenclature-level genotyping results in clinical and research data management and analysis systems

The nomenclatures used to describe HLA and killer-cell immunoglobulin-like receptor (KIR) alleles distinguish unique nucleotide and peptide sequences, and patterns of expression, but are insufficient for describing genotyping results, as description of ambiguities and relations across loci require terminology beyond allele names. The genotype list (GL) String grammar describes genotyping results for genetic systems with defined nomenclatures, like HLA and KIR, documenting what is known and unknown about a given genotyping result. However, the accuracy of a GL String is dependent on the reference database version under which it was generated. Here, we describe the GL string code (GLSC) system, which associates each GL String with meta-data describing the specific reference context in which the GL String was created, and in which it should be interpreted. GLSC is a defined syntax for exchanging GL Strings in the context of a specific gene-family namespace, allele-name code-system, and pertinent reference database version. GLSC allows HLA and KIR genotyping data to be transmitted, parsed and interpreted in the appropriate context, in an unambiguous manner, on modern data-systems, including Health Level 7 Fast Healthcare Interoperability Resource systems. Technical specification for GLSC can be found at https://glstring.org.

Genotype List String 1.1: Extending the Genotype List String grammar for describing HLA and Killer-cell Immunoglobulin-like Receptor genotypes

The Genotype List (GL) String grammar for reporting HLA and Killer-cell Immunoglobulin-like Receptor (KIR) genotypes in a text string was described in 2013. Since this initial description, GL Strings have been used to describe HLA and KIR genotypes for more than 40 million subjects, allowing these data to be recorded, stored and transmitted in an easily parsed, text-based format. After a decade of working with HLA and KIR data in GL String format, with advances in HLA and KIR genotyping technologies that have fostered the generation of full-gene sequence data, the need for an extension of the GL String system has become clear. Here, we introduce the new GL String delimiter "?," which addresses the need to describe ambiguity in assigning a gene sequence to gene paralogs. GL Strings that do not include a "?" delimiter continue to be interpreted as originally described. This extension represents version 1.1 of the GL String grammar.

Snowflake: A deep learning-based human leukocyte antigen matching algorithm considering allele-specific surface accessibility

Histocompatibility in solid-organ transplantation has a strong impact on long-term graft survival. Although recent advances in matching of both B-cell epitopes and T-cell epitopes have improved understanding of allorecognition, the immunogenic determinants are still not fully understood. We hypothesized that HLA solvent accessibility is allele-specific, thus supporting refinement of HLA B-cell epitope prediction. We developed a computational pipeline named Snowflake to calculate solvent accessibility of HLA Class I proteins for deposited HLA crystal structures, supplemented by constructed HLA structures through the AlphaFold protein folding predictor and peptide binding predictions of the APE-Gen docking framework. This dataset trained a four-layer long short-term memory bidirectional recurrent neural network, which in turn inferred solvent accessibility of all known HLA Class I proteins. We extracted 676 HLA Class-I experimental structures from the Protein Data Bank and supplemented it by 37 Class-I alleles for which structures were predicted. For each of the predicted structures, 10 known binding peptides as reported by the Immune Epitope DataBase were rendered into the binding groove. Although HLA Class I proteins predominantly are folded similarly, we found higher variation in root mean square difference of solvent accessibility between experimental structures of different HLAs compared to structures with identical amino acid sequence, suggesting HLA’s solvent accessible surface is protein specific. Hence, residues may be surface-accessible on e.g. HLA-A*02:01, but not on HLA-A*01:01. Mapping these data to antibody-verified epitopes as defined by the HLA Epitope Registry reveals patterns of (1) consistently accessible residues, (2) only subsets of an epitope’s residues being consistently accessible and (3) varying surface accessibility of residues of epitopes. Our data suggest B-cell epitope definitions can be refined by considering allele-specific solvent-accessibility, rather than aggregating HLA protein surface maps by HLA class or locus. To support studies on epitope analyses in organ transplantation, the calculation of donor-allele-specific solvent-accessible amino acid mismatches was implemented as a cloud-based web service.

Challenges for the standardized reporting of NGS HLA genotyping: Surveying gaps between clinical and research laboratories

Next generation sequencing (NGS) is being applied for HLA typing in research and clinical settings. NGS HLA typing has made it feasible to sequence exons, introns and untranslated regions simultaneously, with significantly reduced labor and reagent cost per sample, rapid turnaround time, and improved HLA genotype accuracy. NGS technologies bring challenges for cost-effective computation, data processing and exchange of NGS-based HLA data. To address these challenges, guidelines and specifications such as Genotype List (GL) String, Minimum Information for Reporting Immunogenomic NGS Genotyping (MIRING), and Histoimmunogenetics Markup Language (HML) were proposed to streamline and standardize reporting of HLA genotypes. As part of the 17th International HLA and Immunogenetics Workshop (IHIW), we implemented standards and systems for HLA genotype reporting that included GL String, MIRING and HML, and found that misunderstanding or misinterpretations of these standards led to inconsistencies in the reporting of NGS HLA genotyping results. This may be due in part to a historical lack of centralized data reporting standards in the histocompatibility and immunogenetics community. We have worked with software and database developers, clinicians and scientists to address these issues in a collaborative fashion as part of the Data Standard Hackathons (DaSH) for NGS. Here we report several categories of challenges to the consistent exchange of NGS HLA genotyping data we have observed. We hope to address these challenges in future DaSH for NGS efforts.

Standard reference sequences for submission of HLA genotyping for the 18th International HLA and Immunogenetics Workshop

The International human leukocyte antigen (HLA) and Immunogenetics Workshops (IHIWs) have fostered international collaborations of researchers and experts in the fields of HLA, histocompatibility and immunology. These IHIW collaborations have comprised many projects focused on achieving a variety of specific goals. The international and collaborative nature of these projects necessitates the collection and analysis of complex data generated in multiple laboratories, often using multiple methods of acquisition. Collection and storage of these data in a consistent way adds value to IHIW projects, which can be extended to future work. DNA‐based genotyping data, especially HLA genotyping data, can be transmitted in the form of a Histoimmunogenetics Markup Language (HML) document. HML facilitates clear communication of a genotype and supporting metadata, such as, sequencing platform, laboratory assays, consensus sequence, and interpretation. Sequence information can be reported relative to known reference sequences, which add meaning and context to genotypes. Selecting the correct reference sequence for a given allele sequence is nuanced, and guidelines have emerged through collaborative community efforts such as Data Standards Hackathons. Here, we describe the guidelines established for the selection of reference sequences to be used in transmission of HLA (and MICA/MICB) genotyping data for the 18th IHIW.