Co-evolution of MHC class I and variable NK cell receptors in placental mammals

Genetic susceptibility to sarcoid in Arabian horses: associations with MHC class II and compound MHC class I/KLRA genotypes

Although the Major Histocompatibility Complex (MHC) has been repeatedly associated with susceptibility to equine sarcoid, a disease associated with bovine papillomavirus infection, the role of the MHC in the mechanisms of the disease is not fully understood. The objectives of our work were to analyze associations between polymorphic markers of the MHC genomic subregions and of the Natural Killer Complex (NKC) genomic region and the presence of sarcoid in Arabian horses. Microsatellite loci located in the MHC class I, II and III subregions and two MHC class II genes (DRA, DQA1), along with a set of NKC (KLRA, CLEC subregions) microsatelllite markers were genotyped. Fifteen microsatellites of the standard parentage kit, located outside the MHC and NKC regions, were tested as controls. Standard chi-square and Fisher tests with Bonferroni corrections were used for association analyses. Significant associations of MHC class II and MHC class I_KLRA polymorphic markers with the presence of clinical sarcoid were observed. These findings are consistent with biological theory and indicate a role of MHC class I, class II and KLRA molecules in adaptive as well as in innate immune responses to equine sarcoid. Although limited to Arabian horses, these data point to an as yet unadressed hypothesis regarding the possible roles of NK cells in the pathogenesis of equine sarcoid.

Exploring the genetic mechanisms driving KIR diversification

Killer cell immunoglobulin-like receptors (KIRs) are key modulators of natural killer cell activity, displaying either activating or inhibitory potential upon recognition of major histocompatibility complex (MHC) class I molecules. The genomic organization of KIR genes is complex, involving copy number variation and allelic polymorphism, which is probably due to their coevolution with highly polymorphic MHC ligands. The KIR diversity is reflected by more than 70 similar region configurations encountered in humans, generated through meiotic recombination events. Rhesus macaques happen to display even more diversity, and over 100 distinct configurations were identified in a relatively small cohort of animals. More than half of these region configurations feature hybrid KIR genes, suggesting a more pronounced mode of diversification in macaques. The molecular mechanism facilitating meiotic rearrangements in the KIR region is poorly understood. Examination of 21 rhesus macaque and 14 human KIR region configurations revealed the presence of long terminal repeats and PRDM9 binding motifs associated with recombination hotspots. The variable DNA recognition patterns of PRDM9 could potentially contribute to the differing recombination activities documented for the KIR region in humans and macaques. The diversification process of the KIR repertoire in natural killer cells is fundamentally distinct from the mechanisms generating T and B cell receptor diversity or MHC polymorphisms. This sophisticated recombination machinery preserves the functional integrity by the frequent generation of in-frame KIR genes. A diverse KIR repertoire contributes to the protection of individuals and populations against pathogen evasion and subversion.

WITHDRAWN: MHC-B Diversity and Signs of Respiratory Illness in Wild, East African Chimpanzees ( Pan troglodytes schweinfurthii )

bioRxiv has withdrawn this preprint following a formal investigation by the University of New Mexico Office of Research Integrity and Compliance.

Rhesus Macaque Killer Cell Ig-like Receptor Domain 0 Glycans Impact Surface Expression and Ligand Specificity

Defining the MHC class I ligands of rhesus macaque killer cell Ig-like receptors (KIRs) is fundamental to NK cell biology in this species as a model for infectious diseases and comparative immunogenetics. Several rhesus macaque KIRs belong to a phylogenetically distinct group with a three-amino acid deletion in domain 0 (D0). This deletion results in polymorphic differences in potential N-linked glycosylation (PNG) sites adjacent to a predicted KIR-MHC class I contact site. Whereas most KIRs have two tandem PNG sites in D0 (N36FTN39FT), the KIRs containing the deletion only have a single site in this region (N36FT). To discern the contribution of glycosylation to KIR expression and ligand recognition, we constructed PNG mutants for six lineage II KIR genes that eliminate or create sites for N-glycan addition at these locations. The impact of these mutations on total and surface expression was determined by immunoblotting and flow cytometry. Ligand engagement was assessed by coincubating reporter cell lines bearing chimeric KIR-CD3ζ receptors with target cells expressing individual MHC class I molecules and were corroborated by staining with KIR IgG-Fc fusion proteins. We found that N36FT is glycosylated in KIR with a single site, and at least one site is glycosylated in KIRs with two tandem sites. In general, for rhesus KIRs with a single D0 glycosylation site, that site contributes to surface expression. For KIRs with two tandem sites, the first site can contribute to ligand specificity. This study establishes that D0 glycosylation of rhesus macaque KIRs modulates surface expression and contributes to ligand specificity.

The Primate Major Histocompatibility Complex: An Illustrative Example of Gene Family Evolution

Gene families are groups of evolutionarily-related genes. One large gene family that has experienced rapid evolution is the Major Histocompatibility Complex (MHC), whose proteins serve critical roles in innate and adaptive immunity. Across the ~60 million year history of the primates, some MHC genes have turned over completely, some have changed function, some have converged in function, and others have remained essentially unchanged. Past work has typically focused on identifying MHC alleles within particular species or comparing gene content, but more work is needed to understand the overall evolution of the gene family across species. Thus, despite the immunologic importance of the MHC and its peculiar evolutionary history, we lack a complete picture of MHC evolution in the primates. We readdress this question using sequences from dozens of MHC genes and pseudogenes spanning the entire primate order, building a comprehensive set of gene and allele trees with modern methods. Overall, we find that the Class I gene subfamily is evolving much more quickly than the Class II gene subfamily, with the exception of the Class II MHC-DRB genes. We also pay special attention to the often-ignored pseudogenes, which we use to reconstruct different events in the evolution of the Class I region. We find that despite the shared function of the MHC across species, different species employ different genes, haplotypes, and patterns of variation to achieve a successful immune response. Our trees and extensive literature review represent the most comprehensive look into MHC evolution to date.

Diversity of major histocompatibility complex (MHC) and natural killer cell receptor (NKR) genes and their interactions in domestic horses

The immunogenome is the part of the genome that underlies immune mechanisms and evolves under various selective pressures. Two complex regions of the immunogenome, major histocompatibility complex (MHC) and natural killer cell receptor (NKR) genes, play an important role in the response to selective pressures of pathogens. Their importance is expressed by their genetic polymorphism at the molecular level, and their diversity associated with different types of diseases at the population level. Findings of associations between specific combinations of MHC/NKR haplotypes with different diseases in model species suggest that these gene complexes did not evolve independently. No such associations have been described in horses so far. The aim of the study was to detect associations between MHC and NKR gene/microsatellite haplotypes in three horse breed groups (Camargue, African, and Romanian) by statistical methods; chi-square test, Fisher's exact test, Pearson's goodness-of-fit test and logistic regression. Associations were detected for both MHC/NKR genes and microsatellites; the most significant associations were found between the most variable KLRA3 gene and the EQCA-1 or EQCA-2 genes. This finding supports the assumption that the KLRA3 is an important receptor for MHC I and that interactions of these molecules play important roles in the horse immunity and reproduction. Despite some limitations of the study such as low numbers of horses or lack of knowledge of the selected genes functions, the results were consistent across different statistical methods and remained significant even after overconservative Bonferroni corrections. We therefore consider them biologically plausible.

