HLA class I alleles of archaic origin may have been retained in modern humans because they provide immunity against diseases to which archaic humans had evolved resistance. According to this model, archaic introgressed alleles were somehow distinct from those that evolved in African populations. Here, we show that HLA-B∗73:01, a rare allotype with putative archaic origins, has a relatively rare peptide binding motif with an unusually long-tailed peptide length distribution. We also find that HLA-B∗73:01 combines a restricted and unique peptidome with high-cell surface expression, characteristics that make it well-suited to combat one or a number of closely related pathogens. Furthermore, a crystal structure of HLA-B∗73:01 in complex with KIR2DL2 highlights differences from previously solved structures with HLA-C molecules. These molecular characteristics distinguish HLA-B∗73:01 from other HLA class I alleles previously investigated and may have provided early modern human migrants that inherited this allele with a selective advantage as they colonized Europe and Asia.

Recent in situ multiplexed profiling techniques provide insight into microenvironment formation, maintenance, and transformation through a lens of localized cellular phenotype distribution. In this article, we introduce a method for recovering signatures of microenvironments from such data. We use topic models to identify characteristic cell types overrepresented in neighborhoods that serve as proxies for microenvironment. Furthermore, by assuming spatial coherence among neighboring microenvironments our model limits the number of parameters that need to be learned and permits data-driven decisions about the size of cellular neighborhoods. We apply this method to uncover anatomically known structures in mouse spleen-identifying distinct population of spleen B cells that are defined by their characteristic neighborhoods. Next, we apply the method to a dataset of triple-negative breast cancer tumors from 41 patients to study the structure of tumor-immune boundary. We uncover previously reported tumor-immune microenvironment near the tumor-immune boundary enriched for immune cells with high Indoleamine-pyrrole 2,3-dioxygenase (IDO) and Programmed death ligand 1 (PD-L1) and a novel, immunosuppressed, microenvironment-enriched for cells expressing CD45 and FoxP3.

Genetic diversity in human natural killer (NK) cell receptors is linked to resistance and susceptibility to many diseases. Here, we tested the effect of this diversity on the nanoscale organization of killer cell immunoglobulin-like receptors (KIRs). Using superresolution microscopy, we found that inhibitory KIRs encoded by different genes and alleles were organized differently at the surface of primary human NK cells. KIRs that were found at low abundance assembled into smaller clusters than those formed by KIRs that were more highly abundant, and at low abundance, there was a greater proportion of KIRs in clusters. Upon receptor triggering, a structured interface called the immune synapse assembles, which facilitates signal integration and controls NK cell responses. Here, triggering of low-abundance receptors resulted in less phosphorylation of the downstream phosphatase SHP-1 but more phosphorylation of the adaptor protein Crk than did triggering of high-abundance receptors. In cells with greater KIR abundance, SHP-1 dephosphorylated Crk, which potentiated NK cell spreading during activation. Thus, genetic variation modulates both the abundance and nanoscale organization of inhibitory KIRs. That is, as well as the number of receptors at the cell surface varying with genotype, the way in which these receptors are organized in the membrane also varies. Essentially, a change in the average surface abundance of a protein at the cell surface is a coarse descriptor entwined with changes in local nanoscale clustering. Together, our data indicate that genetic diversity in inhibitory KIRs affects membrane-proximal signaling and, unexpectedly, the formation of activating immune synapses.

The immune system comprises a complex network of specialized cells that protects against infection, eliminates cancerous cells, and regulates tissue repair, thus serving a critical role in homeostasis, health span, and life span. The subterranean-dwelling naked mole-rat (NM-R; Heterocephalus glaber) exhibits prolonged life span relative to its body size, is unusually cancer resistant, and manifests few physiological or molecular changes with advancing age. We therefore hypothesized that the immune system of NM-Rs evolved unique features that confer enhanced cancer immunosurveillance and prevent the age-associated decline in homeostasis. Using single-cell RNA-sequencing (scRNA-seq) we mapped the immune system of the NM-R and compared it to that of the short-lived, cancer-prone mouse. In contrast to the mouse, we find that the NM-R immune system is characterized by a high myeloid-to-lymphoid cell ratio that includes a novel, lipopolysaccharide (LPS)-responsive, granulocyte cell subset. Surprisingly, we also find that NM-Rs lack canonical natural killer (NK) cells. Our comparative genomics analyses support this finding, showing that the NM-R genome lacks an expanded gene family that controls NK cell function in several other species. Furthermore, we reconstructed the evolutionary history that likely led to this genomic state. The NM-R thus challenges our current understanding of mammalian immunity, favoring an atypical, myeloid-biased mode of innate immunosurveillance, which may contribute to its remarkable health span.

