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

The most frequent HLA alleles around the world: A fundamental synopsis

A comprehensive knowledge of human leukocyte antigen (HLA) molecular variation worldwide is essential in human population genetics research and disease association studies and is also indispensable for clinical applications such as allogeneic hematopoietic cell transplantation, where ensuring HLA compatibility between donors and recipients is paramount. Enormous progress has been made in this field thanks to several decades of HLA population studies allowing the development of helpful databases and bioinformatics tools. However, it is still difficult to appraise the global HLA population diversity in a synthetic way. We thus introduce here a novel approach, based on approximately 2000 data sets, to assess this complexity by providing a fundamental synopsis of the most frequent HLA alleles observed in different regions of the world. This new knowledge will be useful not only as a fundamental reference for basic research, but also as an efficient guide for clinicians working in the field of transplantation.

HLA-DR and HLA-DQ Polymorphism Correlation with Sexually Transmitted Infection Caused by Chlamydia trachomatis

Background and Objectives: Chlamydia trachomatis (C. trachomatis) represents one of the most prevalent bacterial sexually transmitted diseases. This study aims to explore the relationship between HLA alleles/genotypes/haplotypes and C. trachomatis infection to better understand high-risk individuals and potential complications. Materials and Methods: This prospective study recruited participants from Transylvania, Romania. Patients with positive NAAT tests for C. trachomatis from cervical/urethral secretion or urine were compared with controls regarding HLA-DR and -DQ alleles. DNA extraction for HLA typing was performed using venous blood samples. Results: Our analysis revealed that the presence of the DRB1*13 allele significantly heightened the likelihood of C. trachomatis infection (p = 0.017). Additionally, we observed that individuals carrying the DRB1*01/DRB1*13 and DQB1*03/DQB1*06 genotype had increased odds of C. trachomatis infection. Upon adjustment, the association between the DRB1*01/DRB1*13 genotype and C. trachomatis remained statistically significant. Conclusions: Our findings underscore the importance of specific HLA alleles and genotypes in influencing susceptibility to C. trachomatis infection. These results highlight the intricate relationship between host genetics and disease susceptibility, offering valuable insights for targeted prevention efforts and personalized healthcare strategies.

High-resolution African HLA resource uncovers HLA-DRB1 expression effects underlying vaccine response

How human genetic variation contributes to vaccine effectiveness in infants is unclear, and data are limited on these relationships in populations with African ancestries. We undertook genetic analyses of vaccine antibody responses in infants from Uganda (n = 1391), Burkina Faso (n = 353) and South Africa (n = 755), identifying associations between human leukocyte antigen (HLA) and antibody response for five of eight tested antigens spanning pertussis, diphtheria and hepatitis B vaccines. In addition, through HLA typing 1,702 individuals from 11 populations of African ancestry derived predominantly from the 1000 Genomes Project, we constructed an imputation resource, fine-mapping class II HLA-DR and DQ associations explaining up to 10% of antibody response variance in our infant cohorts. We observed differences in the genetic architecture of pertussis antibody response between the cohorts with African ancestries and an independent cohort with European ancestry, but found no in silico evidence of differences in HLA peptide binding affinity or breadth. Using immune cell expression quantitative trait loci datasets derived from African-ancestry samples from the 1000 Genomes Project, we found evidence of differential HLA-DRB1 expression correlating with inferred protection from pertussis following vaccination. This work suggests that HLA-DRB1 expression may play a role in vaccine response and should be considered alongside peptide selection to improve vaccine design.

Resolving unknown nucleotides in the IPD-IMGT/HLA database by extended and full-length sequencing of HLA class I and II alleles

In the past, identification of HLA alleles was limited to sequencing the region of the gene coding for the peptide binding groove, resulting in a lack of sequence information in the HLA database, challenging HLA allele assignment software programs. We investigated full-length sequences of 19 HLA class I and 7 HLA class II alleles, and we extended another 47 HLA class I alleles with sequences of 5' and 3' UTR regions that were all not yet available in the IPD-IMGT/HLA database. We resolved 8638 unknown nucleotides in the coding sequence of HLA class I and 2139 of HLA class II. Furthermore, with full-length sequencing of the 26 alleles, more than 90 kb of sequence information was added to the non-coding sequences, whereas extension of the 47 alleles resulted in the addition of 5.5 kb unknown nucleotides to the 5' UTR and > 31.7 kb to the 3' UTR region. With this information, some interesting features were observed, like possible recombination events and lineage evolutionary origins. The continuing increase in the availability of full-length sequences in the HLA database will enable the identification of the evolutionary origin and will help the community to improve the alignment and assignment accuracy of HLA alleles.

Recruiting refugees and migrants as potential hematopoietic stem cell donors to serve patients of comparable ethnicities with rare human leucocyte antigen patterns - The BluStar.NRW project in North Western Germany

Currently, approximately 19 million people with a migration background live in Germany. The majority of those descend from regions where the population has a genetically different distribution of HLA antigens when compared to the HLA frequencies usually found in North Western Europe. In case of severe haematological disorders of these individuals, allogeneic stem cell transplantation may be the treatment of choice. However, finding appropriate histocompatible hematopoietic stem cell donors continues to be a major challenge. If no matching sibling donors are available, there are only few suitable donors with a similar genetic background available in international blood stem cell donor registries. The "BluStar.NRW" project aimed to recruit new blood and hematopoietic stem cell donors with a migration background and to noticeably increase the number of suitable donors for patients within this group. Since December 2017, a total number of 9100 blood and stem cell donors with a migration background were recruited and typed for this project. HLA typing for HLA-A, -B, -C, -DRB1, -DQB1, and -DPB1 was performed by Next Generation Sequencing. We assessed the proportion of rare alleles according to HLA frequency tables, as defined by a frequency of <1:1000. The rare HLA allele frequencies according to HLA frequency tables of the BluStar.NRW cohort were compared with a matched control donor cohort: Rare HLA-A, -B, -C, -DRB1 and -DQB1 alleles occurred three times more frequent than in the control group, but rare HLA-DPB1 alleles occurred more frequently in the control cohort. This difference was highly significant for all HLA alleles (p < 0.0001 for HLA-A, -B, -C, -DRB1, -DPB1; p = 0.0002 for HLA-DQB1). In addition, the distribution of rare alleles differed between the two groups. To date, 29 work-ups were initiated, 12 PBSC, one BM and three DLI were collected so far out of the BluStar.NRW cohort. The apheresis probability is twofold higher (0.18% vs. 0.07%) compared to the control group which clearly shows a serious medical need. However, 13 work-ups were cancelled in the BluStar.NRW donor cohort which represents an almost twice as higher cancellation rate (45% vs. 25%). This single registry analysis with a large sample cohort clearly indicates that hematopoietic stem cell donors with a migration background represent an adequate donor pool to serve patients of comparable ethnicity.

Benchmarking NGS-Based HLA Typing Algorithms

Knowledge of the expected accuracy of HLA typing algorithms is important when choosing between algorithms and when evaluating the HLA typing predictions of an algorithm. This chapter guides the reader through an example benchmarking study that evaluates the performances of four NGS-based HLA typing algorithms as well as outlining factors to consider, when designing and running such a benchmarking study. The code related to this benchmarking workflow can be found at https://github.com/nikolasthuesen/springers-hla-benchmark/ .

Allele Frequency Net Database

