Nature of allelic sequence polymorphism at the KIR3DL3 locus

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.

Remarkably low KIR and HLA diversity in Amerindians reveals signatures of strong purifying selection shaping the centromeric KIR region

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR-HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR-HLA interactions among all described worldwide populations, and that 83-97% of their KIR-HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR-HLA coevolution and its impact on human health and survival.

Quantitative Multiplex Real-Time Reverse Transcriptase-Polymerase Chain Reaction with Fluorescent Probe Detection of Killer Immunoglobulin-Like Receptors, KIR2DL4/3DL3

(1) Background: KIR2DL4/KIR3DL3 are the framework genes present in all KIR haplotypes, with unique expression patterns being present only in women and CD56bright NK cells. KIR genes have a high degree of DNA sequence identity. Consequently, they are one of the most challenging genes for molecular detection—especially regarding expressions; (2) Methods: We developed an effective method to determine KIR3DL3/KIR2DL4 expressions based on a multiplex quantitative real-time Reverse transcription polymerase chain reaction (qRT-PCR )with fluorescent probes using NK92; (3) Results: Standardizations of the singleplex KIR2DL4 and KIR3DL3 were performed to evaluate the sensitivity and specificity for further development of the multiplex assay. The limit of detection was at 500 copies each. There was cross-amplification with the presence of related KIR genes at a level of 5 × 10⁷ copies. This is not biologically significant because this high level of KIR expression has not been found in clinical samples. The multiplex assay was reproducible equivalent to its singleplex (KIR2DL4; R² = 0.995, KIR3DL3; R² = 0.996, but lower sensitivity of 10³ copies). Furthermore, the validation of the developed method on samples of blood donors showed high sensitivity (100%) and specificity (99.9%); (4) Conclusions: The developed method is reliable and highly specific suitable for evaluation of the KIR2DL4/3DL3 mRNA expressions in further applications.

Diversity of Killer Cell Immunoglobulin-Like Receptors and Disease

Natural killer (NK) cells are bone marrow-derived large granular lymphocytes defined by CD3^negCD56^pos and represent 5% to 25% of peripheral blood mononuclear cell fraction of the healthy humans. NK cells have a highly specific and sophisticated target cell recognition receptor system arbitrated by the integration of signals triggered by a multitude of inhibitory and activating receptors. Human NK cells express distinct families of receptors, including (1) killer cell immunoglobulin-like receptors, (2) killer cell lectin-like receptors, (3) leukocyte immunoglobulin-like receptors, and (4) natural cytotoxicity receptors.

Conservation, Extensive Heterozygosity, and Convergence of Signaling Potential All Indicate a Critical Role for KIR3DL3 in Higher Primates

Natural killer (NK) cell functions are modulated by polymorphic killer cell immunoglobulin-like receptors (KIR). Among 13 human KIR genes, which vary by presence and copy number, KIR3DL3 is ubiquitously present in every individual across diverse populations. No ligand or function is known for KIR3DL3, but limited knowledge of expression suggests involvement in reproduction, likely during placentation. With 157 human alleles, KIR3DL3 is also highly polymorphic and we show heterozygosity exceeds that of HLA-B in many populations. The external domains of catarrhine primate KIR3DL3 evolved as a conserved lineage distinct from other KIR. Accordingly, and in contrast to other KIR, we show the focus of natural selection does not correspond exclusively to known ligand binding sites. Instead, a strong signal for diversifying selection occurs in the D1 Ig domain at a site involved in receptor aggregation, which we show is polymorphic in humans worldwide, suggesting differential ability for receptor aggregation. Meanwhile in the cytoplasmic tail, the first of two inhibitory tyrosine motifs (ITIM) is conserved, whereas independent genomic events have mutated the second ITIM of KIR3DL3 alleles in all great apes. Together, these findings suggest that KIR3DL3 binds a conserved ligand, and a function requiring both receptor aggregation and inhibitory signal attenuation. In this model KIR3DL3 resembles other NK cell inhibitory receptors having only one ITIM, which interact with bivalent downstream signaling proteins through dimerization. Due to the extensive conservation across species, selection, and other unusual properties, we consider elucidating the ligand and function of KIR3DL3 to be a pressing question.

Different Selected Mechanisms Attenuated the Inhibitory Interaction of KIR2DL1 with C2+ HLA-C in Two Indigenous Human Populations in Southern Africa

