The combinatorial diversity of KIR and HLA class I allotypes in Peninsular Malaysia

An archaic HLA class I receptor allele diversifies natural killer cell-driven immunity in First Nations peoples of Oceania

Genetic variation in host immunity impacts the disproportionate burden of infectious diseases that can be experienced by First Nations peoples. Polymorphic human leukocyte antigen (HLA) class I and killer cell immunoglobulin-like receptors (KIRs) are key regulators of natural killer (NK) cells, which mediate early infection control. How this variation impacts their responses across populations is unclear. We show that HLA-A∗24:02 became the dominant ligand for inhibitory KIR3DL1 in First Nations peoples across Oceania, through positive natural selection. We identify KIR3DL1∗114, widespread across and unique to Oceania, as an allele lineage derived from archaic humans. KIR3DL1∗114+NK cells from First Nations Australian donors are inhibited through binding HLA-A∗24:02. The KIR3DL1∗114 lineage is defined by phenylalanine at residue 166. Structural and binding studies show phenylalanine 166 forms multiple unique contacts with HLA-peptide complexes, increasing both affinity and specificity. Accordingly, assessing immunogenetic variation and the functional implications for immunity are fundamental toward understanding population-based disease associations.

Human Leukocyte Antigen-G Gene Polymorphism in Peninsular Malaysia: A Preliminary Report

Introduction: Expression of the nonclassical human leukocyte antigen (HLA)-G gene is upregulated in placenta during pregnancy. In other cells, HLA-G is upregulated during parasitic infections and allergic reactions. Polymorphism at the HLA-G gene locus has been reported for many populations, but so far not for any ethnic groups in Malaysia. In this survey, we screened for genetic variation in HLA-G genes from representative Malay, Chinese, and Indian individuals living in Peninsular Malaysia. Materials and Methods: Blood samples were obtained with informed consent, and ethnicity classes were assigned based on self-declared pedigree information. Exons 2, 3, and 4 of the HLA-G gene were amplified by polymerase chain reaction and subjected to Sanger sequencing. Results: The most common genotype in Malays and Indians was found to be HLA-G*01:01:01:01/01:01:01:01 with frequencies of 0.206 and 0.167, respectively, whereas the HLA-G*01:01:03:01/01:01:01:01 genotype was the one most frequently observed in Chinese (0.221). Based on this study, HLA-G*01:01:01:01 (0.427-0.448) is the most frequent HLA-G allele in the all three ethnic groups. In contrast, HLA-G*01:01:02:01 (0.186) was observed as the second most frequent HLA-G allele in Malays and HLA-G*01:04:01 in Chinese and Indians, (0.188-0.198, respectively). Several minor HLA-G alleles were detected at low frequency in Malays, Chinese, or Indians (HLA-G*01:01:05, 01:01:09, 01:04:02, and 01:04:03). These have only rarely, if ever, been reported in other population groups. Subsequent statistical analysis including using principal coordinate data mapping showed the Malays, Chinese, and Indians are distinct but quite closely related to one another as compared with other population groups from across Europe and Africa. Conclusion: The HLA-G population data collected in this study showed that the ancestrally unrelated Malays, Chinese, and Indians are genetically distinct. This new database provides a foundation for further studies to capture HLA-G allelic diversity in uncharacterized populations of Malaysia and for future attempts to identify their roles in disease resistance and susceptibility.

Analysis of KIR and HLA Polymorphism in Chinese Individuals With COVID-19

Killer-cell immunoglobulin-like receptor (KIR) interactions with HLA class I have crucial roles in modulating NK cell function in response to viral infections. To explore the correlation between KIR/HLA and susceptibility to SARS-CoV-2 infection, we analysed polymorphism of KIR genes, haplotypes, HLA allotypes, and the interplay between KIR and HLA in individuals diagnosed with COVID-19. Compared to a population control group, we observed a significantly increased frequency of KIR3DL3*00802 in the COVID-19 group. When encoded by the HLA-B gene, the frequency of HLA-Bw4, a ligand for KIR3DL1, was at lower frequency in the COVID-19 group. Additionally, significantly elevated frequencies of KIR-Bx3, KIR3DL3*00301, 3DL3*048, and C1+HLA-C were identified in the COVID-19 group before multiple test correction, suggesting associations with susceptibility to SARS-CoV-2 infection. Our findings indicate that the KIR3DL3*00802 allele may be a high-risk factor for SARS-CoV-2 infection, while Bw4 encoded by HLA-B gene may confer protective effects against the infection.

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.

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.

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.

KIR genotype and haplotype frequencies in the multi-ethnic population of Malaysia