The allele frequency net database (AFND, http://www.allelefrequencies.net ) is an online web-based repository that contains information on the frequencies of immune-related genes and their corresponding alleles in worldwide human populations. At present, the website contains data from 1784 population samples in more than 14 million individuals from 129 countries on the frequency of genes from different polymorphic regions including data for the human leukocyte antigen (HLA) system. In addition, over the last four years, AFND has also incorporated genotype raw data from 85,000 individuals comprising 215 population samples from 39 countries. Moreover, more population data sets containing next generation sequencing data spanning >3 million individuals have been added. This resource has been widely used in a variety of contexts such as histocompatibility, immunology, epidemiology, pharmacogenetics, epitope prediction algorithms for population coverage in vaccine development, population genetics, among many others. In this chapter, we present an update of the most used searching mechanisms as described in a previous volume and some of the latest developments included in AFND.

Impact of assigning 2-field HLA alleles from real-time PCR on deceased donor assessments and conformance with high resolution alleles

Waitlisted sensitised transplant recipients with HLA allele level antibodies to their own HLA antigen family are disadvantaged by current deficiencies in HLA typing for deceased donors. This is primarily because at time of organ allocation, HLA typing is provided at antigen level whereas solid phase assays provide allele level antibody definition. The gold standard for HLA allele typing is next generation sequencing (NGS), however time limitations with established NGS systems prevent NGS use for deceased donors. Instead, many labs use a real-time PCR (qPCR) antigen level result for deceased donors, which can disadvantage sensitised patients. Here, we compared assigning qPCR 2-field alleles to qPCR antigen level to determine the impact on virtual crossmatch (VXM) and discuss impact on donor-specific antibody (DSA) assignments. 244 consecutive deceased donors were HLA typed to allelic level by qPCR (LinkSeq SABR) and subsequently by NGS (One Lambda Alltype). The impact of qPCR allele assignments on potential DSA identification was investigated, by retrospectively investigating all 3904 VXMs, where recipient DSA assessments were assessed against donor HLA, was performed within the cohort. There was 96.3% concordance between qPCR and NGS for all allele level loci, with HLA-A; DQB1; and DPB1 having best agreement (99.4%, 98.4% and 99.4% respectively). Of the 3904 VXMs with qPCR allele assignment, there were 13 (<1%) occasions where the potential DSA assignment was impacted, with DQA1 having the most impact. Assigning alleles derived from qPCR to define unacceptable antigens for VXMs, can allow improved access to donor offers for sensitised patients by better defining alleles.

Conserved biophysical compatibility among the highly variable germline-encoded regions shapes TCR-MHC interactions

T cells are critically important components of the adaptive immune system primarily responsible for identifying and responding to pathogenic challenges. This recognition of pathogens is driven by the interaction between membrane-bound T cell receptors (TCRs) and antigenic peptides presented on major histocompatibility complex (MHC) molecules. The formation of the TCR-peptide-MHC complex (TCR-pMHC) involves interactions among germline-encoded and hypervariable amino acids. Germline-encoded and hypervariable regions can form contacts critical for complex formation, but only interactions between germline-encoded contacts are likely to be shared across many of all the possible productive TCR-pMHC complexes. Despite this, experimental investigation of these interactions have focused on only a small fraction of the possible interaction space. To address this, we analyzed every possible germline-encoded TCR-MHC contact in humans, thereby generating the first comprehensive characterization of these largely antigen-independent interactions. Our computational analysis suggests that germline-encoded TCR-MHC interactions that are conserved at the sequence level are rare due to the high amino acid diversity of the TCR CDR1 and CDR2 loops, and that such conservation is unlikely to dominate the dynamic protein-protein binding interface. Instead, we propose that binding properties such as the docking orientation are defined by regions of biophysical compatibility between these loops and the MHC surface.

Design of a novel multiepitope vaccine against Chlamydia pneumoniae using the extracellular protein as a target

Chlamydia pneumoniae (C. pneumoniae) infection in humans is universal and causes various respiratory infectious diseases, making a safe and effective preventive vaccine essential. In this study, a multi-epitope vaccine with CTLA-4 extracellular structure was constructed by an immunoinformatics approach. Since MOMP protein is the major extracellular protein in C. pneumoniae and has good immunogenicity and high conservation, we selected the MOMP protein of C. pneumoniae as the antigen target, predicted the T and B cell epitopes of the MOMP protein and then connected the CTLA-4 extracellular structure with the predicted dominant epitopes by various linkers to construct a multi-epitope vaccine. The biochemical characterization of the multi-epitope vaccine showed its immunogenicity and anti-allergic properties. The tertiary structure of this vaccine, along with molecular docking, molecular dynamics simulation, and principal component analysis, showed that the multi-epitope vaccine structure interacted with B7 (B7-1, B7-2) and toll-like receptors (TLR-2, TLR-4). Ultimately, the vaccine was cloned and effectively expressed in silico on an insect baculovirus expression vector (pFastBac1). These analyses showed that the designed vaccine could potentially target antigen-presenting cells and was immune to C. pneumoniae, which provided novel strategies for developing the vaccine.

Use of next-generation sequencing to detect polymorphism of 11 HLA allele loci in the Chinese Han population and variance from other common and well-documented lists

The focus of this study was to analyze polymorphisms in the HLA gene at 11 loci in 4845 Chinese Han populations using next-generation sequencing methods, and to compare common and well-documented (CWD) allelic differences between China and other CWD lists. A total of 44 DPB1 alleles, 13 DPA1 alleles, 20 DQA1 alleles and 19 DRB3/4/5 alleles were detected in this study. About 20%-50% of the CWD alleles in China differ from the American Society for Histocompatibility and Immunogenetics and European Federation for Immunogenetics (EFI) data. The revised list of HLA-CWD alleles in the Han population will provide additional data for the update of the IMGT/HLA database and contribute to a better understanding of hematopoietic stem cell transplantation and organ transplantation.

Distributions of 11-loci HLA alleles typed by amplicon-based next-generation sequencing in South Koreans

The range of HLA typing for successful hematopoietic stem cell transplantation (HSCT) is gradually expanding with the next-generation sequencing (NGS)-based improvement in its quality. However, it is influenced by the allocation of finances and laboratory conditions. HLA-A, -B, -C, -DRB1/3/4/5, -DQA1, -DQB1, -DPA1, and -DPB1 alleles were genotyped at the 3-field level by amplicon-based NGS using MiSeqDx system and compared to our previous study employing long-range PCR and NGS using TruSight HLA v2 kit, in healthy donors from South Korea. Exon 2, exons 2/3, exons 2/3/4 or 5 of 11-loci were amplified by multiplex PCR. The sequence reads of over 53 depth counts were consistently obtained in each sample exon, depending on the target exon determined to match the reference sequence contained in the IPD-IMGT/HLA Database. HLA alleles were investigated by combinations of the determined exons. A total of 18 alleles with a frequency over 10% were found at the 11 HLA loci. Three ambiguities of HLA-A, -C, and -DRB1 were resolved. We observed a total of 26 HLA-A ~ C ~ B and 6 HLA-DRB1 ~ DQA1 ~ DQB1 ~ DPA1 ~ DPB1 haplotypes having significant linkage disequilibrium between alleles at all neighboring HLA loci. This result was compatible with the previous one, using TruSight HLA v2 kit. Advantages are simple and short progress time because one plate is used for each PCR step in one PCR machine and 11-loci HLA typing is possible even if only eight samples. These data suggested that expanded 11-loci HLA typing data by amplicon-based NGS might help perform HSCT.

Autoantigenic properties of the aminoacyl tRNA synthetase family in idiopathic inflammatory myopathies

OBJECTIVES

Autoantibodies are thought to play a key role in the pathogenesis of idiopathic inflammatory myopathies (IIM). However, up to 40% of IIM patients, even those with clinical manifestations of anti-synthetase syndrome (ASSD), test seronegative to known myositis-specific autoantibodies. We hypothesized the existence of new potential autoantigens among human cytoplasmic aminoacyl tRNA synthetases (aaRS) in patients with IIM.

METHODS

Plasma samples from 217 patients with IIM according to 2017 EULAR/ACR criteria, including 50 patients with ASSD, 165 without, and two with unknown ASSD status were identified retrospectively, as well as age and gender-matched sera from 156 population controls, and 219 disease controls. Patients with previously documented ASSD had to test positive for at least one of the five most common anti-aaRS autoantibodies (anti-Jo1, -PL7, -PL12, -EJ, and -OJ) and present with one or more of the following clinical manifestations: interstitial lung disease, myositis, arthritis, Raynaud's phenomenon, fever, or mechanic's hands. Demographics, laboratory, and clinical data of the IIM cohort (ASSD and non-ASSD) were compared. Samples were screened using a multiplex bead array assay for presence of autoantibodies against a panel of 117 recombinant protein variants, representing 33 myositis-related proteins, including all nineteen cytoplasmic aaRS. Prospectively collected clinical data for the IIM cohort were retrieved and compared between groups within the IIM cohort and correlated with the results of the autoantibody screening. Principal component analysis was used to analyze clinical manifestations between ASSD, non-ASSD groups, and individuals with novel anti-aaRS autoantibodies.

RESULTS

We identified reactivity towards 16 aaRS in 72 of the 217 IIM patients. Twelve patients displayed reactivity against nine novel aaRS. The novel autoantibody specificities were detected in four previously seronegative patients for myositis-specific autoantibodies and eight with previously detected myositis-specific autoantibodies. IIM individuals with novel anti-aaRS autoantibodies (n = 12) all had signs of myositis, and they had either muscle weakness and/or muscle enzyme elevation, 2/12 had mechanic's hands, 3/12 had interstitial lung disease, and 2/12 had arthritis. The individuals with novel anti-aaRS and a pathological muscle biopsy all presented widespread up-regulation of major histocompatibility complex class I. The reactivities against novel aaRS could be confirmed in ELISA and western blot. Using the multiplex bead array assay, we could confirm previously known reactivities to four of the most common aaRS (Jo1, PL12, PL7, and EJ (n = 45)) and identified patients positive for anti-Zo, -KS, and -HA (n = 10) that were not previously tested. A low frequency of anti-aaRS autoantibodies was also detected in controls.

CONCLUSION

Our results suggest that most, if not all, cytoplasmic aaRS may become autoantigenic. Autoantibodies against new aaRS may be found in plasma of patients previously classified as seronegative with potential high clinical relevance.

Prognostic characteristics of immune subtypes associated with acute myeloid leukemia and their identification in cell subsets based on single-cell sequencing analysis

Immune genes play an important role in the development and progression of acute myeloid leukemia (AML). However, the role of immune genes in the prognosis and microenvironment of AML remains unclear. In this study, we analyzed 151 AML patients in the TCGA database for relevant immune cell infiltration. AML patients were divided into high and low immune cell infiltration clusters based on ssGSEA results. Immune-related pathways, AML pathways and glucose metabolism pathways were enriched in the high immune cell infiltration cluster. Then we screened the differential immune genes between the two immune cell infiltration clusters. Nine prognostic immune genes were finally identified in the train set by LASSO-Cox regression. We constructed a model in the train set based on the nine prognostic immune genes and validated the predictive capability in the test set. The areas under the ROC curve of the train set and the test set for ROC at 1, 3, 5 years were 0.807, 0.813, 0.815, and 0.731, 0.745, 0.830, respectively. The areas under ROC curve of external validation set in 1, 3, and 5 years were 0.564, 0.619, and 0.614, respectively. People with high risk scores accompanied by high TMB had been detected with the worst prognosis. Single-cell sequencing analysis revealed the expression of prognostic genes in AML cell subsets and pseudo-time analysis described the differentiation trajectory of cell subsets. In conclusion, our results reveal the characteristics of immune microenvironment and cell subsets of AML, while it still needs to be confirmed in larger samples studies. The prognosis model constructed with nine key immune genes can provide a new method to assess the prognosis of AML patients.

Genetic polymorphism of HLA-DRA and alcohol consumption affect hepatitis development in the Korean population

BACKGROUND

Hepatitis is an inflammation of the liver that has several potential causes; however, the genetic association has recently begun to be studied.

OBJECTIVES

Human leukocyte antigen (HLA) is an essential component of the immune response, and in this study, we conducted a correlation analysis to determine whether genetic polymorphisms of HLA and drinking habits affect hepatitis development.

METHODS

Genetic polymorphisms of HLA were investigated using Korean genomic and epidemiological data. A gene association study was performed using PLINK version 1.07. Other statistical analyses and multivariate logistic regression analyses were performed using PASW Statistics version 18.0.

RESULTS

Thirteen single nucleotide polymorphisms (SNPs) in HLA-DRA showed significant statistical correlations with hepatitis. In particular, rs9268645 showed the highest statistical association with hepatitis (P = 3.97 × 10^-5, odds ratio [OR] = 0.72, 95% confidence interval [CI] = 0.61-0.84). In multivariate logistic regression analysis, when considering only genetic factors, the A allele of rs9268644 showed a reduced hepatitis OR of approximately 0.52-fold. However, the group carrying the minor A allele (AA + AC) with alcohol consumption had an approximately 1.58-fold OR of hepatitis compared to that of the group carrying the same allele with no alcohol consumption. This implies that the A allele of rs9268644 has a protective effect on hepatitis by genetic factors and shows sensitivity to alcohol.

CONCLUSIONS

Our results showed that hepatitis is influenced by both genetic and external factors (drinking habits), which can provide new guidelines for the prevention or treatment of hepatitis.

Polymorphic differences within HLA-C alleles contribute to alternatively spliced transcripts lacking exon 5

The HLA genes are amongst the most polymorphic in the human genome. Alternative splicing could add an extra layer of complexity, but has not been studied extensively. Here, we applied an RNA based approach to study the influence of allele polymorphism on alternative splicing of HLA‐C in peripheral blood. RNA was isolated from these peripheral cells, converted into cDNA and amplified specifically for 12 common HLA‐C allele groups. Through subsequent sequencing of HLA‐C, we observed alternative splicing variants of HLA‐C*04 and *16 that resulted in exon 5 skipping and were co‐expressed with the mature transcript. Investigation of intron 4 sequences of HLA‐C*04 and *16 compared with other HLA‐C alleles demonstrated no effect on predicted splice sites and branch point. To further investigate if the unique polymorphic positions in exon 5 of HLA‐C*04 or *16 may facilitate alternative splicing by acting on splicing regulatory elements (SRE), in‐silico splicing analysis was performed. While the HLA‐C*04 specific SNP in exon 5 had no effect on predicted exonic SRE, the HLA‐C*16 specific exon 5 SNP did alter exonic SRE. Our findings provide experimental and theoretical support for the concept that polymorphisms within the HLA‐C alleles influence the alternative splicing of HLA‐C.

A new strategy for systematically classifying HLA alleles into serological specificities

HLA serological specificities were defined by the reactivity of HLA molecules with sets of sera and monoclonal antibodies. Many recently identified alleles defined by molecular typing lack their serotype assignment. We surveyed the literature describing the correlation of the reactivity of serologic reagents with AA residues. 20 - 25 AA residues determining epitopes (DEP) that correlated with 82 WHO serologic specificities were identified for HLA class I loci. Thirteen DEP each located in the beta-1 domains that correlated with 24 WHO serologic specificities were identified for HLA-DRB1 and -DQB1 loci. The designation of possible HLA-DPB1, -DQA1, -DPA1, and additional serological specificities that result from epitopes defined by residues located at both -DQA1 and -DQB1 subunits were also examined. HATS software was developed for automated serotype assignments to HLA alleles in one of the three hierarchical matching criteria: 1) all DEP (FULL); 2) selected DEP specific to each serological specificities (SEROTYPE); 3) one AA mismatch with one or more SEROTYPES (INCOMPLETE). Results were validated by evaluating the alleles whose serotypes do not correspond to the first field of the allele name listed in the HLA dictionary. Additional 85 and 21 DEP patterns that do not correspond to any WHO serologic specificities for common HLA class I and DRB1 alleles were identified, respectively. A comprehensive antibody identification panel would allow for accurate unacceptable antigen listing and compatibility predictions in solid organ transplantations. We propose that antibody-screening panels should include all serologic specificities identified in this study. This article is protected by copyright. All rights reserved.