The functions of human NK cells in defense against pathogens and placental development during reproduction are modulated by interactions of killer cell Ig-like receptors (KIRs) with HLA-A, -B and -C class I ligands. Both receptors and ligands are highly polymorphic and exhibit extensive differences between human populations. Indigenous to southern Africa are the KhoeSan, the most ancient group of modern human populations, who have highest genomic diversity worldwide. We studied two KhoeSan populations, the Nama pastoralists and the ≠Khomani San hunter-gatherers. Comprehensive next-generation sequence analysis of HLA-A, -B, and -C and all KIR genes identified 248 different KIR and 137 HLA class I, which assort into ∼200 haplotypes for each gene family. All 74 Nama and 78 ≠Khomani San studied have different genotypes. Numerous novel KIR alleles were identified, including three arising by intergenic recombination. On average, KhoeSan individuals have seven to eight pairs of interacting KIR and HLA class I ligands, the highest diversity and divergence of polymorphic NK cell receptors and ligands observed to date. In this context of high genetic diversity, both the Nama and the ≠Khomani San have an unusually conserved, centromeric KIR haplotype that has arisen to high frequency and is different in the two KhoeSan populations. Distinguishing these haplotypes are independent mutations in KIR2DL1, which both prevent KIR2DL1 from functioning as an inhibitory receptor for C2⁺ HLA-C. The relatively high frequency of C2⁺ HLA-C in the Nama and the ≠Khomani San appears to have led to natural selection against strong inhibitory C2-specific KIR.

Overview of the killer cell immunoglobulin-like receptor system

Natural killer (NK) cells are more than simple killers and have been implicated in control and clearance of malignant and virally infected cells, regulation of adaptive immune responses, rejection of bone marrow transplants, and autoimmunity and the maintenance of pregnancy. Human NK cells largely use a family of germ-line encoded killer cell immunoglobulin-like receptors (KIR) to respond to the perturbations from self-HLA class I molecules present on infected, malignant, or HLA-disparate fetal or allogenic transplants. Genes encoding KIR receptors and HLA class I ligands are located on different chromosomes, and both feature extraordinary diversity in the number and type of genes. The independent segregation of KIR and HLA gene families produce diversity in the number and type of KIR-HLA gene combinations inherited in individuals, which may determine their immunity and susceptibility to diseases. This chapter provides an overview of NK cells and their unprecedented phenotypic and functional diversity within and between individuals.

Conserved KIR allele-level haplotypes are altered by microvariation in individuals with European ancestry

NK cell immunoglobulin-like receptor (KIR) haplotype-specific DNA fragments were sequenced to identify centromeric and telomeric allele-level haplotype structures and their frequencies from 76 unrelated individuals with European ancestry. Analysis was simplified by redefining the 5' boundary of the centromeric KIR gene cluster to include only exons 7-9 of KIR3DL3. Three consensus allele-level haplotypes were identified for a centromeric gene presence/absence structure designated as Cen-A1. KIR3DL3*00201 (exons 7-9)∼KIR2DL3*001∼KIR2DL1*00302 was the most frequent (37.5%) centromeric structure. Single consensus haplotypes were observed for haplotype structures Cen-B1 and Cen-B2. Six Tel-A1 and two Tel-B1 consensus haplotypes were observed; the most prevalent (23.0%) was KIR2DL4*00102∼KIR3DL1*002∼KIR2DS4*00101∼KIR3DL2*002. A small number of nucleotide substitutions (≤3) in the coding regions of the functional KIR genes created microvariants of the consensus haplotypes. Eight less common haplotype structures were also detected. Four carried hybrid genes formed during gene deletion events, two carried an insertion with a 2DL5/3DP1 fusion gene and two included a very large insertion. These data show that the KIR gene complex is composed of a limited number of conserved allele-level centromeric and telomeric haplotypes that have diversified by mutation, recombination within a locus and unequal crossing over.

Variable NK cell receptors exemplified by human KIR3DL1/S1

Variegated expression of variable NK cell receptors for polymorphic MHC class I broadens the range of an individual's NK cell response and the capacity for populations and species to survive disease epidemics and population bottlenecks. On evolutionary time scales, this component of immunity is exceptionally dynamic, unstable, and short-lived, being dependent on coevolution of ligands and receptors subject to varying, competing selection pressures. Consequently these systems of variable NK cell receptors are largely species specific and have recruited different classes of glycoprotein, even within the primate order of mammals. Such disparity helps to explain substantial differences in NK cell biology between humans and animal models, for which the population genetics is largely ignored. KIR3DL1/S1, which recognizes the Bw4 epitope of HLA-A and -B and is the most extensively studied of the variable NK cell receptors, exemplifies how variation in all possible parameters of function is recruited to diversify the human NK cell response.

