Genetics of Natural Killer Cells in Human Health, Disease, and Survival.

Natural killer (NK) cells have vital functions in human immunity and reproduction. In the innate and adaptive immune responses to infection, particularly by viruses, NK cells respond by secreting inflammatory cytokines and killing infected cells. In reproduction, NK cells are critical for genesis of the placenta, the organ that controls the supply of oxygen and nutrients to the growing fetus. Controlling NK cell functions are interactions of HLA class I with inhibitory NK cell receptors. First evolved was the conserved interaction of HLA-E with CD94:NKG2A; later established were diverse interactions of HLA-A, -B, and -C with killer cell immunoglobulin-like receptors. Characterizing the latter interactions is rapid evolution, which distinguishes human populations and all species of higher primate. Driving this evolution are the different and competing selections imposed by pathogens on NK cell-mediated immunity and by the constraints of human reproduction on NK cell-mediated placentation. Promoting rapid evolution is independent segregation of polymorphic receptors and ligands throughout human populations.

The CD56-CD16+ NK cell subset in chronic infections.

Long-term human diseases can shape the immune system, and natural killer (NK) cells have been documented to differentiate into distinct subsets specifically associated with chronic virus infections. One of these subsets found in large frequencies in HIV-1 are the CD56-CD16+ NK cells, and this population's association with chronic virus infections is the subject of this review. Human NK cells are classically defined by CD56 expression, yet increasing evidence supports the NK cell status of the CD56-CD16+ subset which we discuss herein. We then discuss the evidence linking CD56-CD16+ NK cells to chronic virus infections, and the potential immunological pathways that are altered by long-term infection that could be inducing the population's differentiation. An important aspect of NK cell regulation is their interaction with human leukocyte antigen (HLA) class-I molecules, and we highlight work that indicates both virus and genetic-mediated variations in HLA expression that have been linked to CD56-CD16+ NK cell frequencies. Finally, we offer a perspective on CD56-CD16+ NK cell function, taking into account recent work that implies the subset is comparable to CD56+CD16+ NK cell functionality in antibody-dependent cell cytotoxicity response, and the definition of CD56-CD16+ NK cell subpopulations with varying degranulation capacity against target cells.

Following Transplantation for Acute Myelogenous Leukemia, Donor KIR Cen B02 Better Protects against Relapse than KIR Cen B01.

In the treatment of acute myelogenous leukemia with allogeneic hematopoietic cell transplantation, we previously demonstrated that there is a greater protection from relapse of leukemia when the hematopoietic cell transplantation donor has either the Cen B/B KIR genotype or a genotype having two or more KIR B gene segments. In those earlier analyses, KIR genotyping could only be assessed at the low resolution of gene presence or absence. To give the analysis greater depth, we developed high-resolution KIR sequence-based typing that defines all the KIR alleles and distinguishes the expressed alleles from those that are not expressed. We now describe and analyze high-resolution KIR genotypes for 890 donors of this human transplant cohort. Cen B01 and Cen B02 are the common CenB haplotypes, with Cen B02 having evolved from Cen B01 by deletion of the KIR2DL5, 2DS3/5, 2DP1, and 2DL1 genes. We observed a consistent trend for Cen B02 to provide stronger protection against relapse than Cen B01 This correlation indicates that protection depends on the donor having inhibitory KIR2DL2 and/or activating KIR2DS2, and is enhanced by the donor lacking inhibitory KIR2DL1, 2DL3, and 3DL1. High-resolution KIR typing has allowed us to compare the strength of the interactions between the recipient's HLA class I and the KIR expressed by the donor-derived NK cells and T cells, but no clinically significant interactions were observed. The trend observed between donor Cen B02 and reduced relapse of leukemia points to the value of studying ever larger transplant cohorts.

Adaptive Admixture of HLA Class I Allotypes Enhanced Genetically Determined Strength of Natural Killer Cells in East Asians.