The clinical utility and thresholds of Virtual and Halifaster Flow crossmatches in lung transplantation

INTRODUCTION

Immune sensitization, defined as the presence of alloreactive donor-specific antibodies (DSA), is associated with increased wait-times and inferior transplant outcomes. Identifying pre-transplant DSA with a physical cell-based assay is critical in defining immunological risk. However, improved solid phase antibody detection has provided the potential to forgo this physical assay. Here, we evaluated the association between DSA mean fluorescence intensity (MFI) and the recently introduced Halifaster Flow cytometry crossmatch (FXM) to determine if MFI could predict the outcome of FXM and whether a virtual crossmatch (VXM) would provide an accurate risk assessment.

METHODS

Sera from 134 waitlisted lung patients was retrospectively assessed by Halifaster FXM against lymphocytes preparations from 32 donors, resulting in 265 FXMs. HLA typing was performed to 2-field allelic level and Luminex single antigen beads (SAB) used to identify DSA. The association between FXM and Luminex MFI was calculated using ROC analysis. MFI threshold accuracy was confirmed using a separate validation cohort (174 recipient sera and 34 donors), whereby both virtual crossmatch (VXM) and FXMs were compared.

RESULTS

From the 265 FXM performed, 48 (18%) T-cell (TFXM) and 56 (21%) B-cell (BFXM) were positive. In the evaluation cohort, MFI thresholds of 2000 for HLA-A, B, DRB1 and >4000 for DQB1, were predictive of a positive FXM. The validation cohort of 233 paired FXM and VXM confirmed these MFI thresholds for both TFXM and BFXM with an accuracy of 91.4% and 89.3% respectively.

