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

Unveiling the Significance of HLA and KIR Diversity in Underrepresented Populations

Human leukocyte antigen (HLA) molecules and their relationships with natural killer (NK) cells, specifically through their interaction with killer-cell immunoglobulin-like receptors (KIRs), exhibit robust associations with the outcomes of diverse diseases. Moreover, genetic variations in HLA and KIR immune system genes offer limitless depths of complexity. In recent years, a surge of high-powered genome-wide association studies (GWASs) utilizing single nucleotide polymorphism (SNP) arrays has occurred, significantly advancing our understanding of disease pathogenesis. Additionally, advances in HLA reference panels have enabled higher resolution and more reliable imputation, allowing for finer-grained evaluation of the association between sequence variations and disease risk. However, it is essential to note that the majority of these GWASs have focused primarily on populations of Caucasian and Asian origins, neglecting underrepresented populations in Latin America and Africa. This omission not only leads to disparities in health care access but also restricts our knowledge of novel genetic variants involved in disease pathogenesis within these overlooked populations. Since the KIR and HLA haplotypes prevalent in each population are clearly modelled by the specific environment, the aim of this review is to encourage studies investigating HLA/KIR involvement in infection and autoimmune diseases, reproduction, and transplantation in underrepresented populations.

High-resolution KIR and HLA genotyping in three Chinese ethnic minorities reveals distinct origins

Polymorphism of killer-cell immunoglobulin-like receptors (KIRs) and their HLA class I ligands impacts the effector activity of cytotoxic NK cell and T cell subsets. Therefore, understanding the extent and implications of KIR and HLA class I genetic polymorphism across various populations is important for immunological and medical research. In this study, we conducted a high-resolution investigation of KIR and HLA class I diversity in three distinct Chinese ethnic minority populations. We studied the She, Yugur, and Tajik, and compared them with the Zhejiang Han population (Zhe), which represents the majority Southern Han ethnicity. Our findings revealed that the Tajik population exhibited the most diverse KIR copy number, allele, and haplotype diversity among the four populations. This diversity aligns with their proposed ancestral origin, closely resembling that of Iranian populations, with a relatively higher presence of KIR-B genes, alleles, and haplotypes compared with the other Chinese populations. The Yugur population displayed KIR distributions similar to those of the Tibetans and Southeast Asians, whereas the She population resembled the Zhe and other East Asians, as confirmed by genetic distance analysis of KIR. Additionally, we identified 12.9% of individuals across the three minority populations as having KIR haplotypes characterized by specific gene block insertions or deletions. Genetic analysis based on HLA alleles yielded consistent results, even though there were extensive variations in HLA alleles. The observed variations in KIR interactions, such as higher numbers of 2DL1-C2 interactions in Tajik and Yugur populations and of 2DL3-C1 interactions in the She population, are likely shaped by demographic and evolutionary mechanisms specific to their local environments. Overall, our findings offer valuable insights into the distribution of KIR and HLA diversity among three distinct Chinese ethnic minority populations, which can inform future clinical and population studies.

Determination for KIR genotype and allele copy number via real-time quantitative PCR method

Killer cell immunoglobulin-like receptor (KIR) and human leukocyte antigen (HLA) play crucial roles in regulating NK cell activity. Here, we report a real-time quantitative PCR (qPCR) to genotype all KIR genes and their copy numbers simultaneously. With 18 pairs of locus-specific primers, we identified KIR genes by Ct values and determined KIR copy number using the 2-∆Ct method. Haplotypes were assigned based on KIR gene copy numbers. The real-time qPCR results were consistent with the NGS method, except for one sample with KIR2DL5 discrepancy. qPCR is a multiplex method that can identify KIR copy number, which helps obtain a relatively accurate haplotype structure, facilitating increased KIR research in laboratories where NGS or other high-resolution methods are not available.

High-resolution characterization of KIR genes polymorphism in healthy subjects from the Bulgarian population-A pilot study

