Whole genome comparisons identified introgression from archaic to modern humans. Our analysis of highly polymorphic HLA class I, vital immune system components subject to strong balancing selection, shows how modern humans acquired the HLA-B*73 allele in west Asia through admixture with archaic humans called Denisovans, a likely sister group to the Neandertals. Virtual genotyping of Denisovan and Neandertal genomes identified archaic HLA haplotypes carrying functionally distinctive alleles that have introgressed into modern Eurasian and Oceanian populations. These alleles, of which several encode unique or strong ligands for natural killer cell receptors, now represent more than half the HLA alleles of modern Eurasians and also appear to have been later introduced into Africans. Thus, adaptive introgression of archaic alleles has significantly shaped modern human immune systems.

That NK cell receptors engage fast-evolving MHC class I ligands suggests that they, too, evolve rapidly. To test this hypothesis, the structure and class I specificity of chimpanzee KIR and CD94:NKG2 receptors were determined and compared to their human counterparts. The KIR families are divergent, with only three KIR conserved between chimpanzees and humans. By contrast, CD94:NKG2 receptors are conserved. Whereas receptors for polymorphic class I are divergent, those for nonpolymorphic class I are conserved. Although chimpanzee and human NK cells exhibit identical receptor specificities for MHC-C, they are mediated by nonorthologous KIR. These results demonstrate the rapid evolution of NK cell receptor systems and imply that "catching up" with class I is not the only force driving this evolution.

KIR3DL1 and KIR3DL2 are NK cell receptors for polymorphic HLA-B and -A determinants. The proportion of NK cells that bind anti-KIR3DL1-specific Ab DX9 and their level of binding vary between individuals. To determine whether these differences are due to KIR polymorphism, we assessed KIR3D gene diversity in unrelated individuals and families. Both KIR3DL1 and KIR3DL2 are highly polymorphic genes, with KIR3DS1 segregating like an allele of KIR3DL1. A KIR haplotype lacking KIR3DL1 and KIR3DS1 was defined. The two KIR3DL1 alleles of a heterozygous donor were expressed by different, but overlapping, subsets of NK cell clones. Sequence variation in KIR3DL1 and KIR3DL2 appear distinct; recombination is more evident in KIR3DL1, and point mutation is more evident in KIR3DL2. The KIR3DL1 genotype correlates well with levels of DX9 binding by NK cells, but not with the frequency of DX9-binding cells. Different KIR3DL1 alleles determine high, low, and no binding of DX9 Ab. Consequently, heterozygotes for high and low binding KIR3DL1 alleles have distinct subpopulations of NK cells that bind DX9 at high and low levels, giving characteristic bimodal distributions in flow cytometry. The Z27 Ab gave binding patterns similar to those of DX9. Four KIR3DL1 alleles producing high DX9 binding phenotypes were distinguished from four alleles producing low or no binding phenotypes by substitution at one or more of four positions in the encoded protein: 182 and 283 in the extracellular Ig-like domains, 320 in the transmembrane region, and 373 in the cytoplasmic tail.

The physiological functions of natural killer (NK) cells in human immunity and reproduction depend upon diverse interactions between killer cell immunoglobulin-like receptors (KIRs) and their HLA class I ligands: HLA-A, HLA-B, and HLA-C. The genomic regions containing the KIR and HLA class I genes are unlinked, structurally complex, and highly polymorphic. They are also strongly associated with a wide spectrum of diseases, including infections, autoimmune disorders, cancers, and pregnancy disorders, as well as the efficacy of transplantation and other immunotherapies. To facilitate study of these extraordinary genes, we developed a method that captures, sequences, and analyzes the 13 KIR genes and HLA-A, HLA-B, and HLA-C from genomic DNA. We also devised a bioinformatics pipeline that attributes sequencing reads to specific KIR genes, determines copy number by read depth, and calls high-resolution genotypes for each KIR gene. We validated this method by using DNA from well-characterized cell lines, comparing it to established methods of HLA and KIR genotyping, and determining KIR genotypes from 1000 Genomes sequence data. This identified 116 previously uncharacterized KIR alleles, which were all demonstrated to be authentic by sequencing from source DNA via standard methods. Analysis of just two KIR genes showed that 22% of the 1000 Genomes individuals have a previously uncharacterized allele or a structural variant. The method we describe is suited to the large-scale analyses that are needed for characterizing human populations and defining the precise HLA and KIR factors associated with disease. The methods are applicable to other highly polymorphic genes.

Genomic DNA from a panel of 41 healthy unrelated Japanese individuals was typed for the presence or absence of 16 KIR genes and pseudogenes. Only eight different KIR genotypes were found. Group A haplotypes outnumbered group B haplotypes in frequency by approximately 3:1, with individuals having two group A haplotypes accounting for 56% of the panel. The frequency of A haplotypes in the Japanese is higher than that observed in other populations. Flow cytometric comparison of KIR expression in 19 panel members showed considerable diversity in NK cell repertoire, which was also seen within the group of individuals having two A haplotypes. This diversity is likely due to allelic polymorphism in expressed genes of the A haplotype. In comparison to other populations, the Japanese appear less heterogeneous in KIR genotype as assessed by gene content.