CONCLUSION

A positive VXM, defined with HLA-specific MFI thresholds predicts Halifaster FXM reactivity, and can potentially expedite organ allocation, by minimising the need for the more time-consuming flow cytometry crossmatch. This article is protected by copyright. All rights reserved.

HLA3D: an integrated structure-based computational toolkit for immunotherapy

Motivation

The human major histocompatibility complex (MHC), also known as human leukocyte antigen (HLA), plays an important role in the adaptive immune system by presenting non-self-peptides to T cell receptors. The MHC region has been shown to be associated with a variety of diseases, including autoimmune diseases, organ transplantation and tumours. However, structural analytic tools of HLA are still sparse compared to the number of identified HLA alleles, which hinders the disclosure of its pathogenic mechanism.

Result

To provide an integrative analysis of HLA, we first collected 1296 amino acid sequences, 256 protein data bank structures, 120 000 frequency data of HLA alleles in different populations, 73 000 publications and 39 000 disease-associated single nucleotide polymorphism sites, as well as 212 modelled HLA heterodimer structures. Then, we put forward two new strategies for building up a toolkit for transplantation and tumour immunotherapy, designing risk alignment pipeline and antigenic peptide prediction pipeline by integrating different resources and bioinformatic tools. By integrating 100 000 calculated HLA conformation difference and online tools, risk alignment pipeline provides users with the functions of structural alignment, sequence alignment, residue visualization and risk report generation of mismatched HLA molecules. For tumour antigen prediction, we first predicted 370 000 immunogenic peptides based on the affinity between peptides and MHC to generate the neoantigen catalogue for 11 common tumours. We then designed an antigenic peptide prediction pipeline to provide the functions of mutation prediction, peptide prediction, immunogenicity assessment and docking simulation. We also present a case study of hepatitis B virus mutations associated with liver cancer that demonstrates the high legitimacy of our antigenic peptide prediction process. HLA3D, including different HLA analytic tools and the prediction pipelines, is available at http://www.hla3d.cn/.

Principles of Virtual Crossmatch Testing for Kidney Transplantation

Human leukocyte antigens (HLAs) are the primary determinants of alloimmunity. A crossmatch test is a test that determines the immunologic risk of a recipient with a potential donor by ensuring that there are no transplant-relevant circulating antibodies in the recipient directed against donor antigens. Physical crossmatch (PXM) tests, such as complement-dependent cytotoxicity crossmatch (CDCXM) and flow cytometry crossmatch (FCXM), require mixing of patient serum and donor cells, are labor intensive, and are logistically challenging. Virtual crossmatch (VXM) test assesses immunologic compatibility between recipient and potential donor by analyzing the results of 2 independently done physical laboratory tests—patient anti-HLA antibody and donor HLA typing. The goal of VXM is pretransplant risk stratification—though there is no consensus on whether such risk assessment involves predicting the PXM result or the posttransplant outcome. Although the concept of VXM is not new, the advent of solid-phase assays for detecting circulating antibodies in the recipient directed against individual HLA and DNA-based methods for typing donor HLA specificities at a higher resolution makes the routine use of VXM a reality. Accordingly, VXM may be applied at different scenarios—both for sensitized and nonsensitized patients. Implementation of VXM-based approach has resulted in statistically significant reduction in cold ischemia time without an increase in hyperacute rejection episodes. Though there are considerable challenges, VXM is expected to be used more often in the future, depending on the transplant center’s tolerance of immunologic risk.

HLA-DRB3/4/5 Matching Improves Outcome of Unrelated Hematopoietic Stem Cell Transplantation

The HLA-DRB3/4/5 loci are closely linked to the HLA-DRB1 gene. Mismatches in these loci occur with a frequency of about 8%–12% in otherwise 10/10 HLA-matched transplant pairs. There is preliminary evidence that these disparities may associate with increased acute graft-versus-host disease (GvHD) rates. The aim of this study was to analyze a large cohort of German patients and their donors for HLA-DRB3/4/5 compatibility and to correlate the HLA-DRB3/4/5 matching status with the outcome of unrelated hematopoietic stem cell transplantation (uHSCT). To this end, 3,410 patients and their respective donors were HLA-DRB3/4/5 and HLA-DPB1 typed by amplicon-based next-generation sequencing (NGS). All patients included received their first allogeneic transplant for malignant hematologic diseases between 2000 and 2014. Mismatches in the antigen recognition domain (ARD) of HLA-DRB3/4/5 genes were correlated with clinical outcome. HLA-DRB3/4/5 incompatibility was seen in 12.5% (n = 296) and 17.8% (n = 185) of the 10/10 and 9/10 HLA-matched cases, respectively. HLA-DRB3/4/5 mismatches in the ARD associated with a worse overall survival (OS), as shown in univariate (5-year OS: 46.1% vs. 39.8%, log-rank p = 0.038) and multivariate analyses [hazard ratio (HR) 1.25, 95% CI 1.02–1.54, p = 0.034] in the otherwise 10/10 HLA-matched subgroup. The worse outcome was mainly driven by a significantly higher non-relapse mortality (HR 1.35, 95% CI 1.05–1.73, p = 0.017). In the 9/10 HLA-matched cases, the effect was not statistically significant. Our study results suggest that mismatches within the ARD of HLA-DRB3/4/5 genes significantly impact the outcome of otherwise fully matched uHSCT and support their consideration upon donor selection in the future.