The KIR repertoire of a West African chimpanzee population is characterized by limited gene, allele, and haplotype variation

Introduction

The killer cell immunoglobulin-like receptors (KIR) play a pivotal role in modulating the NK cell responses, for instance, through interaction with major histocompatibility complex (MHC) class I molecules. Both gene systems map to different chromosomes but co-evolved during evolution. The human KIR gene family is characterized by abundant allelic polymorphism and copy number variation. In contrast, our knowledge of the KIR repertoire in chimpanzees is limited to 39 reported alleles, with no available population data. Only three genomic KIR region configurations have been mapped, and seventeen additional ones were deduced by genotyping.

Methods

Previously, we documented that the chimpanzee MHC class I repertoire has been skewed due to an ancient selective sweep. To understand the depth of the sweep, we set out to determine the full-length KIR transcriptome - in our MHC characterized pedigreed West African chimpanzee cohort - using SMRT sequencing (PacBio). In addition, the genomic organization of 14 KIR haplotypes was characterized by applying a Cas9-mediated enrichment approach in concert with long-read sequencing by Oxford Nanopore Technologies.

Results

In the cohort, we discovered 35 undescribed and 15 already recorded Patr-KIR alleles, and a novel hybrid KIR gene. Some KIR transcripts are subject to evolutionary conserved alternative splicing events. A detailed insight on the KIR region dynamics (location and order of genes) was obtained, however, only five new KIR region configurations were detected. The population data allowed to investigate the distribution of the MHC-C1 and C2-epitope specificity of the inhibitory lineage III KIR repertoire, and appears to be skewed towards C2.

Discussion

Although the KIR region is known to evolve fast, as observed in other primate species, our overall conclusion is that the genomic architecture and repertoire in West African chimpanzees exhibit only limited to moderate levels of variation. Hence, the ancient selective sweep that affected the chimpanzee MHC class I region may also have impacted the KIR system.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

Innate receptors modulating adaptive T cell responses: KIR-HLA interactions and T cell-mediated control of chronic viral infections

Killer-cell immunoglobulin-like receptors (KIRs) are mainly expressed on natural killer (NK) cells and are key regulators of innate immune responses. NK cells are the first responders in the face of infection and help promote placentation during pregnancy; the importance of KIRs in these NK-mediated processes is well-established. However, mounting evidence suggests that KIRs also have a prominent and long-lasting effect on the adaptive immune system. Here, we review the evidence for the impact of KIRs on T cell responses with a focus on the clinical significance of this interaction.

The impact of HLA polymorphism on herpesvirus infection and disease

Human Leukocyte Antigens (HLA) are cell surface molecules, central in coordinating innate and adaptive immune responses, that are targets of strong diversifying natural selection by pathogens. Of these pathogens, human herpesviruses have a uniquely ancient relationship with our species, where coevolution likely has reciprocating impact on HLA and viral genomic diversity. Consistent with this notion, genetic variation at multiple HLA loci is strongly associated with modulating immunity to herpesvirus infection. Here, we synthesize published genetic associations of HLA with herpesvirus infection and disease, both from case/control and genome-wide association studies. We analyze genetic associations across the eight human herpesviruses and identify HLA alleles that are associated with diverse herpesvirus-related phenotypes. We find that whereas most HLA genetic associations are virus- or disease-specific, HLA-A*01 and HLA-A*02 allotypes may be more generally associated with immune susceptibility and control, respectively, across multiple herpesviruses. Connecting genetic association data with functional corroboration, we discuss mechanisms by which diverse HLA and cognate receptor allotypes direct variable immune responses during herpesvirus infection and pathogenesis. Together, this review examines the complexity of HLA-herpesvirus interactions driven by differential T cell and Natural Killer cell immune responses.

Comparative genomics of the Leukocyte Receptor Complex in carnivores

Background

The mammalian Leukocyte Receptor Complex (LRC) chromosomal region may contain gene families for the killer cell immunoglobulin-like receptor (KIR) and/or leukocyte immunoglobulin-like receptor (LILR) collections as well as various framing genes. This complex region is well described in humans, mice, and some domestic animals. Although single KIR genes are known in some Carnivora, their complements of LILR genes remain largely unknown due to obstacles in the assembly of regions of high homology in short-read based genomes.

Methods

As part of the analysis of felid immunogenomes, this study focuses on the search for LRC genes in reference genomes and the annotation of LILR genes in Felidae. Chromosome-level genomes based on single-molecule long-read sequencing were preferentially sought and compared to representatives of the Carnivora.

Results

Seven putatively functional LILR genes were found across the Felidae and in the Californian sea lion, four to five genes in Canidae, and four to nine genes in Mustelidae. They form two lineages, as seen in the Bovidae. The ratio of functional genes for activating LILRs to inhibitory LILRs is slightly in favor of inhibitory genes in the Felidae and the Canidae; the reverse is seen in the Californian sea lion. This ratio is even in all of the Mustelidae except the Eurasian otter, which has a predominance of activating LILRs. Various numbers of LILR pseudogenes were identified.

Conclusions

The structure of the LRC is rather conservative in felids and the other Carnivora studied. The LILR sub-region is conserved within the Felidae and has slight differences in the Canidae, but it has taken various evolutionary paths in the Mustelidae. Overall, the process of pseudogenization of LILR genes seems to be more frequent for activating receptors. Phylogenetic analysis found no direct orthologues across the Carnivora which corroborate the rapid evolution of LILRs seen in mammals.

Malaria-driven adaptation of MHC class I in wild bonobo populations

The malaria parasite Plasmodium falciparum causes substantial human mortality, primarily in equatorial Africa. Enriched in affected African populations, the B*53 variant of HLA-B, a cell surface protein that presents peptide antigens to cytotoxic lymphocytes, confers protection against severe malaria. Gorilla, chimpanzee, and bonobo are humans' closest living relatives. These African apes have HLA-B orthologs and are infected by parasites in the same subgenus (Laverania) as P. falciparum, but the consequences of these infections are unclear. Laverania parasites infect bonobos (Pan paniscus) at only one (TL2) of many sites sampled across their range. TL2 spans the Lomami River and has genetically divergent subpopulations of bonobos on each side. Papa-B, the bonobo ortholog of HLA-B, includes variants having a B*53-like (B07) peptide-binding supertype profile. Here we show that B07 Papa-B occur at high frequency in TL2 bonobos and that malaria appears to have independently selected for different B07 alleles in the two subpopulations.

Carcinogenesis resistance in the longest-lived rodent, the naked mole-rat