Killer-cell Ig-like receptor (KIR) genes are inherited as haplotypes. They are expressed by NK cells and linked to outcomes of infectious diseases and pregnancy in humans. Understanding how genotype relates to phenotype is difficult because of the extensive diversity of the KIR family. Indeed, high-resolution KIR genotyping and phenotyping in single NK cells in the context of disease association is lacking. In this article, we describe a new method to separate NK cells expressing allotypes of the KIR2DL1 gene carried by the KIR A haplotype (KIR2DL1A) from those expressing KIR2DL1 alleles carried by the KIR B haplotype (KIR2DL1B). We find that in KIR AB heterozygous individuals, different KIR2DL1 allotypes can be detected in both peripheral blood and uterine NK cells. Using this new method, we demonstrate that both blood and uterine NK cells codominantly express KIR2DL1A and KIR2DL1B allotypes but with a predominance of KIR2DL1A variants, which associate with enhanced NK cell function. In a case-control study of pre-eclampsia, we show that KIR2DL1A, not KIR2DL1B, associates with increased disease risk. This method will facilitate our understanding of how individual KIR2DL1 allelic variants affect NK cell function and contribute to disease risk.

Immunophenotypic differences between closely related human leukocyte antigen (HLA) alleles have been associated with divergent clinical outcomes in infection, autoimmunity, transplantation and drug hypersensitivity. Here we explore the impact of micropolymorphism on peptide antigen presentation by three closely related HLA molecules, HLA-B*57:01, HLA-B*57:03 and HLA-B*58:01, that are differentially associated with the HIV elite controller phenotype and adverse drug reactions. For each allotype, we mine HLA ligand data sets derived from the same parental cell proteome to define qualitative differences in peptide presentation using classical peptide binding motifs and an unbiased statistical approach. The peptide repertoires show marked qualitative overlap, with 982 peptides presented by all allomorphs. However, differences in peptide abundance, HLA-peptide stability, and HLA-bound conformation demonstrate that HLA micropolymorphism impacts more than simply the range of peptide ligands. These differences provide grounds for distinct immune reactivity and insights into the capacity of micropolymorphism to diversify immune outcomes.

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.

HLA-B∗46:01 was formed by an intergenic mini-conversion, between HLA-B∗15:01 and HLA-C∗01:02, in Southeast Asia during the last 50,000 years, and it has since become the most common HLA-B allele in the region. A functional effect of the mini-conversion was introduction of the C1 epitope into HLA-B∗46:01, making it an exceptional HLA-B allotype that is recognized by the C1-specific natural killer (NK) cell receptor KIR2DL3. High-resolution mass spectrometry showed that HLA-B∗46:01 has a low-diversity peptidome that is distinct from those of its parents. A minority (21%) of HLA-B∗46:01 peptides, with common C-terminal characteristics, form ligands for KIR2DL3. The HLA-B∗46:01 peptidome is predicted to be enriched for peptide antigens derived from Mycobacterium leprae. Overall, the results indicate that the distinctive peptidome and functions of HLA-B∗46:01 provide carriers with resistance to leprosy, which drove its rapid rise in frequency in Southeast Asia.

The immune and reproductive functions of human NK cells are regulated by interactions of the C1 and C2 epitopes of HLA-C with C1-specific and C2-specific lineage III killer cell Ig-like receptors (KIR). This rapidly evolving and diverse system of ligands and receptors is restricted to humans and great apes. In this context, the orangutan has particular relevance because it represents an evolutionary intermediate, one having the C1 epitope and corresponding KIR but lacking the C2 epitope. Through a combination of direct sequencing, KIR genotyping, and data mining from the Great Ape Genome Project, we characterized the KIR alleles and haplotypes for panels of 10 Bornean orangutans and 19 Sumatran orangutans. The orangutan KIR haplotypes have between 5 and 10 KIR genes. The seven orangutan lineage III KIR genes all locate to the centromeric region of the KIR locus, whereas their human counterparts also populate the telomeric region. One lineage III KIR gene is Bornean specific, one is Sumatran specific, and five are shared. Of 12 KIR gene-content haplotypes, 5 are Bornean specific, 5 are Sumatran specific, and 2 are shared. The haplotypes have different combinations of genes encoding activating and inhibitory C1 receptors that can be of higher or lower affinity. All haplotypes encode an inhibitory C1 receptor, but only some haplotypes encode an activating C1 receptor. Of 130 KIR alleles, 55 are Bornean specific, 65 are Sumatran specific, and 10 are shared.