The characteristics of allelic polymorphism in killer-immunoglobulin-like receptor framework genes in African Americans

The frequencies of alleles of killer cell immunoglobulin-like receptor genes, KIR3DL3 and KIR3DL2, and the carrier frequency of KIR2DL4 alleles have been determined from a population of African Americans (n = 100) by DNA sequencing of the coding regions. Fifty alleles of KIR3DL3 were observed with the most frequent, KIR3DL3*00901 (13%). KIR3DL2 was also diverse; 32 alleles with KIR3DL2*00103 the most frequent (17%). For KIR2DL4, of the 18 alleles observed, one allele, KIR2DL4*00103, was found in 64 of the 100 individuals. Thirty-six novel alleles encoding a total of 28 unique receptors are described. Pairwise comparisons among all of the alleles at each locus suggest a predominance of synonymous substitutions. The variation at all three framework loci fits a neutral model of evolution.

Primate-specific regulation of natural killer cells

Natural killer (NK) cells are circulating lymphocytes that function in innate immunity and placental reproduction. Regulating both development and function of NK cells is an array of variable and conserved receptors that interact with major histocompatibility complex (MHC) class I molecules. Families of lectin-like and immunoglobulin-like receptors are determined by genes in the natural killer complex (NKC) and leukocyte receptor complex (LRC), respectively. As a consequence of the strong, varying pressures on the immune and reproductive systems, NK cell receptors and their MHC class I ligands evolve rapidly, are highly diverse and exhibit dramatic species-specific differences. The variable, polymorphic family of killer cell immunoglobulin-like receptors (KIR) that regulate human NK cell development and function arose recently, from a single-copy gene during the evolution of simian primates. Our studies of KIR and MHC class I genes in representative species show how these two unlinked but functionally intertwined genetic complexes have co-evolved. In humans, combinations of KIR and HLA class I factors are associated with infectious diseases, including HIV/AIDS, autoimmunity, reproductive success and the outcome of therapeutic transplantation. The extraordinary, and unanticipated, divergence of human NK cell receptors and MHC class I ligands from their mouse counterparts can in part explain the difficulties experienced in finding informative mouse models for human diseases. Non-human primate models have far greater potential, but to realize their promise will first require more complete definition of the genetics and function of KIR and MHC variation in non-human primate species, at a level comparable to that achieved for the human species.

Donor selection for natural killer cell receptor genes leads to superior survival after unrelated transplantation for acute myelogenous leukemia

Killer-cell immunoglobulin-like receptor (KIR) genes form a diverse, immunogenetic system. Group A and B KIR haplotypes have distinctive centromeric (Cen) and telomeric (Tel) gene-content motifs. Aiming to develop a donor selection strategy to improve transplant outcome, we compared the contribution of these motifs to the clinical benefit conferred by B haplotype donors. We KIR genotyped donors from 1409 unrelated transplants for acute myelogenous leukemia (AML; n = 1086) and acute lymphoblastic leukemia (ALL; n = 323). Donor KIR genotype influenced transplantation outcome for AML but not ALL. Compared with A haplotype motifs, centromeric and telomeric B motifs both contributed to relapse protection and improved survival, but Cen-B homozygosity had the strongest independent effect. With Cen-B/B homozygous donors the cumulative incidence of relapse was 15.4% compared with 36.5% for Cen-A/A donors (relative risk of relapse 0.34; 95% confidence interval 0.2-0.57; P < .001). Overall, significantly reduced relapse was achieved with donors having 2 or more B gene-content motifs (relative risk 0.64; 95% confidence interval 0.48-0.86; P = .003) for both HLA-matched and mismatched transplants. KIR genotyping of several best HLA-matched potential unrelated donors should substantially increase the frequency of transplants by using grafts with favorable KIR gene content. Adopting this practice could result in superior disease-free survival for patients with AML.

A small, variable, and irregular killer cell Ig-like receptor locus accompanies the absence of MHC-C and MHC-G in gibbons