Certain mammalian species are resistant to cancer, and a better understanding of how this cancer resistance arises could provide valuable insights for basic cancer research. Recent technological innovations in molecular biology have allowed the study of cancer-resistant mammals, despite the fact that they are not the classical model animals, which are easily studied using genetic approaches. Naked mole-rats (NMRs; Heterocephalus glaber) are the longest-lived rodent, with a maximum lifespan of more than 37 years, and almost never show spontaneous carcinogenesis. NMRs are currently attracting much attention from aging and cancer researchers, and published studies on NMR have continued to increase over the past decade. Cancer development occurs via multiple steps and involves many biological processes. Recent research on the NMR as a model for cancer resistance suggests that they possess various unique carcinogenesis-resistance mechanisms, including efficient DNA repair pathways, cell-autonomous resistance to transformation, and dampened inflammatory response. Here, we summarize the molecular mechanisms of carcinogenesis resistance in NMR, which have been uncovered over the past two decades, and discuss future perspectives.

KIR-HLA gene diversities and susceptibility to lung cancer

Killer-cell immunoglobulin-like receptors (KIR) are essential for acquiring natural killer (NK) cell effector function, which is modulated by a balance between the net input of signals derived from inhibitory and activating receptors through engagement by human leukocyte antigen (HLA) class I ligands. KIR and HLA loci are polygenic and polymorphic and exhibit substantial variation between individuals and populations. We attempted to investigate the contribution of KIR complex and HLA class I ligands to the genetic predisposition to lung cancer in the native population of southern Iran. We genotyped 16 KIR genes for a total of 232 patients with lung cancer and 448 healthy controls (HC), among which 85 patients and 178 HCs were taken into account for evaluating combined KIR-HLA associations. KIR2DL2 and 2DS2 were increased significantly in patients than in controls, individually (OR 1.63, and OR 1.42, respectively) and in combination with HLA-C1 ligands (OR 1.99, and OR 1.93, respectively). KIR3DS1 (OR 0.67) and 2DS1 (OR 0.69) were more likely presented in controls in the absence of their relative ligands. The incidence of CxTx subset was increased in lung cancer patients (OR 1.83), and disease risk strikingly increased by more than fivefold among genotype ID19 carriers (a CxTx genotype that carries 2DL2 in the absence of 2DS2, OR 5.92). We found that genotypes with iKIRs > aKIRs (OR 1.67) were more frequently presented in lung cancer patients. Additionally, patients with lung cancer were more likely to carry the combination of CxTx/2DS2 compared to controls (OR 2.04), and iKIRs > aKIRs genotypes in the presence of 2DL2 (OR 2.05) increased the likelihood of lung cancer development. Here we report new susceptibility factors and the contribution of KIR and HLA-I encoding genes to lung cancer risk, highlighting an array of genetic effects and disease setting which regulates NK cell responsiveness. Our results suggest that inherited KIR genes and HLA-I ligands specifying the educational state of NK cells can modify lung cancer risk.

The natural cytotoxicity receptor (NCR) genes in the family Felidae

Natural killer (NK) cells belong to the innate immune system. The germline-encoded natural killer cell receptors represent activating and inhibitory receptors regulating multiple NK cell activities. The natural cytotoxicity receptors (NCRs) are activating natural cytotoxicity triggering receptors 1, 2 and 3 (NKp46, NKp44, and NKp30), encoded by the genes NCR1, NCR2, and NCR3, respectively. NCRs may be expressed in different cell types engaged in mechanisms of innate and adaptive immunity. The family Felidae, comprising the domestic cat and a wide variety of free-ranging species represents a well-suited model for biomedical and evolutionary studies. We characterized the NCR1, NCR2 and NCR3 genes in a panel of felid species. We confirmed the presence of potentially functional genes NCR1, NCR2 and NCR3 in all species. All three genes are conserved within the family and are similar to other phylogenetically related mammalian families. The NCR1 and NCR2 phylogenetic trees based on both nucleotide and protein sequences corresponded to the current zoological taxonomy, with some exceptions suggesting effects of different selection pressures in some species. Highly conserved NCR3 sequences did not allow a robust phylogenetic analysis. Most interspecific differences both at the nucleotide and protein level were found in NCR2. Within species, the most polymorphic CDS was detected in NCR1. Selection analyses indicated the effects of purifying selection on individual amino acid sites in all three genes. In stray cats, a rather high intraspecific diversity was observed. This article is protected by copyright. All rights reserved.

Immunogenomics of Killer Cell Immunoglobulin-Like Receptor (KIR) and HLA Class I: Coevolution and Consequences for Human Health

Interactions of killer cell immunoglobin-like receptors (KIR) with human leukocyte antigens (HLA) class I regulate effector functions of key cytotoxic cells of innate and adaptive immunity. The extreme diversity of this interaction is genetically determined, having evolved in the ever-changing environment of pathogen exposure. Diversity of KIR and HLA genes is further facilitated by their independent segregation on separate chromosomes. That fetal implantation relies on many of the same types of immune cells as infection control places certain constraints on the evolution of KIR interactions with HLA. Consequently, specific inherited combinations of receptors and ligands may predispose to specific immune-mediated diseases, including autoimmunity. Combinatorial diversity of KIR and HLA class I can also differentiate success rates of immunotherapy directed to these diseases. Progress toward both etiopathology and predicting response to therapy is being achieved through detailed characterization of the extent and consequences of the combinatorial diversity of KIR and HLA. Achieving these goals is more tractable with the development of integrated analyses of molecular evolution, function, and pathology that will establish guidelines for understanding and managing risks. Here, we present what is known about the coevolution of KIR with HLA class I and the impact of their complexity on immune function and homeostasis.

Human herpesvirus diversity is altered in HLA class I binding peptides

Herpesviruses are ubiquitous, genetically diverse DNA viruses, with long-term presence in humans associated with infrequent but significant pathology. Human leukocyte antigen (HLA) class I presents intracellularly derived peptide fragments from infected tissue cells to CD8+ T and natural killer cells, thereby directing antiviral immunity. Allotypes of highly polymorphic HLA class I are distinguished by their peptide binding repertoires. Because this HLA class I variation is a major determinant of herpesvirus disease, we examined if sequence diversity of virus proteins reflects evasion of HLA presentation. Using population genomic data from Epstein–Barr virus (EBV), human cytomegalovirus (HCMV), and Varicella–Zoster virus, we tested whether diversity differed between the regions of herpesvirus proteins that can be recognized, or not, by HLA class I. Herpesviruses exhibit lytic and latent infection stages, with the latter better enabling immune evasion. Whereas HLA binding peptides of lytic proteins are conserved, we found that EBV and HCMV proteins expressed during latency have increased peptide sequence diversity. Similarly, latent, but not lytic, herpesvirus proteins have greater population structure in HLA binding than nonbinding peptides. Finally, we found patterns consistent with EBV adaption to the local HLA environment, with less efficient recognition of EBV isolates by high-frequency HLA class I allotypes. Here, the frequency of CD8+ T cell epitopes inversely correlated with the frequency of HLA class I recognition. Previous analyses have shown that pathogen-mediated natural selection maintains exceptional polymorphism in HLA residues that determine peptide recognition. Here, we show that HLA class I peptide recognition impacts diversity of globally widespread pathogens.

Comparative Genomics of the Major Histocompatibility Complex (MHC) of Felids

