The killer immunoglobulin-like receptor KIR3DL1 in combination with HLA-Bw4 is protective against multiple sclerosis in African Americans

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.

Identification of potential glioma drug resistance target proteins based on ultra-performance liquid chromatography-mass spectrometry differential proteomics

In this study, to screen for candidate markers of temozolomide (TMZ) resistance in glioblastoma, we artificially established TMZ drug-resistant glioblastoma (GBM) cell lines, U251-TMZ and U87-TMZ. In the U251-TMZ and U87-TMZ cell lines, we screened and analyzed differentially expressed proteins using ultra-performance liquid chromatography-mass spectrometry (UPLC-MS) differential proteomics. Compared with the U251 and U87 control cell lines, 95 differential proteins were screened in the U251-TMZ and U87-TMZ cell lines, of which 28 proteins were upregulated and 67 proteins were down-regulated. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of the co-upregulated proteins showed that most of the differentially expressed proteins were located in the cytoplasm and were significantly upregulated in the biological processes related to vesicular transport in the intimal system and inflammatory response mediated by myeloid leukocytes. Seven candidates were identified as potential GBM markers of TMZ resistance. Combined with existing research findings, our study supports that UAP1L1 and BCKDK are promising potential markers of TMZ resistance in GBM. This is important for further understanding the molecular mechanisms that drive the development and enhancement of TMZ resistance.

Is the Development of Ascites in Alcoholic Liver Patients Influenced by Specific KIR/HLA Gene Profiles?

Decompensated cirrhosis is the most common cause of ascites due to hemodynamic and renal alteration by continuous fluid leakage from the hepatic sinusoids and splanchnic capillaries into the interstitial space. Then, fluid leakage exceeds lymphatic return, leading to progressive fluid accumulation directly into the peritoneal cavity. Alcohol consumption is one of the main risks of developing alcoholic cirrhosis (AC), but not all AC patients develop ascites. Avoiding the development of ascites is crucial, given that it deteriorates prognosis and increases the patient mortality patient. The innate immune system plays a crucial role in cirrhosis through natural killer cells, which are abundant in the liver. The aim of this study was to analyze the KIR/HLA-C genetic profile in AC patients with and without ascites to understand this pathology and find predictive clinical susceptibility biomarkers that can help to establish risks and prevent the development of ascites in AC patients. A total of 281 AC patients with and without ascites were analyzed and compared with 319 healthy controls. Genomic DNA was extracted from peripheral blood in all groups. A PCR-SSO assay was performed for KIR/HLA genotyping analysis. A total of 16 activating and inhibitor KIR genes and their corresponding known ligands, epitopes of HLA-C, and their genotypes were analyzed. According to our analysis, C1 epitopes were statistically significantly decreased in AC patients with and without ascites. When comparing AC patients with ascites and healthy controls, a significant decrease in C1 epitope frequency was also observed. A statistically significant decrease was also found when comparing the C1C2 genotype in AC patients without ascites with controls. In conclusion, the absence of KIR2DL2 and KIR3DL1 genes may be a predisposing factor for the development of ascites in AC patients. The KIR2DS2/KIR2DL2 may could be involved in grade I ascites development, and the presence of the C1+ epitope and the homozygous C2C2 genotype may be protective genetic factors against ascites development in AC patients.

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.

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.

Contribution of HLA and KIR Alleles to Systemic Sclerosis Susceptibility and Immunological and Clinical Disease Subtypes

Systemic sclerosis (SSc) is an autoinflammatory, fibrotic condition of unknown aetiology. The presence of detectable autoantibodies against diverse nuclear antigens, as well as strong HLA associations with disease, suggest autoimmune involvement, however the links between endogenous and exogenous risk factors and SSc pathology remain undetermined. We have conducted a genetic analysis of HLA inheritance in two independent and meta-analysed cohorts of 1,465 SSc cases and 13,273 controls, including stratified association analyses in clinical and autoantibody positive subgroups of disease. Additionally, we have used patient genotypes to impute gene dosages across the KIR locus, encoding paired activating and inhibitory lymphocyte receptors for Class I HLA ligands, to conduct the largest analysis of KIR-HLA epistatic interactions in SSc to date. We confirm previous Class II HLA associations with SSc risk and report a new Class I association with haplotype HLA-B*44:03-HLA-C*16:01 at genome-wide significance (GWS). We further report statistically significant HLA associations with clinical and serological subtypes of disease through direct case-case comparison, and report a new association of HLA-DRB1*15:01, previously shown to bind topoisomerase-1 derived peptides, with anti-topoisomerase (ATA) positive disease. Finally, we identify genetic epistasis between KIRs and HLA class I ligands, suggesting genetic modulation of lymphocyte activation may further contribute to an individual’s underlying disease risk. Taken together, these findings support future functional investigation into endogenous immunological and environmental stimuli for disrupted immune tolerance in SSc.