Human natural killer (NK) cells are essential for controlling infection, cancer, and fetal development. NK cell functions are modulated by interactions between polymorphic inhibitory killer cell immunoglobulin-like receptors (KIR) and polymorphic HLA-A, -B, and -C ligands expressed on tissue cells. All HLA-C alleles encode a KIR ligand and contribute to reproduction and immunity. In contrast, only some HLA-A and -B alleles encode KIR ligands and they focus on immunity. By high-resolution analysis of KIR and HLA-A, -B, and -C genes, we show that the Chinese Southern Han (CHS) are significantly enriched for interactions between inhibitory KIR and HLA-A and -B. This enrichment has had substantial input through population admixture with neighboring populations, who contributed HLA class I haplotypes expressing the KIR ligands B*46:01 and B*58:01, which subsequently rose to high frequency by natural selection. Consequently, over 80% of Southern Han HLA haplotypes encode more than one KIR ligand. Complementing the high number of KIR ligands, the CHS KIR locus combines a high frequency of genes expressing potent inhibitory KIR, with a low frequency of those expressing activating KIR. The Southern Han centromeric KIR region encodes strong, conserved, inhibitory HLA-C-specific receptors, and the telomeric region provides a high number and diversity of inhibitory HLA-A and -B-specific receptors. In all these characteristics, the CHS represent other East Asians, whose NK cell repertoires are thus enhanced in quantity, diversity, and effector strength, likely augmenting resistance to endemic viral infections.

Natural LILRB1 D1-D2 variants show frequency differences in populations and bind to HLA class I with various avidities.

Leukocyte immunoglobulin-like receptor B1 (LILRB1) is widely expressed on various immune cells and the engagement of LILRB1 to HLA class I and pathogen-derived proteins can modulate the immune response. In the current study, 108 LILRB1 alleles were identified by screening the LILRB1 locus from the 1000 Genomes Phase 3 database. Forty-six alleles that occurred in three or more individuals encode 28 LILRB1 allotypes, and the inferred LILRB1 allotypes were then grouped into 9 LILRB1 D1-D2 variants for further analysis. We found that variants 1, 2, and 3 represent the three most frequent LILRB1 D1-D2 variants and the nine variants show frequency differences in populations. The binding assay demonstrated that variant 1 bound to HLA class I with the highest avidity, and all tested LILRB1 D1-D2 variants bound to HLA-C with lower avidity than to HLA-A and -B. Locus-specific polymorphisms at positions 183, 189, and 268 in HLA class I and dimorphisms in HLA-A (positions 207 and 253) and in HLA-B (position 194) affect their binding to LILRB1. Notably, the electrostatic interaction plays a critical role in the binding of LILRB1 to HLA class I as revealed by electrostatic analysis and by comparison of different binding avidities caused by polymorphisms at positions 72 and 103 of LILRB1. In this paper, we present a comprehensive study of the population genetics and binding abilities of LILRB1. The data will help us better understand the LILRB1-related diversity of the immune system and lay a foundation for functional studies.

Diversity of KIR, HLA Class I, and Their Interactions in Seven Populations of Sub-Saharan Africans.

HLA class I and KIR sequences were determined for Dogon, Fulani, and Baka populations of western Africa, Mbuti of central Africa, and Datooga, Iraqw, and Hadza of eastern Africa. Study of 162 individuals identified 134 HLA class I alleles (41 HLA-A, 60 HLA-B, and 33 HLA-C). Common to all populations are three HLA-C alleles (C1+C*07:01, C1+C*07:02, and C2+C*06:02) but no HLA-A or -B Unexpectedly, no novel HLA class I was identified in these previously unstudied and anthropologically distinctive populations. In contrast, of 227 KIR detected, 22 are present in all seven populations and 28 are novel. A high diversity of HLA A-C-B haplotypes was observed. In six populations, most haplotypes are represented just once. But in the Hadza, a majority of haplotypes occur more than once, with 2 having high frequencies and 10 having intermediate frequencies. The centromeric (cen) part of the KIR locus exhibits an even balance between cenA and cenB in all seven populations. The telomeric (tel) part has an even balance of telA to telB in East Africa, but this changes across the continent to where telB is vestigial in West Africa. All four KIR ligands (A3/11, Bw4, C1, and C2) are present in six of the populations. HLA haplotypes of the Iraqw and Hadza encode two KIR ligands, whereas the other populations have an even balance between haplotypes encoding one and two KIR ligands. Individuals in these African populations have a mean of 6.8-8.4 different interactions between KIR and HLA class I, compared with 2.9-6.5 for non-Africans.

