HLA and KIR genetic association and NK cells in anti-NMDAR encephalitis.

Introduction: Genetic predisposition to autoimmune encephalitis with antibodies against N-methyl-D-aspartate receptor (NMDAR) is poorly understood. Given the diversity of associated environmental factors (tumors, infections), we hypothesized that human leukocyte antigen (HLA) and killer-cell immunoglobulin-like receptors (KIR), two extremely polymorphic gene complexes key to the immune system, might be relevant for the genetic predisposition to anti-NMDAR encephalitis. Notably, KIR are chiefly expressed by Natural Killer (NK) cells, recognize distinct HLA class I allotypes and play a major role in anti-tumor and anti-infection responses. Methods: We conducted a Genome Wide Association Study (GWAS) with subsequent control-matching using Principal Component Analysis (PCA) and HLA imputation, in a multi-ethnic cohort of anti-NMDAR encephalitis (n=479); KIR and HLA were further sequenced in a large subsample (n=323). PCA-controlled logistic regression was then conducted for carrier frequencies (HLA and KIR) and copy number variation (KIR). HLA-KIR interaction associations were also modeled. Additionally, single cell sequencing was conducted in peripheral blood mononuclear cells from 16 cases and 16 controls, NK cells were sorted and phenotyped. Results: Anti-NMDAR encephalitis showed a weak HLA association with DRB1*01:01~DQA1*01:01~DQB1*05:01 (OR=1.57, 1.51, 1.45; respectively), and DRB1*11:01 (OR=1.60); these effects were stronger in European descendants and in patients without an underlying ovarian teratoma. More interestingly, we found increased copy number variation of KIR2DL5B (OR=1.72), principally due to an overrepresentation of KIR2DL5B*00201. Further, we identified two allele associations in framework genes, KIR2DL4*00103 (25.4% vs. 12.5% in controls, OR=1.98) and KIR3DL3*00302 (5.3% vs. 1.3%, OR=4.44). Notably, the ligands of these KIR2DL4 and KIR3DL3, respectively, HLA-G and HHLA2, are known to act as immune checkpoint with immunosuppressive functions. However, we did not find differences in specific KIR-HLA ligand interactions or HLA-G polymorphisms between cases and controls. Similarly, gene expression of CD56dim or CD56bright NK cells did not differ between cases and controls. Discussion: Our observations for the first time suggest that the HLA-KIR axis might be involved in anti-NMDAR encephalitis. While the genetic risk conferred by the identified polymorphisms appears small, a role of this axis in the pathophysiology of this disease appears highly plausible and should be analyzed in future studies.

MHConstructor: A high-throughput, haplotype-informed solution to the MHC assembly challenge.

The extremely high levels of genetic polymorphism within the human major histocompatibility complex (MHC) limit the usefulness of reference-based alignment methods for sequence assembly. We incorporate a short read de novo assembly algorithm into a workflow for novel application to the MHC. MHConstructor is a containerized pipeline designed for high-throughput, haplotype-informed, reproducible assembly of both whole genome sequencing and target-capture short read data in large, population cohorts. To-date, no other self-contained tool exists for the generation of de novo MHC assemblies from short read data. MHConstructor facilitates wide-spread access to high quality, alignment-free MHC sequence analysis.

Donor KIR2DL1 Allelic Polymorphism Influences Posthematopoietic Progenitor Cell Transplantation Outcomes in the T Cell Depleted and Reduced Intensity Conditioning Setting.

