Extensive variation in gene copy number at the killer immunoglobulin-like receptor locus in humans

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.

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.

Killer-cell immunoglobulin-like receptor assessment algorithms in haemopoietic progenitor cell transplantation: current perspectives and future opportunities

Natural killer cells preferentially target and kill malignant and virally infected cells. Both these properties present compelling clinical utility in the field of haemopoietic progenitor cell transplantation (HPCT), potentially promoting a graft vs leukaemia effect in the absence of graft vs host disease and protecting against cytomegalovirus activation. Killer Ig-like receptors (KIR) play a central role in the cytotoxic action of natural killer cells, providing opportunity for improving transplantation outcomes by prioritising potential donors with optimal characteristics. Numerous algorithms for assessing KIR gene content as part of HPCT donor selection protocols exist, but no single model has been found to be universally applicable in all transplant centres. This review summarises several of the predominant strategies in KIR assessment algorithms, discussing their basic scientific principles, clinical utility and benefits to post-transplant outcomes. Finally, the review will consider how future donor selection protocols could develop towards unifying the concepts of KIR proteomics and genetics for optimising patient care.

Analysis of five deep-sequenced trio-genomes of the Peninsular Malaysia Orang Asli and North Borneo populations

Background

Recent advances in genomic technologies have facilitated genome-wide investigation of human genetic variations. However, most efforts have focused on the major populations, yet trio genomes of indigenous populations from Southeast Asia have been under-investigated.

Results

We analyzed the whole-genome deep sequencing data (~ 30×) of five native trios from Peninsular Malaysia and North Borneo, and characterized the genomic variants, including single nucleotide variants (SNVs), small insertions and deletions (indels) and copy number variants (CNVs). We discovered approximately 6.9 million SNVs, 1.2 million indels, and 9000 CNVs in the 15 samples, of which 2.7% SNVs, 2.3% indels and 22% CNVs were novel, implying the insufficient coverage of population diversity in existing databases. We identified a higher proportion of novel variants in the Orang Asli (OA) samples, i.e., the indigenous people from Peninsular Malaysia, than that of the North Bornean (NB) samples, likely due to more complex demographic history and long-time isolation of the OA groups. We used the pedigree information to identify de novo variants and estimated the autosomal mutation rates to be 0.81 × 10^− 8 – 1.33 × 10^− 8, 1.0 × 10^− 9 – 2.9 × 10^− 9, and ~ 0.001 per site per generation for SNVs, indels, and CNVs, respectively. The trio-genomes also allowed for haplotype phasing with high accuracy, which serves as references to the future genomic studies of OA and NB populations. In addition, high-frequency inherited CNVs specific to OA or NB were identified. One example is a 50-kb duplication in DEFA1B detected only in the Negrito trios, implying plausible effects on host defense against the exposure of diverse microbial in tropical rainforest environment of these hunter-gatherers. The CNVs shared between OA and NB groups were much fewer than those specific to each group. Nevertheless, we identified a 142-kb duplication in AMY1A in all the 15 samples, and this gene is associated with the high-starch diet. Moreover, novel insertions shared with archaic hominids were identified in our samples.

Conclusion

Our study presents a full catalogue of the genome variants of the native Malaysian populations, which is a complement of the genome diversity in Southeast Asians. It implies specific population history of the native inhabitants, and demonstrated the necessity of more genome sequencing efforts on the multi-ethnic native groups of Malaysia and Southeast Asia.

Killer Ig-Like Receptors (KIRs): Their Role in NK Cell Modulation and Developments Leading to Their Clinical Exploitation

Natural killer (NK) cells contribute to the first line of defense against viruses and to the control of tumor growth and metastasis spread. The discovery of HLA class I specific inhibitory receptors, primarily of killer Ig-like receptors (KIRs), and of activating receptors has been fundamental to unravel NK cell function and the molecular mechanisms of tumor cell killing. Stemmed from the seminal discoveries in early '90s, in which Alessandro Moretta was the major actor, an extraordinary amount of research on KIR specificity, genetics, polymorphism, and repertoire has followed. These basic notions on NK cells and their receptors have been successfully translated to clinical applications, primarily to the haploidentical hematopoietic stem cell transplantation to cure otherwise fatal leukemia in patients with no HLA compatible donors. The finding that NK cells may express the PD-1 inhibitory checkpoint, particularly in cancer patients, may allow understanding how anti-PD-1 therapy could function also in case of HLA class I^neg tumors, usually susceptible to NK-mediated killing. This, together with the synergy of therapeutic anti-checkpoint monoclonal antibodies, including those directed against NKG2A or KIRs, emerging in recent or ongoing studies, opened new solid perspectives in cancer therapy.

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