Allele imputation for the killer cell immunoglobulin-like receptor KIR3DL1/S1

Highly polymorphic interaction of KIR3DL1 and KIR3DS1 with HLA class I ligands modulates the effector functions of natural killer (NK) cells and some T cells. This genetically determined diversity affects severity of infections, immune-mediated diseases, and some cancers, and impacts the course of immunotherapies, including transplantation. KIR3DL1 is an inhibitory receptor, and KIR3DS1 is an activating receptor encoded by the KIR3DL1/S1 gene that has more than 200 diverse and divergent alleles. Determination of KIR3DL1/S1 genotypes for medical application is hampered by complex sequence and structural variation, requiring targeted approaches to generate and analyze high-resolution allele data. To overcome these obstacles, we developed and optimized a model for imputing KIR3DL1/S1 alleles at high-resolution from whole-genome SNP data. We designed the model to represent a substantial component of human genetic diversity. Our Global imputation model is effective at genotyping KIR3DL1/S1 alleles with an accuracy ranging from 88% in Africans to 97% in East Asians, with mean specificity of 99% and sensitivity of 95% for alleles >1% frequency. We used the established algorithm of the HIBAG program, in a modification named Pulling Out Natural killer cell Genomics (PONG). Because HIBAG was designed to impute HLA alleles also from whole-genome SNP data, PONG allows combinatorial diversity of KIR3DL1/S1 with HLA-A and -B to be analyzed using complementary techniques on a single data source. The use of PONG thus negates the need for targeted sequencing data in very large-scale association studies where such methods might not be tractable.

Capturing SNP Association across the NK Receptor and HLA Gene Regions in Multiple Sclerosis by Targeted Penalised Regression Models

Conventional genome-wide association studies (GWASs) of complex traits, such as Multiple Sclerosis (MS), are reliant on per-SNP p-values and are therefore heavily burdened by multiple testing correction. Thus, in order to detect more subtle alterations, ever increasing sample sizes are required, while ignoring potentially valuable information that is readily available in existing datasets. To overcome this, we used penalised regression incorporating elastic net with a stability selection method by iterative subsampling to detect the potential interaction of loci with MS risk. Through re-analysis of the ANZgene dataset (1617 cases and 1988 controls) and an IMSGC dataset as a replication cohort (1313 cases and 1458 controls), we identified new association signals for MS predisposition, including SNPs above and below conventional significance thresholds while targeting two natural killer receptor loci and the well-established HLA loci. For example, rs2844482 (98.1% iterations), otherwise ignored by conventional statistics (p = 0.673) in the same dataset, was independently strongly associated with MS in another GWAS that required more than 40 times the number of cases (~45 K). Further comparison of our hits to those present in a large-scale meta-analysis, confirmed that the majority of SNPs identified by the elastic net model reached conventional statistical GWAS thresholds (p < 5 × 10−8) in this much larger dataset. Moreover, we found that gene variants involved in oxidative stress, in addition to innate immunity, were associated with MS. Overall, this study highlights the benefit of using more advanced statistical methods to (re-)analyse subtle genetic variation among loci that have a biological basis for their contribution to disease risk.

Analysis of whole exome sequencing in severe mental illness hints at selection of brain development and immune related genes

Evolutionary trends may underlie some aspects of the risk for common, non-communicable disorders, including psychiatric disease. We analyzed whole exome sequencing data from 80 unique individuals from India coming from families with two or more individuals with severe mental illness. We used Population Branch Statistics (PBS) to identify variants and genes under positive selection and identified 74 genes as candidates for positive selection. Of these, 20 were previously associated with Schizophrenia, Alzheimer’s disease and cognitive abilities in genome wide association studies. We then checked whether any of these 74 genes were involved in common biological pathways or related to specific cellular or molecular functions. We found that immune related pathways and functions related to innate immunity such as antigen binding were over-represented. We also evaluated for the presence of Neanderthal introgressed segments in these genes and found Neanderthal introgression in a single gene out of the 74 candidate genes. However, the introgression pattern indicates the region is unlikely to be the source for selection. Our findings hint at how selection pressures in individuals from families with a history of severe mental illness may diverge from the general population. Further, it also provides insights into the genetic architecture of severe mental illness, such as schizophrenia and its link to immune factors.