This review summarizes the current knowledge on the major histocompatibility complex (MHC) of the family Felidae. This family comprises an important domestic species, the cat, as well as a variety of free-living felids, including several endangered species. As such, the Felidae have the potential to be an informative model for studying different aspects of the biological functions of MHC genes, such as their role in disease mechanisms and adaptation to different environments, as well as the importance of genetic diversity for conservation issues in free-ranging or captive populations. Despite this potential, the current knowledge on the MHC in the family as a whole is fragmentary and based mostly on studies of the domestic cat and selected species of big cats. The overall structure of the domestic cat MHC is similar to other mammalian MHCs following the general scheme "centromere-MHC class I-MHC class III-MHC class II" with some differences in the gene contents. An unambiguously defined orthologue of the non-classical class I HLA-E gene has not been identified so far and the class II DQ and DP genes are missing or pseudogenized, respectively. A comparison with available genomes of other felids showed a generally high level of structural and sequence conservation of the MHC region. Very little and fragmentary information on in vitro and/or in vivo biological functions of felid MHC genes is available. So far, no association studies have indicated effects of MHC genetic diversity on a particular disease. No information is available on the role of MHC class I molecules in interactions with Natural Killer (NK) cell receptors or on the putative evolutionary interactions (co-evolution) of the underlying genes. A comparison of complex genomic regions encoding NK cell receptors (the Leukocyte Receptor Complex, LRC and the Natural Killer Cell Complex, NKC) in the available felid genomes showed a higher variability in the NKC compared to the LRC and the MHC regions. Studies of the genetic diversity of domestic cat populations and/or specific breeds have focused mainly on DRB genes. Not surprisingly, higher levels of MHC diversity were observed in stray cats compared to pure breeds, as evaluated by DRB sequencing as well as by MHC-linked microsatellite typing. Immunogenetic analysis in wild felids has only been performed on MHC class I and II loci in tigers, Namibian leopards and cheetahs. This information is important as part of current conservation tasks to assess the adaptive potential of endangered wild species at the human-wildlife interface, which will be essential for preserving biodiversity in a functional ecosystem.

Highly-contiguous bovine genomes underpin accurate functional analyses and updated nomenclature of MHC class I

The major histocompatibility complex (MHC) class I region of cattle is both highly polymorphic and, unlike many species, highly variable in gene content between haplotypes. Cattle MHC class I alleles were historically grouped by sequence similarity in the more conserved 3' end of the coding sequence to form phylogenetic allele groups. This has formed the basis of current cattle MHC class I nomenclature. We presently describe and compare five fully assembled MHC class I haplotypes using the latest cattle and yak genome assemblies. Of the five previously described "pseudogenes" in the cattle MHC class I region, Pseudogene 3 is putatively functional in all haplotypes and Pseudogene 6 and Pseudogene 7 are putatively functional in some haplotypes. This was reinforced by evidence of transcription. Based on full gene sequences as well as 3' coding sequence, we identified distinct subgroups of BoLA-3 and BoLA-6 that represent distinct genetic loci. We further examined allele-specific expression using transcriptomic data revealing that certain alleles are consistently weakly expressed compared to others. These observations will help to inform further studies into how MHC class I region variability influences T cell and natural killer cell functions in cattle.

The Genomic Organization of the LILR Region Remained Largely Conserved Throughout Primate Evolution: Implications for Health And Disease

The genes of the leukocyte immunoglobulin-like receptor (LILR) family map to the leukocyte receptor complex (LRC) on chromosome 19, and consist of both activating and inhibiting entities. These receptors are often involved in regulating immune responses, and are considered to play a role in health and disease. The human LILR region and evolutionary equivalents in some rodent and bird species have been thoroughly characterized. In non-human primates, the LILR region is annotated, but a thorough comparison between humans and non-human primates has not yet been documented. Therefore, it was decided to undertake a comprehensive comparison of the human and non-human primate LILR region at the genomic level. During primate evolution the organization of the LILR region remained largely conserved. One major exception, however, is provided by the common marmoset, a New World monkey species, which seems to feature a substantial contraction of the number of LILR genes in both the centromeric and the telomeric region. Furthermore, genomic analysis revealed that the killer-cell immunoglobulin-like receptor gene KIR3DX1, which maps in the LILR region, features one copy in humans and great ape species. A second copy, which might have been introduced by a duplication event, was observed in the lesser apes, and in Old and New World monkey species. The highly conserved gene organization allowed us to standardize the LILR gene nomenclature for non-human primate species, and implies that most of the receptors encoded by these genes likely fulfill highly preserved functions.

Rapid Characterization of Complex Killer Cell Immunoglobulin-Like Receptor (KIR) Regions Using Cas9 Enrichment and Nanopore Sequencing

Long-read sequencing approaches have considerably improved the quality and contiguity of genome assemblies. Such platforms bear the potential to resolve even extremely complex regions, such as multigenic immune families and repetitive stretches of DNA. Deep sequencing coverage, however, is required to overcome low nucleotide accuracy, especially in regions with high homopolymer density, copy number variation, and sequence similarity, such as the MHC and KIR gene clusters of the immune system. Therefore, we have adapted a targeted enrichment protocol in combination with long-read sequencing to efficiently annotate complex KIR gene regions. Using Cas9 endonuclease activity, segments of the KIR gene cluster were enriched and sequenced on an Oxford Nanopore Technologies platform. This provided sufficient coverage to accurately resolve and phase highly complex KIR haplotypes. Our strategy eliminates PCR-induced amplification errors, facilitates rapid characterization of large and complex multigenic regions, including its epigenetic footprint, and is applicable in multiple species, even in the absence of a reference genome.

Analysis of Genomic DNA from Medieval Plague Victims Suggests Long-Term Effect of Yersinia pestis on Human Immunity Genes

Pathogens and associated outbreaks of infectious disease exert selective pressure on human populations, and any changes in allele frequencies that result may be especially evident for genes involved in immunity. In this regard, the 1346-1353 Yersinia pestis-caused Black Death pandemic, with continued plague outbreaks spanning several hundred years, is one of the most devastating recorded in human history. To investigate the potential impact of Y. pestis on human immunity genes, we extracted DNA from 36 plague victims buried in a mass grave in Ellwangen, Germany in the 16th century. We targeted 488 immune-related genes, including HLA, using a novel in-solution hybridization capture approach. In comparison with 50 modern native inhabitants of Ellwangen, we find differences in allele frequencies for variants of the innate immunity proteins Ficolin-2 and NLRP14 at sites involved in determining specificity. We also observed that HLA-DRB1*13 is more than twice as frequent in the modern population, whereas HLA-B alleles encoding an isoleucine at position 80 (I-80+), HLA C*06:02 and HLA-DPB1 alleles encoding histidine at position 9 are half as frequent in the modern population. Simulations show that natural selection has likely driven these allele frequency changes. Thus, our data suggest that allele frequencies of HLA genes involved in innate and adaptive immunity responsible for extracellular and intracellular responses to pathogenic bacteria, such as Y. pestis, could have been affected by the historical epidemics that occurred in Europe.

ELF3 activated by a superenhancer and an autoregulatory feedback loop is required for high-level HLA-C expression on extravillous trophoblasts