The killer-cell immunoglobulin-like receptor (KIR) genes regulate natural killer cell activity, influencing predisposition to immune mediated disease, and affecting hematopoietic stem cell transplantation (HSCT) outcome. Owing to the complexity of the KIR locus, with extensive gene copy number variation (CNV) and allelic diversity, high-resolution characterization of KIR has so far been applied only to relatively small cohorts. Here, we present a comprehensive high-throughput KIR genotyping approach based on next generation sequencing. Through PCR amplification of specific exons, our approach delivers both copy numbers of the individual genes and allelic information for every KIR gene. Ten-fold replicate analysis of a set of 190 samples revealed a precision of 99.9%. Genotyping of an independent set of 360 samples resulted in an accuracy of more than 99% taking into account consistent copy number prediction. We applied the workflow to genotype 1.8 million stem cell donor registry samples. We report on the observed KIR allele diversity and relative abundance of alleles based on a subset of more than 300,000 samples. Furthermore, we identified more than 2,000 previously unreported KIR variants repeatedly in independent samples, underscoring the large diversity of the KIR region that awaits discovery. This cost-efficient high-resolution KIR genotyping approach is now applied to samples of volunteers registering as potential donors for HSCT. This will facilitate the utilization of KIR as additional selection criterion to improve unrelated donor stem cell transplantation outcome. In addition, the approach may serve studies requiring high-resolution KIR genotyping, like population genetics and disease association studies.

CD8dim but not CD8bright cells positive to CD56 dominantly express KIR and are cytotoxic during visceral leishmaniasis

This study reports a structural and functional heterogeneity of CD8⁺CD56⁺NKT cells, which usually decrease quantitatively during visceral leishmaniasis. Based on fluorescence intensity of CD8 receptors on CD56⁺NKT cells, two populations of CD8⁺CD56⁺NKT cells have been identified. These cells were recognized as CD8^dimCD56⁺NKT and CD8^brightCD56⁺NKT cells. We further analyzed the functional nature of CD8^dim and CD8^bright positive CD56⁺NKT cells. In comparison to CD8^brightCD56⁺NKT cells, a significantly higher percentage of CD8^dimCD56⁺NKT cells expressed KIR during VL. The percentage of CD8^dimCD56⁺NKT cells expressing KIR was found 4 fold higher in VL as compared to healthy subjects. But, the difference was insignificant in case of CD8^brightCD56⁺NKT cells. CD8⁺CD56⁺NKT cells release granzyme B to kill the infected cells. A categorical difference was also observed in the function of CD8^dimCD56⁺NKT and CD8^brightCD56⁺NKT cells during visceral leishmaniasis. The percentage of granzyme B expressing CD8^dimCD56⁺NKT cells was 2.83 fold higher in VL compared to healthy subjects. But, there was no significant difference in granzyme B expressing CD8^brightCD56⁺NKT cells in samples from healthy and VL subjects. However, within VL subject, the percentage of granzyme B expressing CD8^dimCD56⁺NKT cells was 5.7 fold higher in comparison to CD8^brightCD56⁺NKT cells. This study concludes that CD8^dimCD56⁺NKT cells are more cytotoxic than CD8^brightCD56⁺NKT cells during VL.

Killer Immunoglobulin-Like Receptor Allele Determination Using Next-Generation Sequencing Technology

The impact of natural killer (NK) cell alloreactivity on hematopoietic stem cell transplantation (HSCT) outcome is still debated due to the complexity of graft parameters, HLA class I environment, the nature of killer cell immunoglobulin-like receptor (KIR)/KIR ligand genetic combinations studied, and KIR⁺ NK cell repertoire size. KIR genes are known to be polymorphic in terms of gene content, copy number variation, and number of alleles. These allelic polymorphisms may impact both the phenotype and function of KIR⁺ NK cells. We, therefore, speculate that polymorphisms may alter donor KIR⁺ NK cell phenotype/function thus modulating post-HSCT KIR⁺ NK cell alloreactivity. To investigate KIR allele polymorphisms of all KIR genes, we developed a next-generation sequencing (NGS) technology on a MiSeq platform. To ensure the reliability and specificity of our method, genomic DNA from well-characterized cell lines were used; high-resolution KIR typing results obtained were then compared to those previously reported. Two different bioinformatic pipelines were used allowing the attribution of sequencing reads to specific KIR genes and the assignment of KIR alleles for each KIR gene. Our results demonstrated successful long-range KIR gene amplifications of all reference samples using intergenic KIR primers. The alignment of reads to the human genome reference (hg19) using BiRD pipeline or visualization of data using Profiler software demonstrated that all KIR genes were completely sequenced with a sufficient read depth (mean 317× for all loci) and a high percentage of mapping (mean 93% for all loci). Comparison of high-resolution KIR typing obtained to those published data using exome capture resulted in a reported concordance rate of 95% for centromeric and telomeric KIR genes. Overall, our results suggest that NGS can be used to investigate the broad KIR allelic polymorphism. Hence, these data improve our knowledge, not only on KIR⁺ NK cell alloreactivity in HSCT but also on the role of KIR⁺ NK cell populations in control of viral infections and diseases.