An immunogenetic view of COVID-19

Meeting the challenges brought by the COVID-19 pandemic requires an interdisciplinary approach. In this context, integrating knowledge of immune function with an understanding of how genetic variation influences the nature of immunity is a key challenge. Immunogenetics can help explain the heterogeneity of susceptibility and protection to the viral infection and disease progression. Here, we review the knowledge developed so far, discussing fundamental genes for triggering the innate and adaptive immune responses associated with a viral infection, especially with the SARS-CoV-2 mechanisms. We emphasize the role of the HLA and KIR genes, discussing what has been uncovered about their role in COVID-19 and addressing methodological challenges of studying these genes. Finally, we comment on questions that arise when studying admixed populations, highlighting the case of Brazil. We argue that the interplay between immunology and an understanding of genetic associations can provide an important contribution to our knowledge of COVID-19.

High Resolution Haplotype Analyses of Classical HLA Genes in Families With Multiple Sclerosis Highlights the Role of HLA-DP Alleles in Disease Susceptibility

Multiple sclerosis (MS) susceptibility shows strong genetic associations with HLA alleles and haplotypes. We genotyped 11 HLA genes in 477 non-Hispanic European MS patients and their 954 unaffected parents using a validated next-generation sequencing (NGS) methodology. HLA haplotypes were assigned unequivocally by tracing HLA allele transmissions. We explored HLA haplotype/allele associations with MS using the genotypic transmission disequilibrium test (gTDT) and multiallelic TDT (mTDT). We also conducted a case-control (CC) study with all patients and 2029 healthy unrelated ethnically matched controls. We performed separate analyses of 54 extended multi-case families by reviewing transmission of haplotype blocks. The haplotype fragment including DRB5*01:01:01~DRB1*15:01:01:01 was significantly associated with predisposition (gTDT: p < 2.20e-16; mTDT: p =1.61e-07; CC: p < 2.22e-16) as reported previously. A second risk allele, DPB1*104:01 (gTDT: p = 3.69e-03; mTDT: p = 2.99e-03; CC: p = 1.00e-02), independent from the haplotype bearing DRB1*15:01 was newly identified. The allele DRB1*01:01:01 showed significant protection (gTDT: p = 8.68e-06; mTDT: p = 4.50e-03; CC: p = 1.96e-06). Two DQB1 alleles, DQB1*03:01 (gTDT: p = 2.86e-03; mTDT: p = 5.56e-02; CC: p = 4.08e-05) and DQB1*03:03 (gTDT: p = 1.17e-02; mTDT: p = 1.16e-02; CC: p = 1.21e-02), defined at two-field level also showed protective effects. The HLA class I block, A*02:01:01:01~C*03:04:01:01~B*40:01:02 (gTDT: p = 5.86e-03; mTDT: p = 3.65e-02; CC: p = 9.69e-03) and the alleles B*27:05 (gTDT: p = 6.28e-04; mTDT: p = 2.15e-03; CC: p = 1.47e-02) and B*38:01 (gTDT: p = 3.20e-03; mTDT: p = 6.14e-03; CC: p = 1.70e-02) showed moderately protective effects independently from each other and from the class II associated factors. By comparing statistical significance of 11 HLA loci and 19 haplotype segments with both untruncated and two-field allele names, we precisely mapped MS candidate alleles/haplotypes while eliminating false signals resulting from ‘hitchhiking’ alleles. We assessed genetic burden for the HLA allele/haplotype identified in this study. This family-based study including the highest-resolution of HLA alleles proved to be powerful and efficient for precise identification of HLA genotypes associated with both, susceptibility and protection to development of MS.

High-Resolution Characterization of KIR Genes in a Large North American Cohort Reveals Novel Details of Structural and Sequence Diversity