Critical Protein-Protein Interactions Determine the Biological Activity of Elk-1, a Master Regulator of Stimulus-Induced Gene Transcription.

Elk-1 is a transcription factor that binds together with a dimer of the serum response factor (SRF) to the serum-response element (SRE), a genetic element that connects cellular stimulation with gene transcription. Elk-1 plays an important role in the regulation of cellular proliferation and apoptosis, thymocyte development, glucose homeostasis and brain function. The biological function of Elk-1 relies essentially on the interaction with other proteins. Elk-1 binds to SRF and generates a functional ternary complex that is required to activate SRE-mediated gene transcription. Elk-1 is kept in an inactive state under basal conditions via binding of a SUMO-histone deacetylase complex. Phosphorylation by extracellular signal-regulated protein kinase, c-Jun N-terminal protein kinase or p38 upregulates the transcriptional activity of Elk-1, mediated by binding to the mediator of RNA polymerase II transcription (Mediator) and the transcriptional coactivator p300. Strong and extended phosphorylation of Elk-1 attenuates Mediator and p300 recruitment and allows the binding of the mSin3A-histone deacetylase corepressor complex. The subsequent dephosphorylation of Elk-1, catalyzed by the protein phosphatase calcineurin, facilitates the re-SUMOylation of Elk-1, transforming Elk-1 back to a transcriptionally inactive state. Thus, numerous protein-protein interactions control the activation cycle of Elk-1 and are essential for its biological function.

Nomenclature report 2019: major histocompatibility complex genes and alleles of Great and Small Ape and Old and New World monkey species.

The major histocompatibility complex (MHC) is central to the innate and adaptive immune responses of jawed vertebrates. Characteristic of the MHC are high gene density, gene copy number variation, and allelic polymorphism. Because apes and monkeys are the closest living relatives of humans, the MHCs of these non-human primates (NHP) are studied in depth in the context of evolution, biomedicine, and conservation biology. The Immuno Polymorphism Database (IPD)-MHC NHP Database (IPD-MHC NHP), which curates MHC data of great and small apes, as well as Old and New World monkeys, has been upgraded. The curators of the database are responsible for providing official designations for newly discovered alleles. This nomenclature report updates the 2012 report, and summarizes important nomenclature issues and relevant novel features of the IPD-MHC NHP Database.

Nomenclature report for killer-cell immunoglobulin-like receptors (KIR) in macaque species: new genes/alleles, renaming recombinant entities and IPD-NHKIR updates.

The Killer-cell Immunoglobulin-like Receptors (KIR) are encoded by a diverse group of genes, which are characterized by allelic polymorphism, gene duplications, and recombinations, which may generate recombinant entities. The number of reported macaque KIR sequences is steadily increasing, and these data illustrate a gene system that may match or exceed the complexity of the human KIR cluster. This report lists the names of quality controlled and annotated KIR genes/alleles with all the relevant references for two different macaque species: rhesus and cynomolgus macaques. Numerous recombinant KIR genes in these species necessitate a revision of some of the earlier-published nomenclature guidelines. In addition, this report summarizes the latest information on the Immuno Polymorphism Database (IPD)-NHKIR Database, which contains annotated KIR sequences from four non-human primate species.

Correction to: Nomenclature report 2019: major histocompatibility complex genes and alleles of great and small ape and old and new world monkey species.

The original version of this article contained a spelling error in the Acknowledgments regarding the name of the funding organisation supporting GM and JAH. UKRI-BBSCR should have been UKRI-BBSRC, as is now indicated correctly below.