Self-mediated positive selection of T cells sets an obstacle to the recognition of nonself

Adaptive immune recognition is mediated by the binding of peptide-human leukocyte antigen complexes by T cells. Positive selection of T cells in the thymus is a fundamental step in the generation of a responding T cell repertoire: only those T cells survive that recognize human peptides presented on the surface of cortical thymic epithelial cells. We propose that while this step is essential for optimal immune function, the process results in a defective T cell repertoire because it is mediated by self-peptides. To test our hypothesis, we focused on amino acid motifs of peptides in contact with T cell receptors. We found that motifs rarely or not found in the human proteome are unlikely to be recognized by the immune system just like the ones that are not expressed in cortical thymic epithelial cells or not presented on their surface. Peptides carrying such motifs were especially dissimilar to human proteins. Importantly, we present our main findings on two independent T cell activation datasets and directly demonstrate the absence of naïve T cells in the repertoire of healthy individuals. We also show that T cell cross-reactivity is unable to compensate for the absence of positively selected T cells. Additionally, we show that the proposed mechanism could influence the risk for different infectious diseases. In sum, our results suggest a side effect of T cell positive selection, which could explain the nonresponsiveness to many nonself peptides and could improve the understanding of adaptive immune recognition.

Like Wings of a Bird: Functional Divergence and Complementarity between HLA-A and HLA-B Molecules

Human leukocyte antigen (HLA) genes are among the most polymorphic of our genome, as a likely consequence of balancing selection related to their central role in adaptive immunity. HLA-A and HLA-B genes were recently suggested to evolve through a model of joint divergent asymmetric selection conferring all human populations, including those with severe loss of diversity, an equivalent immune potential. However, the mechanisms by which these two genes might undergo joint evolution while displaying very distinct allelic profiles in populations are still unknown. To address this issue, we carried out extensive data analyses (among which factorial correspondence analysis and linear modeling) on 2,909 common and rare HLA-A, HLA-B, and HLA-C alleles and 200,000 simulated pathogenic peptides by taking into account sequence variation, predicted peptide-binding affinity and HLA allele frequencies in 123 populations worldwide. Our results show that HLA-A and HLA-B (but not HLA-C) molecules maintain considerable functional divergence in almost all populations, which likely plays an instrumental role in their immune defense. We also provide robust evidence of functional complementarity between HLA-A and HLA-B molecules, which display asymmetric relationships in terms of amino acid diversity at both inter- and intraprotein levels and in terms of promiscuous or fastidious peptide-binding specificities. Like two wings of a flying bird, the functional complementarity of HLA-A and HLA-B is a perfect example, in our genome, of duplicated genes sharing their capacity of assuming common vital functions while being submitted to complex and sometimes distinct environmental pressures.

MHCVision: estimation of global and local false discovery rate for MHC class I peptide binding prediction

Motivation

MHC-peptide binding prediction has been widely used for understanding the immune response of individuals or populations, each carrying different MHC molecules as well as for the development of immunotherapeutics. The results from MHC-peptide binding prediction tools are mostly reported as a predicted binding affinity (IC₅₀) and the percentile rank score, and global thresholds e.g. IC₅₀ value < 500 nM or percentile rank < 2% are generally recommended for distinguishing binding peptides from non-binding peptides. However, it is difficult to evaluate statistically the probability of an individual peptide binding prediction to be true or false solely considering predicted scores. Therefore, statistics describing the overall global false discovery rate (FDR) and local FDR, also called posterior error probability (PEP) are required to give statistical context to the natively produced scores.

Result

We have developed an algorithm and code implementation, called MHCVision, for estimation of FDR and PEP values for the predicted results of MHC-peptide binding prediction from the NetMHCpan tool. MHCVision performs parameter estimation using a modified expectation maximization framework for a two-component beta mixture model, representing the distribution of true and false scores of the predicted dataset. We can then estimate the PEP of an individual peptide’s predicted score, and conversely the probability that it is true. We demonstrate that the use of global FDR and PEP estimation can provide a better trade-off between sensitivity and precision over using currently recommended thresholds from tools.

Availability and implementation

https://github.com/PGB-LIV/MHCVision.

Supplementary information

Supplementary data are available at Bioinformatics online.

T-Cell Receptor Mimic Antibodies for Cancer Immunotherapy

Antibody-based immunotherapies show clinical effectiveness in various cancer types. However, the target repertoire is limited to surface or soluble antigens, which are a relatively small percentage of the cancer proteome. Most proteins of the human proteome are intracellular. Short peptides from intracellular targets can be presented by MHC class I (MHC-I) molecules on cell surface, making them potential targets for cancer immunotherapy. Antibodies can be developed to target these peptide/MHC complexes, similar to the recognition of such complexes by the T-cell receptor (TCR). These antibodies are referred to as T-cell receptor mimic (TCRm) or TCR-like antibodies. Ongoing preclinical and clinical studies will help us understand their mechanisms of action and selection of target epitopes for immunotherapy. The present review will summarize and discuss the selection of intracellular antigens, production of the peptide/MHC complexes, isolation of TCRm antibodies for therapeutic applications, limitations of TCRm antibodies, and possible ways to advance TCRm antibody-based approaches into the clinic.

Allele and haplotype frequencies of human leukocyte antigen-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, -DPA1, and -DPB1 by next generation sequencing-based typing in Koreans in South Korea

Allele frequencies and haplotype frequencies of HLA-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, -DPA1, and -DPB1 have been rarely reported in South Koreans using unambiguous, phase-resolved next generation DNA sequencing. In this study, HLA typing of 11 loci in 173 healthy South Koreans were performed using next generation DNA sequencing with long-range PCR, TruSight^® HLA v2 kit, Illumina MiSeqDx platform system, and Assign^™ for TruSight^™ HLA software. Haplotype frequencies were calculated using the PyPop software. Direct counting methods were used to investigate the association with DRB1 for samples with only one copy of a particular secondary DRB locus. We compared these allele types with the ambiguous allele combinations of the IPD-IMGT/HLA database. We identified 20, 40, 26, 31, 19, 16, 4, and 16 alleles of HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQA1, HLA-DQB1, HLA-DPA1, and HLA-DPB1, respectively. The number of HLA-DRB3/4/5 alleles was 4, 5, and 3, respectively. The haplotype frequencies of most common haplotypes were as follows: A*33:03:01-B*44:03:01-C*14:03-DRB1*13:02:01-DQB1*06:04:01-DPB1*04:01:01 (2.89%), A*33:03:01-B*44:03:01-C*14:03 (4.91%), DRB1*08:03:02-DQA1*01:03:01-DQB1*06:01:01-DPA1*02:02:02-DPB1*05:01:01 (5.41%), DRB1*04:05:01-DRB4*01:03:01 (12.72%), DQA1*01:03:01-DQB1*06:01:01 (13.01%), and DPA1*02:02:02-DPB1*05:01:01 (30.83%). In samples with only one copy of a specific secondary DRB locus, we examined its association with DRB1. We, thus, resolved 10 allele ambiguities in HLA-B, -C (each exon 2+3), -DRB1, -DQB1, -DQA1, and -DPB1 (each exon 2) of the IPD-IMGT/HLA database. Korean population was geographically close to Japanese and Han Chinese populations in the genetic distances by multidimensional scaling (MDS) plots. The information obtained by HLA typing of the 11 extended loci by next generation sequencing may be useful for more exact diagnostic tests on various transplantations and the genetic population relationship studies in South Koreans.