Although killer-cell immunoglobulin-like receptor (KIR) gene content has been widely studied in health and disease, with the advancement of next-generation sequencing (NGS) technology the high-resolution characterization of this complex gene region has become achievable. KIR allele-level diversity has lately been described across human populations. The present study aimed to analyze for the first time the allele-level polymorphism of nine KIR genes in 155 healthy, unrelated individuals from the Bulgarian population by applying NGS. The highest degree of polymorphism was detected for the KIR3DL3 gene with 40 observed alleles at five-digit resolution in total, 22 of which were common. On the other hand, the KIR3DS1 gene was found to have the lowest degree of polymorphism among the studied KIR genes with one common allele: KIR3DS1*01301 (31.6%). To better understand KIR allelic associations and patterns in Bulgarians, we have estimated the pairwise linkage disequilibrium (LD) for the 10 KIR loci, where KIR2DL3*00501 allele was found in strong LD with KIR2DL1*00101 (D' = 1.00, R2  = 0.742). This is the first study investigating KIR polymorphism at the allele level in a population from the South-East European region. Considering the effect of the populationally shaped KIR allelic polymorphism on NK cell function, this data could lead to a better understanding of the genetic heterogeneity of this region and can be carried into clinical practice by improvement of the strategies taken for NK-mediated diseases.

Exceptional diversity of KIR and HLA class I in Egypt

Genetically determined variation of killer cell immunoglobulin like receptors (KIR) and their HLA class I ligands affects multiple aspects of human health. Their extreme diversity is generated through complex interplay of natural selection for pathogen resistance and reproductive health, combined with demographic structure and dispersal. Despite significant importance to multiple health conditions of differential effect across populations, the nature and extent of immunogenetic diversity is under-studied for many geographic regions. Here, we describe the first high-resolution analysis of KIR and HLA class I combinatorial diversity in Northern Africa. Analysis of 125 healthy unrelated individuals from Cairo in Egypt yielded 186 KIR alleles arranged in 146 distinct centromeric and 79 distinct telomeric haplotypes. The most frequent haplotypes observed were KIR-A, encoding two inhibitory receptors specific for HLA-C, two that are specific for HLA-A and -B, and no activating receptors. Together with 141 alleles of HLA class I, 75 of which encode a KIR ligand, we identified a mean of six distinct interacting pairs of inhibitory KIR and HLA allotypes per individual. We additionally characterize 16 KIR alleles newly identified in the study population. Our findings place Egyptians as one of the most highly diverse populations worldwide, with important implications for transplant matching and studies of immune-mediated diseases.

HLA-Bw4 in association with KIR3DL1 favors natural killer cell-mediated protection against severe COVID-19

Replicating SARS-CoV-2 has been shown to degrade HLA class I on target cells to evade the cytotoxic T-cell (CTL) response. HLA-I downregulation can be sensed by NK cells to unleash killer cell immunoglobulin-like receptor (KIR)-mediated self-inhibition by the cognate HLA-I ligands. Here, we investigated the impact of HLA and KIR genotypes and HLA-KIR combinations on COVID-19 outcome. We found that the peptide affinities of HLA alleles were not correlated with COVID-19 severity. The predicted poor binders for SARS-CoV-2 peptides belong to HLA-B subtypes that encode KIR ligands, including Bw4 and C1 (introduced by B*46:01), which have a small F pocket and cannot accommodate SARS-CoV-2 CTL epitopes. However, HLA-Bw4 weak binders were beneficial for COVID-19 outcome, and individuals lacking the HLA-Bw4 motif were at higher risk for serious illness from COVID-19. The presence of the HLA-Bw4 and KIR3DL1 combination had a 58.8% lower risk of developing severe COVID-19 (OR = 0.412, 95% CI = 0.187-0.904, p = 0.02). This suggests that HLA-Bw4 alleles that impair their ability to load SARS-CoV-2 peptides will become targets for NK-mediated destruction. Thus, we proposed that the synergistic responsiveness of CTLs and NK cells can efficiently control SARS-CoV-2 infection and replication, and NK-cell-mediated anti-SARS-CoV-2 immune responses being mostly involved in severe infection when the level of ORF8 is high enough to degrade HLA-I. The HLA-Bw4/KIR3DL1 genotype may be particularly important for East Asians undergoing COVID-19 who are enriched in HLA-Bw4-inhibitory KIR interactions and carry a high frequency of HLA-Bw4 alleles that bind poorly to coronavirus peptides.

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.

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.

Panobinostat enhances NK cell cytotoxicity in soft tissue sarcoma

Sarcoma is a rare and heterogeneous class of mesenchymal malignancies with poor prognosis. Panobinostat (LBH589) as one of histone deacetylase (HDAC) inhibitors, has demonstrated anti-tumor activity in patients with sarcoma, but its mechanisms remains unclear. Here, we found that LBH589 alone inhibited the proliferation and colony formation of soft tissue sarcoma(STS) cell lines. Transcriptome analysis showed that treatment with LBH589 augmented the NK cell mediated cytotoxicity. Quantitative real-time PCR and flow cytometric analysis (FACS) further confirmed that LBH589 increased the expression of NKG2D ligands MICA/MICB. Mechanistically, LBH589 activated the Wnt/β-catenin pathway by upregulating the histone acetylation in β-catenin promoter. In vitro co-culture experiments and in vivo animal experiments showed that LBH589 increased the cytotoxicity of natural killer (NK) cells while Wnt/β-catenin inhibitor decreased the effects. Our findings suggests that LBH589 facilitates the anti-tumor effect of NK cells, highlights LBH589 an effective assistance drug in NK cell-based immunotherapies.

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.

Donor with HLA-C2 is associated with acute rejection following liver transplantation in southern Chinese

Apart from presenting peptides to T cells, class I HLA molecules serve as ligands for KIRs and regulate the response of NK cells. The role played by HLA and KIR in the acute rejection (AR) following liver transplantation has been controversial. In this retrospective study, we assessed the influence of class I HLA alleles, HLA matching between donor-recipient pairs, recipient KIR and donor HLA ligands on AR following liver transplantation in southern Chinese. 143 recipients and 78 donors obtained from a single transplant center were included in the study cohort. 33 recipients with histologically confirmed AR were observed. We found that the incidence of AR did not correlate with donor or recipient class I HLA alleles and HLA matching. Neither recipient KIR gene nor the KIR genotype was associated with AR, moreover, high-resolution genotyping of 14 functional KIR genes of recipients showed that no KIR allele was independently associated with AR. However, the frequency of HLA-C2⁺ donor significantly increased in AR group compared with NAR group (52.9% vs. 24.6%, P = 0.03). In the presence of HLA-C2 by the donor allograft, AR was more frequently observed in recipients with normal expressed KIR2DS4 (43.8% vs. 15.0%, P = 0.03). Donor with HLA-C2 is therefore a major determinant of AR, which can confer risk effect in liver transplantation. Our findings can provide valuable clues for better understanding pathogenesis of AR and have important clinical implications in liver transplantation for Chinese. This article is protected by copyright. All rights reserved.

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.

COVID-19 Pandemic: Escape of Pathogenic Variants and MHC Evolution

We propose a new hypothesis that explains the maintenance and evolution of MHC polymorphism. It is based on two phenomena: the constitution of the repertoire of naive T lymphocytes and the evolution of the pathogen and its impact on the immune memory of T lymphocytes. Concerning the latter, pathogen evolution will have a different impact on reinfection depending on the MHC allomorph. If a mutation occurs in a given region, in the case of MHC allotypes, which do not recognize the peptide in this region, the mutation will have no impact on the memory repertoire. In the case where the MHC allomorph binds to the ancestral peptides and not to the mutated peptide, that individual will have a higher chance of being reinfected. This difference in fitness will lead to a variation of the allele frequency in the next generation. Data from the SARS-CoV-2 pandemic already support a significant part of this hypothesis and following up on these data may enable it to be confirmed. This hypothesis could explain why some individuals after vaccination respond less well than others to variants and leads to predict the probability of reinfection after a first infection depending upon the variant and the HLA allomorph.

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.

Remarkably low KIR and HLA diversity in Amerindians reveals signatures of strong purifying selection shaping the centromeric KIR region

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR-HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR-HLA interactions among all described worldwide populations, and that 83-97% of their KIR-HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR-HLA coevolution and its impact on human health and survival.

High-Resolution Analysis Identifies High Frequency of KIR-A Haplotypes and Inhibitory Interactions of KIR With HLA Class I in Zhejiang Han

Killer cell immunoglobulin-like receptors (KIR) interact with human leukocyte antigen (HLA) class I molecules, modulating critical NK cell functions in the maintenance of human health. Characterizing the distribution and characteristics of KIR and HLA allotype diversity across defined human populations is thus essential for understanding the multiple associations with disease, and for directing therapies. In this study of 176 Zhejiang Han individuals from Southeastern China, we describe diversity of the highly polymorphic KIR and HLA class I genes at high resolution. KIR-A haplotypes, which carry four inhibitory receptors specific for HLA-A, B or C, are known to associate with protection from infection and some cancers. We show the Chinese Southern Han from Zhejiang are characterized by a high frequency of KIR-A haplotypes and a high frequency of C1 KIR ligands. Accordingly, interactions of inhibitory KIR2DL3 with C1+HLA are more frequent in Zhejiang Han than populations outside East Asia. Zhejiang Han exhibit greater diversity of inhibitory than activating KIR, with three-domain inhibitory KIR exhibiting the greatest degree of polymorphism. As distinguished by gene copy number and allele content, 54 centromeric and 37 telomeric haplotypes were observed. We observed 6% of the population to have KIR haplotypes containing large-scale duplications or deletions that include complete genes. A unique truncated haplotype containing only KIR2DL4 in the telomeric region was also identified. An additional feature is the high frequency of HLA-B*46:01, which may have arisen due to selection pressure from infectious disease. This study will provide further insight into the role of KIR and HLA polymorphism in disease susceptibility of Zhejiang Chinese.