Sequences of 95 human MHC haplotypes reveal extreme coding variation in genes other than highly polymorphic HLA class I and II.

The most polymorphic part of the human genome, the MHC, encodes over 160 proteins of diverse function. Half of them, including the HLA class I and II genes, are directly involved in immune responses. Consequently, the MHC region strongly associates with numerous diseases and clinical therapies. Notoriously, the MHC region has been intractable to high-throughput analysis at complete sequence resolution, and current reference haplotypes are inadequate for large-scale studies. To address these challenges, we developed a method that specifically captures and sequences the 4.8-Mbp MHC region from genomic DNA. For 95 MHC homozygous cell lines we assembled, de novo, a set of high-fidelity contigs and a sequence scaffold, representing a mean 98% of the target region. Included are six alternative MHC reference sequences of the human genome that we completed and refined. Characterization of the sequence and structural diversity of the MHC region shows the approach accurately determines the sequences of the highly polymorphic HLA class I and HLA class II genes and the complex structural diversity of complement factor C4A/C4B It has also uncovered extensive and unexpected diversity in other MHC genes; an example is MUC22, which encodes a lung mucin and exhibits more coding sequence alleles than any HLA class I or II gene studied here. More than 60% of the coding sequence alleles analyzed were previously uncharacterized. We have created a substantial database of robust reference MHC haplotype sequences that will enable future population scale studies of this complicated and clinically important region of the human genome.

Human NK Cells Downregulate Zap70 and Syk in Response to Prolonged Activation or DNA Damage.

The extent of NK cell activity during the innate immune response affects downstream immune functions and, ultimately, the outcome of infectious or malignant disease. However, the mechanisms that terminate human NK cell responses have yet to be defined. When activation receptors expressed on NK cell surfaces bind to ligands on diseased cells, they initiate a signal that is propagated by a number of intracellular kinases, including Zap70 and Syk, eventually leading to NK cell activation. We assayed Zap70 and Syk content in NK cells from healthy human donors and identified a subset of NK cells with unusually low levels of these two kinases. We found that this Zap70lowSyklow subset consisted of NK cells expressing a range of surface markers, including CD56hi and CD56low NK cells. Upon in vitro stimulation with target cells, Zap70lowSyklow NK cells failed to produce IFN-γ and lysed target cells at one third the capacity of Zap70hiSykhi NK cells. We determined two independent in vitro conditions that induce the Zap70lowSyklow phenotype in NK cells: continuous stimulation with activation beads and DNA damage. The expression of inhibitory receptors, including NKG2A and inhibitory killer Ig-like receptors (KIRs), was negatively correlated with the Zap70lowSyklow phenotype. Moreover, expression of multiple KIRs reduced the likelihood of Zap70 downregulation during continuous activation, regardless of whether NK cells had been educated through KIR-HLA interactions in vivo. Our findings show that human NK cells are able to terminate their functional activity without the aid of other immune cells through the downregulation of activation kinases.

Distinguishing functional polymorphism from random variation in the sequences of >10,000 HLA-A, -B and -C alleles.

HLA class I glycoproteins contain the functional sites that bind peptide antigens and engage lymphocyte receptors. Recently, clinical application of sequence-based HLA typing has uncovered an unprecedented number of novel HLA class I alleles. Here we define the nature and extent of the variation in 3,489 HLA-A, 4,356 HLA-B and 3,111 HLA-C alleles. This analysis required development of suites of methods, having general applicability, for comparing and analyzing large numbers of homologous sequences. At least three amino-acid substitutions are present at every position in the polymorphic α1 and α2 domains of HLA-A, -B and -C. A minority of positions have an incidence >1% for the 'second' most frequent nucleotide, comprising 70 positions in HLA-A, 85 in HLA-B and 54 in HLA-C. The majority of these positions have three or four alternative nucleotides. These positions were subject to positive selection and correspond to binding sites for peptides and receptors. Most alleles of HLA class I (>80%) are very rare, often identified in one person or family, and they differ by point mutation from older, more common alleles. These alleles with single nucleotide polymorphisms reflect the germ-line mutation rate. Their frequency predicts the human population harbors 8-9 million HLA class I variants. The common alleles of human populations comprise 42 core alleles, which represent all selected polymorphism, and recombinants that have assorted this polymorphism.