The KIR (killer-cell immunoglobulin-like receptor) region is characterized by structural variation and high sequence similarity among genes, imposing technical difficulties for analysis. We undertook the most comprehensive study to date of KIR genetic diversity in a large population sample, applying next-generation sequencing in 2,130 United States European-descendant individuals. Data were analyzed using our custom bioinformatics pipeline specifically designed to address technical obstacles in determining KIR genotypes. Precise gene copy number determination allowed us to identify a set of uncommon gene-content KIR haplotypes accounting for 5.2% of structural variation. In this cohort, KIR2DL4 is the framework gene that most varies in copy number (6.5% of all individuals). We identified phased high-resolution alleles in large multi-locus insertions and also likely founder haplotypes from which they were deleted. Additionally, we observed 250 alleles at 5-digit resolution, of which 90 have frequencies ≥1%. We found sequence patterns that were consistent with the presence of novel alleles in 398 (18.7%) individuals and contextualized multiple orphan dbSNPs within the KIR complex. We also identified a novel KIR2DL1 variant, Pro151Arg, and demonstrated by molecular dynamics that this substitution is predicted to affect interaction with HLA-C. No previous studies have fully explored the full range of structural and sequence variation of KIR as we present here. We demonstrate that pairing high-throughput sequencing with state-of-art computational tools in a large cohort permits exploration of all aspects of KIR variation including determination of population-level haplotype diversity, improving understanding of the KIR system, and providing an important reference for future studies.

The role of KIR positive NK cells in diseases and its importance in clinical intervention

Natural killer (NK) cells are essential for the elimination of the transformed and cancerous cells. Killer cell immunoglobulin-like receptors (KIRs) which expressed by T and NK cells, are key regulator of NK cell function. The KIR and their ligands, MHC class I (HLA-A, B and C) molecules, are highly polymorphic and their related genes are located on 19 q13.4 and 6 q21.3 chromosomes, respectively. It is clear that particular interaction between the KIRs and their related ligands can influence on the prevalence, progression and outcome of several diseases, like complications of pregnancy, viral infection, autoimmune diseases, and hematological malignancies. The mechanisms of immune signaling in particular NK cells involvement in causing pathological conditions are not completely understood yet. Therefore, better understanding of the molecular mechanism of KIR-MHC class I interaction could facilitate the treatment strategy of diseases. The present review focused on the main characteristics and functional details of various KIR and their combination with related ligands in diseases and also highlights ongoing efforts to manipulate the key checkpoints in NK cell-based immunotherapy.

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.

The Roles of Uterine Natural Killer (NK) Cells and KIR/HLA-C Combination in the Development of Preeclampsia: A Systematic Review

Preeclampsia (PE) is termed as a systemic disease that involves multiple organs; however, the exact etiology is still quite unclear. It is believed that the poor remodeling of uterine spiral arteries triggers PE, thereby causing failed placentation and producing inflammatory factors. The decline of blood flow results in lowering the nutrients and oxygen received by the fetus and augmenting the placental pressure in PE. Decidual immune cells, especially uterine natural killer (uNK) cells, are involved in the process of placentation. Decidual NK (dNK) cells significantly contribute to the vascular remodeling through the secretion of cytokines and angiogenic mediators in normal placental development. The abnormal activation of NK cells in both the peripheral blood and the decidua was counted among the causes leading to PE. The correlation existing between maternal killer cell immunoglobulin-like receptor (KIR) and HLA-C in trophoblast cells constitutes a robust evidence for the genetic etiology of PE. The combinations of the two kinds of gene systems, together with the KIR genotype in the mother and the HLA-C group in her fetus, are likely to exactly decide the pregnancy outcome. The women, who have the inappropriate match of KIR/HLA-C, are likely to be prone to the augmented risk of PE. However, the combinations of KIR/HLA-C in PE undergo ethnic changes. The extensive prospective research works in Europe, Asia, and Africa are required for providing more findings in PE patients.

KIR3DL1/S1 Allotypes Contribute Differentially to the Development of Behçet Disease