Introduction: KIR2DS5 is an activating human NK cell receptor of lineage III KIR. These include both inhibitory KIR2DL1, 2 and 3 and activating KIR2DS1 that recognize either the C1 or C2 epitope of HLA-C. In Europeans KIR2DS5 is essentially monomorphic, with KIR2DS5*002 being predominant. Pioneering investigations showed that KIR2DS5*002 has activating potential, but cannot recognize HLA-A, -B, or -C. Subsequent studies have shown that KIR2DS5 is highly polymorphic in Africans, and that KIR2DS5*006 protects pregnant Ugandan women from preeclampsia. Because inhibitory C2-specific KIR2DL1 correlates with preeclampsia, whereas activating C2-specific KIR2DS1 protects, this association pointed to KIR2DS5*006 being an activating C2-specific receptor. To test this hypothesis we made KIR-Fc fusion proteins from all ten KIR2DS5 allotypes and tested their binding to a representative set of HLA-A, -B and -C allotypes. Results: Six African-specific KIR2DS5 bound to C2(+)HLA-C but not to other HLA class I. Their avidity for C2 is similar to 20% that of C2-specific KIR2DL1 and similar to 40% that of C2-specific KIR2DS1. Among the African C2 receptors is KIR2DS5*006, which protected a cohort of pregnant Ugandans from pre-eclampsia. Three African KIR2DS5 allotypes and KIR2DS5*002, bound no HLA-A, -B or -C. As a group the C2-binding KIR2DS5 allotypes protect against pre-eclampsia compared to the nonbinding KIR2DS5 allotypes. Natural substitutions that contribute to loss or reduction of C2 receptor function are at positions 127, 158, and 176 in the D2 domain. Conclusions: KIR2DS5*005 has the KIR2DS5 consensus sequence, is the only allele found at both centromeric and telomeric locations of KIR2DS5, and is likely the common ancestor of all KIR2DS5 alleles. That KIR2DS5*005 has C2 receptor activity, points to KIR2DS5*002, and other allotypes lacking C2 receptor function, being products of attenuation, a characteristic feature of most KIR B haplotype genes. Alleles encoding attenuated and active KIR2DS5 are present in both centromeric and telomeric locations.

Natural killer (NK) cells are fast-acting and versatile lymphocytes that are critical effectors of innate immunity, adaptive immunity, and placental development. Controlling NK cell function are the interactions between killer-cell immunoglobulin-like receptors (KIRs) and their HLA-A, HLA-B and HLA-C ligands. Due to the extensive polymorphism of both KIR and HLA class I, these interactions are highly diversified and specific combinations correlate with protection or susceptibility to a range of infectious, autoimmune, and reproductive disorders. Evolutionary, genetic, and functional studies are consistent with the interactions between KIR and HLA-C being the dominant control mechanism of human NK cells. In addition to their recognition of the C1 and C2 epitopes, increasing evidence points to KIR having a previously unrecognized selectivity for the peptide presented by HLA-C. This selectivity appears to be a conserved feature of activating KIR and may partly explain the slow progress made in identifying their HLA class I ligands. The peptide selectivity of KIR allows NK cells to respond, not only to changes in the surface expression of HLA-C, but also to the more subtle changes in the HLA-C peptidome, such as occur during viral infection and malignant transformation. Here, we review recent advances in understanding of human-specific KIR evolution and how the inhibitory and activating HLA-C receptors allow NK cells to respond to healthy cells, diseased cells, and the semi-allogeneic cells of the fetus.

KIR2DP1 is an inactive member of the human lineage III KIR family, which includes all HLA-C-specific receptor genes. The lethal, and only, defect in KIR2DP1 is a nucleotide deletion in codon 88. Fixed in modern humans, the deletion is also in archaic human genomes. KIR2DP1 is polymorphic, with dimorphism at specificity-determining position 44. By repairing the deletion, we resurrected 11 alleles of KIR2DP1(F), the functional antecedent of KIR2DP1. We demonstrate how K44-KIR2DP1(F) with lysine 44 recognized C1(+) HLA-C, whereas T44-KIR2DP1(F) recognized C2(+)HLA-C. Dimorphisms at 12 other KIR2DP1(F) residues modulate receptor avidity or signaling. KIR2DP1 and KIR2DL1 are neighbors in the centromeric KIR region and are in tight linkage disequilibrium. Like KIR2DL1, KIR2DP1 contributed to CenA and CenB KIR haplotype differences. Encoded on CenA, C1-specific K44-KIR2DP1(F) were stronger receptors than the attenuated C2-specific T44-KIR2DP1(F) encoded on CenB. The last common ancestor of humans and chimpanzees had diverse lineage III KIR that passed on to chimpanzees but not to humans. Early humans inherited activating KIR2DS4 and an inhibitory lineage III KIR, likely encoding a C1-specific receptor. The latter spawned the modern family of HLA-C receptors. KIR2DP1(F) has properties consistent with KIR2DP1(F) having been the founder gene. The first KIR2DP1(F) alleles encoded K44-C1 receptors; subsequently KIR2DP1(F) alleles encoding T44-C2 receptors evolved. The emergence of dedicated KIR2DL2/3 and KIR2DL1 genes encoding C1 and C2 receptors, respectively, could have led to obsolescence of KIR2DP1(F). Alternatively, pathogen subversion caused its demise. Preservation of KIR2DP1(F) functional polymorphism was a side effect of fixation of the deletion in KIR2DP1(F) by micro gene conversion.