KIR Donor Selection: Feasibility in Identifying better Donors.

We previously reported that acute myelogenous leukemia (AML) transplants using killer cell immunoglobulin-type receptor (KIR) B haplotype better or best (≥2 B activating gene loci ± Cen B/B) unrelated donors (URDs) yield less relapse and better survival. In this prospective trial we evaluated 535 AML searches from 14 participating centers with centralized donor KIR genotyping for donor selection. This represented 3% to 48% of all AML searches (median 20%) per center, totaling 3 to 172 patients (median 22) per center. Donor KIR genotype was reported at a median of 14 days after request (≤26 days for 76% of searches). In 535 searches, 2080 donors were requested for KIR genotyping (mean 4.3 per search); and a median of 1.8 (range, 0 to 4.5) per search were KIR typed. Choosing more donors for confirmatory HLA and KIR haplotype identification enriched the likelihood of finding KIR better or best donors. The search process identified a mean of 30% KIR better or best donors; the success ranged from 24% to 38% in the 11 centers enrolling ≥8 patients. More donors requested for KIR genotyping increased the likelihood of identifying KIR better or best haplotype donors. Of the 247 transplants, 9.3% used KIR best, 19% used KIR better, and 48% used KIR neutral donors while 24% used a non-KIR-tested donor. KIR genotyping did not delay transplantation. The time from search to transplant was identical for transplants using a KIR-genotyped versus a non-KIR-genotyped donor. Prospective evaluation can rapidly identify KIR favorable genotype donors, but choosing more donors per search would substantially increase the likelihood of having a KIR best or better donor available for transplantation. Transplant centers and donor registries must both commit extra effort to incorporate new characteristics (beyond HLA, age, and parity) into improved donor selection. Deliberate efforts to present additional genetic factors for donor selection will require novel procedures.

Donor Killer Cell Immunoglobulin-Like Receptor Genotype Does Not Improve Graft-versus-Leukemia Responses in Chronic Lymphocytic Leukemia after Unrelated Donor Transplant: A Center for International Blood and Marrow Transplant Research Analysis.

Allogeneic hematopoietic cell transplantation (alloHCT) remains the sole curative therapy for patients with chronic lymphocytic leukemia (CLL), leading to 40% to 45% long-term survival. The impact of donor killer immunoglobulin-like receptor (KIR) genotype on outcomes of unrelated donor (URD) alloHCT for CLL is unknown. We examined 573 adult URD CLL recipient pairs. KIR genotype (presence/absence) was determined for each donor, and comprehensive modeling of interactions with recipient HLA class I loci (KIR ligands) was used to evaluate their effect on relapse and survival. Recipients had a median age of 56 years, and most were not in remission (65%). Both 8/8 HLA-matched (81%) or 7/8 HLA matched grafts (19%) were studied. Factors associated with improved overall survival (OS) were reduced-intensity conditioning (hazard ratio [HR] of death, .76) and good performance status (HR, .46), whereas alloHCT in nonremission (HR, 1.96) and mismatched donors (HR, 2.01) increased mortality. No models demonstrated a relationship between donor KIR genotype and transplant outcomes. Cox regression models comparing donors with A/A versus B/x KIR haplotypes and those with KIR gene content scores of 0 versus 1 versus ≥2 yielded similar rates of nonrelapse mortality, relapse, acute graft-versus-host disease (GVHD), and chronic GVHD and the same progression-free survival and OS. Relapse risk was not different for grafts from donors with KIR3DL1 transplanted into HLA C1/1 versus C2 recipients. This large analysis failed to demonstrate an association between URD KIR genotype and transplant outcome for patients with CLL, and thus KIR genotyping should not be used as a donor selection criterion in this setting.