Cataloguing experimentally confirmed 80.7 kb-long ACKR1 haplotypes from the 1000 Genomes Project database

Background

Clinically effective and safe genotyping relies on correct reference sequences, often represented by haplotypes. The 1000 Genomes Project recorded individual genotypes across 26 different populations and, using computerized genotype phasing, reported haplotype data. In contrast, we identified long reference sequences by analyzing the homozygous genomic regions in this online database, a concept that has rarely been reported since next generation sequencing data became available.

Study design and methods

Phased genotype data for a 80.6 kb region of chromosome 1 was downloaded for all 2,504 unrelated individuals of the 1000 Genome Project Phase 3 cohort. The data was centered on the ACKR1 gene and bordered by the CADM3 and FCER1A genes. Individuals with heterozygosity at a single site or with complete homozygosity allowed unambiguous assignment of an ACKR1 haplotype. A computer algorithm was developed for extracting these haplotypes from the 1000 Genome Project in an automated fashion. A manual analysis validated the data extracted by the algorithm.

Results

We confirmed 902 ACKR1 haplotypes of varying lengths, the longest at 80,584 nucleotides and shortest at 1,901 nucleotides. The combined length of haplotype sequences comprised 19,895,388 nucleotides with a median of 16,014 nucleotides. Based on our approach, all haplotypes can be considered experimentally confirmed and not affected by the known errors of computerized genotype phasing.

Conclusions

Tracts of homozygosity can provide definitive reference sequences for any gene. They are particularly useful when observed in unrelated individuals of large scale sequence databases. As a proof of principle, we explored the 1000 Genomes Project database for ACKR1 gene data and mined long haplotypes. These haplotypes are useful for high throughput analysis with next generation sequencing. Our approach is scalable, using automated bioinformatics tools, and can be applied to any gene.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04169-6.

High-throughput sequencing defines donor and recipient HLA B-cell epitope frequencies for prospective matching in transplantation

Compatibility for human leukocyte antigen (HLA) genes between transplant donors and recipients improves graft survival but prospective matching is rarely performed due to the vast heterogeneity of this gene complex. To reduce complexity, we have combined next-generation sequencing and in silico mapping to determine transplant population frequencies and matching probabilities of 150 antibody-binding eplets across all 11 classical HLA genes in 2000 ethnically heterogeneous renal patients and donors. We show that eplets are more common and uniformly distributed between donors and recipients than the respective HLA isoforms. Simulations of targeted eplet matching shows that a high degree of overall compatibility, and perfect identity at the clinically important HLA class II loci, can be obtained within a patient waiting list of approximately 250 subjects. Internal epitope-based allocation is thus feasible for most major renal transplant programs, while regional or national sharing may be required for other solid organs.

Tran et al. combine high throughput sequencing, structural biology and computational simulation to determine the HLA allele and antibody-defined epitope frequencies in renal transplant patients and donors. These results demonstrate the feasibility of HLA epitope matching using data from a national transplantation program.

HLA-DPB1 allele frequencies in the West Midlands region of the United Kingdom: A critical evaluation against the common, intermediate and well-documented allele catalogues CWD 2.0.0, EFI CWD and CIWD 3.0.0

The Birmingham H&I laboratory performed HLA typing on 456 potential deceased solid organ donors in the UK between 2014 and 2016. Accurate DPB1 typing is essential for determining HLA compatibility in transplantation, thus we report HLA-DPB1 for potential deceased solid organ donors. To correctly interpret HLA typing data, laboratories must understand both international and local HLA allele frequencies. In this analysis, we determined HLA-DPB1 allele and genotype frequencies for these 456 donors. HLA-DPB1 diversity was evaluated against the common and well-documented (CWD) alleles 2.0.0 catalogue, the European Federation for Immunogenetics (EFI) CWD catalogue and the common, intermediate and well-documented (CIWD) 3.0.0 catalogue. Additionally, we determined which alleles are common in our local deceased donor population. We observed 27 HLA-DPB1 alleles with DPB1*04:01 being the most frequently observed (allele frequency = 0.4342). All alleles detected locally were present in CIWD 3.0.0, however, DPB1*124:01 and *135:01 were not present in CWD 2.0.0 and DPB1*104:01 and *135:01 were not present in EFI CWD. Twenty of 27 DPB1 alleles identified were defined as locally common and also listed as common in CIWD 3.0.0 representing 62.5% of common alleles in the subset of CIWD 3.0.0 from individuals of a European geographic, ancestral or ethnic background. The alleles HLA-DPB1*16:01 and *20:01 are locally common but not listed as common in EFI CWD and DPB1*104:01 is not listed as common in CWD 2.0.0 catalogue. Our analysis showed that the detected alleles and locally common alleles within our population were aligned with the CIWD 3.0.0 catalogue.

HLA Variants and Inhibitor Development in Hemophilia A: A Retrospective Case-Controlled Study Using the ATHNdataset

In hemophilia A (HA) patients, F8 gene-defects as genetic risk-factors for developing inhibitors to Factor VIII have been extensively studied. Here we provide estimates of inhibitor-risk associated with the patient's Human Leukocyte Antigen (HLA). We used next generation sequencing for high-resolution HLA Class II typing of 997 HA patients. Using inhibitor prevalence reports from the My Life Our Future (MLOF) research repository, we calculated Odds Ratios (OR) for inhibitor development in a multivariate model considering HLA-DRB1/3/4/5, HLA-DPB1, HLA-DQB1, race, F8 pathogenic variant type, and age. Participants with 1 HLA variant (DPB1*02:02) had developed inhibitors at a higher rate while participants with 2 HLA variants (DRB1*04:07; DRB1*11:04) had developed inhibitors at a lower rate. Additionally, patients with missense variants had developed inhibitors at a lower rate and participants with large structural changes (>50 bp) had developed inhibitors at a higher rate (both compared to Intron 22 inversion). Using a cohort of participants with a distribution of HLA-DRB1 alleles comparable to that in the North American population we show that the HLA repertoire of a HA patient can be a risk-factor for inhibitor development.

High Frequency of Shared Clonotypes in Human T Cell Receptor Repertoires

The collection of T cell receptors (TCRs) generated by somatic recombination is large but unknown. We generate large TCR repertoire datasets as a resource to facilitate detailed studies of the role of TCR clonotypes and repertoires in health and disease. We estimate the size of individual human recombined and expressed TCRs by sequence analysis and determine the extent of sharing between individual repertoires. Our experiments reveal that each blood sample contains between 5 million and 21 million TCR clonotypes. Three individuals share 8% of TCRβ- or 11% of TCRα-chain clonotypes. Sorting by T cell phenotypes in four individuals shows that 5% of naive CD4+ and 3.5% of naive CD8+ subsets share their TCRβ clonotypes, whereas memory CD4+ and CD8+ subsets share 2.3% and 0.4% of their clonotypes, respectively. We identify the sequences of these shared TCR clonotypes that are of interest for studies of human T cell biology.