Behçet disease is a chronic, relapsing-remitting autoinflammatory syndrome with a strong HLA-B*51 association. In this paper, we describe a human cohort of 267 individuals with Behçet disease and 445 matched controls from a tertiary referral center in the U.K. HLA-B*51 was confirmed as a genetic risk factor in this group (p = 0.0006, Bonferroni-Dunn correction for multiple testing [Pc] = 0.0192, odds ratio [OR] 1.92, 95% confidence interval [CI] 1.33-2.76). KIR3DL1/S1 allele-level analysis indicated that low-expressing KIR3DL1/S1 alleles in combination with KIR3DS1 increased the risk of developing Behçet disease (KIR3DL1^LOW/KIR3DS1: p = 0.0004, Pc = 0.0040, OR 2.47, 95% CI 1.43-4.25), whereas high-expressing KIR3DL1/S1 alleles in combination with a null-expressing KIR3DL1 reduced the risk of disease (KIR3DL1^HIGH/KIR3DL1^NULL: p = 0.0035, Pc = 0.0350, OR 0.53, 95% CI 0.33-0.87). Behçet disease can manifest as a purely mucocutaneous disease or can involve other organ systems such as the eyes. In the U.K. cohort studied in this study, KIR3DL1^LOW/KIR3DS1 increased the risk of ophthalmic disease (p = 1.2 × 10^-5, OR 3.92, 95% CI 2.06-7.47), whereas KIR3DL1^HIGH/KIR3DL1^NULL reduced the risk of having purely mucocutaneous disease (p = 0.0048, OR 0.45, 95% CI 0.25-0.81). To our knowledge, this is the first analysis of KIR3DL1/S1 allelic variation in Behçet disease and may provide insight into the pathogenic role of HLA-B*51 and its interaction with KIR3DL1/S1.

Autoantibodies to killer cell immunoglobulin-like receptor 3DL1 in patients with systemic lupus erythematosus

A genetic variant of the killer immunoglobulin-like receptor 3DL1 (KIR3DL1) has been found in patients with systemic lupus erythematosus (SLE). Herein, we investigated the presence of autoantibodies to KIR3DL1 in a cohort of patients with SLE. We tested sera from 28 patients with SLE, 11 patients with rheumatoid arthritis (RA) and 17 healthy control subjects for anti-KIR3DL1 activity by an enzyme-linked immunosorbent assay (ELISA) using recombinant KIR3DL1-enhanced green fluorescent protein (EGFP) and EGFP proteins. Anti-KIR3DL1 antibodies were detected in 22 (79%) of the 28 patients with SLE, whereas they were present in only three (27%) of the 11 patients with RA examined. Notably, 10 (91%) of the 11 samples from patients with SLE prior to therapy had anti-KIR3DL1 antibodies. None of the samples from healthy donors were positive for the antibodies. Here, we report the presence of anti-KIR3DL1 antibodies in the sera of patients with SLE for the first time. Anti-KIR3DL1 autoantibodies may be involved in the pathogenesis of autoimmune diseases.

Allele-Level KIR Genotyping of More Than a Million Samples: Workflow, Algorithm, and Observations

The killer-cell immunoglobulin-like receptor (KIR) genes regulate natural killer cell activity, influencing predisposition to immune mediated disease, and affecting hematopoietic stem cell transplantation (HSCT) outcome. Owing to the complexity of the KIR locus, with extensive gene copy number variation (CNV) and allelic diversity, high-resolution characterization of KIR has so far been applied only to relatively small cohorts. Here, we present a comprehensive high-throughput KIR genotyping approach based on next generation sequencing. Through PCR amplification of specific exons, our approach delivers both copy numbers of the individual genes and allelic information for every KIR gene. Ten-fold replicate analysis of a set of 190 samples revealed a precision of 99.9%. Genotyping of an independent set of 360 samples resulted in an accuracy of more than 99% taking into account consistent copy number prediction. We applied the workflow to genotype 1.8 million stem cell donor registry samples. We report on the observed KIR allele diversity and relative abundance of alleles based on a subset of more than 300,000 samples. Furthermore, we identified more than 2,000 previously unreported KIR variants repeatedly in independent samples, underscoring the large diversity of the KIR region that awaits discovery. This cost-efficient high-resolution KIR genotyping approach is now applied to samples of volunteers registering as potential donors for HSCT. This will facilitate the utilization of KIR as additional selection criterion to improve unrelated donor stem cell transplantation outcome. In addition, the approach may serve studies requiring high-resolution KIR genotyping, like population genetics and disease association studies.

KIR2DL4-HLAG interaction at human NK cell-oligodendrocyte interfaces regulates IFN-γ-mediated effects