HLA-C arose during evolution of pregnancy in the great apes 10 to 15 million years ago. It has a dual function on placental extravillous trophoblasts (EVTs) as it contributes to both tolerance and immunity at the maternal-fetal interface. The mode of its regulation is of considerable interest in connection with the biology of pregnancy and pregnancy abnormalities. First-trimester primary EVTs in which HLA-C is highly expressed, as well as JEG3, an EVT model cell line, were employed. Single-cell RNA-seq data and quantitative PCR identified high expression of the transcription factor ELF3 in those cells. Chromatin immunoprecipitation (ChIP)-PCR confirmed that both ELF3 and MED1 bound to the proximal HLA-C promoter region. However, binding of RFX5 to this region was absent or severely reduced, and the adjacent HLA-B locus remained closed. Expression of HLA-C was inhibited by ELF3 small interfering RNAs (siRNAs) and by wrenchnolol treatment. Wrenchnolol is a cell-permeable synthetic organic molecule that mimics ELF3 and is relatively specific for binding to ELF3's coactivator, MED23, as our data also showed in JEG3. Moreover, the ELF3 gene is regulated by a superenhancer that spans more than 5 Mb, identified by assay for transposase-accessible chromatin using sequencing (ATAC-seq), as well as by its sensitivity to (+)-JQ1 (inhibitor of BRD4). ELF3 bound to its own promoter, thus creating an autoregulatory feedback loop that establishes expression of ELF3 and HLA-C in trophoblasts. Wrenchnolol blocked binding of MED23 to ELF3, thus disrupting the positive-feedback loop that drives ELF3 expression, with down-regulation of HLA-C expression as a consequence.

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.

The Unusual Immune System of the Naked Mole-Rat

The immune system plays a critical role in host defense to pathogens, tissue homeostasis, cancer development, and several aging-associated chronic inflammatory diseases. The naked mole-rat (Heterocephalus glaber) is a subterranean rodent with both extraordinary longevity and cancer-resistant phenotypes. Unlike the immune system of standard laboratory rodents, that of the naked mole-rat features a higher myeloid-to-lymphoid ratio, lacks natural killer cells, has higher pro-inflammatory cytokine production in macrophages, and exhibits a novel LPS-responsive neutrophil subset that highly expresses several antimicrobials. Given these unusual features, the potential involvement of the naked mole-rat's immune system in their longevity and cancer-resistance remains enigmatic. In this chapter, we summarize the current knowledge of the immune system in the naked mole-rat, including the immune cell repertoire, the primary and secondary lymphoid organs, and the inflammatory responses to the pathogenic stimulation such as bacterial toxins. We compare these findings to published studies of the other subterranean rodents and discuss how the environmental factors in which they have evolved may have influenced their immune function.

NK cells in the tumor microenvironment: Prognostic and theranostic impact. Recent advances and trends

NK cells orchestrate the tumor destruction and control metastasis in a coordinated way with other immune cells of the tumor microenvironment. However, NK cell infiltration in the tumor microenvironment is limited, and tumor cells have developed numerous mechanisms to escape NK cell attack. As a result, NK cells that have been able to infiltrate the tumors are exhausted, and metabolically and functionally impaired. Depending this impairment the prognostic and theranostic values of NK cells differ depending on the studies, the type of cancer, the stage of tumor and the nature of the tumor microenvironment. Extensive studies have been done to investigate different strategies to improve the NK cell function, and nowadays, a battery of therapeutic tools are being tested, with promising results.

Killer Cell Immunoglobulin-like Receptor Variants Are Associated with Protection from Symptoms Associated with More Severe Course in Parkinson Disease

Immune dysfunction plays a role in the development of Parkinson disease (PD). NK cells regulate immune functions and are modulated by killer cell immunoglobulin-like receptors (KIR). KIR are expressed on the surface of NK cells and interact with HLA class I ligands on the surface of all nucleated cells. We investigated KIR-allelic polymorphism to interrogate the role of NK cells in PD. We sequenced KIR genes from 1314 PD patients and 1978 controls using next-generation methods and identified KIR genotypes using custom bioinformatics. We examined associations of KIR with PD susceptibility and disease features, including age at disease onset and clinical symptoms. We identified two KIR3DL1 alleles encoding highly expressed inhibitory receptors associated with protection from PD clinical features in the presence of their cognate ligand: KIR3DL1*015/HLA-Bw4 from rigidity (p _c = 0.02, odds ratio [OR] = 0.39, 95% confidence interval [CI] 0.23-0.69) and KIR3DL1*002/HLA-Bw4i from gait difficulties (p _c = 0.05, OR = 0.62, 95% CI 0.44-0.88), as well as composite symptoms associated with more severe disease. We also developed a KIR3DL1/HLA interaction strength metric and found that weak KIR3DL1/HLA interactions were associated with rigidity (p_c = 0.05, OR = 9.73, 95% CI 2.13-172.5). Highly expressed KIR3DL1 variants protect against more debilitating symptoms of PD, strongly implying a role of NK cells in PD progression and manifestation.

Variations in killer-cell immunoglobulin-like receptor and human leukocyte antigen genes and immunity to malaria

Malaria is one of the deadliest infectious diseases in the world. Immune responses to Plasmodium falciparum malaria vary among individuals and between populations. Human genetic variation in immune system genes is likely to play a role in this heterogeneity. Natural killer (NK) cells produce inflammatory cytokines in response to malaria infection, kill intraerythrocytic Plasmodium falciparum parasites by cytolysis, and participate in the initiation and development of adaptive immune responses to plasmodial infection. These functions are modulated by interactions between killer-cell immunoglobulin-like receptors (KIRs) and human leukocyte antigens (HLAs). Therefore, variations in KIR and HLA genes can have a direct impact on NK cell functions. Understanding the role of KIRs and HLAs in immunity to malaria can help to better characterize antimalarial immune responses. In this review, we summarize the different KIRs and HLAs associated with immunity to malaria thus far.

Allele frequency net database (AFND) 2020 update: gold-standard data classification, open access genotype data and new query tools

The Allele Frequency Net Database (AFND, www.allelefrequencies.net) provides the scientific community with a freely available repository for the storage of frequency data (alleles, genes, haplotypes and genotypes) related to human leukocyte antigens (HLA), killer-cell immunoglobulin-like receptors (KIR), major histocompatibility complex Class I chain related genes (MIC) and a number of cytokine gene polymorphisms in worldwide populations. In the last five years, AFND has become more popular in terms of clinical and scientific usage, with a recent increase in genotyping data as a necessary component of Short Population Report article submissions to another scientific journal. In addition, we have developed a user-friendly desktop application for HLA and KIR genotype/population data submissions. We have also focused on classification of existing and new data into ‘gold–silver–bronze’ criteria, allowing users to filter and query depending on their needs. Moreover, we have also continued to expand other features, for example focussed on HLA associations with adverse drug reactions. At present, AFND contains >1600 populations from >10 million healthy individuals, making AFND a valuable resource for the analysis of some of the most polymorphic regions in the human genome.

Comparative genetics of the major histocompatibility complex in humans and nonhuman primates

