2.7 million samples genotyped for HLA by next generation sequencing: lessons learned

Association of human leucocyte antigen loci with vaccine-induced immune thrombotic thrombocytopenia: Potential role of the interaction between platelet factor 4-derived peptides and MHC-II

No risk factors have been identified for vaccine-induced immune thrombotic thrombocytopenia (VITT) so far. The aim of this study was to identify human leucocyte antigen (HLA) alleles potentially associated with VITT susceptibility. Specific HLA class II alleles were detected with significantly higher frequency in VITT patients compared with Italian controls: DPB1*17:01, DQA1*05:01, and DRB1*11:04. In silico analysis revealed increased affinity of DRB1*11:04 for a platelet factor 4 (PF4)-derived peptide, ITSLEVIKA, that contains two amino acids present in the specific binding site of anti-PF4 antibodies from VITT patients. Our findings show for the first time a genetic predisposition to developing anti-PF4 antibodies in response to Ad-vector vaccines.

The 18th International HLA & Immunogenetics workshop project report: Creating fully representative MHC reference haplotypes

A fundamental endeavor of the International Histocompatibility and Immunogenetics Workshop (IHIW) was assembling a collection of DNA samples homozygous through the MHC genomic region. This collection proved invaluable for assay development in the histocompatibility and immunogenetics field, for generating the human reference genome, and furthered our understanding of MHC diversity. Defined by their HLA-A, -B, -C and -DRB1 alleles, the combined frequency of the haplotypes from these individuals is ~20% in Europe. Thus, a significant proportion of MHC haplotypes, both common and rare throughout the world, and including many associated with disease, are not yet represented. In this workshop component, we are collecting the next generation of MHC -homozygous samples, to expand, diversify and modernize this critical community resource that has been foundational to the field. We asked laboratories worldwide to identify samples homozygous through all HLA class I and/or HLA class II genes, or through whole-genome SNP genotyping or sequencing, to have extensive homozygosity tracts within the MHC region. The focus is non-Europeans or those having HLA haplotypes less common in Europeans. Through this effort, we have obtained samples from 537 individuals representing 294 distinct haplotypes, as determined by their HLA class I and II alleles, and an additional 50 haplotypes distinct in HLA class I or II alleles. Although we have expanded the diversity, many populations remain underrepresented, particularly from Africa, and we encourage further participation. The data will serve as a resource for investigators seeking to characterize variation across the MHC genomic region for disease and population studies.

Body odours as putative chemosignals in the father-child relationship: New insights on paternal olfactory kin recognition and preference from infancy to adolescence

Children's body odours are effective chemical cues in the parent-child relationship. Mothers can recognize the odour of their child and prefer this odour over that of unfamiliar children. This effect is mediated by genetic similarity and developmental stage and is therefore suited to promote parental care at pre-pubertal stage, while facilitating incest avoidance at (post-)pubertal stage. The present study tested whether similar mechanisms apply to fathers. Therefore n = 56 fathers evaluated body odour samples of their own and of unfamiliar children in varying genetic and developmental stages. Genetic status was determined by human leucocyte antigen (HLA) profiling, developmental status by standardized assessment of pubertal status and steroid hormone concentration (estradiol, testosterone). Similar to mothers, fathers identified their own child's body odour above chance and preferred that odour. The paternal preference did not relate to HLA similarity but decreased with increasing age of the child. The decline was associated with higher pubertal stages in daughters only, which supports the hypothesis of odour-mediated incest prevention in opposite-sex parent-child dyads.

Providing hematopoietic stem cell products from unrelated donors to the world: DKMS donor centers and DKMS Registry

Allogeneic hematopoietic stem cell (HSC) transplantation is a curative therapy for many severe blood diseases. As many patients have no suitable family donor, large unrelated donor registries and donor centers have been established in many countries, along with an international system for the provision of unrelated donor HSC products. As an essential part of this system, DKMS operates donor centers in 7 countries with a total of 12.2 million donors and over 114,000 donations so far, and a multinational donor registry. In 2022, DKMS donors contributed 57.5% of all cross-border donations worldwide. In this review, we describe the international system for the provision of unrelated donor HSC products as well as tasks and responsibilities of donor registries and donor centers. We also discuss relevant aspects of DKMS donor centers, namely donor file composition, matching and donation probabilities and actual donations, and the unique multinational approach of the DKMS Registry.

High population frequencies of MICA copy number variations originate from independent recombination events

MICA is a stress-induced ligand of the NKG2D receptor that stimulates NK and T cell responses and was identified as a key determinant of anti-tumor immunity. The MICA gene is located inside the MHC complex and is in strong linkage disequilibrium with HLA-B. While an HLA-B*48-linked MICA deletion-haplotype was previously described in Asian populations, little is known about other MICA copy number variations. Here, we report the genotyping of more than two million individuals revealing high frequencies of MICA duplications (1%) and MICA deletions (0.4%). Their prevalence differs between ethnic groups and can rise to 2.8% (Croatia) and 9.2% (Mexico), respectively. Targeted sequencing of more than 70 samples indicates that these copy number variations originate from independent nonallelic homologous recombination events between segmental duplications upstream of MICA and MICB. Overall, our data warrant further investigation of disease associations and consideration of MICA copy number data in oncological study protocols.

Novel alleles in the era of next-generation sequencing-based HLA typing calls for standardization and policy

Next-Generation Sequencing (NGS) has transformed clinical histocompatibility laboratories through its capacity to provide accurate, high-throughput, high-resolution typing of Human Leukocyte Antigen (HLA) genes, which is critical for transplant safety and success. As this technology becomes widely used for clinical genotyping, histocompatibility laboratories now have an increased capability to identify novel HLA alleles that previously would not be detected using traditional genotyping methods. Standard guidelines for the clinical verification and reporting of novelties in the era of NGS are greatly needed. Here, we describe the experience of a clinical histocompatibility laboratory's use of NGS for HLA genotyping and its management of novel alleles detected in an ethnically-diverse population of British Columbia, Canada. Over a period of 18 months, 3,450 clinical samples collected for the purpose of solid organ or hematopoietic stem cell transplantation were sequenced using NGS. Overall, 29 unique novel alleles were identified at a rate of ∼1.6 per month. The majority of novelties (52%) were detected in the alpha chains of class II (HLA-DQA1 and -DPA1). Novelties were found in all 11 HLA classical genes except for HLA-DRB3, -DRB4, and -DQB1. All novelties were single nucleotide polymorphisms, where more than half led to an amino acid change, and one resulted in a premature stop codon. Missense mutations were evaluated for changes in their amino acid properties to assess the potential effect on the novel HLA protein. All novelties identified were confirmed independently at another accredited HLA laboratory using a different NGS assay and platform to ensure validity in the reporting of novelties. The novel alleles were submitted to the Immuno Polymorphism Database-Immunogenetics/HLA (IPD-IMGT/HLA) for official allele name designation and inclusion in future database releases. A nationwide survey involving all Canadian HLA laboratories confirmed the common occurrence of novel allele detection but identified a wide variability in the assessment and reporting of novelties. In summary, a considerable proportion of novel alleles were identified in routine clinical testing. We propose a framework for the standardization of policies on the clinical management of novel alleles and inclusion in proficiency testing programs in the era of NGS-based HLA genotyping.