Deciphering the killer-cell immunoglobulin-like receptor system at super-resolution for natural killer and T-cell biology

Killer-cell immunoglobulin-like receptors (KIRs) are components of two fundamental biological systems essential for human health and survival. First, they contribute to host immune responses, both innate and adaptive, through their expression by natural killer cells and T cells. Second, KIR play a key role in regulating placentation, and hence reproductive success. Analogous to the diversity of their human leucocyte antigen class I ligands, KIR are extremely polymorphic. In this review, we describe recent developments, fuelled by methodological advances, that are helping to decipher the KIR system in terms of haplotypes, polymorphisms, expression patterns and their ligand interactions. These developments are delivering deeper insight into the relevance of KIR in immune system function, evolution and disease.

Estimating KIR Haplotype Frequencies on a Cohort of 10,000 Individuals: A Comprehensive Study on Population Variations, Typing Resolutions, and Reference Haplotypes

The killer cell immunoglobulin-like receptors (KIR) mediate human natural killer (NK) cell cytotoxicity via activating or inhibiting signals. Although informative and functional haplotype patterns have been reported, most genotyping has been performed at resolutions that are structurally ambiguous. In order to leverage structural information given low-resolution genotypes, we performed experiments to quantify the effects of population variations, reference haplotypes, and genotyping resolutions on population-level haplotype frequency estimations as well as predictions of individual haplotypes. We genotyped 10,157 unrelated individuals in 5 populations (518 African American[AFA], 258 Asian or Pacific Islander[API], 8,245 European[EUR], 1,073 Hispanic[HIS], and 63 Native American[NAM]) for KIR gene presence/absence (PA), and additionally half of the AFA samples for KIR gene copy number variation (CNV). A custom EM algorithm was used to estimate haplotype frequencies for each population by interpretation in the context of three sets of reference haplotypes. The algorithm also assigns each individual the haplotype pairs of maximum likelihood. Generally, our haplotype frequency estimates agree with similar previous publications to within <5% difference for all haplotypes. The exception is that estimates for NAM from the U.S. showed higher frequency association of cB02 with tA01 (+14%) instead of tB01 (-8.5%) compared to a previous study of NAM from south of the U.S. The higher-resolution CNV genotyping on the AFA samples allowed unambiguous haplotype-pair assignments for the majority of individuals, resulting in a 22% higher median typing resolution score (TRS), which measures likelihood of self-match in the context of population-specific haplo- and geno-types. The use of TRS to quantify reduced ambiguity with CNV data clearly revealed the few individuals with ambiguous genotypes as outliers. It is observed that typing resolution and reference haplotype set influence haplotype frequency estimates. For example, PA resolution may be used with reference haplotype sets up to the point where certain haplotypes are gene-content subsets of others. At that point, CNV must be used for all genes.

Novel insights in the genomic organization and hotspots of recombination in the human KIR locus through analysis of intergenic regions

The Killer Immunoglobulin-like Receptor (KIR) proteins constitute a family of highly homologous surface receptors involved in the regulation of the innate cytotoxicity of natural killer (NK) cells. Within the human genome, 17 KIR genes are present, many of which show large variation across the population owing to the high number of allelic variants and copy number variation (CNV). KIR genotyping and CNV determination were used to map the KIR locus in a large cohort of >400 Caucasian individuals. Gene order and structure was determined by sequence-specific polymerase chain reaction of the intergenic regions. In this way, we could show that KIR3DL1 and KIR2DS4 gene variants are linked and that--contrary to current views--the gene KIR2DS5 is only present in the telomeric half of the KIR locus. Our study revealed novel insights in the highly organized distribution of KIR genes. Novel recombination hotspots were identified that contribute to the diversity of KIR gene distribution in the Caucasian population. Next-generation sequencing of the KIR intergenic regions allowed for a detailed single-nucleotide polymorphism analysis, which demonstrated several gene-specific as well as haplotype-specific nucleotides for a more accurate genotyping of this notoriously complex gene cluster.