The major histocompatibility complex (MHC) is one of the most gene-dense regions of the mammalian genome. Multiple genes within the human MHC (HLA) show extensive polymorphism, and currently, more than 26,000 alleles divided over 39 different genes are known. Nonhuman primate (NHP) species are grouped into great and lesser apes and Old and New World monkeys, and their MHC is studied mostly because of their important role as animal models in preclinical research or in connection with conservation biology purposes. The evolutionary equivalents of many of the HLA genes are present in NHP species, and these genes may also show abundant levels of polymorphism. This review is intended to provide a comprehensive comparison relating to the organization and polymorphism of human and NHP MHC regions.

KIR Variation in Iranians Combines High Haplotype and Allotype Diversity With an Abundance of Functional Inhibitory Receptors

Natural killer (NK) cells are innate lymphocytes that eliminate infected and transformed cells. They discriminate healthy from diseased tissue through killer cell Ig-like receptor (KIR) recognition of HLA class I ligands. Directly impacting NK cell function, KIR polymorphism associates with infection control and multiple autoimmune and pregnancy syndromes. Here we analyze KIR diversity of 241 individuals from five groups of Iranians. These five populations represent Baloch, Kurd, and Lur, together comprising 15% of the ethnically diverse Iranian population. We identified 159 KIR alleles, including 11 not previously characterized. We also identified 170 centromeric and 94 telomeric haplotypes, and 15 different KIR haplotypes carrying either a deletion or duplication encompassing one or more complete KIR genes. As expected, comparing our data with those representing major worldwide populations revealed the greatest similarity between Iranians and Europeans. Despite this similarity we observed higher frequencies of KIR3DL1^*001 in Iran than any other population, and the highest frequency of HLA-B^*51, a Bw4-containing allotype that acts as a strong educator of KIR3DL1^*001⁺ NK cells. Compared to Europeans, the Iranians we studied also have a reduced frequency of 3DL1^*004, which encodes an allotype that is not expressed at the NK cell surface. Concurrent with the resulting high frequency of strong viable interactions between inhibitory KIR and polymorphic HLA class I, the majority of KIR-A haplotypes characterized do not express a functional activating receptor. By contrast, the most frequent KIR-B haplotype in Iran expresses only one functional inhibitory KIR and the maximum number of activating KIR. This first complete, high-resolution, characterization of the KIR locus of Iranians will form a valuable reference for future clinical and population studies.

Regulation of NK-Cell Function by HLA Class II

Natural Killer (NK) cells were initially described as part of the innate immune system and characterized by their ability to lyse malignant and virus-infected cells. The cytolytic function of NK cells is tightly controlled by activating and inhibitory receptors expressed on the cell surface. Ligands that interact with a variety of NK-cell receptors include the human leukocyte antigen (HLA) molecules, and the regulation of NK-cell function by HLA class I molecules is well-established. Earlier studies also suggested a role of HLA class II molecules in regulating NK cell activity; yet, interactions between HLA class II molecules and NK cell receptors have not been well-characterized. We recently identified a subset of HLA-DP molecules that can serve as ligands for the natural cytotoxicity receptor NKp44 and activate NK cells. This novel receptor-ligand interaction provides a potential mechanism to explain the strong associations of HLA-DP molecules with HBV infection outcomes, graft-vs.-host disease and inflammatory bowel disease. Furthermore, it adds a new mechanism for NK-cell crosstalk with immune cells expressing HLA class II molecules. In this perspective article, we discuss the potential implications of NK cell receptor interactions with HLA class II molecules for the regulation of immune responses.

In vitro education of human natural killer cells by KIR3DL1

Using NK cells isolated from individuals who lack the Bw4 epitope on HLA-B, Pugh et al reveal that KIR3DL1⁺ NK cells can be educated in vitro by co-culturing them with target cells that display the missing epitope.

During development, NK cells are “educated” to respond aggressively to cells with low surface expression of HLA class I, a hallmark of malignant and infected cells. The mechanism of education involves interactions between inhibitory killer immunoglobulin–like receptors (KIRs) and specific HLA epitopes, but the details of this process are unknown. Because of the genetic diversity of HLA class I genes, most people have NK cells that are incompletely educated, representing an untapped source of human immunity. We demonstrate how mature peripheral KIR3DL1⁺ human NK cells can be educated in vitro. To accomplish this, we trained NK cells expressing the inhibitory KIR3DL1 receptor by co-culturing them with target cells that expressed its ligand, Bw4⁺HLA-B. After this training, KIR3DL1⁺ NK cells increased their inflammatory and lytic responses toward target cells lacking Bw4⁺HLA-B, as though they had been educated in vivo. By varying the conditions of this basic protocol, we provide mechanistic and translational insights into the process NK cell education.

Single-cell transcriptomics of the naked mole-rat reveals unexpected features of mammalian immunity

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.

Natural Killer Cells Offer Differential Protection From Leukemia in Chinese Southern Han

Interactions of human natural killer (NK) cell inhibitory receptors with polymorphic HLA-A, -B and -C molecules educate NK cells for immune surveillance against tumor cells. The KIR A haplotype encodes a distinctive set of HLA-specific NK cell inhibiting receptors having strong influence on immunity. We observed higher frequency of KIR A homozygosity among 745 healthy Chinese Southern Han than 836 adult patients representing three types of leukemia: ALL (OR = 0.68, 95% CI = 0.52-0.89, p = 0.004), AML (OR = 0.76, 95% CI = 0.59-0.98, p = 0.034), and CML (OR = 0.72 95% CI = 0.51-1.0, ns). We observed the same trend for NHL (OR = 0.47 95% CI = 0.26-0.88 p = 0.017). For ALL, the protective effect of the KIR AA genotype was greater in the presence of KIR ligands C1 (Pc = 0.01) and Bw4 (Pc = 0.001), which are tightly linked in East Asians. By contrast, the C2 ligand strengthened protection from CML (Pc = 0.004). NK cells isolated from KIR AA individuals were significantly more cytotoxic toward leukemic cells than those from other KIR genotypes (p < 0.0001). These data suggest KIR allotypes encoded by East Asian KIR A haplotypes are strongly inhibitory, arming NK cells to respond to leukemogenic cells having altered HLA expression. Thus, the study of populations with distinct KIR and HLA distributions enlightens understanding of immune mechanisms that significantly impact leukemia pathogenesis.

Natural Killer Cell Receptor Genes in Camels: Another Mammalian Model