Individuals heterozygous for HLA haplotypes have increased NK cell diversity because of distinct education pathways.

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.

Shaping natural killer (NK) cell functions in human immunity and reproduction are diverse killer cell immunoglobulin-like receptors (KIRs) that recognize polymorphic MHC class I determinants. A survey of placental mammals suggests that KIRs serve as variable NK cell receptors only in certain primates and artiodactyls. Divergence of the functional and variable KIRs in primates and artiodactyls predates placental reproduction. Among artiodactyls, cattle but not pigs have diverse KIRs. Catarrhine (humans, apes, and Old World monkeys) and platyrrhine (New World monkeys) primates, but not prosimians, have diverse KIRs. Platyrrhine and catarrhine systems of KIR and MHC class I are highly diverged, but within the catarrhines, a stepwise co-evolution of MHC class I and KIR is discerned. In Old World monkeys, diversification focuses on MHC-A and MHC-B and their cognate lineage II KIR. With evolution of C1-bearing MHC-C from MHC-B, as informed by orangutan, the focus changes to MHC-C and its cognate lineage III KIR. Evolution of C2 from C1 and fixation of MHC-C drove further elaboration of MHC-C-specific KIR, as exemplified by chimpanzee. In humans, the evolutionary trajectory changes again. Emerging from reorganization of the KIR locus and selective attenuation of KIR avidity for MHC class I are the functionally distinctive KIR A and KIR B haplotypes.

The human killer cell Ig-like receptor (KIR) locus comprises two groups of KIR haplotypes, termed A and B. These are present in all human populations but with different relative frequencies, suggesting they have different functional properties that underlie their balancing selection. We studied the genomic organization and functional properties of the alleles of the inhibitory and activating HLA-C receptors encoded by KIR haplotypes. Because every HLA-C allotype functions as a ligand for KIR, the interactions between KIR and HLA-C dominate the HLA class I-mediated regulation of human NK cells. The C2 epitope is recognized by inhibitory KIR2DL1 and activating KIR2DS1, whereas the C1 epitope is recognized by inhibitory KIR2DL2 and KIR2DL3. This study shows that the KIR2DL1, KIR2DS1, and KIR2DL2/3 alleles form distinctive phylogenetic clades that associate with specific KIR haplotypes. KIR A haplotypes are characterized by KIR2DL1 alleles that encode strong inhibitory C2 receptors and KIR2DL3 alleles encoding weak inhibitory C1 receptors. In striking contrast, KIR B haplotypes are characterized by KIR2DL1 alleles that encode weak inhibitory C2 receptors and KIR2DL2 alleles encoding strong inhibitory C1 receptors. The wide-ranging properties of KIR allotypes arise from substitutions throughout the KIR molecule. Such substitutions can influence cell surface expression, as well as the avidity and specificity for HLA-C ligands. Consistent with the crucial role of inhibitory HLA-C receptors in self-recognition, as well as NK cell education and response, most KIR haplotypes have both a functional C1 and C2 receptor, despite the considerable variation that occurs in ligand recognition and surface expression.

Chimpanzees have orthologs of the six fixed, functional human MHC class I genes. But, in addition, the chimpanzee has a seventh functional gene, Patr-AL, which is not polymorphic but contributes substantially to population diversity by its presence on only 50% of MHC haplotypes. The ancestral AL gene emerged long before the separation of human and chimpanzee ancestors and then subsequently and specifically lost function during human evolution, but was maintained in chimpanzees. Patr-AL is an alloantigen that participates in negative and positive selection of the T cell repertoire. The three-dimensional structure and the peptide-binding repertoire of Patr-AL and HLA-A*02 are surprisingly similar. In contrast, the expression of these two molecules is very different, as shown using specific mAbs and polyclonal Abs made against Patr-AL. Peripheral blood cells and B cell lines express low levels of Patr-AL at the cell surface. Higher levels are seen for 221-cell transfectants expressing Patr-AL, but in these cells a large majority of Patr-AL molecules are retained in the early compartments of the secretory pathway: mainly the endoplasmic reticulum, but also cis-Golgi. Replacing the cytoplasmic tail of Patr-AL with that of HLA-A*02 increased the cell-surface expression of Patr-AL substantially. Four substitutions distinguish the Patr-AL and HLA-A*02 cytoplasmic tails. Systematic mutagenesis showed that each substitution contributes changes in cell-surface expression. The combination of residues present in Patr-AL appears unique, but each individual residue is present in other primate MHC class I molecules, notably MHC-E, the most ancient of the functional human MHC class I molecules.