Paul J. Norman

Synergistic Polymorphism at Two Positions Distal to the Ligand-Binding Site Makes KIR2DL2 a Stronger Receptor for HLA-C Than KIR2DL3

Interactions between HLA-C ligands and inhibitory killer cell Ig-like receptors (KIR) control the development and response of human NK cells. This regulatory mechanism is usually described by mutually exclusive interactions of KIR2DL1 with C2 having lysine 80, and KIR2DL2/3 with C1 having asparagine 80. Consistent with this simple rule, we found from functional analysis and binding assays to 93 HLA-A, HLA-B, and HLA-C isoforms that KIR2DL1*003 bound all C2, and only C2, allotypes. The allotypically related KIR2DL2*001 and KIR2DL3*001 interacted with all C1, but they violated the simple rule through interactions with several C2 allotypes, notably Cw*0501 and Cw*0202, and two HLA-B allotypes (B*4601 and B*7301) that share polymorphisms with HLA-C. Although the specificities of the “cross-reactions” were similar for KIR2DL2*001 and KIR2DL3*001, they were stronger for KIR2DL2*001, as were the reactions with C1. Mutagenesis explored the avidity difference between KIR2DL2*001 and KIR2DL3*001. Recombinant mutants mapped the difference to the Ig-like domains, where site-directed mutagenesis showed that the combination, but not the individual substitutions, of arginine for proline 16 in D1 and cysteine for arginine 148 in D2 made KIR2DL2*001 a stronger receptor than KIR2DL3*001. Neither residue 16 or 148 is part of, or near to, the ligand-binding site. Instead, their juxtaposition near the flexible hinge between D1 and D2 suggests that their polymorphisms affect the ligand-binding site by changing the hinge angle and the relative orientation of the two domains. This study demonstrates how allelic polymorphism at sites distal to the ligand-binding site of KIR2DL2/3 has diversified this receptor’s interactions with HLA-C.

The Shaping of Modern Human Immune Systems by Multiregional Admixture with Archaic Humans