5 years DKMS Chile: approach, results and impact of the first unrelated stem cell donor center in Chile

Introduction

Allogeneic hematopoietic stem cell transplantation (allo-HSCT) is performed worldwide to treat blood cancer and other life-threatening blood disorders. As successful transplantation requires an HLA-compatible donor, unrelated donor centers and registries have been established worldwide to identify donors for patients without a family match. Ethnic minorities are underrepresented in large donor registries. Matching probabilities are higher when donors and patients share the same ethnic background, making it desirable to increase the diversity of the global donor pool by recruiting donors in new regions. Here, we report the establishment and the first 5 years of operation of the first unrelated stem cell donor center in Chile, a high-income country in South America with a population of over 19 million.

Methods

We used online and in-person donor recruitment practices through patient appeals and donor drives in companies, universities, the armed forces, and public services. After confirmatory typing donors were subjected to medical work-up and cleared for donation.

Results

We recruited almost 170,000 donors in 5 years. There were 1,488 requests received for confirmatory typing and donor availability checks, of which 333 resulted in medical work-up, leading to 194 stem cell collections. Products were shipped to Chile (48.5%) and abroad. Even when the COVID-19 pandemic challenged our activities, the number of donors recruited and shipped stem cell products remained steady. In Chile there was an almost 8-fold increase in unrelated donor transplantation activity from 16 procedures in 2016-2018 to 124 procedures in 2019-2021, mainly for pediatric patients following the center's establishment. We estimate that 49.6% of Chilean patients would find at least one matched unrelated donor in the global DKMS donor pool.

Discussion

Establishing a DKMS donor center in Chile has significantly increased donor availability for Chilean patients and contributed to an increase of unrelated donor stem cell transplant activity.

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.

Tutorial: a statistical genetics guide to identifying HLA alleles driving complex disease