Soto et al. examine the extent to which five healthy adults share their T cell receptor (TCR) repertoire. Using sequencing and bioinformatics, they show a high prevalence of shared clonotypes even considering different T cell phenotypes. Possible functions for some clonotypes are inferred based on homology with TCRs in GenBank.

Coexistence of inhibitory and activating killer-cell immunoglobulin-like receptors to the same cognate HLA-C2 and Bw4 ligands confer breast cancer risk

Human leukocyte antigen (HLA) class I-specific killer-cell immunoglobulin-like receptors (KIR) regulate natural killer (NK) cell function in eliminating malignancy. Breast cancer (BC) patients exhibit reduced NK-cytotoxicity in peripheral blood. To test the hypothesis that certain KIR-HLA combinations impairing NK-cytotoxicity predispose to BC risk, we analyzed KIR and HLA polymorphisms in 162 women with BC and 278 controls. KIR-Bx genotypes increased significantly in BC than controls (83.3% vs. 71.9%, OR 1.95), and the increase was more pronounced in advanced-cancer (OR 5.3). No difference was observed with inhibitory KIR (iKIR) and HLA-ligand combinations. The activating KIR (aKIR) and HLA-ligand combinations, 2DS1 + C2 (OR 2.98) and 3DS1 + Bw4 (OR 2.6), were significantly increased in advanced-BC. All patients with advanced-cancer carrying 2DS1 + C2 or 3DS1 + Bw4 also have their iKIR counterparts 2DL1 and 3DL1, respectively. Contrarily, the 2DL1 + C2 and 3DL1 + Bw4 pairs without their aKIR counterparts are significantly higher in controls. These data suggest that NK cells expressing iKIR to the cognate HLA-ligands in the absence of putative aKIR counterpart are instrumental in antitumor response. These data provide a new framework for improving the utility of genetic risk scores for individualized surveillance.

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.

BSHI guideline: HLA matching and donor selection for haematopoietic progenitor cell transplantation

A review of the British Society for Histocompatibility and Immunogenetics (BSHI) Guideline 'HLA matching and donor selection for haematopoietic progenitor cell transplantation' published in 2016 was undertaken by a BSHI appointed writing committee. Literature searches were performed and the data extracted were presented as recommendations according to the GRADE nomenclature.

New challenges, new opportunities: Next generation sequencing and its place in the advancement of HLA typing

The Human Leukocyte Antigen (HLA) system has a critical role in immunorecognition, transplantation, and disease association. Early typing techniques provided the foundation for genotyping methods that revealed HLA as one of the most complex, polymorphic regions of the human genome. Next Generation Sequencing (NGS), the latest molecular technology introduced in clinical tissue typing laboratories, has demonstrated advantages over other established methods. NGS offers high-resolution sequencing of entire genes in time frames and price points considered unthinkable just a few years ago, contributing a wealth of data informing histocompatibility assessment and standards of clinical care. Although the NGS platforms share a high-throughput massively parallel processing model, differing chemistries provide specific strengths and weaknesses. Research-oriented Third Generation Sequencing and related advances in bioengineering continue to broaden the future of NGS in clinical settings. These diverse applications have demanded equally innovative strategies for data management and computational bioinformatics to support and analyze the unprecedented volume and complexity of data generated by NGS. We discuss some of the challenges and opportunities associated with NGS technologies, providing a comprehensive picture of the historical developments that paved the way for the NGS revolution, its current state and future possibilities for HLA typing.

Allogeneic hematopoietic stem cell transplant recipients in Spain: Human leukocyte antigen characteristics and diversity by high-resolution analysis

There are many studies on the polymorphism of the HLA system in healthy donor populations, such as registries of unrelated bone marrow donors. Investigations on the characterization of the HLA complex in hematopoietic stem cell transplant (HSCT) patients, however, are scarce, at least in the Spanish population. This study presents a large-scale analysis of allelic diversity and HLA distribution at a high-resolution level in 2886 patients undergoing HSCT in Spanish centres of the "Grupo Español de Trasplante Hematopoyético y Terapia Celular" during a period of 11 years. Allelic diversity analysis identified 67 HLA-A, 133 HLA-B, 60 HLA-C, 63 HLA-DRB1, 24 HLA-DQB1 and 27 HLA-DPB1 different alleles. Rare alleles were detected among which 33 alleles had not been reported in the European catalog of common and well-documented HLA alleles. Regarding the distribution of five genes-haplotypes, it was observed that the five most frequent extended haplotypes found in our population were between the most common in other Spanish populations, both in patients and in healthy subjects. However, some particular haplotypes were also detected. Bilocus associations HLA-C ~ B and -DRB1 ~ DQB1 were analyzed in order to predict the probability of finding 10/10 matched donors in registries. We found HLA-B alleles showing a great diversity of combinations with HLA-C alleles and unusual associations involving a negative predicting factor. In the field of adoptive therapies, our work supports the necessity to expand further research of TCR-engineered cells, adoptive transfer of virus-specific T-cells and vaccines to target HLA alleles other than A*02:01. HLA alleles such as A*01:01, A*03:01, A*24:02, B*44:03, B*07:02 or B*51:01, might be considered new targets due to its high frequency in our population.

Distributions of HLA-A, -B, and -DRB1 alleles typed by amplicon-based next generation sequencing in Korean volunteer donors for unrelated hematopoietic stem cell transplantation

HLA genes play a pivotal role for successful hematopoietic stem cell transplantation (HSCT). There is an increasing need for sophisticated screening of donor HLA genotypes for unrelated HSCT. Next generation sequencing (NGS) has emerged as an alternative for classical Sanger sequence for HLA typing. In this study, HLA-A, -B, and -DRB1 alleles were genotyped at the allelic (6-digit) level using MiSeqDx in 26,202 volunteers from the Korean Network for Organ Sharing. Exon 2 and 3 of HLA-A and -B and exon 2 of HLA-DRB1 were amplified by polymerase chain reaction (PCR) and each allele was determined by matching the targeted exons and the reference sequence consisting of the IPD-IMGT/HLA Database. Seventy alleles of HLA-A, 102 alleles of HLA-B, and 69 alleles of HLA-DRB1 were identified. According to common and well-documented catalogs, 34 alleles in HLA-A, 61 in HLA-B, and 45 in HLA-DRB1 locus were common alleles, and 12, 14, and 11 kinds, were well-documented alleles, respectively. Thirteen novel alleles including 3 alleles in HLA-A, 8 alleles in HLA-B, and 2 alleles in HLA-DRB1 loci were found. Ten haplotypes with a frequency of more than 1.0% accounted for 22.4% of the total haplotype frequencies. Cis/trans ambiguities of HLA-A and -B loci by combination of exons 2 and 3 were analyzed to be 0.17% of 3 and 3.95% of 22 genotypes, respectively. This information on rare and novel alleles found by accurate HLA typing with NGS may be helpful for unrelated HSCT among Koreans.

A snapshot of human leukocyte antigen (HLA) diversity using data from the Allele Frequency Net Database