Viral defense and embryo implantation mechanisms have been shaped by contributions from Neandertal and Denisovan genes. Whole genome comparisons identified introgression from archaic to modern humans. Our analysis of highly polymorphic human leukocyte antigen (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.

Hunter-gatherer genomic diversity suggests a southern African origin for modern humans

Africa is inferred to be the continent of origin for all modern human populations, but the details of human prehistory and evolution in Africa remain largely obscure owing to the complex histories of hundreds of distinct populations. We present data for more than 580,000 SNPs for several hunter-gatherer populations: the Hadza and Sandawe of Tanzania, and the ≠Khomani Bushmen of South Africa, including speakers of the nearly extinct N|u language. We find that African hunter-gatherer populations today remain highly differentiated, encompassing major components of variation that are not found in other African populations. Hunter-gatherer populations also tend to have the lowest levels of genome-wide linkage disequilibrium among 27 African populations. We analyzed geographic patterns of linkage disequilibrium and population differentiation, as measured by FST, in Africa. The observed patterns are consistent with an origin of modern humans in southern Africa rather than eastern Africa, as is generally assumed. Additionally, genetic variation in African hunter-gatherer populations has been significantly affected by interaction with farmers and herders over the past 5,000 y, through both severe population bottlenecks and sex-biased migration. However, African hunter-gatherer populations continue to maintain the highest levels of genetic diversity in the world.

Genetic and Environmental Determinants of Human NK Cell Diversity Revealed by Mass Cytometry

Both genetics and environment contribute to human NK cell diversity. NK Cell Nature Versus Nurture Natural killer (NK) cells were first discovered because of their ability to kill tumor cells without any previous exposure. However, this population is actually quite heterogeneous: Different subgroups of NK cells express different combinations of activating and inhibiting receptors that govern their specificity. Now, Horowitz et al. use mass cytometry to examine NK cell diversity in humans. The authors examined 35 parameters simultaneously in 5 sets of monozygotic twins as well as 12 unrelated donors. They found up to 30,000 phenotypic NK cell populations in a given individual. What’s more, by comparing the twins versus unrelated donors, they determined that although genetics primarily determined inhibitory receptor expression, activating receptors were controlled by the environment. These data suggest that inhibitory receptors may contribute more to NK cell self-tolerance, whereas activating receptors may guide response to pathogens and tumors. Natural killer (NK) cells play critical roles in immune defense and reproduction, yet remain the most poorly understood major lymphocyte population. Because their activation is controlled by a variety of combinatorially expressed activating and inhibitory receptors, NK cell diversity and function are closely linked. To provide an unprecedented understanding of NK cell repertoire diversity, we used mass cytometry to simultaneously analyze 37 parameters, including 28 NK cell receptors, on peripheral blood NK cells from 5 sets of monozygotic twins and 12 unrelated donors of defined human leukocyte antigen (HLA) and killer cell immunoglobulin-like receptor (KIR) genotype. This analysis revealed a remarkable degree of NK cell diversity, with an estimated 6000 to 30,000 phenotypic populations within an individual and >100,000 phenotypes in the donor panel. Genetics largely determined inhibitory receptor expression, whereas activation receptor expression was heavily environmentally influenced. Therefore, NK cells may maintain self-tolerance through strictly regulated expression of inhibitory receptors while using adaptable expression patterns of activating and costimulatory receptors to respond to pathogens and tumors. These findings further suggest the possibility that discrete NK cell subpopulations could be harnessed for immunotherapeutic strategies in the settings of infection, reproduction, and transplantation.

Natural killer cell immunoglobulin-like receptor (KIR) locus profiles in African and South Asian populations

Natural killer (NK) and some T cells express killer cell immunoglobulin-like receptors (KIRs), which interact with HLA class I expressed by target cells and consequently regulate cytolytic activity. The number of KIR loci can vary and so a range of genetic profiles is observed. We have determined the KIR genetic profiles from one African (n = 62) and two South Asian (n = 108, n = 78) populations. Several of the KIRs are present at significantly different frequencies between the two major ethnic groups (eg KIR2DS4 gene frequency 0.82 African, 0.47 S Asian. Pc < 1 × 10−6) and this is due to uneven distribution of two KIR haplotype families ‘A’ and ‘B’. All three populations described here displayed a greater degree of diversity of KIR genetic profiles than other populations investigated, which indicates further complexity of underlying haplotypes; in this respect we describe two individuals who appear homozygous for a large deletion including the previously ubiquitous 2DL4. We have also reanalysed three populations that we studied previously, for the presence of a KIR which is now known to be an indicator of the ‘B’ haplotype. South Asians had the highest overall frequencies of all KIR loci characteristic of ‘B’ haplotypes (Pc < 0.0001 to < 0.004). Furthermore, gene frequency independent deviances in the linkage disequilibrium were apparent between populations.

Genomic evidence for the Pleistocene and recent population history of Native Americans

Genetic history of Native Americans Several theories have been put forth as to the origin and timing of when Native American ancestors entered the Americas. To clarify this controversy, Raghavan et al. examined the genomic variation among ancient and modern individuals from Asia and the Americas. There is no evidence for multiple waves of entry or recurrent gene flow with Asians in northern populations. The earliest migrations occurred no earlier than 23,000 years ago from Siberian ancestors. Amerindians and Athabascans originated from a single population, splitting approximately 13,000 years ago. Science, this issue 10.1126/science.aab3884 Genetic variation within ancient and extant Native American populations informs on their migration into the Americas. INTRODUCTION The consensus view on the peopling of the Americas is that ancestors of modern Native Americans entered the Americas from Siberia via the Bering Land Bridge and that this occurred at least ~14.6 thousand years ago (ka). However, the number and timing of migrations into the Americas remain controversial, with conflicting interpretations based on anatomical and genetic evidence. RATIONALE In this study, we address four major unresolved issues regarding the Pleistocene and recent population history of Native Americans: (i) the timing of their divergence from their ancestral group, (ii) the number of migrations into the Americas, (iii) whether there was ~15,000 years of isolation of ancestral Native Americans in Beringia (Beringian Incubation Model), and (iv) whether there was post-Pleistocene survival of relict populations in the Americas related to Australo-Melanesians, as suggested by apparent differences in cranial morphologies between some early (“Paleoamerican”) remains and those of more recent Native Americans. We generated 31 high-coverage modern genomes from the Americas, Siberia, and Oceania; 23 ancient genomic sequences from the Americas dating between ~0.2 and 6 ka; and SNP chip genotype data from 79 present-day individuals belonging to 28 populations from the Americas and Siberia. The above data sets were analyzed together with published modern and ancient genomic data from worldwide populations, after masking some present-day Native Americans for recent European admixture. RESULTS Using three different methods, we determined the divergence time for all Native Americans (Athabascans and Amerindians) from their Siberian ancestors to be ~20 ka, and no earlier than ~23 ka. Furthermore, we dated the divergence between Athabascans (northern Native American branch, together with northern North American Amerindians) and southern North Americans and South and Central Americans (southern Native American branch) to be ~13 ka. Similar divergence times from East Asian populations and a divergence time between the two branches that is close in age to the earliest well-established archaeological sites in the Americas suggest that the split between the branches occurred within the Americas. We additionally found that several sequenced Holocene individuals from the Americas are related to present-day populations from the same geographical regions, implying genetic continuity of ancient and modern populations in some parts of the Americas over at least the past 8500 years. Moreover, our results suggest that there has been gene flow between some Native Americans from both North and South America and groups related to East Asians and Australo-Melanesians, the latter possibly through an East Asian route that might have included ancestors of modern Aleutian Islanders. Last, using both genomic and morphometric analyses, we found that historical Native American groups such as the Pericúes and Fuego-Patagonians were not “relicts” of Paleoamericans, and hence, our results do not support an early migration of populations directly related to Australo-Melanesians into the Americas. CONCLUSION Our results provide an upper bound of ~23 ka on the initial divergence of ancestral Native Americans from their East Asian ancestors, followed by a short isolation period of no more than ~8000 years, and subsequent entrance and spread across the Americas. The data presented are consistent with a single-migration model for all Native Americans, with later gene flow from sources related to East Asians and, indirectly, Australo-Melanesians. The single wave diversified ~13 ka, likely within the Americas, giving rise to the northern and southern branches of present-day Native Americans. Population history of present-day Native Americans. The ancestors of all Native Americans entered the Americas as a single migration wave from Siberia (purple) no earlier than ~23 ka, separate from the Inuit (green), and diversified into “northern” and “southern” Native American branches ~13 ka. There is evidence of post-divergence gene flow between some Native Americans and groups related to East Asians/Inuit and Australo-Melanesians (yellow). How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative “Paleoamerican” relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.

KIR2DS4 is a product of gene conversion with KIR3DL2 that introduced specificity for HLA-A*11 while diminishing avidity for HLA-C

Human killer cell immunoglobulin-like receptors (KIRs) are distinguished by expansion of activating KIR2DS, whose ligands and functions remain poorly understood. The oldest, most prevalent KIR2DS is KIR2DS4, which is represented by a variable balance between “full-length” and “deleted” forms. We find that full-length 2DS4 is a human histocompatibility leukocyte antigen (HLA) class I receptor that binds specifically to subsets of C1+ and C2+ HLA-C and to HLA-A*11, whereas deleted 2DS4 is nonfunctional. Activation of 2DS4+ NKL cells was achieved with A*1102 as ligand, which differs from A*1101 by unique substitution of lysine 19 for glutamate, but not with A*1101 or HLA-C. Distinguishing KIR2DS4 from other KIR2DS is the proline–valine motif at positions 71–72, which is shared with KIR3DL2 and was introduced by gene conversion before separation of the human and chimpanzee lineages. Site-directed swap mutagenesis shows that these two residues are largely responsible for the unique HLA class I specificity of KIR2DS4. Determination of the crystallographic structure of KIR2DS4 shows two major differences from KIR2DL: displacement of contact loop L2 and altered bonding potential because of the substitutions at positions 71 and 72. Correlation between the worldwide distributions of functional KIR2DS4 and HLA-A*11 points to the physiological importance of their mutual interaction.

Different patterns of evolution in the centromeric and telomeric regions of group A and B haplotypes of the human killer cell Ig-like receptor locus.

The fast evolving human KIR gene family encodes variable lymphocyte receptors specific for polymorphic HLA class I determinants. Nucleotide sequences for 24 representative human KIR haplotypes were determined. With three previously defined haplotypes, this gave a set of 12 group A and 15 group B haplotypes for assessment of KIR variation. The seven gene-content haplotypes are all combinations of four centromeric and two telomeric motifs. 2DL5, 2DS5 and 2DS3 can be present in centromeric and telomeric locations. With one exception, haplotypes having identical gene content differed in their combinations of KIR alleles. Sequence diversity varied between haplotype groups and between centromeric and telomeric halves of the KIR locus. The most variable A haplotype genes are in the telomeric half, whereas the most variable genes characterizing B haplotypes are in the centromeric half. Of the highly polymorphic genes, only the 3DL3 framework gene exhibits a similar diversity when carried by A and B haplotypes. Phylogenetic analysis and divergence time estimates, point to the centromeric gene-content motifs that distinguish A and B haplotypes having emerged ∼6 million years ago, contemporaneously with the separation of human and chimpanzee ancestors. In contrast, the telomeric motifs that distinguish A and B haplotypes emerged more recently, ∼1.7 million years ago, before the emergence of Homo sapiens. Thus the centromeric and telomeric motifs that typify A and B haplotypes have likely been present throughout human evolution. The results suggest the common ancestor of A and B haplotypes combined a B-like centromeric region with an A-like telomeric region.

Reconstructing the population genetic history of the Caribbean.

The Caribbean basin is home to some of the most complex interactions in recent history among previously diverged human populations. Here, we investigate the population genetic history of this region by characterizing patterns of genome-wide variation among 330 individuals from three of the Greater Antilles (Cuba, Puerto Rico, Hispaniola), two mainland (Honduras, Colombia), and three Native South American (Yukpa, Bari, and Warao) populations. We combine these data with a unique database of genomic variation in over 3,000 individuals from diverse European, African, and Native American populations. We use local ancestry inference and tract length distributions to test different demographic scenarios for the pre- and post-colonial history of the region. We develop a novel ancestry-specific PCA (ASPCA) method to reconstruct the sub-continental origin of Native American, European, and African haplotypes from admixed genomes. We find that the most likely source of the indigenous ancestry in Caribbean islanders is a Native South American component shared among inland Amazonian tribes, Central America, and the Yucatan peninsula, suggesting extensive gene flow across the Caribbean in pre-Columbian times. We find evidence of two pulses of African migration. The first pulse—which today is reflected by shorter, older ancestry tracts—consists of a genetic component more similar to coastal West African regions involved in early stages of the trans-Atlantic slave trade. The second pulse—reflected by longer, younger tracts—is more similar to present-day West-Central African populations, supporting historical records of later transatlantic deportation. Surprisingly, we also identify a Latino-specific European component that has significantly diverged from its parental Iberian source populations, presumably as a result of small European founder population size. We demonstrate that the ancestral components in admixed genomes can be traced back to distinct sub-continental source populations with far greater resolution than previously thought, even when limited pre-Columbian Caribbean haplotypes have survived.

Human-specific evolution of killer cell immunoglobulin-like receptor recognition of major histocompatibility complex class I molecules

In placental mammals, natural killer (NK) cells are a population of lymphocytes that make unique contributions to immune defence and reproduction, functions essential for survival of individuals, populations and species. Modulating these functions are conserved and variable NK-cell receptors that recognize epitopes of major histocompatibility complex (MHC) class I molecules. In humans, for example, recognition of human leucocyte antigen (HLA)-E by the CD94:NKG2A receptor is conserved, whereas recognition of HLA-A, B and C by the killer cell immunoglobulin-like receptors (KIRs) is diversified. Competing demands of the immune and reproductive systems, and of T-cell and NK-cell immunity—combined with the segregation on different chromosomes of variable NK-cell receptors and their MHC class I ligands—drive an unusually rapid evolution that has resulted in unprecedented levels of species specificity, as first appreciated from comparison of mice and humans. Counterparts to human KIR are present only in simian primates. Observed in these species is the coevolution of KIR and the four MHC class I epitopes to which human KIR recognition is restricted. Unique to hominids is the emergence of the MHC-C locus as a supplier of specialized and superior ligands for KIR. This evolutionary trend is most highly elaborated in the chimpanzee. Unique to the human KIR locus are two groups of KIR haplotypes that are present in all human populations and subject to balancing selection. Group A KIR haplotypes resemble chimpanzee KIR haplotypes and are enriched for genes encoding KIR that bind HLA class I, whereas group B KIR haplotypes are enriched for genes encoding receptors with diminished capacity to bind HLA class I. Correlating with their balance in human populations, B haplotypes favour reproductive success, whereas A haplotypes favour successful immune defence. Evolution of the B KIR haplotypes is thus unique to the human species.