The majority of established KIR clinical assessment algorithms used for donor selection for hematopoietic progenitor cell transplantation (HPCT) evaluate gene content (presence/absence) of the KIR gene complex. In comparison, relatively little is known about the impact of KIR allelic polymorphism. By analyzing donors of T cell depleted (TcD) reduced intensity conditioning (RIC) HPCT, this study investigated the influence on post-transplant outcome of 2 polymorphic residues of the inhibitory KIR2DL1. The aim of this study was to expand upon existing research into the influence of KIR2DL1 allelic polymorphism upon post-transplant outcome. The effects of allele groups upon transplant outcomes were investigated within a patient cohort using a defined treatment protocol of RIC with TcD. Using phylogenetic data, KIR2DL1 allelic polymorphism was categorized into groups on the basis of variation within codons 114 and 245 (positive or negative for the following groups: KIR2DL1*002/001g, KIR2DL1*003, KIR2DL1*004g) and the identification of null alleles. The influence of these KIR2DL1 allele groups in hematopoietic progenitor cell transplantation (HPCT) donors was assessed in the post-transplant data of 86 acute myelogenous leukemia patients receiving RIC TcD HPCT at a single center. KIR2DL1 allele groups in the donor significantly impacted upon 5-year post-transplant outcomes in RIC TcD HPCT. Donor KIR2DL1*003 presented the greatest influence upon post-transplant outcomes, with KIR2DL1*003 positive donors severely reducing 5-year post-transplant overall survival (OS) compared to those receiving a transplant from a KIR2DL1*003 negative donor (KIR2DL1*003 pos versus neg: 27.0% versus 60.0%, P = .008, pc = 0.024) and disease-free survival (DFS) (KIR2DL1*003 pos versus neg: 23.5% versus 60.0%, P = .004, pc = 0.012), and increasing 5-year relapse incidence (KIR2DL1*003 pos versus neg: 63.9% versus 27.2%, P = .009, pc = 0.027). KIR2DL1*003 homozygous and KIR2DL1*003 heterozygous grafts did not present significantly different post-transplant outcomes. Donors possessing the KIR2DL1*002/001 allele group were found to significantly improve post-transplant outcomes, with donors positive for the KIR2DL1*004 allele group presenting a trend towards improvement. KIR2DL1*002/001 allele group (KIR2DL1*002/001g) positive donors improved 5-year OS (KIR2DL1*002/001g pos versus neg: 56.4% versus 27.2%, P = .009, pc = 0.024) and DFS (KIR2DL1*002/001g pos versus neg: 53.8% versus 25.5%, P = .018, pc = 0.036). KIR2DL1*004 allele group (KIR2DL1*004g) positive donors trended towards improving 5-year OS (KIR2DL1*004g pos versus neg: 53.3% versus 35.5%, P = .097, pc = 0.097) and DFS (KIR2DL1*004g pos versus neg: 50.0% versus 33.9%, P = .121, pc = 0.121), and reducing relapse incidence (KIR2DL1*004g pos versus neg: 33.1% versus 54.0%, P = .079, pc = 0.152). The presented findings suggest donor selection algorithms for TcD RIC HPCT should consider avoiding KIR2DL1*003 positive donors, where possible, and contributes to the mounting evidence that KIR assessment in donor selection algorithms should reflect the conditioning regime protocol used.

HLA-DPB1*13:01 associates with enhanced, and KIR2DS4*001 with diminished protection from developing severe COVID-19.

Extreme polymorphism of HLA and killer-cell immunoglobulin-like receptors (KIR) differentiates immune responses across individuals. Additional to T cell receptor interactions, subsets of HLA class I act as ligands for inhibitory and activating KIR, allowing natural killer (NK) cells to detect and kill infected cells. We investigated the impact of HLA and KIR polymorphism on the severity of COVID-19. High resolution HLA class I and II and KIR genotypes were determined from 403 non-hospitalized and 1575 hospitalized SARS-CoV-2 infected patients from Italy collected in 2020. We observed that possession of the activating KIR2DS4*001 allotype is associated with severe disease, requiring hospitalization (OR = 1.48, 95% CI 1.20-1.85, pc = 0.017), and this effect is greater in individuals homozygous for KIR2DS4*001 (OR = 3.74, 95% CI 1.75-9.29, pc = 0.003). We also observed the HLA class II allotype, HLA-DPB1*13:01 protects SARS-CoV-2 infected patients from severe disease (OR = 0.49, 95% CI 0.33-0.74, pc = 0.019). These association analyses were replicated using logistic regression with sex and age as covariates. Autoantibodies against IFN-α associated with COVID-19 severity were detected in 26% of 156 hospitalized patients tested. HLA-C*08:02 was more frequent in patients with IFN-α autoantibodies than those without, and KIR3DL1*01502 was only present in patients lacking IFN-α antibodies. These findings suggest that KIR and HLA polymorphism is integral in determining the clinical outcome following SARS-CoV-2 infection, by influencing the course both of innate and adaptive immunity.

Complete sequences of six major histocompatibility complex haplotypes, including all the major MHC class II structures.

Accurate and comprehensive immunogenetic reference panels are key to the successful implementation of population-scale immunogenomics. The 5Mbp Major Histocompatibility Complex (MHC) is the most polymorphic region of the human genome and associated with multiple immune-mediated diseases, transplant matching and therapy responses. Analysis of MHC genetic variation is severely complicated by complex patterns of sequence variation, linkage disequilibrium and a lack of fully resolved MHC reference haplotypes, increasing the risk of spurious findings on analyzing this medically important region. Integrating Illumina, ultra-long Nanopore, and PacBio HiFi sequencing as well as bespoke bioinformatics, we completed five of the alternative MHC reference haplotypes of the current (GRCh38/hg38) build of the human reference genome and added one other. The six assembled MHC haplotypes encompass the DR1 and DR4 haplotype structures in addition to the previously completed DR2 and DR3, as well as six distinct classes of the structurally variable C4 region. Analysis of the assembled haplotypes showed that MHC class II sequence structures, including repeat element positions, are generally conserved within the DR haplotype supergroups, and that sequence diversity peaks in three regions around HLA-A, HLA-B+C, and the HLA class II genes. Demonstrating the potential for improved short-read analysis, the number of proper read pairs recruited to the MHC was found to be increased by 0.06%-0.49% in a 1000 Genomes Project read remapping experiment with seven diverse samples. Furthermore, the assembled haplotypes can serve as references for the community and provide the basis of a structurally accurate genotyping graph of the complete MHC region.

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.