Allele-Level KIR Genotyping of More Than a Million Samples: Workflow, Algorithm, and Observations

Donor KIR genotype based outcome prediction after allogeneic stem cell transplantation: no land in sight

Optimizing natural killer (NK) cell alloreactivity could further improve outcome after allogeneic hematopoietic cell transplantation (alloHCT). The donor's Killer-cell Immunoglobulin-like Receptor (KIR) genotype may provide important information in this regard. In the past decade, different models have been proposed aiming at maximizing NK cell activation by activating KIR-ligand interactions or minimizing inhibitory KIR-ligand interactions. Alternative classifications intended predicting outcome after alloHCT by donor KIR-haplotypes. In the present study, we aimed at validating proposed models and exploring more classification approaches. To this end, we analyzed samples stored at the Collaborative Biobank from HLA-compatible unrelated stem cell donors who had donated for patients with acute myeloid leukemia (AML) or myelodysplastic neoplasm (MDS) and whose outcome data had been reported to EBMT or CIBMTR. The donor KIR genotype was determined by high resolution amplicon-based next generation sequencing. We analyzed data from 5,017 transplants. The median patient age at alloHCT was 56 years. Patients were transplanted for AML between 2013 and 2018. Donor-recipient pairs were matched for HLA-A, -B, -C, -DRB1, and -DQB1 (79%) or had single HLA mismatches. Myeloablative conditioning was given to 56% of patients. Fifty-two percent of patients received anti-thymocyte-globulin-based graft-versus-host disease prophylaxis, 32% calcineurin-inhibitor-based prophylaxis, and 7% post-transplant cyclophosphamide-based prophylaxis. We tested several previously reported classifications in multivariable regression analyses but could not confirm outcome associations. Exploratory analyses in 1,939 patients (39%) who were transplanted from donors with homozygous centromeric (cen) or telomeric (tel) A or B motifs, showed that the donor cen B/B-tel A/A diplotype was associated with a trend to better event-free survival (HR 0.84, p=.08) and reduced risk of non-relapse mortality (NRM) (HR 0.65, p=.01). When we further dissected the contribution of B subtypes, we found that only the cen B01/B01-telA/A diplotype was associated with a reduced risk of relapse (HR 0.40, p=.04) while all subtype combinations contributed to a reduced risk of NRM. This exploratory finding has to be validated in an independent data set. In summary, the existing body of evidence is not (yet) consistent enough to recommend use of donor KIR genotype information for donor selection in routine clinical practice.

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.

Haplotype reconstruction for genetically complex regions with ambiguous genotype calls: Illustration by the KIR gene region

Advances in DNA sequencing technologies have enabled genotyping of complex genetic regions exhibiting copy number variation and high allelic diversity, yet it is impossible to derive exact genotypes in all cases, often resulting in ambiguous genotype calls, that is, partially missing data. An example of such a gene region is the killer-cell immunoglobulin-like receptor (KIR) genes. These genes are of special interest in the context of allogeneic hematopoietic stem cell transplantation. For such complex gene regions, current haplotype reconstruction methods are not feasible as they cannot cope with the complexity of the data. We present an expectation-maximization (EM)-algorithm to estimate haplotype frequencies (HTFs) which deals with the missing data components, and takes into account linkage disequilibrium (LD) between genes. To cope with the exponential increase in the number of haplotypes as genes are added, we add three components to a standard EM-algorithm implementation. First, reconstruction is performed iteratively, adding one gene at a time. Second, after each step, haplotypes with frequencies below a threshold are collapsed in a rare haplotype group. Third, the HTF of the rare haplotype group is profiled in subsequent iterations to improve estimates. A simulation study evaluates the effect of combining information of multiple genes on the estimates of these frequencies. We show that estimated HTFs are approximately unbiased. Our simulation study shows that the EM-algorithm is able to combine information from multiple genes when LD is high, whereas increased ambiguity levels increase bias. Linear regression models based on this EM, show that a large number of haplotypes can be problematic for unbiased effect size estimation and that models need to be sparse. In a real data analysis of KIR genotypes, we compare HTFs to those obtained in an independent study. Our new EM-algorithm-based method is the first to account for the full genetic architecture of complex gene regions, such as the KIR gene region. This algorithm can handle the numerous observed ambiguities, and allows for the collapsing of haplotypes to perform implicit dimension reduction. Combining information from multiple genes improves haplotype reconstruction.

High-resolution characterization of KIR genes polymorphism in healthy subjects from the Bulgarian population-A pilot study