Due to production of special homodimeric heavy chain antibodies, somatic hypermutation of their T-cell receptor genes and unusually low diversity of their major histocompatibility complex genes, camels represent an important model for immunogenetic studies. Here, we analyzed genes encoding selected natural killer cell receptors with a special focus on genes encoding receptors for major histocompatibility complex (MHC) class I ligands in the two domestic camel species, Camelus dromedarius and Camelus bactrianus. Based on the dromedary genome assembly CamDro2, we characterized the genetic contents, organization, and variability of two complex genomic regions, the leukocyte receptor complex and the natural killer complex, along with the natural cytotoxicity receptor genes NCR1, NCR2, and NCR3. The genomic organization of the natural killer complex region of camels differs from cattle, the phylogenetically most closely related species. With its minimal set of KLR genes, it resembles this complex in the domestic pig. Similarly, the leukocyte receptor complex of camels is strikingly different from its cattle counterpart. With KIR pseudogenes and few LILR genes, it seems to be simpler than in the pig. The syntenies and protein sequences of the NCR1, NCR2, and NCR3 genes in the dromedary suggest that they could be human orthologues. However, only NCR1 and NCR2 have a structure of functional genes, while NCR3 appears to be a pseudogene. High sequence similarities between the two camel species as well as with the alpaca Vicugna pacos were observed. The polymorphism in all genes analyzed seems to be generally low, similar to the rest of the camel genomes. This first report on natural killer cell receptor genes in camelids adds new data to our understanding of specificities of the camel immune system and its functions, extends our genetic knowledge of the innate immune variation in dromedaries and Bactrian camels, and contributes to studies of natural killer cell receptors evolution in mammals.

The novel KIR2DL1 allele, KIR2DL1*037, defined in the cell line SPO010 (IHW9036)

The novel KIR2DL1*037 allele discovered and characterised by single molecule real-time (SMRT) DNA sequencing.

Human Natural Killer Receptors, Co-Receptors, and Their Ligands

In the last 20 years, the study of human natural killer (NK) cells has moved from the first molecular characterizations of very few receptor molecules to the identification of a plethora of receptors displaying surprisingly divergent functions. We have contributed to the description of inhibitory receptors and their signaling pathways, important in fine regulation in many cell types, but unknown until their discovery in the NK cells. Inhibitory function is central to regulating NK-mediated cytolysis, with different molecular structures evolving during speciation to assure its persistence. More recently, it has become possible to characterize the NK triggering receptors mediating natural cytotoxicity, unveiling the existence of a network of cellular interactions between effectors of both natural and adaptive immunity. This unit reviews the contemporary history of molecular studies of receptors and ligands involved in NK cell function, characterizing the ligands of the triggering receptor and the mechanisms for finely regulating their expression in pathogen-infected or tumor cells. © 2018 by John Wiley & Sons, Inc.

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.

The effect of LILRB1 but not LILRA3 gene polymorphism in immunopathology of ankylosing spondylitis-A parallel to KIR genes

LILR and KIR receptors recognize HLA-B27 and may influence immune response in ankylosing spondylitis (AS) development. Purpose of the study was to analyse LILRB1/LILRA3 polymorphisms in AS. We observed a possible protective effect of the T allele of LILRB1 rs1061680:T>C and no association with insertion/deletion polymorphisms of LILRA3 with AS.

Two to Tango: Co-evolution of Hominid Natural Killer Cell Receptors and MHC

Natural killer (NK) cells have diverse roles in hominid immunity and reproduction. Modulating these functions are the interactions between major histocompatibility complex (MHC) class I molecules that are ligands for two NK cell surface receptor types. Diverse killer cell immunoglobulin-like receptors (KIR) bind specific motifs encoded within the polymorphic MHC class I cell surface glycoproteins, while, in more conserved interactions, CD94:NKG2A receptors recognize MHC-E with bound peptides derived from MHC class I leader sequences. The hominid lineage presents a choreographed co-evolution of KIR with their MHC class I ligands. MHC-A, -B, and -C are present in all great apes with species-specific haplotypic variation in gene content. The Bw4 epitope recognized by lineage II KIR is restricted to MHC-B but also present on some gorilla and human MHC-A. Common to great apes, but rare in humans, are MHC-B possessing a C1 epitope recognized by lineage III KIR. MHC-C arose from duplication of MHC-B and is fixed in all great apes except orangutan, where it exists on approximately 50% of haplotypes and all allotypes are C1-bearing. Recent study showed that gorillas possess yet another intermediate MHC organization compared to humans. Like orangutans, but unlike the Pan-Homo species, duplication of MHC-B occurred. However, MHC-C is fixed, and the MHC-C C2 epitope (absent in orangutans) emerges. The evolution of MHC-C drove expansion of its cognate lineage III KIR. Recently, position −21 of the MHC-B leader sequence has been shown to be critical in determining NK cell educational outcome. In humans, methionine (−21M) results in CD94:NKG2A-focused education whereas threonine (−21T) produces KIR-focused education. This is another dynamic position among hominids. Orangutans have exclusively −21M, consistent with their intermediate stage in lineage III KIR-focused evolution. Gorillas have both −21M and −21T, like humans, but they are unequally encoded by their duplicated B genes. Chimpanzees have near-fixed −21T, indicative of KIR-focused NK education. Harmonious with this observation, chimpanzee KIR exhibit strong binding and, compared to humans, smaller differences between binding levels of activating and inhibitory KIR. Consistent between these MHC-NK cell receptor systems over the course of hominid evolution is the evolution of polymorphism favoring the more novel and dynamic KIR system.

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.

Inhibitory Receptor Crosslinking Quantitatively Dampens Calcium Flux Induced by Activating Receptor Triggering in NK Cells

Natural killer (NK) cell function is regulated by a balance between activating and inhibitory receptors, but the details of this receptor interplay are not extensively understood. We developed a flow cytometry-based assay system in which Ca²⁺ flux downstream of antibody-mediated activating receptor triggering was studied in the presence or absence of inhibitory receptor co-crosslinking. We show that the inhibitory influence on activating receptor-induced Ca²⁺ flux is quantitatively regulated, both on murine and human NK cells. Furthermore, both activating and inhibitory receptors operate in an additive way, suggesting that a fine-tuned balance between activating and inhibitory receptors regulate proximal NK cell signaling. We also demonstrate that murine NK cell expression of H2D^d lowered the capacity of Ly49A to deliver inhibitory signals after antibody crosslinking, suggesting that the cis interaction between H2D^d and Ly49A reduce the signaling capacity of Ly49A in this setting. Finally, we show that priming of NK cells by IL-15 rapidly augments Ca²⁺ flux after activating receptor signaling without attenuating the potential of inhibitory receptors to reduce Ca²⁺ flux. Our data shed new light on NK cell inhibition and raises new questions for further studies.

Exploiting NK Cell Surveillance Pathways for Cancer Therapy

Natural killer (NK) cells can evoke potent anti-tumour activity. This function is largely mediated through a battery of specialised cell-surface receptors which probe the tissue microenvironment for changes in surface and secretory phenotypes that may alert to the presence of infection or malignancy. These receptors have the potential to arouse the robust cytotoxic and cytokine-secreting functions of NK cells and so must be tightly regulated to prevent autoimmunity. However, such functions also hold great promise for clinical intervention. In this review, we highlight some of the latest breakthroughs in fundamental NK cell receptor biology that have illuminated our understanding of the molecular strategies NK cells employ to perceive malignant cells from normal healthy cells. Moreover, we highlight how these sophisticated tumour recognition strategies are being harnessed for cancer immunotherapies in the clinic.

Cost-effective and fast KIR gene-content genotyping by multiplex melting curve analysis