The killer cell Ig-like receptors (KIRs) of NK cells recognize MHC class I ligands and function in placental reproduction and immune defense against pathogens. During the evolution of monkeys, great apes, and humans, an ancestral KIR3DL gene expanded to become a diverse and rapidly evolving gene family of four KIR lineages. Characterizing the KIR locus are three framework regions, defining two intervals of variable gene content. By analysis of four KIR haplotypes from two species of gibbon, we find that the smaller apes do not conform to these rules. Although diverse and irregular in structure, the gibbon haplotypes are unusually small, containing only two to five functional genes. Comparison with the predicted ancestral hominoid KIR haplotype indicates that modern gibbon KIR haplotypes were formed by a series of deletion events, which created new hybrid genes as well as eliminating ancestral genes. Of the three framework regions, only KIR3DL3 (lineage V), defining the 5' end of the KIR locus, is present and intact on all gibbon KIR haplotypes. KIR2DL4 (lineage I) defining the central framework region has been a major target for elimination or inactivation, correlating with the absence of its putative ligand, MHC-G, in gibbons. Similarly, the MHC-C-driven expansion of lineage III KIR genes in great apes has not occurred in gibbons because they lack MHC-C. Our results indicate that the selective forces shaping the size and organization of the gibbon KIR locus differed from those acting upon the KIR of other hominoid species.

Methods for assessing gene content diversity of KIR with examples from a global set of populations

A number of statistical methods are widely used to describe allelic variation at specific genetic loci and its implication on the evolutionary history of these loci. Although the methods were developed primarily to study allelic variation at loci that are virtually always present in the genome, they are often applied to data of gene content variation (i.e., presence/absence of multiple homologous genes) at the killer cell immunoglobulin-like receptor (KIR) gene cluster. In this paper, we discuss methodological issues involved in the analysis of gene content variation data in the KIR region and also its covariation with polymorphism at the human leukocyte antigen class I loci, which encode ligands for KIR. A comparison of several statistical methods and measures (gene frequency, haplotype frequency, and linkage disequilibrium estimation) using the Centre d'Etude du Polymorphisme Humain data will be provided using KIR haplotypes that have been determined by segregation analysis, noting the strengths and weaknesses of the methods when only the presence/absence data is considered. Finally, application of these methods to a set of globally distributed populations is described (see Single et al., Nat Genet 39:1114-1119, 2007) in order to illustrate the challenges faced when inferring the joint effects of natural selection and demographic history on these immune-related genes.

KIR2DL5 alleles mark certain combination of activating KIR genes

Killer cell Ig-like receptors (KIR) control the immune response of NK cells and some T cells to infections and tumors. KIR genes evolve rapidly and are variable between individuals in their number, type and sequence. Here, we determined the nature of KIR2DL5 gene polymorphism in four ethnic groups using direct DNA sequencing method. Nine new sequences were discovered. Within the panel of 248 KIR2DL5-positive individuals, 14 KIR2DL5-sequences differing in coding regions were observed. They differed at only seven amino acid positions, and such limited polymorphism is consistent with its conserved nature throughout the hominoid lineage. Ethnic deviation was seen in the distribution of KIR2DL5A, KIR2DL5B and their alleles. African Americans had more KIR2DL5 alleles than other populations indicating that more polymorphisms are yet to be discovered in Africans. Linkage between KIR2DL5-alleles and certain activating-KIR genes were observed, but frequency of these linked clusters differed substantially between populations. Consequently, KIR2DL5 alleles can be used as markers to predict the activating-KIR gene content. Typing system distinguishing A*001 and B*002 alleles can serve as a powerful screening test to assess the content of most variable activating-KIR genes that have been implicated in human disease and in the outcome of hematopoietic stem cell transplantation.

Seventeen novel alleles add to the already extensive KIR3DL3 diversity

Exons 2-9 of KIR3DL3 alleles were characterized by genomic DNA sequencing in two families and in 78 bone marrow transplant samples. Several strategies were used to isolate single alleles for characterization and to resolve alternative allele combinations. We describe 17 novel 3DL3 alleles carried by 30 individuals. Compared with the most closely matched alleles, the new alleles differ by from one to three nucleotides and from zero to three amino acids. The majority of the substitutions were shared with other 3DL3 alleles although three novel polymorphic codons, 151 (exon 4), 327 (exon 7) and 352 (exon 9), are described. Of the 36 different 3DL3 alleles detected in the transplant population, the three most common alleles accounting for 47% of the total were KIR3DL3*00101 (13.5%), KIR3DL3*003 (21.2%) and KIR3DL3*00901 (12.2%).

Facilitation of KIR genotyping by a PCR-SSP method that amplifies short DNA fragments