The human leukocyte antigen (HLA) locus is associated with more complex diseases than any other locus in the human genome. In many diseases, HLA explains more heritability than all other known loci combined. In silico HLA imputation methods enable rapid and accurate estimation of HLA alleles in the millions of individuals that are already genotyped on microarrays. HLA imputation has been used to define causal variation in autoimmune diseases, such as type I diabetes, and in human immunodeficiency virus infection control. However, there are few guidelines on performing HLA imputation, association testing, and fine mapping. Here, we present a comprehensive tutorial to impute HLA alleles from genotype data. We provide detailed guidance on performing standard quality control measures for input genotyping data and describe options to impute HLA alleles and amino acids either locally or using the web-based Michigan Imputation Server, which hosts a multi-ancestry HLA imputation reference panel. We also offer best practice recommendations to conduct association tests to define the alleles, amino acids, and haplotypes that affect human traits. Along with the pipeline, we provide a step-by-step online guide with scripts and available software ( https://github.com/immunogenomics/HLA_analyses_tutorial ). This tutorial will be broadly applicable to large-scale genotyping data and will contribute to defining the role of HLA in human diseases across global populations.

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.

HLA allele and haplotype frequencies of registered stem cell donors in Chile

Patients in need of hematopoietic stem cell transplantation often rely on unrelated stem cell donors matched in certain human leukocyte antigen (HLA) genes. Donor search is complicated by the extensive allelic variability of the HLA system. Therefore, large registries of potential donors are maintained in many countries worldwide. Population-specific HLA characteristics determine the registry benefits for patients and also the need for further regional donor recruitment. In this work, we analyzed HLA allele and haplotype frequencies of donors of DKMS Chile, the first Chilean donor registry, with self-assessed "non-Indigenous" (n=92,788) and "Mapuche" (n=1,993) ancestry. We identified HLA alleles that were distinctly more abundant in the Chilean subpopulations than in worldwide reference populations, four of them particularly characteristic for the Mapuche subpopulation, namely B*39:09g, B*35:09, DRB1*04:07g, and DRB1*16:02g. Both population subsamples carried haplotypes of both Native American and European origin at high frequencies, reflecting Chile's complex history of admixture and immigration. Matching probability analysis revealed limited benefits for Chilean patients (both non-Indigenous and Mapuche) from donor registries of non-Chilean donors, thus indicating a need for ongoing significant donor recruitment efforts in Chile.

Validation of a genotyping technique for a surrogate marker of HLA-B∗58:01 for allopurinol-induced Stevens-Johnson syndrome and toxic epidermal necrolysis in the Japanese population

Stevens-Johnson syndrome and toxic epidermal necrolysis (SJS/TEN) are rare but severe cutaneous adverse drug reactions. Certain human leukocyte antigen (HLA) types have been associated with SJS/TEN onset, e.g., HLA-B^∗58:01 with allopurinol-induced SJS/TEN, but HLA typing is time-consuming and expensive; thus, it is not commonly used in clinical situations. In the previous work, we demonstrated that the single-nucleotide polymorphisms (SNP) rs9263726 was in absolute linkage disequilibrium with HLA-B^∗58:01 in the Japanese population, and can be used as a surrogate marker for the HLA. Here, we developed a new genotyping method for the surrogate SNP using the single-stranded tag hybridization chromatographic printed-array strip (STH-PAS) technique and performed an analytical validation. The results of genotyping rs9263726 using STH-PAS correlated well with those obtained using the TaqMan SNP Genotyping Assay for 15 HLA-B^∗58:01-positive and 13 HLA-B^∗58:01-negative patients (analytical sensitivity and specificity were both 100%). Additionally, at least 1.11 ng of genomic DNA was sufficient to digitally and manually detect positive signals on the strip. Robustness studies showed that the annealing temperature (66 °C) was the most important condition related to reliable results. Collectively, we developed an STH-PAS method that can rapidly and easily detect rs9263726 for predicting SJS/TEN onset.

Simultaneous genotyping for human platelet antigen systems and HLA-A and HLA-B loci by targeted next-generation sequencing

In order to treat the alloimmunization platelet transfusion refractoriness (PTR), human leukocyte antigen (HLA)-type and/or human platelet antigen (HPA)-type matched platelets between donors and patients are usually used. Therefore, genotyping of HLA-A and HLA-B loci, as well as HPA systems, for donors and patients, is of great significance. However, there is a rare report of genotyping for HLA-A and HLA-B loci as well as HPA systems at the same time. In this study, a high-throughput method for simultaneous genotyping of HLA-A and HLA-B loci, as well as HPA genotyping, was developed. A RNA capture probe panel was designed covering all exon sequences of the GP1BA, GP1BB, ITGA2, CD109, ITGB3, and ITGA2B genes and HLA-A and HLA-B loci. The HLA-A, HLA-B, and 34 HPA systems were genotyped using a targeted next-generation sequencing (NGS) method. The genotypes of the HLA-A and HLA-B loci, as well as the HPA, were assigned based on the nucleotides in the polymorphism sites. Using the NGS method, 204 unrelated blood specimens were successfully genotyped for all 34 HPA systems as well as HLA-A and HLA-B loci. The accuracy of the NGS method was 100%. Only HPA-2, HPA-3, HPA-5, HPA-6w, HPA-15, and HPA-21w showed polymorphism with frequencies of 0.9412, 0.6863, 0.9853, 0.9779, 0.4314, and 0.9951 for a allele, respectively. Thirty-two single nucleotide variants (SNVs) were detected. Of them, 12 SNVs can lead to amino acid change. HLA-A*11:01 and HLA-B*46:01 are the most common alleles for HLA-A and HLA-B loci. A targeted next-generation sequencing method for simultaneously genotyping HPA systems and HLA-A and HLA-B loci was first established, which could be used to create a database of HLA-typed and/or HPA-typed unrelated donors.

A Human and Rhesus Macaque Interferon-Stimulated Gene Screen Shows That Over-Expression of ARHGEF3/XPLN Inhibits Replication of Hepatitis C Virus and Other Flavivirids

Natural hepatitis C virus (HCV) infection is restricted to humans, whereas other primates such as rhesus macaques are non-permissive for infection. To identify human and rhesus macaque genes that differ or share the ability to inhibit HCV replication, we conducted a medium-throughput screen of lentivirus-expressed host genes that disrupt replication of HCV subgenomic replicon RNA expressing secreted Gaussia luciferase. A combined total of >800 interferon-stimulated genes (ISGs) were screened. Our findings confirmed established anti-HCV ISGs, such as IRF1, PKR and DDX60. Novel species−specific inhibitors were also identified and independently validated. Using a cell-based system that recapitulates productive HCV infection, we identified that over-expression of the ‘Rho Guanine Nucleotide Exchange Factor 3’ gene (ARHGEF3) from both species inhibits full-length virus replication. Additionally, replication of two mosquito-borne flaviviruses, yellow fever virus (YFV) and Zika virus (ZIKV), were also reduced in cell lines over-expressing ARHGEF3 compared to controls. In conclusion, we ascribe novel antiviral activity to the cellular gene ARHGEF3 that inhibits replication of HCV and other important human viral pathogens belonging to the Flaviviridae, and which is conserved between humans and rhesus macaques.

Host KIR/HLA-C Genotypes Determine HIV-Mediated Changes of the NK Cell Repertoire and Are Associated With Vpu Sequence Variations Impacting Downmodulation of HLA-C

NK cells play a pivotal role in viral immunity, utilizing a large array of activating and inhibitory receptors to identify and eliminate virus-infected cells. Killer-cell immunoglobulin-like receptors (KIRs) represent a highly polymorphic receptor family, regulating NK cell activity and determining the ability to recognize target cells. Human leukocyte antigen (HLA) class I molecules serve as the primary ligand for KIRs. Herein, HLA-C stands out as being the dominant ligand for the majority of KIRs. Accumulating evidence indicated that interactions between HLA-C and its inhibitory KIR2DL receptors (KIR2DL1/L2/L3) can drive HIV-1-mediated immune evasion and thus may contribute to the intrinsic control of HIV-1 infection. Of particular interest in this context is the recent observation that HIV-1 is able to adapt to host HLA-C genotypes through Vpu-mediated downmodulation of HLA-C. However, our understanding of the complex interplay between KIR/HLA immunogenetics, NK cell-mediated immune pressure and HIV-1 immune escape is still limited. Therefore, we investigated the impact of specific KIR/HLA-C combinations on the NK cell receptor repertoire and HIV-1 Vpu protein sequence variations of 122 viremic, untreated HIV-1⁺ individuals. Compared to 60 HIV-1^- controls, HIV-1 infection was associated with significant changes within the NK cell receptor repertoire, including reduced percentages of NK cells expressing NKG2A, CD8, and KIR2DS4. In contrast, the NKG2C⁺ and KIR3DL2⁺ NK cell sub-populations from HIV-1⁺ individuals was enlarged compared to HIV-1^- controls. Stratification along KIR/HLA-C genotypes revealed a genotype-dependent expansion of KIR2DL1⁺ NK cells that was ultimately associated with increased binding affinities between KIR2DL1 and HLA-C allotypes. Lastly, our data hinted to a preferential selection of Vpu sequence variants that were associated with HLA-C downmodulation in individuals with high KIR2DL/HLA-C binding affinities. Altogether, our study provides evidence that HIV-1-associated changes in the KIR repertoire of NK cells are to some extent predetermined by host KIR2DL/HLA-C genotypes. Furthermore, analysis of Vpu sequence polymorphisms indicates that differential KIR2DL/HLA-C binding affinities may serve as an additional mechanism how host genetics impact immune evasion by HIV-1.

A novel swab storage gel is superior to dry swab DNA collection, and enables long-range high resolution NGS HLA typing from buccal cell samples

HLA sequence-based DNA typing (SBT) by long-range PCR amplification (LR PCR) and next-generation sequencing (NGS) is a high-throughput DNA sequencing method (LR-NGS-SBT) for the efficient and sensitive detection of novel and null HLA alleles to the field-4 level of allelic resolution without phase ambiguity. However, the accuracy and reliability of the HLA typing results using buccal cells (BCs) and saliva as genetic source materials for the LR-NGS-SBT method are dependent largely on the quality of the extracted genomic DNA (gDNA) because a large degree of gDNA fragmentation can result in insufficient PCR amplification with the incorrect assignment of HLA alleles due to allele dropouts. In this study, we developed a new cost-efficient swab storage gel (SSG) for wet swab collection of BCs (BC-SSG) and evaluated its usefulness by performing different DNA analytical parameters including LR-NGS-SBT to compare the quality and quantity of gDNA extracted from BCs (in SSG or air dried), blood and saliva of 30 subjects. The BC-SSG samples after 5 days of storage revealed qualitative and quantitative DNA values equivalent to that of blood and/or saliva and better than swabs that were only air-dried (BC-nSSG). Moreover, all the gDNA extracted from blood, saliva and BC-SSG samples were HLA-typed successfully to an equivalent total of 408 alleles for each sample type. Therefore, the application of BC-SSG collection media for LR-NGS-SBT has benefits over BC dried samples (dry swabs) such as reducing retesting and the number of untestable BC samples due to insufficient DNA amplification. This article is protected by copyright. All rights reserved.

Personalized HLA typing leads to the discovery of novel HLA alleles and tumor‐specific HLA variants

Accurate and full‐length typing of the HLA region is important in many clinical and research settings. With the advent of next generation sequencing (NGS), several HLA typing algorithms have been developed, including many that are applicable to whole exome sequencing (WES). However, most of these solutions operate by providing the closest‐matched HLA allele among the known alleles in IPD‐IMGT/HLA Database. These database‐matching approaches have demonstrated very high performance when typing well characterized HLA alleles. However, as they rely on the completeness of the HLA database, they are not optimal for detecting novel or less well characterized alleles. Furthermore, the database‐matching approaches are also not adequate in the context of cancer, where a comprehensive characterization of somatic HLA variation and expression patterns of a tumor's HLA locus may guide therapy and clinical outcome, because of the pivotal role HLA alleles play in tumor antigen recognition and immune escape. Here, we describe a personalized HLA typing approach applied to WES data that leverages the strengths of database‐matching approaches while simultaneously allowing for the discovery of novel HLA alleles and tumor‐specific HLA variants, through the systematic integration of germline and somatic variant calling. We applied this approach on WES from 10 metastatic melanoma patients and validated the HLA typing results using HLA targeted NGS sequencing from patients where at least one HLA germline candidate was detected on Class I HLA. Targeted NGS sequencing confirmed 100% performance for the 1st and 2nd fields. In total, five out of the six detected HLA germline variants were because of Class I ambiguities at the third or fourth fields, and their detection recovered the correct HLA allele genotype. The sixth germline variant let to the formal discovery of a novel Class I allele. Finally, we demonstrated a substantially improved somatic variant detection accuracy in HLA alleles with a 91% of success rate in simulated experiments. The approach described here may allow the field to genotype more accurately using WES data, leading to the discovery of novel HLA alleles and help characterize the relationship between somatic variation in the HLA region and immunosurveillance.

Graphical user interface for the haplotype frequency estimation software Hapl-o-Mat

Population-specific human leukocyte antigen (HLA) haplotype frequencies are an essential basis of advanced algorithms for donor selection in unrelated hematopoietic stem cell transplantation. In 2016, we introduced Hapl-o-Mat, a versatile tool for haplotype frequency estimation based on an expectation-maximization algorithm (https://github.com/DKMS/hapl-o-Mat). Hapl-o-Mat is specifically tailored to the analysis of HLA genes and able to cope with the heterogeneous genotyping data usually found in donor registries. To make Hapl-o-Mat accessible to a wider range of users, we designed a graphical user interface module that considerably facilitates the interaction with the application (https://github.com/DKMS/hapl-o-Mat_GUI). We further provide a precompiled version of Hapl-o-Mat that can be used on Windows personal computers without dependency on additional software libraries (https://github.com/DKMS/hapl-o-Mat_WinBin). We are confident that these new, user-oriented features will encourage more researchers to apply Hapl-o-Mat to their data, thereby increasing knowledge and public availability of population-specific HLA haplotype frequencies.

Noninvasive Determination of CMV Serostatus From Dried Buccal Swab Samples: Assay Development, Validation, and Application to 1.2 Million Samples

BACKGROUND

Buccal swab sampling constitutes an attractive noninvasive alternative to blood drawings for antibody serostatus assays. Here we describe a method to determine the cytomegalovirus immunoglobulin G (CMV IgG) serostatus from dried buccal swab samples.

METHODS

Upon solubilization, CMV IgG is determined by an ELISA assay specifically adapted to cope with low IgG concentrations. The derived CMV titer is normalized against the total protein concentration to adjust for incorrectly or less efficiently sampled buccal swabs. Assay parameters were optimized on a set of 713 samples.

RESULTS

Validation with 1784 samples revealed distinct results for > 80% of samples with 98.6% specificity and 99.1% sensitivity. Based on the analysis of 1.2 million samples we derived age- and sex-stratified CMV prevalence statistics for Germany, Poland, United Kingdom, and Chile. To confirm accuracy of the assay in routine operation, the CMV status of 6518 donors was reassessed by independent laboratories based on conventional blood samples revealing 96.9% specificity and 97.4% sensitivity.

CONCLUSIONS

The assay accurately delivers the CMV IgG serostatus from dried buccal swab samples for > 80% of the participants. Thereby it provides a noninvasive alternative to plasma-based CMV monitoring for nondiagnostic purposes such as hematopoietic stem cell transplantation donor screening or population studies.

Estimating HLA haplotype frequencies from homozygous individuals - A Technical Report

We estimated HLA haplotype frequencies based on individuals homozygous for 4, 5 or 6 loci. Validation of our approach using a sample of over 3.4 million German individuals was successful. Compared to an expectation‐maximization algorithm, the errors were larger. However, our approach allows the unequivocal detection of rare haplotypes.

Next-generation sequencing and clinical histocompatibility testing

Histocompatibility testing is essential for donor identification and risk assessment in solid organ and hematopoietic stem cell transplant. Additionally, it is useful for identifying donor specific alleles for monitoring donor specific antibodies in post-transplant patients. Next-generation sequence (NGS) based human leukocyte antigen (HLA) typing has improved many aspects of histocompatibility testing in hematopoietic stem cell and solid organ transplant. HLA disease association testing and research has also benefited from the advent of NGS technologies. In this review we discuss the current impact and future applications of NGS typing on clinical histocompatibility testing for transplant and non-transplant purposes.

Relevance of Polymorphic KIR and HLA Class I Genes in NK-Cell-Based Immunotherapies for Adult Leukemic Patients

Simple Summary

Immunotherapies are promising approaches to curing different acute leukemias. Natural killer (NK) cells are lymphocytes that are efficient in the elimination of leukemic cells. NK-cell-based immunotherapies are particularly attractive, but the landscape of the heterogeneity of NK cells must be deciphered. This review provides an overview of the polymorphic KIR and HLA class I genes that modulate the NK cell repertoire and how these markers can improve the outcomes of patients with acute leukemia. A better knowledge of these genetic markers that are linked to NK cell subsets that are efficient against hematological diseases will optimize hematopoietic stem-cell donor selection and NK immunotherapy design.

Abstract

Since the mid-1990s, the biology and functions of natural killer (NK) cells have been deeply investigated in healthy individuals and in people with diseases. These effector cells play a particularly crucial role after allogeneic hematopoietic stem-cell transplantation (HSCT) through their graft-versus-leukemia (GvL) effect, which is mainly mediated through polymorphic killer-cell immunoglobulin-like receptors (KIRs) and their cognates, HLA class I ligands. In this review, we present how KIRs and HLA class I ligands modulate the structural formation and the functional education of NK cells. In particular, we decipher the current knowledge about the extent of KIR and HLA class I gene polymorphisms, as well as their expression, interaction, and functional impact on the KIR⁺ NK cell repertoire in a physiological context and in a leukemic context. In addition, we present the impact of NK cell alloreactivity on the outcomes of HSCT in adult patients with acute leukemia, as well as a description of genetic models of KIRs and NK cell reconstitution, with a focus on emergent T-cell-repleted haplo-identical HSCT using cyclosphosphamide post-grafting (haplo-PTCy). Then, we document how the immunogenetics of KIR/HLA and the immunobiology of NK cells could improve the relapse incidence after haplo-PTCy. Ultimately, we review the emerging NK-cell-based immunotherapies for leukemic patients in addition to HSCT.

Individual HLA-A, -B, -C, and -DRB1 Genotypes Are No Major Factors Which Determine COVID-19 Severity

HLA molecules are key restrictive elements to present intracellular antigens at the crossroads of an effective T-cell response against SARS-CoV-2. To determine the impact of the HLA genotype on the severity of SARS-CoV-2 courses, we investigated data from 6,919 infected individuals. HLA-A, -B, and -DRB1 allotypes grouped into HLA supertypes by functional or predicted structural similarities of the peptide-binding grooves did not predict COVID-19 severity. Further, we did not observe a heterozygote advantage or a benefit from HLA diplotypes with more divergent physicochemical peptide-binding properties. Finally, numbers of in silico predicted viral T-cell epitopes did not correlate with the severity of SARS-CoV-2 infections. These findings suggest that the HLA genotype is no major factor determining COVID-19 severity. Moreover, our data suggest that the spike glycoprotein alone may allow for abundant T-cell epitopes to mount robust T-cell responses not limited by the HLA genotype.

Impact of Previously Unrecognized HLA Mismatches Using Ultrahigh Resolution Typing in Unrelated Donor Hematopoietic Cell Transplantation

PURPOSE

Ultrahigh resolution (UHR) HLA matching is reported to result in better outcomes following unrelated donor hematopoietic cell transplantation, improving survival and reducing post-transplant complications. However, most studies included relatively small numbers of patients. Here we report the findings from a large, multicenter validation study.

METHODS

UHR HLA typing was available on 5,140 conventionally 10 out of 10 HLA-matched patients with malignant disease transplanted between 2008 and 2017.

RESULTS

After UHR HLA typing, 82% of pairs remained 10 out of 10 UHR-matched; 12.3% of patients were 12 out of 12 UHR HLA-matched. Compared with 12 out of 12 UHR-matched patients, probabilities of grade 2-4 acute graft-versus-host disease (aGVHD) were significantly increased with UHR mismatches (overall P = .0019) and in those patients who were HLA-DPB1 T-cell epitope permissively mismatched or nonpermissively mismatched (overall P = .0011). In the T-cell-depleted subset, the degree of UHR HLA mismatch was only associated with increased transplant-related mortality (TRM) (overall P = .0068). In the T-cell-replete subset, UHR HLA matching was associated with a lower probability of aGVHD (overall P = .0020); 12 out of 12 UHR matching was associated with reduced TRM risk when compared with HLA-DPB1 T-cell epitope permissively mismatched patients, whereas nonpermissive mismatching resulted in a greater risk (overall P = .0003).

CONCLUSION

This study did not confirm that UHR 12 out of 12 HLA matching increases the probability of overall survival but does demonstrate that aGVHD risk, and in certain settings TRM, is lowest in UHR HLA-matched pairs and thus warrants consideration when multiple 10 out of 10 HLA-matched donors of equivalent age are available.

HLA-E typing of more than 2.5 million potential hematopoietic stem cell donors: Methods and population-specific allele frequencies

We developed a cost-efficient workflow for genotyping HLA-E by NGS and applied it for genotyping more than 2.5 million potential stem cell donors. The data obtained were used to determine HLA-E allele frequency distributions for 104 populations. Our results confirm the known dominance of the alleles E*01:01 and E*01:03, which have a combined frequency of more than 0.99 in 97 of the 104 populations. E*01:01 is more frequent in Africa and the western part of South America, E*01:03 in Southeast and East Asia. E*01:03 shows a pronounced regional substructure at the high-resolution level with E*01:03:01G being particularly common in a large connected region extending from Turkey to China, E*01:03:02G in Northwestern Europe and E*01:03:03 in Central and Eastern Europe as well as Central Asia. The presented results are relevant both as a basis for further population genetics studies and for optimizing stem cell donor searches.

DR2S: an integrated algorithm providing reference-grade haplotype sequences from heterozygous samples

Background

High resolution HLA genotyping of donors and recipients is a crucially important prerequisite for haematopoetic stem-cell transplantation and relies heavily on the quality and completeness of immunogenetic reference sequence databases of allelic variation.

Results

Here, we report on DR2S, an R package that leverages the strengths of two sequencing technologies—the accuracy of next-generation sequencing with the read length of third-generation sequencing technologies like PacBio’s SMRT sequencing or ONT’s nanopore sequencing—to reconstruct fully-phased high-quality full-length haplotype sequences. Although optimised for HLA and KIR genes, DR2S is applicable to all loci with known reference sequences provided that full-length sequencing data is available for analysis. In addition, DR2S integrates supporting tools for easy visualisation and quality control of the reconstructed haplotype to ensure suitability for submission to public allele databases.

Conclusions

DR2S is a largely automated workflow designed to create high-quality fully-phased reference allele sequences for highly polymorphic gene regions such as HLA or KIR. It has been used by biologists to successfully characterise and submit more than 500 HLA alleles and more than 500 KIR alleles to the IPD-IMGT/HLA and IPD-KIR databases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04153-0.

Concurrent typing of over 4000 samples by long-range PCR amplicon-based NGS and rSSO revealed the need to verify NGS typing for HLA allelic dropouts

Hematopoietic stem cell transplantation (HSCT) from HLA-matched donors significantly decreases the risks of graft-rejection and graft-versus-host disease. Long-range PCR- amplicon-based next-generation sequencing (NGS) is increasingly used as a standalone method in clinical laboratories to determine HLA compatibility for HSCT and solid-organ transplantation. We hypothesized that an allelic dropout is a frequent event in the long-range PCR amplicon-based NGS HLA typing method. To test the hypothesis, we typed 4,006 samples concurrently using a commercially available long-range PCR amplicon-based NGS-typing and short exon-specific amplicon-based reverse sequence-specific oligonucleotide (rSSO) methods. The concordance between the NGS and rSSO typing results was 100% at HLA-A, -B, -C, -DRB1, -DRB3, -DRB5, -DQA1, DPA1 loci. However, 4.5% of the samples (179/4006) showed allelic-dropouts at one of the other three loci: HLA-DRB4 (3.9%), HLA-DPB1 (0.4%), and HLA-DQB1*(0.15%). The allelic-dropouts are not associated with specific haplotypes, and some dropouts can be reagent lot-specific. Although DRB1-DRB3/4/5-DQB1 linkages help to diagnose these allelic-dropouts in some cases, the rSSO typing was crucial to identify the dropouts in DQB1 and DPB1 loci. These results uncover the critical limitations of using long-range PCR amplicon-based NGS as a standalone method in clinical histocompatibility laboratories and advocate the need for strategies to diagnose and resolve allelic-dropouts.

Challenges in the application of NGS in the clinical laboratory

Next-generation sequencing (NGS), also known as massively parallel sequencing, has revolutionized genomic research. The current advances in NGS technology make it possible to provide high resolution, high throughput HLA typing in clinical laboratories. The focus of this review is on the recent development and implementation of NGS in clinical laboratories. Here, we examine the critical role of NGS technologies in clinical immunology for HLA genotyping. Two major NGS platforms (Illumina and Ion Torrent) are characterized including NGS library preparation, data analysis, and validation. Challenges of NGS implementation in the clinical laboratory are also discussed, including sequencing error rate, bioinformatics, result interpretation, analytic sensitivity, as well as large data storage. This review aims to promote the broader applications of NGS technology in clinical laboratories and advocate for the novel applications of NGS to drive future research.

Stem cell donor registry activities during the COVID-19 pandemic: a field report by DKMS

The COVID-19 pandemic has serious implications also for patients with other diseases. Here, we describe the effects of the pandemic on unrelated hematopoietic stem cell donation and transplantation from the perspective of DKMS, a large international donor registry. Especially, we cover the development of PBSC and bone marrow collection figures, donor management including Health and Availability Check (HAC), transport and cryopreservation of stem cell products, donor recruitment and business continuity measures. The total number of stem cell products provided declined by around 15% during the crisis with a particularly strong decrease in bone marrow products. We modified donor management processes to ensure donor and product safety. HAC instead of confirmatory typing was helpful especially in countries with strict lockdowns. New transport modes were developed so that stem cell products could be safely delivered despite COVID-19-related travel restrictions. Cryopreservation of stem cell products became the new temporary standard during the pandemic to minimize risks related to transport logistics and donor availability. However, many products from unrelated donors will never be transfused. DKMS discontinued public offline donor recruitment, leading to a 40% decline in new donors during the crisis. Most DKMS employees worked from home to ensure business continuity during the crisis.

HLA RNA Sequencing With Unique Molecular Identifiers Reveals High Allele-Specific Variability in mRNA Expression

The HLA gene complex is the most important single genetic factor in susceptibility to most diseases with autoimmune or autoinflammatory origin and in transplantation matching. Most studies have focused on the vast allelic variation in these genes; only a few studies have explored differences in the expression levels of HLA alleles. In this study, we quantified mRNA expression levels of HLA class I and II genes from peripheral blood samples of 50 healthy individuals. The gene- and allele-specific mRNA expression was assessed using unique molecular identifiers, which enabled PCR bias removal and calculation of the number of original mRNA transcripts. We identified differences in mRNA expression between different HLA genes and alleles. Our results suggest that HLA alleles are differentially expressed and these differences in expression levels are quantifiable using RNA sequencing technology. Our method provides novel insights into HLA research, and it can be applied to quantify expression differences of HLA alleles in various tissues and to evaluate the role of this type of variation in transplantation matching and susceptibility to autoimmune diseases.

HLA-A, -B, -C, -DRB1 allele and haplotype frequencies of the Korean population and performance characteristics of HLA typing by next-generation sequencing

INTRODUCTION

Human leukocyte antigen (HLA) identification at the allelic level is important for haematopoietic stem cell transplantation (HSCT). Next-generation sequencing (NGS) resolves ambiguous alleles by determining the phase of the polymorphisms. The aim of this study was to validate the software for HLA-SBT (sequence-based typing), assess Korean allele frequency, and characterise the performance of NGS-HLA typing.

METHODS

From the 2009 to 2016 registry, 1293 unrelated healthy donors with a complete dataset of previously characterised HLA-A, -B, -C, and -DRB1 loci were selected and assessed for frequency, haplotype inference, and relative linkage disequilibrium. For performance characteristics of NGS-HLA, alleles included in 1293 cases and ambiguous or alleles assigned as new by SBT-HLA software, or unassigned alleles were included. A total of 91 and 41 quality control samples resulted in 1056 alleles (132 samples × 4 loci × 2 diploid) for analysis. The GenDx NGSgo kit was used for NGS-HLA typing using the Illumina MiSeq platform.

RESULTS

A panel of 132 samples covered 231 alleles, including 53 HLA-A, 80 HLA-B, 43 HLA-C, and 55 HLA-DRB1 by HLA-SBT typing. Comparison of SBT-HLA and NGS-HLA typing showed 99.7% (1053/1056) concordance and discrepant cases were resolved by manual evaluation. Typing by NGS resulted in 67 HLA-A, 112 HLA-B, 71 HLA-C, and 72 HLA-DRB1 alleles. A total of 132 ambiguous, 4 new, and 1 unassigned alleles by HLA-SBT were resolved by NGS-HLA typing.

CONCLUSIONS

NGS-HLA typing provided robust and conclusive results without ambiguities, and its implementation could support HSCT in clinical settings.

State-of-the-art genome inference in the human MHC

The Major Histocompatibility Complex (MHC) on the short arm of chromosome 6 is associated with more diseases than any other region of the genome; it encodes the antigen-presenting Human Leukocyte Antigen (HLA) proteins and is one of the key immunogenetic regions of the genome. Accurate genome inference and interpretation of MHC association signals have traditionally been hampered by the region's uniquely complex features, such as high levels of polymorphism; inter-gene sequence homologies; structural variation; and long-range haplotype structures. Recent algorithmic and technological advances have, however, significantly increased the accessibility of genetic variation in the MHC; these developments include (i) accurate SNP-based HLA type imputation; (ii) genome graph approaches for variation-aware genome inference from next-generation sequencing data; (iii) long-read-based diploid de novo assembly of the MHC; (iv) cost-effective targeted MHC sequencing methods. Applied to hundreds of thousands of samples over the last years, these technologies have already enabled significant biological discoveries, for example in the field of autoimmune disease genetics. Remaining challenges concern the development of integrated methods that leverage haplotype-resolved de novo assembly of the MHC for the development of improved MHC genotyping methods for short reads and the construction of improved reference panels for SNP-based imputation. Improved genome inference in the MHC can crucially contribute to an improved genetic and functional understanding of many immune-related phenotypes and diseases.

High-resolution HLA allele and haplotype frequencies of the Saudi Arabian population based on 45,457 individuals and corresponding stem cell donor matching probabilities

We estimated HLA allele and haplotype frequencies of the Saudi Arabian population from a sample of 45,457 registered stem cell donors. The most frequent HLA alleles were A*02:01g (18.5%), C*06:02g (16.1%), B*51:01g (14.1%), DRB1*07:01g (16.2%), DQB1*02:01g (30.5%), and DPB1*04:01g (33.6%). The most frequent 5-locus haplotypes were A*02:05g~C*06:02g~B*50:01g~DRB1*07:01g~DQB1*02:01g (1.73%), A*02:01g~C*06:02g~B*50:01g~DRB1*07:01g~DQB1*02:01g (1.66%), and A*26:01g~C*07:02g~B*08:01g~DRB1*03:01g~DQB1*02:01g (1.38%). Furthermore, we used the calculated haplotype frequencies to estimate stem cell donor matching probabilities for Saudi Arabian donor and patient populations under various matching requirements. These results are relevant for strategic donor registry planning in the Kingdom of Saudi Arabia.

Marriage does not relate to major histocompatibility complex: a genetic analysis based on 3691 couples

Optimization of chances for healthy offspring is thought to be one of the factors driving mate choice and compatibility of the major histocompatibility complex (MHC) is assumed to determine the offspring's fitness. While humans have been claimed to be able to perceive information of MHC compatibility via the olfactory channel, it remains unknown whether humans use such information for mate choice. By investigation of 3691 married couples, we observed that the high polymorphism of MHC leads to a low chance for homozygous offspring. MHC similarity between couples did not differ from chance, we hence observed no MHC effect in married couples. Hormonal contraception at the time of relationship initiation had no significant effect towards enhanced similarity. A low variety of alleles within a postcode area led to a higher likelihood of homozygous offspring. Based on this data, we conclude that there is no pattern of MHC dis-assortative mating in a genetically diverse Western society. We discuss the question of olfactory mate preference, in-group mating bias and the high polymorphism as potential explanations.

High Metabolic Function and Resilience of NKG2A-Educated NK Cells

Natural killer (NK) cells are an important component of the innate immune system for the control of intracellular pathogens and cancer cells. NK cells demonstrate heterogeneous expression of inhibitory surface receptors. Signaling through these various receptors during NK cell development promotes functionality, referred to as NK cell education. Here we investigated the impact of education on NK cell metabolism through functional assessment of critical metabolic pathways and calcium signaling. Educated NK cells had an increased uptake of the metabolic substrates 2-NBDG, a fluorescent glucose analog, and BODIPY FL C₁₆, a fluorescent palmitate, compared to uneducated NK cells. Comparison of NK cells educated via KIRs or NKG2A showed that NKG2A-educated NK cells were the main contributor to these differences in uptake of metabolites, and that NKG2A-educated NK cells were functionally more resilient in response to metabolic blockade of oxidative phosphorylation. Furthermore, NKG2A-educated NK cells exhibited higher peak calcium concentration following stimulation, indicating stronger signaling events taking place in these educated NK cells. These results demonstrate that cellular metabolism plays an important role in the functional differences observed between educated and uneducated NK cells, and show that NKG2A-educated NK cells remain more functionally competent than KIR-educated NK cells when oxidative phosphorylation is restricted. Understanding metabolic programming during NK cell education may unveil future targets to manipulate NK cell function for use in clinical settings, such as cancer therapies.

Association between the HLA-DRB1*03:01-DQB1*02:01 haplotype and PF4/heparin antibodies

HLA-DRB1*03:01-DQB1*02:01 haplotype patients develop higher anti–platelet factor 4/heparin (PF4/H) levels following IV heparin exposure. HLA-restricted antigen presentation may play a role in the formation of anti-PF4/H antibodies.

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.

Body odours as a chemosignal in the mother-child relationship: new insights based on an human leucocyte antigen-genotyped family cohort

Mothers are able to identify the body odour (BO) of their own child and prefer this smell above other BOs. It has hence been assumed that the infantile BO functions as a chemosignal promoting targeted parental care. We tested this hypothesis and examined whether children's BOs signal genetic similarity and developmental status to mothers. In addition, we assessed whether BOs facilitate inbreeding avoidance (Westermarck effect). In a cross-sectional design, N = 164 mothers participated with their biological children (N = 226 children, aged 0-18 years) and evaluated BO probes of their own and four other, sex-matched children. Those varied in age and in genetic similarity, which was assessed by human leucocyte antigen profiling. The study showed not only that mothers identified and preferred their own child's BO, but also that genetic similarity and developmental status are transcribed in BOs. Accordingly, maternal preference of their own child's odour changes throughout development. Our data partly supported the Westermarck effect: mothers' preference of pubertal boys' BOs was negatively related to testosterone for the own son, but not for unfamiliar children. This article is part of the Theo Murphy meeting issue 'Olfactory communication in humans'.

A long road/read to rapid high-resolution HLA typing: The nanopore perspective

Next-generation sequencing (NGS) has been widely adopted for clinical HLA typing and advanced immunogenetics researches. Current methodologies still face challenges in resolving cis-trans ambiguity involving distant variant positions, and the turnaround time is affected by testing volume and batching. Nanopore sequencing may become a promising addition to the existing options for HLA typing. The technology delivered by the MinION sequencer of Oxford Nanopore Technologies (ONT) can record the ionic current changes during the translocation of DNA/RNA strands through transmembrane pores and translate the signals to sequence reads. It features simple and flexible library preparations, long sequencing reads, portable and affordable sequencing devices, and rapid, real-time sequencing. However, the error rate of the sequencing reads is high and remains a hurdle for its broad application. This review article will provide a brief overview of this technology and then focus on the opportunities and challenges of using nanopore sequencing for high-resolution HLA typing and immunogenetics research.

High-Resolution HLA Typing of HLA-A, -B, -C, -DRB1, and -DQB1 in Kinh Vietnamese by Using Next-Generation Sequencing

Human leukocyte antigen (HLA) genotyping displays the particular characteristics of HLA alleles and haplotype frequencies in each population. Although it is considered the current gold standard for HLA typing, high-resolution sequence-based HLA typing is currently unavailable in Kinh Vietnamese populations. In this study, high-resolution sequence-based HLA typing (3-field) was performed using an amplicon-based next-generation sequencing platform to identify the HLA-A, -B, -C, -DRB1, and -DQB1 alleles of 101 unrelated healthy Kinh Vietnamese individuals from southern Vietnam. A total of 28 HLA-A, 41 HLA-B, 21 HLA-C, 26 HLA-DRB1, and 25 HLA-DQB1 alleles were identified. The most frequently occurring HLA alleles were A^∗11:01:01, B^∗15:02:01, C^∗07:02:01, DRB1^∗12:02:01, and DQB1^∗03:01:01. Haplotype calculation showed that A^∗29:01:01∼B^∗07:05:01, DRB1^∗12:02:01∼DQB1^∗3:01:01, A^∗29:01:01∼C^∗15:05:02∼B^∗07:05:01, A^∗33:03:01∼B^∗58:01:01∼DRB1^∗03:01:01, and A^∗29:01:01∼C^∗15:05:02∼B^∗07:05:01∼DRB1^∗10:01:01∼DQB1^∗05:01:01 were the most common haplotypes in the southern Kinh Vietnamese population. Allele distribution and haplotype analyses demonstrated that the Vietnamese population shares HLA features with South-East Asians but retains unique characteristics. Data from this study will be potentially applicable in medicine and anthropology.

Performance of a multiplexed amplicon-based next-generation sequencing assay for HLA typing

BACKGROUND

Next-generation sequencing (NGS) has enabled efficient high-resolution typing of human leukocyte antigen (HLA) genes with minimal ambiguity. Most commercially available assays amplify individual or subgroup of HLA genes by long-range PCR followed by library preparation and sequencing. The AllType assay simplifies the workflow by amplifying 11 transplant-relevant HLA genes in one PCR reaction. Here, we report the performance of this unique workflow evaluated using 218 genetically diverse samples.

METHODS

Five whole genes (HLA-A/B/C/DQA1/DPA1) and six near-whole genes (HLA-DRB1/DRB345/DQB1/DPB1; excluding exon 1 and part of intron 1) were amplified in a multiplexed, long-range PCR. Manual library preparation was performed per manufacturer's protocol, followed by template preparation and chip loading on the Ion Chef, and sequencing on the Ion S5 sequencer. Pre-specified rules for quality control and repeat testing were followed; technologists were blinded to the reference results. The concordance between AllType and reference results was determined at 2-field resolution. We also describe the ranges of input DNA and library concentrations, read number per sample and per locus, and key health metrics in relation to typing results.

RESULTS

The concordance rates were 98.6%, 99.8% and 99.9% at the sample (n = 218), genotype (n = 1688), and allele (n = 3376) levels, respectively. Three genotypes were discordant, all of which shared the same G group typing results with the reference. Most ambiguous genotypes (116 out of 144, 80.6%) were due to the lack of exon 1 and intron 1 coverage for HLA-DRB1/DRB345/DQB1/DPB1 genes. A broad range of input DNA concentrations and library concentrations were tolerated. Per sample read numbers were adequate for accurate genotyping. Per locus read numbers showed some inter-lot variations, and a trend toward improved inter-locus balance was observed with later lots of reagents.

CONCLUSION

The AllType assay on the Ion Chef/Ion S5 platform offers a robust and efficient workflow for clinical HLA typing at the 2-field resolution. The multiplex PCR strategy simplifies the laboratory procedure without compromising the typing accuracy.

Immunogenetics in stem cell donor registry work: The DKMS example (Part 2)

DKMS is a leading stem cell donor registry with more than 9 million donors. Donor registry activities share many touch points with topics from immunogenetics or population genetics. In this two-part review article, we deal with these aspects of donor registry work by using the example of DKMS. In the second part of the review, we focus on donor typing of non-HLA genes, the impact of donor age, gender and CMV serostatus on donation probabilities, the identification of novel HLA, KIR and MIC alleles by high-throughput donor typing, the activities of the Collaborative Biobank and pharmacogenetics in the donor registry context.

Estimation of German KIR Allele Group Haplotype Frequencies

The impact of the highly polymorphic Killer-cell immunoglobulin-like receptor (KIR) gene cluster on the outcome of hematopoietic stem cell transplantation (HCST) is subject of current research. To further understand the involvement of this gene family into Natural Killer (NK) cell-mediated graft-versus-leukemia reactions, knowledge of haplotype structures, and allelic linkage is of importance. In this analysis, we estimate population-specific KIR haplotype frequencies at allele group resolution in a cohort of n = 458 German families. We addressed the polymorphism of the KIR gene complex and phasing ambiguities by a combined approach. Haplotype inference within first-degree family relations allowed us to limit the number of possible diplotypes. Structural restriction to a pattern set of 92 previously described KIR copy number haplotypes further reduced ambiguities. KIR haplotype frequency estimation was finally accomplished by means of an expectation-maximization algorithm. Applying a resolution threshold of ½ n, we were able to identify a set of 551 KIR allele group haplotypes, representing 21 KIR copy number haplotypes. The haplotype frequencies allow studying linkage disequilibrium in two-locus as well as in multi-locus analyses. Our study reveals associations between KIR haplotype structures and allele group frequencies, thereby broadening our understanding of the KIR gene complex.

High-Throughput MICA/B Genotyping of Over Two Million Samples: Workflow and Allele Frequencies

MICA and MICB are ligands of the NKG2D receptor and thereby influence NK and T cell activity. MICA/B gene polymorphisms, expression levels and the amount of soluble MICA/B in the serum have been linked to autoimmune diseases, infections, and cancer. In hematopoietic stem cell transplantation, MICA matching between donor and patient has been correlated with reduced acute and chronic graft-vs.-host disease and improved survival. Hence, we developed an extremely cost-efficient high-throughput workflow for genotyping MICA/B for newly registered potential stem cell donors. Since mid-2017, we have genotyped over two million samples using NGS amplicon sequencing for MICA/B exons 2–5. In donors of German origin, MICA*008 is the most common MICA allele with a frequency of 42.3%. It is followed by MICA*002 (11.7%) and MICA*009 (8.8%). The three most common MICB alleles are MICB*005 (43.9%), MICB*004 (21.7%), and MICB*002 (18.9%). In general, MICB is less diverse than MICA and only 6 alleles, instead of 15, account for a cumulative allele frequency of 99.5%. In 0.5% of the samples we observed at least one allele of MICA or MICB which has so far not been reported to the IPD/IMGT-HLA database. By providing MICA/B typed voluntary donors, clinicians become empowered to include MICA/B into their donor selection process to further improve unrelated hematopoietic stem cell transplantation.

Immunogenetics in stem cell donor registry work: The DKMS example (Part 1)

Currently, stem cell donor registries include more than 35 million potential donors worldwide to provide HLA-matched stem cell products for patients in need of an unrelated donor transplant. DKMS is a leading stem cell donor registry with more than 9 million donors from Germany, Poland, the United States, the United Kingdom, India and Chile. DKMS donors have donated hematopoietic stem cells more than 80,000 times. Many aspects of donor registry work are closely related to topics from immunogenetics or population genetics. In this two-part review article, we describe, analyse and discuss these areas of donor registry work by using the example of DKMS. Part 1 of the review gives a general overview on DKMS and includes typical donor registry activities with special focus on the HLA system: high-throughput HLA typing of potential stem cell donors, HLA haplotype frequencies and resulting matching probabilities, and donor file optimization with regard to HLA diversity.

Apparent Lack of BRAF V600E Derived HLA Class I Presented Neoantigens Hampers Neoplastic Cell Targeting by CD8+ T Cells in Langerhans Cell Histiocytosis

Langerhans Cell Histiocytosis (LCH) is a neoplastic disorder of hematopoietic origin characterized by inflammatory lesions containing clonal histiocytes (LCH-cells) intermixed with various immune cells, including T cells. In 50-60% of LCH-patients, the somatic BRAF ^V600E driver mutation, which is common in many cancers, is detected in these LCH-cells in an otherwise quiet genomic landscape. Non-synonymous mutations like BRAF ^V600E can be a source of neoantigens capable of eliciting effective antitumor CD8⁺ T cell responses. This requires neopeptides to be stably presented by Human Leukocyte Antigen (HLA) class I molecules and sufficient numbers of CD8⁺ T cells at tumor sites. Here, we demonstrate substantial heterogeneity in CD8⁺ T cell density in n = 101 LCH-lesions, with BRAF ^V600E mutated lesions displaying significantly lower CD8⁺ T cell:CD1a⁺ LCH-cell ratios (p = 0.01) than BRAF wildtype lesions. Because LCH-lesional CD8⁺ T cell density had no significant impact on event-free survival, we investigated whether the intracellularly expressed BRAF ^V600E protein is degraded into neopeptides that are naturally processed and presented by cell surface HLA class I molecules. Epitope prediction tools revealed a single HLA class I binding BRAF ^V600E derived neopeptide (KIGDFGLATEK), which indeed displayed strong to intermediate binding capacity to HLA-A^*03:01 and HLA-A^*11:01 in an in vitro peptide-HLA binding assay. Mass spectrometry-based targeted peptidomics was used to investigate the presence of this neopeptide in HLA class I presented peptides isolated from several BRAF ^V600E expressing cell lines with various HLA genotypes. While the HLA-A^*02:01 binding BRAF wildtype peptide KIGDFGLATV was traced in peptides isolated from all five cell lines expressing this HLA subtype, KIGDFGLATEK was not detected in the HLA class I peptidomes of two distinct BRAF ^V600E transduced cell lines with confirmed expression of HLA-A^*03:01 or HLA-A^*11:01. These data indicate that the in silico predicted HLA class I binding and proteasome-generated neopeptides derived from the BRAF ^V600E protein are not presented by HLA class I molecules. Given that the BRAF ^V600E mutation is highly prevalent in chemotherapy refractory LCH-patients who may qualify for immunotherapy, this study therefore questions the efficacy of immune checkpoint inhibitor therapy in LCH.

Next-Generation Sequencing Technologies in Blood Group Typing

Sequencing of the human genome has led to the definition of the genes for most of the relevant blood group systems, and the polymorphisms responsible for most of the clinically relevant blood group antigens are characterized. Molecular blood group typing is used in situations where erythrocytes are not available or where serological testing was inconclusive or not possible due to the lack of antisera. Also, molecular testing may be more cost-effective in certain situations. Molecular typing approaches are mostly based on either PCR with specific primers, DNA hybridization, or DNA sequencing. Particularly the transition of sequencing techniques from Sanger-based sequencing to next-generation sequencing (NGS) technologies has led to exciting new possibilities in blood group genotyping. We describe briefly the currently available NGS platforms and their specifications, depict the genetic background of blood group polymorphisms, and discuss applications for NGS approaches in immunohematology. As an example, we delineate a protocol for large-scale donor blood group screening established and in use at our institution. Furthermore, we discuss technical challenges and limitations as well as the prospect for future developments, including long-read sequencing technologies.