Killer cell immunoglobulin-like receptor (KIR) genes encode cell surface molecules that recognize HLA molecules and modulate the activity of natural killer (NK) cells. KIR genes exhibit presence and absence polymorphism, which generates a variety of gene-content haplotypes in worldwide populations. KIR gene-content variation is implicated in many diseases and is also important for placentation and transplantation. Because of the complexity of KIR polymorphism, variation in this family is still mostly studied at the gene-content level, even with the advent of next-generation sequencing (NGS) methods. Gene-content determination is generally expensive and/or time-consuming. To overcome these difficulties, we developed a method based on multiplex polymerase chain reaction with specific sequence primers (PCR-SSP) followed by melting curve analysis that allows cost-effective, precise and fast generation of results. Our method was 100% concordant with a gel-based method and 99.9% concordant with presence and absence determination by NGS. The limit of detection for accurate typing was 30 ng of DNA (0.42 μM) with 260/230 and 260/280 ratios as low as 0.19 and of 0.44. In addition, we developed a user-friendly Java-based computational application called killerPeak that interprets the raw data generated by Viia7 or QuantStudio 7 quantitative PCR machines and reliably exports the final genotyping results in spreadsheet file format. The combination of a reliable method that requires low amount of DNA with an automated interpretation of results allows scaling the KIR genotyping in large cohorts with reduced turnaround time.

Fetal human leukocyte antigen-C and maternal killer-cell immunoglobulin-like receptors in cases of severe preeclampsia

INTRODUCTION

The pathogenesis of preeclampsia may involve inadequate trophoblast invasion caused by excessive inhibition of uterine natural killer cells (uNK) by extravillous trophoblast cells (EVT). This may be the result of a combination of maternal killer-cell immunoglobin-like receptor (KIR) AA genotype and fetal human leukocyte antigen-C2 (HLA-C2) genotype. A few studies have reported a significantly increased frequency of the maternal KIR AA/fetal HLA-C2 combination in cases of preeclampsia compared to controls.

METHODS

Study subjects were 259 cases of severe preeclampsia/eclampsia and 259 matched pregnant women without preeclampsia or eclampsia. All pregnancies were singleton pregnancies, and mothers were preferentially primigravidae. Blood samples from women and their newborns were obtained from the Danish National Birth Cohort (DNBC) and the Danish Neonatal Screening Biobank. Significant differences in the frequencies of KIR AA and HLA-C2 between cases and controls were investigated.

RESULTS

No significant difference was observed between cases and controls in the frequency of maternal KIR AA (OR = 0.86, 95%CI = 0.60-1.23, P = 0.41), neither when the fetus carried an HLA-C2 allele (OR = 0.85, 95%CI = 0.52-1.38, P = 0.51), nor when the fetus carried an HLA-C2 allele more than its mother (OR = 0.75, 95%CI = 0.34-1.64, P = 0.47).

CONCLUSION

The Results show no influence of HLA-C/KIR genetic variation on the risk of severe preeclampsia, contrary to what some previous studies have observed. An explanation could be that severe preeclampsia represents a separate pathological entity compared to mild preeclampsia.

Polymicrobial sepsis influences NK-cell-mediated immunity by diminishing NK-cell-intrinsic receptor-mediated effector responses to viral ligands or infections

The sepsis-induced cytokine storm leads to severe lymphopenia and reduced effector capacity of remaining/surviving cells. This results in a prolonged state of immunoparalysis, that contributes to enhanced morbidity/mortality of sepsis survivors upon secondary infection. The impact of sepsis on several lymphoid subsets has been characterized, yet its impact on NK-cells remains underappreciated-despite their critical role in controlling infection(s). Here, we observed numerical loss of NK-cells in multiple tissues after cecal-ligation-and-puncture (CLP)-induced sepsis. To elucidate the sepsis-induced lesions in surviving NK-cells, transcriptional profiles were evaluated and indicated changes consistent with impaired effector functionality. A corresponding deficit in NK-cell capacity to produce effector molecules following secondary infection and/or cytokine stimulation (IL-12,IL-18) further suggested a sepsis-induced NK-cell intrinsic impairment. To specifically probe NK-cell receptor-mediated function, the activating Ly49H receptor, that recognizes the murine cytomegalovirus (MCMV) m157 protein, served as a model receptor. Although relative expression of Ly49H receptor did not change, the number of Ly49H+ NK-cells in CLP hosts was reduced leading to impaired in vivo cytotoxicity and the capacity of NK-cells (on per-cell basis) to perform Ly49H-mediated degranulation, killing, and effector molecule production in vitro was also severely reduced. Mechanistically, Ly49H adaptor protein (DAP12) activation and clustering, assessed by TIRF microscopy, was compromised. This was further associated with diminished AKT phosphorylation and capacity to flux calcium following receptor stimulation. Importantly, DAP12 overexpression in NK-cells restored Ly49H/D receptors-mediated effector functions in CLP hosts. Finally, as a consequence of sepsis-dependent numerical and functional lesions in Ly49H+ NK-cells, host capacity to control MCMV infection was significantly impaired. Importantly, IL-2 complex (IL-2c) therapy after CLP improved numbers but not a function of NK-cells leading to enhanced immunity to MCMV challenge. Thus, the sepsis-induced immunoparalysis state includes numerical and NK-cell-intrinsic functional impairments, an instructive notion for future studies aimed in restoring NK-cell immunity in sepsis survivors.

The unique evolution of the pig LRC, a single KIR but expansion of LILR and a novel Ig receptor family

The leukocyte receptor complex (LRC) encodes numerous immunoglobulin (Ig)-like receptors involved in innate immunity. These include the killer-cell Ig-like receptors (KIR) and the leukocyte Ig-like receptors (LILR) which can be polymorphic and vary greatly in number between species. Using the recent long-read genome assembly, Sscrofa11.1, we have characterized the porcine LRC on chromosome 6. We identified a ~ 197-kb region containing numerous LILR genes that were missing in previous assemblies. Out of 17 such LILR genes and fragments, six encode functional proteins, of which three are inhibitory and three are activating, while the majority of pseudogenes had the potential to encode activating receptors. Elsewhere in the LRC, between FCAR and GP6, we identified a novel gene that encodes two Ig-like domains and a long inhibitory intracellular tail. Comparison with two other porcine assemblies revealed a second, nearly identical, non-functional gene encoding a short intracellular tail with ambiguous function. These novel genes were found in a diverse range of mammalian species, including a pseudogene in humans, and typically consist of a single long-tailed receptor and a variable number of short-tailed receptors. Using porcine transcriptome data, both the novel inhibitory gene and the LILR were highly expressed in peripheral blood, while the single KIR gene, KIR2DL1, was either very poorly expressed or not at all. These observations are a prerequisite for improved understanding of immune cell functions in the pig and other species.

Nomenclature for the KIR of non-human species

The increasing number of Killer Immunoglobulin-like Receptor (KIR) sequences available for non-human primate species and cattle has prompted development of a centralized database, guidelines for a standardized nomenclature, and minimum requirements for database submission. The guidelines and nomenclature are based on those used for human KIR and incorporate modifications made for inclusion of non-human species in the companion IPD-NHKIR database. Included in this first release are the rhesus macaque (Macaca mulatta), chimpanzee (Pan troglodytes), orangutan (Pongo abelii and Pongo pygmaeus), and cattle (Bos taurus).