Detection of killer-cell immunoglobulin-like receptors (KIR) genes by polymerase chain reaction with sequence-specific primers (PCR-SSP) led in 1997 to the discovery that human genomes diverge largely in the KIR they encode. While only a few KIR genes are conserved in all humans, most individuals lack several those genes, which tend to associate in diverse haplotypic combinations. The PCR-SSP technique, updated to detect the more recently identified KIR genes and alleles, is still used widely to analyze the diversity of human populations, and to study the influence of KIR-gene variability on human health. Several published PCR-SSP methods for KIR genotyping, although simple and robust, have the drawback of relying on the amplification of DNA fragments spanning 0.5-2.0 kbp, which tends to fail in low-quality DNAs. Valuable collections of DNAs often include such poor quality samples, which lead to loss of data and resources. Even worse, undetected falsely negative or positive reactions may result in erroneous gene frequencies and in odd gene combinations. To address those problems, we have redesigned our previously published KIR genotyping method so that it produces short amplicons (less than 200 bp for most genes). This modification minimizes amplification failures, thus conferring greater consistency and reliability to KIR genotyping. In addition, the new PCR-SSP method detects recently described alleles of several KIR genes, and allows for discrimination between the major structural variants of KIR2DS4 and KIR3DP1 without increasing the number of reactions.

Comprehensive approach to high-resolution KIR typing

Multiple studies suggest that prospective KIR typing may be beneficial for the outcome of bone marrow transplants, but to date no practical high-resolution KIR typing system has been developed. Here we propose a comprehensive high-resolution typing approach that provides allele level KIR typing. Based on the low-resolution typing obtained by SSO, the 14 KIR loci are divided in groups according to the level of polymorphism in exons coding for extracellular Ig-like domains and cytoplasmic tail. The first group is typed by sequence-specific oligonucleotide only; the second is typed by sequence-based typing (SBT) based on the amplification of a fragment coding the Ig-like domains; and the third is typed by SBT based on amplification of a fragment coding the cytoplasmic tail. SBT for the fourth group includes both the Ig-like and cytoplasmic domains. Because of a considerable number of polymorphisms scattered throughout all nine KIR exons, SBT results may still produce a number of ambiguities, which can be resolved by sequence-specific primers. This combined high-resolution approach was applied to the complete KIR typing of 205 Caucasian hematopoietic stem cell donors in support of the National Marrow Donor Program High-resolution KIR Typing Pilot Project. High-resolution typing of several KIR loci produced numerous novel alleles, whereas some loci demonstrated very limited polymorphism. Several of the novel alleles appeared in more than four donors, suggesting that these alleles are not rare. Our results showed that the comprehensive KIR typing approach presented here provides the balance of high-resolution typing and cost effectiveness.

Amazonian Amerindians exhibit high variability of KIR profiles

Natural killer cell immunoglobulin-like receptors (KIRs) mediate cell lysis through the recognition of human leukocyte antigen class I complexes in target cells, playing an important role in innate immune response. In this context, disease-based selective pressures could be relevant, leaving signatures detected by population studies. However, most population studies on KIR variability have focused on Europe and Asia, while Americas, Oceania, and Africa remain poorly studied. The aim of this study was to analyze the variability of KIR genes in Amerindian tribes from the Amazon region to infer about their evolutionary history. KIR profiles were estimated in 40 individuals from six Amazonian Amerindian tribes using single specific primer polymerase chain reaction. Twenty-five different profiles were identified, and surprisingly, the haplogroup A frequency was the lowest observed in human populations (16%). Results showed also that KIR variability was higher in this group in contrast to Venezuelan Amerindians. Principal components analysis evidenced that Amerindians formed a separated group from other worldwide populations and showed a higher intraethnic differentiation in comparison to other ethnic groups. Such pattern may reflect small effective size and intense genetic drift. However, because of the role of KIR in immune response, selective pressures cannot be entirely ruled out.

Diversity of the killer cell immunoglobulin-like receptor gene KIR2DS4 in the Chinese population

Human killer cell immunoglobulin-like receptors are a subfamily of the immunoglobulin superfamily, which map to the leukocyte receptor complex on chromosome 19. Here, we established polymerase chain reaction-sequence-based typing (PCR-SBT) procedures to identify alleles of the KIR2DS4 gene. The method was designed around the specific amplification of exons 4-5 of the KIR2DS4 gene. Genomic DNA from 105 healthy unrelated Chinese Han individuals were typed for the KIR2DS4 alleles. Each sample was assigned the KIR2DS4 alleles combination, consistent with the pairwise combinations of sequences of all the known KIR2DS4 alleles. We observed eleven different genotypes and four KIR2DS4 alleles in the population, with the KIR2DS4*00101 having the highest frequency, 0.576, and also confirmed the new KIR2DS4*007 allele. Our data demonstrated that the established PCR-SBT method for the KIR2DS4 allele typing was reliable.