Killer cell immunoglobulin-like receptors (KIR) genotype and haplotype frequencies have been reported to vary distinctly between populations, which in turn contributes to variation in the alloreactivity of natural killer (NK) cells. Utilizing the diverse KIR genes to identify suitable transplant donors would prove challenging in multi-ethnic countries, even more in resource-limited countries where KIR genotyping has not been established. In this study, we determined the KIR genotypes from 124 unrelated Malaysians consisting of the Malays, Chinese, Indians, and aboriginal people through polymerase chain reaction sequence-specific primer (PCR-SSP) genotyping and employing an expectation-maximization (EM) algorithm to assign haplotypes based on pre-established reference haplotypes. A total of 27 distinct KIR haplotypes were discerned with higher frequencies of haplotype A (55.2%) than haplotype B (44.8%). The most frequent haplotypes were cA01:tA01 (55.2%), cB01:tB01 (18.1%), and cB02:tA01 (13.3%), while the least frequent haplotypes were cB03:tB01 (1.2%), cB04:tB03 (0.4%), and cB03:tA01 (0.4%). Several haplotypes were identified to be unique to a specific ethnic group. The genotype with the highest frequency was genotype AB (71.8%), followed by AA (19.4%), and BB (8.9%). The Indians exhibited the lowest genotype AA but the highest genotype BB, whereas genotype BB was absent in the aboriginal people. Despite the limitations, the genotype and haplotypes in the Malaysian population were successfully highlighted. The identification of ethnic-specific KIR genotypes and haplotypes provides the first step to utilizing KIR in identifying suitable transplant donors to further improve the transplant outcome in the Malaysian population.

Distribution of 22 Single Nucleotide Polymorphisms in 13 Cytokine Genes in Malays, Chinese, and Indians in Peninsular Malaysia

Introduction: Cytokines are cell signaling glycoproteins that are particularly important in immunity and inflammatory responses. Therefore, variations, such as single nucleotide polymorphisms (SNPs), in genes encoding for cytokines may have important consequences for their roles in health. Materials and Methods: A total of 222 unrelated, healthy, and un-admixed Malays (n = 97), Chinese (n = 77), and Indians (n = 48) with a median age of 30 years old (range 21-50) were typed for 22 cytokine gene SNPs: IL-1α -889 T/C, IL-1β (-511 T/C, +3962 T/C), IL-1R pst1 1970 T/C, IL-1RA mspa1 11100 T/C, IL-4Rα +1902 G/A, IL-12 - 1188 C/A, IFN-γ +874 A/T, TGF-β (cdn 10 C/T, cdn 25 G/C), TNF-α (-308 A/G, -238 A/G) IL-2 (+166 G/T, -330 T/G), IL-4 (-1098 T/G, -590 T/C, -33 T/C), IL-6 (-174 C/G, nt565 G/A), and IL-10 (-1082 G/A, -819 C/T, -592 A/C). This involved using well-established polymerase chain reaction procedures with sequence-specific primers and restriction fragment length polymorphism methods. Results: The majority of the screened cytokine gene SNPs are polymorphic in all three ethnicities. Exceptions include TGF-β cdn 25 (G/C), IL-1β +3962 (T/C), and TNF-α -238 (A/G), which were all observed to be monomorphic in Malays, Chinese and Indians. Many of the analyzed cytokine gene SNP genotypes deviated from Hardy-Weinberg equilibrium and the three ethnic study groups were all well-separated from reference Asian, African and European populations in a principal component analysis plot. Conclusion: We successfully typed 22 SNPs in 13 cytokine genes from genetic material collected from unrelated and un-admixed Malay, Chinese and Indian individuals in Peninsular Malaysia. These new cytokine gene population datasets reveal interesting contrasts with other populations.

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.

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 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.

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.

Pregnancy induced TMA in severe preeclampsia results from complement-mediated thromboinflammation

Preeclampsia is a multifactorial vascular disease unique to human pregnancy. While genetic and antiangiogenic factors are important contributors to preeclampsia susceptibility, recent studies have shown that dysregulation and/or over-activation of the complement system has an integral role in disease etiology. Furthermore, the role of the coagulation cascade may be underappreciated in the development of the disease. Traditionally, for research purposes, the pool of preeclampsia cases has been divided into non-severe and severe disease depending on the onset and severity of the symptoms. However, of particular interest are a small but important minority of cases that present with symptoms likening to those of hemolysis, elevated liver enzymes and low platelets syndrome, atypical hemolytic uremic syndrome, or thrombotic thrombocytopenic purpura, all thrombotic microangiopathy (TMA) diseases, with the hallmark mechanisms of endothelial dysfunction and aberrant activation of complement and coagulation cascades. We therefore propose a third class, severe TMA-like preeclampsia to be included in the categorization of preeclampsia patients. Identifying these patients would target research, diagnostic differentiation, and novel treatment options to the subclass of patients with life-threatening disease that are most likely to benefit from next-generation drug development.