Defining KIR and HLA Class I Genotypes at Highest Resolution via High-Throughput Sequencing.

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.

Resurrecting KIR2DP1: A Key Intermediate in the Evolution of Human Inhibitory NK Cell Receptors That Recognize HLA-C.

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 KIR2DP1F , the functional antecedent of KIR2DP1 We demonstrate how K44-KIR2DP1F with lysine 44 recognized C1+HLA-C, whereas T44-KIR2DP1F recognized C2+HLA-C. Dimorphisms at 12 other KIR2DP1F 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-KIR2DP1F were stronger receptors than the attenuated C2-specific T44-KIR2DP1F 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. KIR2DP1F has properties consistent with KIR2DP1F having been the founder gene. The first KIR2DP1F alleles encoded K44-C1 receptors; subsequently KIR2DP1F 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 KIR2DP1F Alternatively, pathogen subversion caused its demise. Preservation of KIR2DP1F functional polymorphism was a side effect of fixation of the deletion in KIR2DP1F by micro gene conversion.

Two Orangutan Species Have Evolved Different KIR Alleles and Haplotypes.

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.

Bonobos Maintain Immune System Diversity with Three Functional Types of MHC-B.

Fast-evolving MHC class I polymorphism serves to diversify NK cell and CD8 T cell responses in individuals, families, and populations. Because only chimpanzee and bonobo have strict orthologs of all HLA class I, their study gives unique perspectives on the human condition. We defined polymorphism of Papa-B, the bonobo ortholog of HLA-B, for six wild bonobo populations. Sequences for Papa-B exon 2 and 3 were determined from the genomic DNA in 255 fecal samples, minimally representing 110 individuals. Twenty-two Papa-B alleles were defined, each encoding a different Papa-B protein. No Papa-B is identical to any chimpanzee Patr-B, human HLA-B, or gorilla Gogo-B. Phylogenetic analysis identified a clade of MHC-B, defined by residues 45-74 of the α1 domain, which is broadly conserved among bonobo, chimpanzee, and gorilla. Bonobo populations have 3-14 Papa-B allotypes. Three Papa-B are in all populations, and they are each of a different functional type: allotypes having the Bw4 epitope recognized by killer cell Ig-like receptors of NK cells, allotypes having the C1 epitope also recognized by killer cell Ig-like receptors, and allotypes having neither epitope. For population Malebo, these three Papa-B are the only Papa-B allotypes. Although small in number, their sequence divergence is such that the nucleotide diversity (mean proportional distance) of Papa-B in Malebo is greater than in the other populations and is also greater than expected for random combinations of three Papa-B Overall, Papa-B has substantially less diversity than Patr-B in chimpanzee subspecies and HLA-B in indigenous human populations, consistent with bonobo having experienced narrower population bottlenecks.

IPD-MHC 2.0: an improved inter-species database for the study of the major histocompatibility complex.

The IPD-MHC Database project (http://www.ebi.ac.uk/ipd/mhc/) collects and expertly curates sequences of the major histocompatibility complex from non-human species and provides the infrastructure and tools to enable accurate analysis. Since the first release of the database in 2003, IPD-MHC has grown and currently hosts a number of specific sections, with more than 7000 alleles from 70 species, including non-human primates, canines, felines, equids, ovids, suids, bovins, salmonids and murids. These sequences are expertly curated and made publicly available through an open access website. The IPD-MHC Database is a key resource in its field, and this has led to an average of 1500 unique visitors and more than 5000 viewed pages per month. As the database has grown in size and complexity, it has created a number of challenges in maintaining and organizing information, particularly the need to standardize nomenclature and taxonomic classification, while incorporating new allele submissions. Here, we describe the latest database release, the IPD-MHC 2.0 and discuss planned developments. This release incorporates sequence updates and new tools that enhance database queries and improve the submission procedure by utilizing common tools that are able to handle the varied requirements of each MHC-group.