Although killer-cell immunoglobulin-like receptor (KIR) gene content has been widely studied in health and disease, with the advancement of next-generation sequencing (NGS) technology the high-resolution characterization of this complex gene region has become achievable. KIR allele-level diversity has lately been described across human populations. The present study aimed to analyze for the first time the allele-level polymorphism of nine KIR genes in 155 healthy, unrelated individuals from the Bulgarian population by applying NGS. The highest degree of polymorphism was detected for the KIR3DL3 gene with 40 observed alleles at five-digit resolution in total, 22 of which were common. On the other hand, the KIR3DS1 gene was found to have the lowest degree of polymorphism among the studied KIR genes with one common allele: KIR3DS1*01301 (31.6%). To better understand KIR allelic associations and patterns in Bulgarians, we have estimated the pairwise linkage disequilibrium (LD) for the 10 KIR loci, where KIR2DL3*00501 allele was found in strong LD with KIR2DL1*00101 (D' = 1.00, R2  = 0.742). This is the first study investigating KIR polymorphism at the allele level in a population from the South-East European region. Considering the effect of the populationally shaped KIR allelic polymorphism on NK cell function, this data could lead to a better understanding of the genetic heterogeneity of this region and can be carried into clinical practice by improvement of the strategies taken for NK-mediated diseases.

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.

Killer Cell Immunoglobulin-like Receptor Genotypes and Reproductive Outcomes in a Group of Infertile Women: A Romanian Study

A growing body of evidence suggests that endometrial immune disorders may be responsible for endometrial dysfunctions that can lead to gynecological and obstetrical pathology. The aim of this study was to explore the potential relationship between different killer cell immunoglobulin-like receptor (KIR) genotypes and reproductive outcomes. We conducted a prospective cohort study that included 104 infertile patients undergoing an in vitro fertilization procedure. All participants underwent clinical and ultrasound examination, genetic evaluation (KIR genotyping), endometrial washing fluid sampling for cytokine determination, endometrial tissue sampling for histologic assessment and hysteroscopic evaluation. Our analysis showed statistically significant lower levels of uterine cytokines TNF-α (p = 0.001) and IL-1beta (p = 0.000) in the KIR AA genotype group as compared to KIR AB and BB among study participants with chronic endometritis. The study results suggest that the KIR AA genotype population subgroups may be more susceptible to developing endometrial disorders such as chronic endometritis. The changes in the behavior of NK cells seem to be subtle and expressed as an altered regulatory pattern.

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.

KIR genotype and haplotype frequencies in the multi-ethnic population of Malaysia

Killer cell immunoglobulin-like receptors (KIR) genotype and haplotype frequencies have been reported to vary distinctly between populations, which in turn contributes to variation in the alloreactivity of natural killer (NK) cells. Utilizing the diverse KIR genes to identify suitable transplant donors would prove challenging in multi-ethnic countries, even more in resource-limited countries where KIR genotyping has not been established. In this study, we determined the KIR genotypes from 124 unrelated Malaysians consisting of the Malays, Chinese, Indians, and aboriginal people through polymerase chain reaction sequence-specific primer (PCR-SSP) genotyping and employing an expectation-maximization (EM) algorithm to assign haplotypes based on pre-established reference haplotypes. A total of 27 distinct KIR haplotypes were discerned with higher frequencies of haplotype A (55.2%) than haplotype B (44.8%). The most frequent haplotypes were cA01:tA01 (55.2%), cB01:tB01 (18.1%), and cB02:tA01 (13.3%), while the least frequent haplotypes were cB03:tB01 (1.2%), cB04:tB03 (0.4%), and cB03:tA01 (0.4%). Several haplotypes were identified to be unique to a specific ethnic group. The genotype with the highest frequency was genotype AB (71.8%), followed by AA (19.4%), and BB (8.9%). The Indians exhibited the lowest genotype AA but the highest genotype BB, whereas genotype BB was absent in the aboriginal people. Despite the limitations, the genotype and haplotypes in the Malaysian population were successfully highlighted. The identification of ethnic-specific KIR genotypes and haplotypes provides the first step to utilizing KIR in identifying suitable transplant donors to further improve the transplant outcome in the Malaysian population.

Comparison of NK alloreactivity prediction models based on KIR-MHC interactions in haematopoeitic stem cell transplantation

The biological processes underlying NK cell alloreactivity in haematopoietic stem cell transplantation (HSCT) remain unclear. Many different models to predict NK alloreactivity through KIR and MHC genotyping exist, raising ambiguities in its utility and application for clinicians. We assessed 27 predictive models, broadly divided into six categories of alloreactivity prediction: ligand-ligand, receptor-ligand, educational, KIR haplotype-based, KIR matching and KIR allelic polymorphism. The models were applied to 78 NGS-typed donor/recipient pairs undergoing allogeneic HSCT in genoidentical (n=43) or haploidentical (n=35) matchings. Correlations between different predictive models differed widely, suggesting that the choice of the model in predicting NK alloreactivity matters. For example, two broadly used models, educational and receptor-ligand, led to opposing predictions especially in the genoidentical cohort. Correlations also depended on the matching fashion, suggesting that this parameter should also be taken into account in the choice of the scoring strategy. The number of centromeric B-motifs was the only model strongly correlated with the incidence of acute graft-versus-host disease in our set of patients in both the genoidentical and the haploidentical cohorts, suggesting that KIR-based alloreactivity, not MHC mismatches, are responsible for it. To our best knowledge, this paper is the first to experimentally compare NK alloreactivity prediction models within a cohort of genoidentical and haploidentical donor-recipient pairs. This study helps to resolve current discrepancies in KIR-based alloreactivity predictions and highlights the need for deeper consideration of the models used in clinical studies as well as in medical practice.

Analysis of KIR gene variants in The Cancer Genome Atlas and UK Biobank using KIRCLE

BACKGROUND

Natural killer (NK) cells represent a critical component of the innate immune system's response against cancer and viral infections, among other diseases. To distinguish healthy host cells from infected or tumor cells, killer immunoglobulin receptors (KIR) on NK cells bind and recognize Human Leukocyte Antigen (HLA) complexes on their target cells. However, NK cells exhibit great diversity in their mechanism of activation, and the outcomes of their activation are not yet understood fully. Just like the HLAs they bind, KIR receptors exhibit high allelic diversity in the human population. Here we provide a method to identify KIR allele variants from whole exome sequencing data and uncover novel associations between these variants and various molecular and clinical correlates.

RESULTS

In order to better understand KIRs, we have developed KIRCLE, a novel method for genotyping individual KIR genes from whole exome sequencing data, and used it to analyze approximately sixty-thousand patient samples in The Cancer Genome Atlas (TCGA) and UK Biobank. We were able to assess population frequencies for different KIR alleles and demonstrate that, similar to HLA alleles, individuals' KIR alleles correlate strongly with their ethnicities. In addition, we observed associations between different KIR alleles and HLA alleles, including HLA-B*53 with KIR3DL2*013 (Fisher's exact FDR = 7.64e-51). Finally, we showcased statistically significant associations between KIR alleles and various clinical correlates, including peptic ulcer disease (Fisher's exact FDR = 0.0429) and age of onset of atopy (Mann-Whitney U FDR = 0.0751).

CONCLUSIONS

We show that KIRCLE is able to infer KIR variants accurately and consistently, and we demonstrate its utility using data from approximately sixty-thousand individuals from TCGA and UK Biobank to discover novel molecular and clinical correlations with KIR germline variants. Peptic ulcer disease and atopy are just two diseases in which NK cells may play a role beyond their "classical" realm of anti-tumor and anti-viral responses. This tool may be used both as a benchmark for future KIR-variant-inference algorithms, and to better understand the immunogenomics of and disease processes involving KIRs.

Skeletal muscle provides the immunological micro-milieu for specific plasma cells in anti-synthetase syndrome-associated myositis

Anti-synthetase syndrome (ASyS)-associated myositis is a major subgroup of the idiopathic inflammatory myopathies (IIM) and is characterized by disease chronicity with musculoskeletal, dermatological and pulmonary manifestations. One of eight autoantibodies against the aminoacyl-transferase RNA synthetases (ARS) is detectable in the serum of affected patients. However, disease-specific therapeutic approaches have not yet been established.To obtain a deeper understanding of the underlying pathogenesis and to identify putative therapeutic targets, we comparatively investigated the most common forms of ASyS associated with anti-PL-7, anti-PL-12 and anti-Jo-1. Our cohort consisted of 80 ASyS patients as well as healthy controls (n = 40), diseased controls (n = 40) and non-diseased controls (n = 20). We detected a reduced extent of necrosis and regeneration in muscle biopsies from PL-12⁺ patients compared to Jo-1⁺ patients, while PL-7⁺ patients had higher capillary dropout in biopsies of skeletal muscle. Aside from these subtle alterations, no significant differences between ASyS subgroups were observed. Interestingly, a tissue-specific subpopulation of CD138⁺ plasma cells and CXCL12⁺/CXCL13⁺CD20⁺ B cells common to ASyS myositis were identified. These cells were localized in the endomysium associated with alkaline phosphatase⁺ activated mesenchymal fibroblasts and CD68⁺MHC-II⁺CD169⁺ macrophages. An MHC-I⁺ and MHC-II⁺ MxA negative type II interferon-driven milieu of myofiber activation, topographically restricted to the perifascicular area and the adjacent perimysium, as well as perimysial clusters of T follicular helper cells defined an extra-medullary immunological niche for plasma cells and activated B cells. Consistent with this, proteomic analyses of muscle tissues from ASyS patients demonstrated alterations in antigen processing and presentation. In-depth immunological analyses of peripheral blood supported a B-cell/plasma-cell-driven pathology with a shift towards immature B cells, an increase of B-cell-related cytokines and chemokines, and activation of the complement system. We hypothesize that a B-cell-driven pathology with the presence and persistence of a specific subtype of plasma cells in the skeletal muscle is crucially involved in the self-perpetuating chronicity of ASyS myositis. This work provides the conceptual framework for the application of plasma-cell-targeting therapies in ASyS myositis.

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.

Immunogenomics of Killer Cell Immunoglobulin-Like Receptor (KIR) and HLA Class I: Coevolution and Consequences for Human Health

Interactions of killer cell immunoglobin-like receptors (KIR) with human leukocyte antigens (HLA) class I regulate effector functions of key cytotoxic cells of innate and adaptive immunity. The extreme diversity of this interaction is genetically determined, having evolved in the ever-changing environment of pathogen exposure. Diversity of KIR and HLA genes is further facilitated by their independent segregation on separate chromosomes. That fetal implantation relies on many of the same types of immune cells as infection control places certain constraints on the evolution of KIR interactions with HLA. Consequently, specific inherited combinations of receptors and ligands may predispose to specific immune-mediated diseases, including autoimmunity. Combinatorial diversity of KIR and HLA class I can also differentiate success rates of immunotherapy directed to these diseases. Progress toward both etiopathology and predicting response to therapy is being achieved through detailed characterization of the extent and consequences of the combinatorial diversity of KIR and HLA. Achieving these goals is more tractable with the development of integrated analyses of molecular evolution, function, and pathology that will establish guidelines for understanding and managing risks. Here, we present what is known about the coevolution of KIR with HLA class I and the impact of their complexity on immune function and homeostasis.

Immune determinants of endometrial receptivity: a biological perspective

Immune cells are essential for endometrial receptivity to embryo implantation and early placental development. They exert tissue-remodeling and immune regulatory roles-acting to promote epithelial attachment competence, regulate the differentiation of decidual cells, remodel the uterine vasculature, control and resolve inflammatory activation, and suppress destructive immunity to paternally inherited alloantigens. From a biological perspective, the endometrial immune response exerts a form of "quality control"-it promotes implantation success when conditions are favorable but constrains receptivity when physiological circumstances are not ideal. Women with recurrent implantation failure and recurrent miscarriage may exhibit altered numbers or disturbed function of certain uterine immune cell populations-most notably uterine natural killer cells and regulatory T cells. Preclinical and animal studies indicate that deficiencies or aberrant activation states in these cells can be causal in the pathophysiological mechanisms of infertility. Immune cells are, therefore, targets for diagnostic evaluation and therapeutic intervention. However, current diagnostic tests are overly simplistic and have limited clinical utility. To be more informative, they need to account for the full complexity and reflect the range of perturbations that can occur in uterine immune cell phenotypes and networks. Moreover, safe and effective interventions to modulate these cells are in their infancy, and personalized approaches matched to specific diagnostic criteria will be needed. Here we summarize current biological understanding and identify knowledge gaps to be resolved before the promise of therapies to target the uterine immune response can be fully realized.

Allelic Polymorphisms of Killer Immunoglobulin-Like Receptor Genes in Malay and Orang Asli Populations of Peninsular Malaysia

Next-generation DNA sequencing (NGS) technology advancements provide new insight into the level of variation in killer immunoglobulin-like receptor (KIR) genes. High resolution allele genotyping of seven KIR genes was conducted among 94 unrelated Malay and Orang Asli (OA) individuals of Peninsular Malaysia. A manual bioinformatics analysis is performed and optimised by Sanger sequencing method. The Malays expressed a total of 22 alleles, as compared to only 15 alleles in the OA population. In total, 12 centromeric and 9 telomeric allelic haplotypes were identified in the Malays, whereas 8 centromeric and 5 telomeric allelic haplotypes were identified in the OA. The KIR2DL1, KIR2DL3, and KIR2DS4 genes exhibited a high degree of variation and balanced distribution in the Malay and OA populations. On the other hand, KIR2DL4, KIR3DL1, KIR3DL2 and KIR3DL3 genes exhibited a high degree of conservation, with less number of alleles identified and the dominance of a single allele at high frequency. High-resolution KIR allele genotyping has revealed unique sequence variations and allelic haplotypes between individuals and populations. The distributions of KIR alleles and haplotypes are useful for genetic population studies and serve as a baseline for future transplantation matching and disease association research.

Donor–Recipient Non-HLA Variants, Mismatches and Renal Allograft Outcomes: Evolving Paradigms

Despite significant improvement in the rates of acute allograft rejection, proportionate improvements in kidney allograft longevity have not been realized, and are a source of intense research efforts. Emerging translational data and natural history studies suggest a role for anti-donor immune mechanisms in a majority of cases of allograft loss without patient death, even when overt evidence of acute rejection is not identified. At the level of the donor and recipient genome, differences in highly polymorphic HLA genes are routinely evaluated between donor and recipient pairs as part of organ allocation process, and utilized for patient-tailored induction and maintenance immunosuppression. However, a growing body of data have characterized specific variants in donor and recipient genes, outside of HLA loci, that induce phenotypic changes in donor organs or the recipient immune system, impacting transplant outcomes. Newer mechanisms for “mismatches” in these non-HLA loci have also been proposed during donor–recipient genome interactions with transplantation. Here, we review important recent data evaluating the role of non-HLA genetic loci and genome-wide donor-recipient mismatches in kidney allograft outcomes.

Allele imputation for the killer cell immunoglobulin-like receptor KIR3DL1/S1

Highly polymorphic interaction of KIR3DL1 and KIR3DS1 with HLA class I ligands modulates the effector functions of natural killer (NK) cells and some T cells. This genetically determined diversity affects severity of infections, immune-mediated diseases, and some cancers, and impacts the course of immunotherapies, including transplantation. KIR3DL1 is an inhibitory receptor, and KIR3DS1 is an activating receptor encoded by the KIR3DL1/S1 gene that has more than 200 diverse and divergent alleles. Determination of KIR3DL1/S1 genotypes for medical application is hampered by complex sequence and structural variation, requiring targeted approaches to generate and analyze high-resolution allele data. To overcome these obstacles, we developed and optimized a model for imputing KIR3DL1/S1 alleles at high-resolution from whole-genome SNP data. We designed the model to represent a substantial component of human genetic diversity. Our Global imputation model is effective at genotyping KIR3DL1/S1 alleles with an accuracy ranging from 88% in Africans to 97% in East Asians, with mean specificity of 99% and sensitivity of 95% for alleles >1% frequency. We used the established algorithm of the HIBAG program, in a modification named Pulling Out Natural killer cell Genomics (PONG). Because HIBAG was designed to impute HLA alleles also from whole-genome SNP data, PONG allows combinatorial diversity of KIR3DL1/S1 with HLA-A and -B to be analyzed using complementary techniques on a single data source. The use of PONG thus negates the need for targeted sequencing data in very large-scale association studies where such methods might not be tractable.

High-resolution human KIR genotyping

Killer immunoglobulin-like receptors (KIR) regulate the function of natural killer cells through interactions with various ligands on the surface of cells, thereby determining whether natural killer (NK) cells are to be activated or inhibited from killing the cell being interrogated. The genes encoding these proteins display extensive variation through variable gene content, copy number and allele polymorphism. The combination of KIR genes and their ligands is implicated in various clinical settings including haematopoietic stem cell and solid organ transplant and infectious disease progression. The determination of KIR genes has been used as a factor in the selection of optimal stem cell donors with haplotype variations in recipient and donor giving differential clinical outcomes. Methods to determine KIR genes have primarily involved ascertaining the presence or absence of genes in an individual. With the more recent introduction of massively parallel clonal next-generation sequencing and single molecule very long read length third-generation sequencing, high-resolution determination of KIR alleles has become feasible. Determining the extent and functional impact of allele variation has the potential to lead to further optimisation of clinical outcomes as well as a deeper understanding of the functional properties of the receptors and their interactions with ligands. This review summarizes recently published high-resolution KIR genotyping methods and considers the various advantages and disadvantages of the approaches taken. In addition the application of allele level genotyping in the setting of transplantation and infectious disease control is discussed.

KIR gene content imputation from single-nucleotide polymorphisms in the Finnish population

The killer cell immunoglobulin-like receptor (KIR) gene cluster on chromosome 19 encodes cell surface glycoproteins that bind class I human leukocyte antigen (HLA) molecules as well as some other ligands. Through regulation of natural killer (NK) cell activity KIRs participate in tumour surveillance and clearing viral infections. KIR gene gene copy number variation associates with the outcome of transplantations and susceptibility to immune-mediated diseases. Inferring KIR gene content from genetic variant data is therefore desirable for immunogenetic analysis, particularly in the context of growing biobank genome data collections that rely on genotyping by microarray. Here we describe a stand-alone and freely available gene content imputation for 12 KIR genes. The models were trained using 807 Finnish biobank samples genotyped for 5900 KIR-region SNPs and analysed for KIR gene content with targeted sequencing. Cross-validation results demonstrate a high mean overall accuracy of 98.5% (95% CI [97.0-99.2]%) which compares favourably with previous methods including short-read sequencing based approaches.

High-throughput Interpretation of Killer-cell Immunoglobulin-like Receptor Short-read Sequencing Data with PING

The killer-cell immunoglobulin-like receptor (KIR) complex on chromosome 19 encodes receptors that modulate the activity of natural killer cells, and variation in these genes has been linked to infectious and autoimmune disease, as well as having bearing on pregnancy and transplant outcomes. The medical relevance and high variability of KIR genes makes short-read sequencing an attractive technology for interrogating the region, providing a high-throughput, high-fidelity sequencing method that is cost-effective. However, because this gene complex is characterized by extensive nucleotide polymorphism, structural variation including gene fusions and deletions, and a high level of homology between genes, its interrogation at high resolution has been thwarted by bioinformatic challenges, with most studies limited to examining presence or absence of specific genes. Here, we present the PING (Pushing Immunogenetics to the Next Generation) pipeline, which incorporates empirical data, novel alignment strategies and a custom alignment processing workflow to enable high-throughput KIR sequence analysis from short-read data. PING provides KIR gene copy number classification functionality for all KIR genes through use of a comprehensive alignment reference. The gene copy number determined per individual enables an innovative genotype determination workflow using genotype-matched references. Together, these methods address the challenges imposed by the structural complexity and overall homology of the KIR complex. To determine copy number and genotype determination accuracy, we applied PING to European and African validation cohorts and a synthetic dataset. PING demonstrated exceptional copy number determination performance across all datasets and robust genotype determination performance. Finally, an investigation into discordant genotypes for the synthetic dataset provides insight into misaligned reads, advancing our understanding in interpretation of short-read sequencing data in complex genomic regions. PING promises to support a new era of studies of KIR polymorphism, delivering high-resolution KIR genotypes that are highly accurate, enabling high-quality, high-throughput KIR genotyping for disease and population studies.

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.

HLA Class I Molecules as Immune Checkpoints for NK Cell Alloreactivity and Anti-Viral Immunity in Kidney Transplantation

Natural killer (NK) cells are innate lymphocytes that can kill diseased- or virally-infected cells, mediate antibody dependent cytotoxicity and produce type I immune-associated cytokines upon activation. NK cells also contribute to the allo-immune response upon kidney transplantation either by promoting allograft rejection through lysis of cells of the transplanted organ or by promoting alloreactive T cells. In addition, they protect against viral infections upon transplantation which may be especially relevant in patients receiving high dose immune suppression. NK cell activation is tightly regulated through the integrated balance of signaling via inhibitory- and activating receptors. HLA class I molecules are critical regulators of NK cell activation through the interaction with inhibitory- as well as activating NK cell receptors, hence, HLA molecules act as critical immune checkpoints for NK cells. In the current review, we evaluate how NK cell alloreactivity and anti-viral immunity are regulated by NK cell receptors belonging to the KIR family and interacting with classical HLA class I molecules, or by NKG2A/C and LILRB1/KIR2DL4 engaging non-classical HLA-E or -G. In addition, we provide an overview of the methods to determine genetic variation in these receptors and their HLA ligands.

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.

High-Resolution Characterization of KIR Genes in a Large North American Cohort Reveals Novel Details of Structural and Sequence Diversity

The KIR (killer-cell immunoglobulin-like receptor) region is characterized by structural variation and high sequence similarity among genes, imposing technical difficulties for analysis. We undertook the most comprehensive study to date of KIR genetic diversity in a large population sample, applying next-generation sequencing in 2,130 United States European-descendant individuals. Data were analyzed using our custom bioinformatics pipeline specifically designed to address technical obstacles in determining KIR genotypes. Precise gene copy number determination allowed us to identify a set of uncommon gene-content KIR haplotypes accounting for 5.2% of structural variation. In this cohort, KIR2DL4 is the framework gene that most varies in copy number (6.5% of all individuals). We identified phased high-resolution alleles in large multi-locus insertions and also likely founder haplotypes from which they were deleted. Additionally, we observed 250 alleles at 5-digit resolution, of which 90 have frequencies ≥1%. We found sequence patterns that were consistent with the presence of novel alleles in 398 (18.7%) individuals and contextualized multiple orphan dbSNPs within the KIR complex. We also identified a novel KIR2DL1 variant, Pro151Arg, and demonstrated by molecular dynamics that this substitution is predicted to affect interaction with HLA-C. No previous studies have fully explored the full range of structural and sequence variation of KIR as we present here. We demonstrate that pairing high-throughput sequencing with state-of-art computational tools in a large cohort permits exploration of all aspects of KIR variation including determination of population-level haplotype diversity, improving understanding of the KIR system, and providing an important reference for future studies.

High-Resolution Analysis Identifies High Frequency of KIR-A Haplotypes and Inhibitory Interactions of KIR With HLA Class I in Zhejiang Han

Killer cell immunoglobulin-like receptors (KIR) interact with human leukocyte antigen (HLA) class I molecules, modulating critical NK cell functions in the maintenance of human health. Characterizing the distribution and characteristics of KIR and HLA allotype diversity across defined human populations is thus essential for understanding the multiple associations with disease, and for directing therapies. In this study of 176 Zhejiang Han individuals from Southeastern China, we describe diversity of the highly polymorphic KIR and HLA class I genes at high resolution. KIR-A haplotypes, which carry four inhibitory receptors specific for HLA-A, B or C, are known to associate with protection from infection and some cancers. We show the Chinese Southern Han from Zhejiang are characterized by a high frequency of KIR-A haplotypes and a high frequency of C1 KIR ligands. Accordingly, interactions of inhibitory KIR2DL3 with C1+HLA are more frequent in Zhejiang Han than populations outside East Asia. Zhejiang Han exhibit greater diversity of inhibitory than activating KIR, with three-domain inhibitory KIR exhibiting the greatest degree of polymorphism. As distinguished by gene copy number and allele content, 54 centromeric and 37 telomeric haplotypes were observed. We observed 6% of the population to have KIR haplotypes containing large-scale duplications or deletions that include complete genes. A unique truncated haplotype containing only KIR2DL4 in the telomeric region was also identified. An additional feature is the high frequency of HLA-B*46:01, which may have arisen due to selection pressure from infectious disease. This study will provide further insight into the role of KIR and HLA polymorphism in disease susceptibility of Zhejiang Chinese.

KIR2DL4 genetic diversity in a Brazilian population sample: implications for transcription regulation and protein diversity in samples with different ancestry backgrounds

KIR2DL4 is an important immune modulator expressed in natural killer cells; HLA-G is its main ligand. We have characterized the KIR2DL4 genetic diversity by considering the promoter, all exons, and all introns in a highly admixed Brazilian population sample and by using massively parallel sequencing. We introduce a molecular method to amplify and to sequence the complete KIR2DL4 gene. To avoid the mapping bias and genotype errors commonly observed in gene families, we have developed and validated a bioinformatic pipeline designed to minimize these errors and applied it to survey the variability of 220 individuals from the State of São Paulo, southeastern Brazil. We have also compared the KIR2DL4 genetic diversity in the Brazilian cohort with the diversity previously reported by the 1000Genomes consortium. KIR2DL4 presents high linkage disequilibrium throughout the gene, with coding sequences associated with specific promoters. There are few but divergent promoter haplotypes. We have also detected many new KIR2DL4 sequences, all bearing nucleotide exchanges in introns and encoding previously described proteins. Exons 3 and 4, which encode the external domains, are the most variable. The ancestry background influences the KIR2DL4 allele frequencies and must be considered for association studies regarding KIR2DL4.

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.

Immunogenicity Challenges Associated with Subcutaneous Delivery of Therapeutic Proteins

The subcutaneous route of administration has provided convenient and non-inferior delivery of therapeutic proteins compared to intravenous infusion, but there is potential for enhanced immunogenicity toward subcutaneously administered proteins in a subset of patients. Unwanted anti-drug antibody response toward proteins or monoclonal antibodies upon repeated administration is shown to impact the pharmacokinetics and efficacy of multiple biologics. Unique immunogenicity challenges of the subcutaneous route have been realized through various preclinical and clinical examples, although subcutaneous delivery has often demonstrated comparable immunogenicity to intravenous administration. Beyond route of administration as a treatment-related factor of immunogenicity, certain product-related risk factors are particularly relevant to subcutaneously administered proteins. This review attempts to provide an overview of the mechanism of immune response toward proteins administered subcutaneously (subcutaneous proteins) and comments on product-related risk factors related to protein structure and stability, dosage form, and aggregation. A two-wave mechanism of antigen presentation in the immune response toward subcutaneous proteins is described, and interaction with dynamic antigen-presenting cells possessing high antigen processing efficiency and migratory activity may drive immunogenicity. Mitigation strategies for immunogenicity are discussed, including those in general use clinically and those currently in development. Mechanistic insights along with consideration of risk factors involved inspire theoretical strategies to provide antigen-specific, long-lasting effects for maintaining the safety and efficacy of therapeutic proteins.

Haplotype Motif-Based Models for KIR-Genotype Informed Selection of Hematopoietic Cell Donors Fail to Predict Outcome of Patients With Myelodysplastic Syndromes or Secondary Acute Myeloid Leukemia

Results from registry studies suggest that harnessing Natural Killer (NK) cell reactivity mediated through Killer cell Immunoglobulin-like Receptors (KIR) could reduce the risk of relapse after allogeneic Hematopoietic Cell Transplantation (HCT). Several competing models have been developed to classify donors as KIR-advantageous or disadvantageous. Basically, these models differ by grouping donors based on distinct KIR–KIR–ligand combinations or by haplotype motif assignment. This study aimed to validate different models for unrelated donor selection for patients with Myelodysplatic Syndromes (MDS) or secondary Acute Myeloid Leukemia (sAML). In a joint retrospective study of the European Society for Blood and Marrow Transplantation (EBMT) and the Center for International Blood and Marrow Transplant Research (CIBMTR) registry data from 1704 patients with secondary AML or MDS were analysed. The cohort consisted mainly of older patients (median age 61 years) with high risk disease who had received chemotherapy-based reduced intensity conditioning and anti-thymocyte globulin prior to allogeneic HCT from well-matched unrelated stem cell donors. The impact of the predictors on Overall Survival (OS) and relapse incidence was tested in Cox regression models adjusted for patient age, a modified disease risk index, performance status, donor age, HLA-match, sex-match, CMV-match, conditioning intensity, type of T-cell depletion and graft type. KIR genes were typed using high-resolution amplicon-based next generation sequencing. In univariable and multivariable analyses none of the models predicted OS and the risk of relapse consistently. Our results do not support the hypothesis that optimizing NK-mediated alloreactivity is possible by KIR-genotype informed selection of HLA-matched unrelated donors. However, in the context of allogeneic transplantation, NK-cell biology is complex and only partly understood. KIR-genes are highly diverse and current assignment of haplotype motifs based on the presence or absence of selected KIR genes is over-simplistic. As a consequence, further research is highly warranted and should integrate cutting edge knowledge on KIR genetics, and NK-cell biology into future studies focused on homogeneous groups of patients and treatment modalities.

Centromeric KIR AA Individuals Harbor Particular KIR Alleles Conferring Beneficial NK Cell Features with Implications in Haplo-Identical Hematopoietic Stem Cell Transplantation

We have recently shown a broad disparity of Natural Killer (NK) cell responses against leukemia highlighting good and bad responders resting on the Killer cell Immunoglobulin-like Receptors (KIR) and HLA genetics. In this study, we deeply studied KIR2D allele expression, HLA-C recognition and functional effect on NK cells in 108 blood donors in combining high-resolution KIR allele typing and multicolor flow cytometry. The KIR2DL1*003 allotype is associated with centromeric (cen) AA motif and confers the highest NK cell frequency, expression level and strength of KIR/HLA-C interactions compared to the KIR2DL1*002 and KIR2DL1*004 allotypes respectively associated with cenAB and BB motifs. KIR2DL2*001 and *003 allotypes negatively affect the frequency of KIR2DL1⁺ and KIR2DL3⁺ NK cells. Altogether, our data suggest that cenAA individuals display more efficient KIR2DL alleles (L1*003 and L3*001) to mount a consistent frequency of KIR2DL⁺ NK cells and to confer an effective NK cell responsiveness. The transposition of our in vitro observations in the T-replete haplo-identical HSCT context led us to observe that cenAA HSC grafts limit significantly the incidence of relapse in patients with myeloid diseases after T-replete haplo-identical HSCT. As NK cells are crucial in HSCT reconstitution, one could expect that the consideration of KIR2DL1/2/3 allelic polymorphism could help to refine scores used for HSC donor selection.

Dynamic Changes of Inhibitory Killer-Immunoglobulin-Like Receptors on NK Cells after Allogeneic Hematopoietic Stem Cell Transplantation: An Initial Study

Killer-immunoglobulin-like receptors (KIRs) are critical natural killer (NK) cell regulators. The expression of KIRs is a dynamic process influenced by many factors. Their ligands-HLA(Human Leukocyte Antigen) class I molecules-are expressed on all nucleated cells that keep NK cells under control. In hematopoietic stem cell transplantation (HSCT), NK cells play an essential role in relapse protection. In the presented pilot study, we characterized the dynamic expression of inhibitory KIRS (iKIRs), which protect cells against untoward lysis, in donors and patients during the first three months after HSCT using flow cytometry. The expression of all iKIRs was highly variable and sometimes correlated with patients' clinical presentation and therapy regiment. Cyclophosphamide (Cy) in the graft-versus-host disease (GvHD) prevention protocol downregulated KIR2DL1 to just 25% of the original donor value, and the FEAM (Fludarabine + Etoposid + Ara-C + Melphalan) conditioning protocol reduced KIR2DL3. In lymphoid neoplasms, there was a slightly increased KIR2DL3 expression compared to myeloid malignancies. Additionally, we showed that the ex vivo activation of NK cells did not alter the level of iKIRs. Our study shows the influence of pre- and post-transplantation protocols on iKIR expression on the surface of NK cells and the importance of monitoring their cell surface.

External validation of models for KIR2DS1/KIR3DL1-informed selection of hematopoietic cell donors fails

Several studies suggest that harnessing natural killer (NK) cell reactivity mediated through killer cell immunoglobulin-like receptors (KIRs) could reduce the risk of relapse after allogeneic hematopoietic cell transplantation. Based on one promising model, information on KIR2DS1 and KIR3DL1 and their cognate ligands can be used to classify donors as KIR-advantageous or KIR-disadvantageous. This study was aimed at externally validating this model in unrelated donor hematopoietic cell transplantation. The impact of the predictor on overall survival (OS) and relapse incidence was tested in a Cox regression model adjusted for patient age, a modified disease risk index, Karnofsky performance status, donor age, HLA match, sex match, cytomegalovirus match, conditioning intensity, type of T-cell depletion, and graft type. Data from 2222 patients with acute myeloid leukemia or myelodysplastic syndrome were analyzed. KIR genes were typed by using high-resolution amplicon-based next-generation sequencing. In univariable analyses and subgroup analyses, OS and the cumulative incidence of relapse of patients with a KIR-advantageous donor were comparable to patients with a KIR-disadvantageous donor. The adjusted hazard ratio from the multivariable Cox regression model was 0.99 (Wald test, P = .93) for OS and 1.04 (Wald test, P = .78) for relapse incidence. We also tested the impact of activating donor KIR2DS1 and inhibition by KIR3DL1 separately but found no significant impact on OS and the risk of relapse. Thus, our study shows that the proposed model does not universally predict NK-mediated disease control. Deeper knowledge of NK-mediated alloreactivity is necessary to predict its contribution to graft-versus-leukemia reactions and to eventually use KIR genotype information for donor selection.

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.

In silico tools for accurate HLA and KIR inference from clinical sequencing data empower immunogenetics on individual-patient and population scales

Immunogenetic variation in humans is important in research, clinical diagnosis and increasingly a target for therapeutic intervention. Two highly polymorphic loci play critical roles, namely the human leukocyte antigen (HLA) system, which is the human version of the major histocompatibility complex (MHC), and the Killer-cell immunoglobulin-like receptors (KIR) that are relevant for responses of natural killer (NK) and some subsets of T cells. Their accurate classification has typically required the use of dedicated biological specimens and a combination of in vitro and in silico efforts. Increased availability of next generation sequencing data has led to the development of ancillary computational solutions. Here, we report an evaluation of recently published algorithms to computationally infer complex immunogenetic variation in the form of HLA alleles and KIR haplotypes from whole-genome or whole-exome sequencing data. For both HLA allele and KIR gene typing, we identified tools that yielded >97% overall accuracy for four-digit HLA types, and >99% overall accuracy for KIR gene presence, suggesting the readiness of in silico solutions for use in clinical and high-throughput research settings.

Downregulation of HLA-I by the molluscum contagiosum virus mc080 impacts NK-cell recognition and promotes CD8+ T-cell evasion

Molluscum contagiosum virus (MCV) is a common cause of benign skin lesions in young children and currently the only endemic human poxvirus. Following the infection of primary keratinocytes in the epidermis, MCV induces the proliferation of infected cells and this results in the production of wart-like growths. Full productive infection is observed only after the infected cells differentiate. During this prolonged replication cycle the virus must avoid elimination by the host immune system. We therefore sought to investigate the function of the two major histocompatibility complex class-I-related genes encoded by the MCV genes mc033 and mc080. Following insertion into a replication-deficient adenovirus vector, codon-optimized versions of mc033 and mc080 were expressed as endoglycosidase-sensitive glycoproteins that localized primarily in the endoplasmic reticulum. MC080, but not MC033, downregulated cell-surface expression of endogenous classical human leucocyte antigen (HLA) class I and non-classical HLA-E by a transporter associated with antigen processing (TAP)-independent mechanism. MC080 exhibited a capacity to inhibit or activate NK cells in autologous assays in a donor-specific manner. MC080 consistently inhibited antigen-specific T cells being activated by peptide-pulsed targets. We therefore propose that MC080 acts to promote evasion of HLA-I-restricted cytotoxic T cells.

SPECTRUM AND FREQUENCY OF NK CELL RECEPTOR GENES AMONG CYSTIC FIBROSIS PATIENTS

Aim – to establish and analyze the spectrum of KIR genes in people with a confirmed diagnosis of Cystic fibrosis (CF), homozygote of F508del mutation of the СFTR gene for understanding the genetic predisposition of congenital immunity key part functioning during CF. Materials and Methods. Examined 48 people with a confirmed diagnosis of CF, homozygotes of the F508del mutation of the CFTR gene, and 104 practically healthy people without the F508del mutation of the CFTR gene from the control group. The following molecular genetic methods were used: DNA extraction from peripheral blood cells, KIR genotyping by PCR-SSP for the presence or absence of the 14 KIR genes (KIR2DL1, 2DL2, 2DL3, 2DL4, 2DL5, 2DS1, 2DS2, 2DS3, 2DS4, 2DS5, 3DL1, 3DL2, 3DL3, 3DS1). Results – molecular genetic studies of KIR-genes repertoire in the group of cystic fibrosis patients showed a decrease in the frequency of genes, responsible for activating NK cells receptors. Of the five examined NK cell activation genes, one gene was completely absent, namely 2DS4, and another (2DS1) was detected in only 3 of 48 patients examined, which was 6.25 %, and this figure is significantly lower in comparison with the control group (c2=4.801, p<0.05). Regarding the genes of NK-cell inhibitory receptors, all investigated genes were detected in the study group (8 in general). By detection frequency, they mostly correspond to the control group, with the exception of the 2DL3 gene, found in patients with CF with a significantly lower frequency (c2=11.97, p<0.005). Conclusion – for the first time in the group of patients with CF, a study was performed on the frequency and spectrum of KIR-genes, responsible for NK cell receptors. Reducing the frequency of activation NK cell receptor genes in patients with CF can lead to a weakening of congenital immunity and the severity of infectious processes during CF

TypeLoader2: Automated submission of novel HLA and killer-cell immunoglobulin-like receptor alleles in full length

The Immuno Polymorphism Database (IPD) databases provide global, curated repositories for information regarding polymorphisms of genes of the immune system, thereby generating immense value for the research and clinical communities. The advent of high‐throughput genotyping in immunogenetics has led to dramatically growing numbers of heretofore unknown HLA and lately also killer‐cell immunoglobulin‐like receptor (KIR) alleles, which are to be curated and deposited in the IPD‐IMGT/HLA and IPD‐KIR databases, respectively. It is highly desirable that these novel alleles are characterised and submitted in full length, and that known alleles are extended to cover the complete gene sequence.

However, the manual annotation and submission of sequences to European Molecular Biology Laboratory's European Nucleotide Archive and the IPD‐IMGT/HLA and IPD‐KIR databases is time‐consuming and error‐prone. Here, we report the substantial extension of the HLA allele submission tool TypeLoader, which now also supports the annotation and submission of KIR alleles. To enable a more widespread use of this tool, we have made it available as a stand‐alone application that can easily be installed on standard Windows or Linux computers. Furthermore, an internal SQLite database was added to store a wide range of metadata about each allele. This allows TypeLoader2 to be used as a lab's central information platform for the annotation, curation and submission of full‐length HLA and KIR allele sequences. The software is freely available from GitHub (https://github.com/DKMS-LSL/typeloader).

We hope that the increased convenience and scope of TypeLoader2 will foster the submission of more full‐length sequences to the IPD‐IMGT/HLA and IPD‐KIR databases, ultimately promoting the use of full‐length sequencing for genotyping both HLA and KIR.

KIR Variation in Iranians Combines High Haplotype and Allotype Diversity With an Abundance of Functional Inhibitory Receptors

Natural killer (NK) cells are innate lymphocytes that eliminate infected and transformed cells. They discriminate healthy from diseased tissue through killer cell Ig-like receptor (KIR) recognition of HLA class I ligands. Directly impacting NK cell function, KIR polymorphism associates with infection control and multiple autoimmune and pregnancy syndromes. Here we analyze KIR diversity of 241 individuals from five groups of Iranians. These five populations represent Baloch, Kurd, and Lur, together comprising 15% of the ethnically diverse Iranian population. We identified 159 KIR alleles, including 11 not previously characterized. We also identified 170 centromeric and 94 telomeric haplotypes, and 15 different KIR haplotypes carrying either a deletion or duplication encompassing one or more complete KIR genes. As expected, comparing our data with those representing major worldwide populations revealed the greatest similarity between Iranians and Europeans. Despite this similarity we observed higher frequencies of KIR3DL1^*001 in Iran than any other population, and the highest frequency of HLA-B^*51, a Bw4-containing allotype that acts as a strong educator of KIR3DL1^*001⁺ NK cells. Compared to Europeans, the Iranians we studied also have a reduced frequency of 3DL1^*004, which encodes an allotype that is not expressed at the NK cell surface. Concurrent with the resulting high frequency of strong viable interactions between inhibitory KIR and polymorphic HLA class I, the majority of KIR-A haplotypes characterized do not express a functional activating receptor. By contrast, the most frequent KIR-B haplotype in Iran expresses only one functional inhibitory KIR and the maximum number of activating KIR. This first complete, high-resolution, characterization of the KIR locus of Iranians will form a valuable reference for future clinical and population studies.

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.

Haplotype-Based Analysis of KIR-Gene Profiles in a South European Population-Distribution of Standard and Variant Haplotypes, and Identification of Novel Recombinant Structures

Inhibitory Killer-cell Immunoglobulin-like Receptors (KIR) specific for HLA class I molecules enable human natural killer cells to monitor altered antigen presentation in pathogen-infected and tumor cells. KIR genes display extensive copy-number variation and allelic polymorphism. They organize in a series of variable arrangements, designated KIR haplotypes, which derive from duplications of ancestral genes and sequence diversification through point mutation and unequal crossing-over events. Genomic studies have established the organization of multiple KIR haplotypes—many of them are fixed in most human populations, whereas variants of those have less certain distributions. Whilst KIR-gene diversity of many populations and ethnicities has been explored superficially (frequencies of individual genes and presence/absence profiles), less abundant are in-depth analyses of how such diversity emerges from KIR-haplotype structures. We characterize here the genetic diversity of KIR in a sample of 414 Spanish individuals. Using a parsimonious approach, we manage to explain all 38 observed KIR-gene profiles by homo- or heterozygous combinations of six fixed centromeric and telomeric motifs; of six variant gene arrangements characterized previously by us and others; and of two novel haplotypes never detected before in Caucasoids. Associated to the latter haplotypes, we also identified the novel transcribed KIR2DL5B^*0020202 allele, and a chimeric KIR2DS2/KIR2DL3 gene (designated KIR2DL3^*033) that challenges current criteria for classification and nomenclature of KIR genes and haplotypes.

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.

NK Cell-Based Immunotherapies in Cancer

With the development of technologies that can transform immune cells into therapeutic modalities, immunotherapy has remarkably changed the current paradigm of cancer treatment in recent years. NK cells are components of the innate immune system that act as key regulators and exhibit a potent tumor cytolytic function. Unlike T cells, NK cells exhibit tumor cytotoxicity by recognizing non-self, without deliberate immunization or activation. Currently, researchers have developed various approaches to improve the number and anti-tumor function of NK cells. These approaches include the use of cytokines and Abs to stimulate the efficacy of NK cell function, adoptive transfer of autologous or allogeneic ex vivo expanded NK cells, establishment of homogeneous NK cell lines using the NK cells of patients with cancer or healthy donors, derivation of NK cells from induced pluripotent stem cells (iPSCs), and modification of NK cells with cutting-edge genetic engineering technologies to generate chimeric Ag receptor (CAR)-NK cells. Such NK cell-based immunotherapies are currently reported as being promising anti-tumor strategies that have shown enhanced functional specificity in several clinical trials investigating malignant tumors. Here, we summarize the recent advances in NK cell-based cancer immunotherapies that have focused on providing improved function through the use of the latest genetic engineering technologies. We also discuss the different types of NK cells developed for cancer immunotherapy and present the clinical trials being conducted to test their safety and efficacy.

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.

The IPD Project: a centralised resource for the study of polymorphism in genes of the immune system

The Immuno Polymorphism Database (IPD), https://www.ebi.ac.uk/ipd/, is a set of specialist databases that enable the study of polymorphic genes which function as part of the vertebrate immune system. The major focus is on the hyperpolymorphic major histocompatibility complex (MHC) genes and the killer-cell immunoglobulin-like receptor (KIR) genes, by providing the official repository and primary source of sequence data. Databases are centred around humans as well as animals important for food security, for companionship and as disease models. The IPD project works with specialist groups or nomenclature committees who provide and manually curate individual sections before they are submitted for online publication. To reflect the recent advance of allele sequencing technologies and the increasing demands of novel tools for the analysis of genomic variation, the IPD project is undergoing a progressive redesign and reorganisation. In this review, recent updates and future developments are discussed, with a focus on the core concepts to better future-proof the project.

Selection of matched unrelated donors moving forward: from HLA allele counting to functional matching

Matched unrelated donors (URD) are the most frequent source of stem cells for allogeneic hematopoietic cell transplantation (HCT) to date, with HCT performed mainly under conventional immunosuppression by methotrexate and cyclosporine. In this setting, every single allelic donor-recipient mismatch for HLA-A, -B, -C, -DRB1 (8/8), but not for HLA-DQB1, -DPB1, has a significant negative effect on overall survival (OS). When several 8/8 HLA-matched URD are available, donor age is the most important factor impacting OS. Moving forward from the traditional way of counting the number of donor-recipient HLA allele mismatches to biology-driven algorithms for functional matching has led to the unraveling of an association between permissive, low-risk HLA-DPB1 mismatches and improved outcome after URD HCT for malignant disease but not for nonmalignant disease. Functional HLA matching might prove to have increasing importance for URD selection in the era of new immunosuppressive regimens that have the potential to substantially reshuffle the role of HLA mismatches in URD HCT.

Bioinformatics Strategies, Challenges, and Opportunities for Next Generation Sequencing-Based HLA Genotyping

The advent of next generation sequencing (NGS) has altered the face of genotyping the human leukocyte antigen (HLA) system in clinical, stem cell donor registry, and research contexts. NGS has led to a dramatically increased sequencing throughput at high accuracy, while being more time and cost efficient than precursor technologies. This has led to a broader and deeper profiling of the key genes in the human immunogenetic make-up. The rapid evolution of sequencing technologies is evidenced by the development of varied short-read sequencing platforms with differing read lengths and sequencing capacities to long-read sequencing platforms capable of profiling full genes without fragmentation. Concomitantly, there has been development of a diverse set of computational analyses and software tools developed to deal with the various strengths and limitations of the sequencing data generated by the different sequencing platforms. This review surveys the different modalities involved in generating NGS HLA profiling sequence data. It systematically describes various computational approaches that have been developed to achieve HLA genotyping to different degrees of resolution. At each stage, this review enumerates the drawbacks and advantages of each of the platforms and analysis approaches, thus providing a comprehensive picture of the current state of HLA genotyping technologies.

Novel Approach to Cell Surface Discrimination Between KIR2DL1 Subtypes and KIR2DS1 Identifies Hierarchies in NK Repertoire, Education, and Tolerance

Cumulative activating and inhibitory receptor signaling controls the functional output of individual natural killer (NK) cells. Investigation of how competing signals impact response, however, has been hampered by the lack of available antibodies capable of distinguishing inhibitory and activating receptors with highly homologous ectodomains. Utilizing a novel combination of monoclonal antibodies with specificity for discrete inhibitory KIR2DL1 and activating KIR2DS1 allotypes found among 230 healthy donors, we investigated allele-specific receptor expression and function driven by KIR2DL1 and KIR2DS1 alleles. We found that co-expression of the HLA-C2 ligand diminishes KIR2DL1, but not KIR2DS1, cell surface staining, but does not impact the respective frequencies of KIR2DL1- and KIR2DS1-expressing cells within the NK repertoire. We can distinguish by flow cytometry NK cell populations expressing the most common KIR2DL1-C245 allotypes from those expressing the most common KIR2DL1-R245 allotypes, and we show that the informative differential binding anti-KIR2DL1/S1 clone 1127B is determined by amino acid residue T154. Although both KIR2DL1-C245 and KIR2DL1-R245 subtypes can be co-expressed in the same cell, NK cells preferentially express one or the other. Cells expressing KIR2DL1-C245 exhibited a lower KIR2DL1 cell surface density and lower missing-self reactivity in comparison to cells expressing KIR2DL1-R245. We found no difference, however, in sensitivity to inhibition or cell surface stability between the two KIR2DL1 isoforms, and both demonstrated similar expansion among NKG2C+ KIR2DL1+ NK cells in HCMV-seropositive healthy individuals. In addition to cell surface density of KIR2DL1, copy number of cognate HLA-C2 hierarchically impacted the effector capacity of both KIR2DL1+ cells and the tolerization of KIR2DS1+ NK cells. HLA-C2 tolerization of KIR2DS1+ NK cells could be overridden, however, by education via co-expressed self-specific inhibitory receptors, such as the heterodimer CD94/NKG2A. Our results demonstrate that effector function of NK cells expressing KIR2DL1 or KIR2DS1 is highly influenced by genetic variability and is calibrated by co-expression of additional NK receptors and cognate HLA-C2 ligands. These findings define the molecular conditions under which NK cells are activated or inhibited, potentially informing selection of donors for adoptive NK therapies.