The extensive allelic variability observed in several genes related to the immune response and its significance in different areas including transplantation, disease association studies, diversity in human populations, among many others, has led the scientific community to analyse these variants among individuals. Serving as an electronic data warehouse, the Allele Frequency Net Database (AFND, http://www.allelefrequencies.net) contains data on the frequency of immune related genes and their corresponding alleles from more than 1700 worldwide population samples covering more than ten million unrelated individuals. The collection of population data sets available in AFND encompasses different polymorphic regions including the highly-polymorphic human leukocyte antigen (HLA) system for which more than 1200 populations are available. In this article, we provide an insight of the high diversity found in the HLA region by examining population data sets stored in AFND, as well as a description of the available data sets for further analyses.

Correlation of the two most frequent HLA haplotypes in the Italian population to the differential regional incidence of Covid-19

BACKGROUND

Understanding how HLA polymorphisms may affect both susceptibility, course and severity of Covid-19 infection could help both at the clinical level to identify individuals at higher risk from the disease and at the epidemiological one to explain the differences in the epidemic trend among countries or even within a specific country. Covid-19 disease in Italy showed a peculiar geographical distribution from the northern most affected regions to the southern ones only slightly touched.

METHODS

In this study we analysed the regional frequencies for the most common Italian haplotypes from the Italian Bone Marrow Donor Registry (HLA-A, -B, -C and -DRB1 at four-digit level). Then we performed Pearson correlation analyses among regional haplotypes estimated frequency in the population and Covid-19 incidence and mortality.

RESULTS

In this study we found that the two most frequent HLA haplotypes in the Italian population, HLA-A*:01:01g-B*08:01 g-C*07:01g-DRB1*03:01g and HLA-A*02.01g-B*18.01g-C*07.01g-DRB1*11.04g, had a regional distribution overlapping that of Covid-19 and showed respectively a positive (suggestive of susceptibility) and negative (suggestive of protection) significant correlation with both Covid-19 incidence and mortality.

CONCLUSIONS

Based on these results, in order to define such HLA haplotypes as a factor effectively associated to the disease susceptibility, the creation of national networks that can collect patients' samples from all regions for HLA typing should be highly encouraged.

Improvement in HLA-C typing by a new sequence-specific oligonucleotides kit

PCR-sequence-specific oligonucleotide (PCR-SSO) is commonly used for HLA-typing. We recently replaced LabType SSO HD kit for HLA-C typing with the XR kit (OneLambda, Inc.). We used 137 patients' samples representing unique most likely HLA-C allele combinations to compare these kits' performance. The XR kit decreased the number of allele ambiguities with a median of 26.1% (first to third quartiles 20% to 36.2%, range 0% to 96.1%). The XR kit did not resolve all the common, intermediate and well-documented HLA allele ambiguities, which may be important in the clinical setting. The XR kit eliminated 23.6% of null allele ambiguities. In conclusion, the HLA-C XR kit provides a moderate yet valuable improvement of HLA-typing resolution in comparison with the HLA-C HD kit.

Frequencies and haplotype associations of non-expressed HLA alleles in ethnically diverse populations on the National Marrow Donor Program's Be The Match Registry

HLA allele matching is critical to successful bone marrow transplantation between a patient and donor. Non-functional HLA alleles, so called 'null alleles', are not well described within a large population of well HLA-typed ethnically diverse individuals despite their impact on donor selection. A retrospective analysis was performed on 833,789 unrelated donors (URDs) in the National Marrow Donor Program's Be The Match Registry® typed for HLA-A, -B, -C, -DRB1, -DQB1, and -DPB1 by next-generation DNA sequencing. Results showed that null alleles occur in low frequency (2.30E-04) compared to expressed alleles. Their overall frequency ranged from 6.00E-07 to 9.25E-04 with a median of 1.20E-06. The expected allele associations were commonly observed for HLA-A*24:09N, HLA-B*51:11N, and HLA-C*04:09N; however, associations outside of the expected were also observed. Notably, 82% of the National Marrow Donor Program Registry URDs carrying HLA-A*24:11N showed a different HLA-C allele association, HLA-C*05:01, compared to the allele described by prior published work characterizing German donor populations, HLA-C*04:01. The frequencies of these observed null alleles and linkage disequilibrium information could be invaluable and helpful in guiding the HLA testing decisions.

Concordance between predicted HLA type using next generation sequencing data generated for non-HLA purposes and clinical HLA type

We explored the feasibility of obtaining accurate HLA type using pre-existing NGS data not generated for HLA purposes. 83 exomes and 500 targeted NGS pharmacogenomic panels were analyzed using Omixon HLA Explore, OptiType, and/or HLA-Genotyper software. Results were compared against clinical HLA genotyping. 765 (94.2%) Omixon and 769 (94.7%) HLA-Genotyper of 812 germline allele calls across class I/II loci and 402 (99.5%) of 404 OptiType class I calls were concordant to the second field (i.e. HLA-A*02:01). An additional 19 (2.3%) Omixon, 39 (4.8%) HLA-Genotyper, and 2 (0.5%) OptiType allele calls were first field concordant (i.e. HLA-A*02). Using Omixon, four alleles (0.4%) were discordant and 24 (3.0%) failed to call, while 4 alleles (0.4%) were discordant using HLA-Genotyper. Tumor exomes were also evaluated and were 85.4%, 91.6%, and 100% concordant (Omixon and HLA-Genotyper with 96 alleles tested, and Optitype with 48 class I alleles, respectively). The 15 exomes and 500 pharmacogenomic panels were 100% concordant for each pharmacogenomic allele tested. This work has broad implications spanning future clinical care (pharmacogenomics, tumor response to immunotherapy, autoimmunity, etc.) and research applications.

Common, intermediate and well-documented HLA alleles in world populations: CIWD version 3.0.0

A catalog of common, intermediate and well-documented (CIWD) HLA-A, -B, -C, -DRB1, -DRB3, -DRB4, -DRB5, -DQB1 and -DPB1 alleles has been compiled from over 8 million individuals using data from 20 unrelated hematopoietic stem cell volunteer donor registries. Individuals are divided into seven geographic/ancestral/ethnic groups and data are summarized for each group and for the total population. P (two-field) and G group assignments are divided into one of four frequency categories: common (≥1 in 10 000), intermediate (≥1 in 100 000), well-documented (≥5 occurrences) or not-CIWD. Overall 26% of alleles in IPD-IMGT/HLA version 3.31.0 at P group resolution fall into the three CIWD categories. The two-field catalog includes 18% (n = 545) common, 17% (n = 513) intermediate, and 65% (n = 1997) well-documented alleles. Full-field allele frequency data are provided but are limited in value by the variations in resolution used by the registries. A recommended CIWD list is based on the most frequent category in the total or any of the seven geographic/ancestral/ethnic groups. Data are also provided so users can compile a catalog specific to the population groups that they serve. Comparisons are made to three previous CWD reports representing more limited population groups. This catalog, CIWD version 3.0.0, is a step closer to the collection of global HLA frequencies and to a clearer view of HLA diversity in the human population as a whole.

Synthesis of site-specific antibody-drug conjugates by ADP-ribosyl cyclases

Most of the current antibody-drug conjugates (ADCs) in clinic are heterogeneous mixtures. To produce homogeneous ADCs, established procedures often require multiple steps or long reaction times. The introduced mutations or foreign sequences may cause high immunogenicity. Here, we explore a new concept of transforming CD38 enzymatic activity into a facile approach for generating site-specific ADCs. This was achieved through coupling bifunctional antibody-CD38 fusion proteins with designer dinucleotide-based covalent inhibitors with stably attached payloads. The resulting adenosine diphosphate-ribosyl cyclase-enabled ADC (ARC-ADC) with a drug-to-antibody ratio of 2 could be rapidly generated through single-step conjugation. The generated ARC-ADC targeting human epidermal growth factor receptor 2 (HER2) displays excellent stability and potency against HER2-positive breast cancer both in vitro and in vivo. This proof-of-concept study demonstrates a new strategy for production of site-specific ADCs. It may provide a general approach for the development of a novel class of ADCs with potentially enhanced properties.