Interactions between germline-encoded natural killer (NK) cell receptors and their respective ligands on tumorigenic or virus-infected cells determine NK cell cytotoxic activity and/or cytokine secretion. NK cell cytokine responses can be augmented in and can potentially contribute to multiple sclerosis (MS), an inflammatory disease of the central nervous system focused upon the oligodendrocytes (OLs). To investigate mechanisms by which NK cells may contribute to MS pathogenesis, we developed an in vitro human model of OL-NK cell interaction. We found that activated, but not resting human NK cells form conjugates with, and mediate cytotoxicity against, human oligodendrocytes. NK cells, when in conjugate with OLs, rapidly synthesize and polarize IFN-γ toward the OLs. IFN-γ  is capable of reducing myelin oligodendrocyte and myelin associated glycoproteins (MOG and MAG) content. This activity is independent of MHC class-I mediated inhibition via KIR2DL1, but dependent upon the interaction between NK cell-expressed KIR2DL4 and its oligodendrocyte-expressed ligand, HLA-G. NK cells from patients with MS express higher levels of IFN-γ following conjugation to OLs, more actively promote in vitro reduction of MOG and MAG and have higher frequencies of the KIR2DL4 positive population. These data collectively suggest a mechanism by which NK cells can promote pathogenic effects upon OLs.

The immunogenetics of neurological disease

Genes encoding antigen-presenting molecules within the human major histocompatibility complex (MHC) account for the highest component of genetic risk for many neurological diseases, such as multiple sclerosis, neuromyelitis optica, Parkinson's disease, Alzheimer's disease, schizophrenia, myasthenia gravis and amyotrophic lateral sclerosis. Myriad genetic, immunological and environmental factors may contribute to an individual's susceptibility to neurological disease. Here, we review and discuss the decades long research on the influence of genetic variation at the MHC locus and the role of immunogenetic killer cell immunoglobulin-like receptor (KIR) loci in neurological diseases, including multiple sclerosis, neuromyelitis optica, Parkinson's disease, Alzheimer's disease, schizophrenia, myasthenia gravis and amyotrophic lateral sclerosis. The findings of immunogenetic association studies are consistent with a polygenic model of inheritance in the heterogeneous and multifactorial nature of complex traits in various neurological diseases. Future investigation is highly recommended to evaluate both coding and non-coding variation in immunogenetic loci using high-throughput high-resolution next-generation sequencing technologies in diverse ethnic groups to fully appreciate their role in neurological diseases.

High resolution HLA analysis reveals independent class I haplotypes and amino-acid motifs protective for multiple sclerosis

We investigated association between HLA class I and class II alleles and haplotypes, and KIR loci and their HLA class I ligands, with multiple sclerosis (MS) in 412 European-American MS patients and 419 ethnically-matched controls, using next generation sequencing. The DRB1*15:01~DQB1*06:02 haplotype was highly predisposing (odds ratio (OR) = 3.98; 95% confidence interval (CI) = 3−5.31; p-value (p) = 2.22E−16), as was DRB1*03:01~DQB1*02:01 (OR = 1.63; CI = 1.19–2.24; p = 1.41E−03). Hardy-Weinberg (HW) analysis in MS patients revealed a significant DRB1*03:01~DQB1*02:01 homozyote excess (15 observed, 8.6 expected; p = 0.016). The OR for this genotype (5.27; CI = 1.47–28.52; p = 0.0036) suggests a recessive MS risk model. Controls displayed no HW deviations. The C*03:04~B*40:01 haplotype (OR = 0.27; CI = 0.14–0.51; p = 6.76E−06) was highly protective for MS, especially in haplotypes with A*02:01 (OR = 0.15; CI = 0.04–0.45; p = 6.51E−05). By itself, A*02:01 is moderately protective, (OR = 0.69; CI = 0.54–0.87; p = 1.46E−03), and haplotypes of A*02:01 with the HLA-B Thr80 Bw4 variant (Bw4T) more so (OR = 0.53; CI = 0.35–0.78; p = 7.55E−04). Protective associations with the Bw4 KIR ligand resulted from linkage disequilibrium (LD) with DRB1*15:01, but the Bw4T variant was protective (OR = 0.64; CI = 0.49–0.82; p = 3.37E−04) independent of LD with DRB1*15:01. The Bw4I variant was not associated with MS. Overall, we find specific class I HLA polymorphisms to be protective for MS, independent of the strong predisposition conferred by DRB1*15:01.

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.