A deep look at KIR–HLA in Amerindians: Comprehensive meta-analysis reveals limited diversity of KIR haplotypes

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

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.

Amerindian genetic ancestry as a risk factor for tuberculosis in an amazonian population

In recent years, the incidence of tuberculosis (TB) has declined worldwide, although this disease still occurs at relatively high rates in Amerindian populations. This suggests that the genetic ancestry of Amerindians may be an important factor in the development of infections, and may account for at least some of the variation in infection rates in the different populations. The present study investigated the potential influence of Amerindian genetic ancestry on susceptibility to tuberculosis in an Amazon population. The study included 280 patients diagnosed with tuberculosis and 138 asymptomatic hospital employees with no history of TB, but who were in contact with bacterially active TB patients. Ancestry analysis was run on a set of 61 Ancestry-Informative Markers to estimate European, African, and Amerindian genetic ancestry using STRUCTURE v2.2. The TB group had significantly higher Amerindian ancestry in comparison with the control group, and significantly lower European ancestry. Amerindian ancestry in the 20–60% range was found to be the principal risk factor for increased susceptibility to TB. The results of the study indicate that Amerindian ancestry is an important risk factor for susceptibility to TB in the admixed population of the Brazilian Amazon region.

Fluctuating and Geographically Specific Selection Characterize Rapid Evolution of the Human KIR Region

The killer-cell immunoglobulin-like receptor (KIR) region comprises a fast-evolving family of genes that encode receptors for natural killer (NK) cells and have crucial role in host defense. Evolution of KIR was examined in the context of the human genome. Gene-content diversity and single nucleotide polymorphisms (SNP) in the KIR genes and flanking regions were compared to >660,000 genome-wide SNPs in over 800 individuals from 52 populations of the human genome diversity panel (HGDP). KIR allelic diversity was further examined using next generation sequencing in a subset of 56 individuals. We identified the SNP rs587560 located in KIR3DL3 as a marker of KIR2DL2 and KIR2DL3 and, consequently, Cen A and Cen B haplotypes. We also show that combinations of two KIR2DL4 SNPs (rs35656676 and rs592645) distinguish KIR3DL1 from KIR3DS1 and also define the major KIR3DL1 high- and low-expressing alleles lineages. Comparing the diversity of the SNPs within the KIR region to remainder of the genome, we observed a high diversity for the centromeric KIR region consistent with balancing selection (p < 0.01); in contrast, centromeric KIR diversity is significantly reduced in East Asian populations (p < 0.01), indicating purifying selection. By analyzing SNP haplotypes in a region spanning ~500 kb that includes the KIR cluster, we observed evidence of strong positive selection in Africa for high-expressing KIR3DL1 alleles, favored over the low-expressing alleles (p < 0.01). In sharp contrast, the strong positive selection (p < 0.01) that we also observed in the telomeric KIR region in Oceanic populations tracked with a high frequency of KIR3DS1. In addition, we demonstrated that worldwide frequency of high-expression KIR3DL1 alleles was correlated with virus with virus (r = 0.64, p < 10⁻⁶) and protozoa (r = 0.69, p < 10⁻⁶) loads, which points to selection globally on KIR3DL1 high-expressing alleles attributable to pathogen exposure.

KIR-HLA distribution in a Vietnamese population from Hanoi

The KIR (killer cell immunoglobulin-like receptors) gene family codifies a group of receptors that recognize human leukocyte antigens (HLA) and modulate natural killer (NK) cells response. Genetic diversity of KIR genes and HLA ligands has not yet been deeply investigated in South East Asia. Here, we characterized KIR gene presence and absence polymorphism of 14 KIR genes and two pseudogenes, as well as the frequencies of the ligands HLA-Bw4, HLA-C1 and HLA-C2 in a Vietnamese population from Hanoi (n = 140). Genotyping was performed by polymerase chain reaction with specific sequence primers (PCR-SSP). We compared KIR frequencies and performed principal component analysis with 43 worldwide populations of different ancestries. KIR carrier frequencies in Vietnamese were similar to those reported for Thai and Chinese Han, but differed significantly from other geographically close populations such as Japanese and South Korean. This similarity was also observed in KIR gene-content genotypes and is in accordance with the origin from Southern China and Thailand proposed for the Vietnamese population. The frequencies of HLA ligands observed in Vietnamese did not differ from those reported for other East-Asian populations (p > .05). Studies regarding KIR-HLA in populations are of prime importance to understand their evolution, function and role in diseases.

Distribution of KIR genes and KIR2DS4 gene variants in two Mexican Mestizo populations

Killer immunoglobulin-like receptors (KIR) are transmembrane proteins that regulate NK and T cell subsets by recognizing HLA-I molecules as ligands. The KIR gene family consists of 16 genes, located at chromosome 19q13.4. KIR gene frequencies vary among populations. In Mexico, HLA and genetic ancestry studies show that Mestizo populations have different genetic backgrounds based on admixture with European, African, and Asian ancestry. This study aimed to evaluate the frequencies of KIR genes and genotypes in Guerrero and Jalisco, two Mexican Mestizo populations located in the south and the west of the country, respectively, and to compare these frequencies with those of other populations. KIR genotyping was performed by SSP-PCR. We observed that KIR gene frequencies were similar in both populations. There were 24 genotypes observed in Guerrero, 38 genotypes observed in Jalisco, 15 genotypes shared in both populations and 32 genotypes unique to one population or the other. In 10 individuals, nine novel genotypes were identified. KIR2DS4 gene variants showed significant differences: The KIR2DS4full gene was more common in Guerrero (p<0.0001), and the KIR2DS4del variant was more common in Jalisco (p<0.05). Differences in KIR2DS4 gene variants and genotypic profiles could be influenced by the genetic admixture in both regions.

Colorectal cancer and the KIR genes in the human genome: A meta-analysis

Colorectal cancer is one of the most common types of inflammation-based cancers and is occurred due to growth and spread of cancer cells in colon and/or rectum. Previously genetic association of cell cycle genes, both proto-oncogenes and the tumor suppressors has been proved. But there were few studies about association of immune related genes such as killer-cell immunoglobulin-like receptors (KIRs). Thus we intend to perform a meta-analysis to find the association of different genes of KIR and susceptibility to be affected by colorectal cancer. The overall population of the four studies investigated in our meta-analysis was 953 individuals (470 individuals with colorectal cancer and 483 individuals in control groups). After the analyses, we concluded that colorectal cancer is affected by KIR2DS5 and also there were no protecting gene. This result shows the inflammatory basis of this cancer. In other words, in contrast to leukemia and blood cancers, colorectal cancers seem to be affected by hyper activity of natural killer-cells (NKs). Whys and therefore of this paradox, is suggested to be investigated further.

Estimating KIR Haplotype Frequencies on a Cohort of 10,000 Individuals: A Comprehensive Study on Population Variations, Typing Resolutions, and Reference Haplotypes

The killer cell immunoglobulin-like receptors (KIR) mediate human natural killer (NK) cell cytotoxicity via activating or inhibiting signals. Although informative and functional haplotype patterns have been reported, most genotyping has been performed at resolutions that are structurally ambiguous. In order to leverage structural information given low-resolution genotypes, we performed experiments to quantify the effects of population variations, reference haplotypes, and genotyping resolutions on population-level haplotype frequency estimations as well as predictions of individual haplotypes. We genotyped 10,157 unrelated individuals in 5 populations (518 African American[AFA], 258 Asian or Pacific Islander[API], 8,245 European[EUR], 1,073 Hispanic[HIS], and 63 Native American[NAM]) for KIR gene presence/absence (PA), and additionally half of the AFA samples for KIR gene copy number variation (CNV). A custom EM algorithm was used to estimate haplotype frequencies for each population by interpretation in the context of three sets of reference haplotypes. The algorithm also assigns each individual the haplotype pairs of maximum likelihood. Generally, our haplotype frequency estimates agree with similar previous publications to within <5% difference for all haplotypes. The exception is that estimates for NAM from the U.S. showed higher frequency association of cB02 with tA01 (+14%) instead of tB01 (-8.5%) compared to a previous study of NAM from south of the U.S. The higher-resolution CNV genotyping on the AFA samples allowed unambiguous haplotype-pair assignments for the majority of individuals, resulting in a 22% higher median typing resolution score (TRS), which measures likelihood of self-match in the context of population-specific haplo- and geno-types. The use of TRS to quantify reduced ambiguity with CNV data clearly revealed the few individuals with ambiguous genotypes as outliers. It is observed that typing resolution and reference haplotype set influence haplotype frequency estimates. For example, PA resolution may be used with reference haplotype sets up to the point where certain haplotypes are gene-content subsets of others. At that point, CNV must be used for all genes.

The Impact of KIR Polymorphism on the Risk of Developing Cancer: Not as Strong as Imagined?

The polymorphism of killer cell immunoglobulin-like receptors (KIR) has been associated with several diseases, including infection, autoimmunity and cancer. KIR molecules are a family of receptors expressed on the surface of natural killer cells (NK), frontline defense of innate immunity against microorganisms and neoplastic cells. Some studies have shown conflicting results concerning the role that KIR polymorphism plays in tumor susceptibility, particularly in leukemia and lymphoma. Interestingly, the presence of HLA ligands is sometimes strongly associated with several types of cancer and apparently is not related with their interaction with KIR. This manuscript briefly reviews the uncommon polymorphism of KIR and critically summarizes the recent findings with regards of the importance of KIR variation for cancer susceptibility.

KIR and HLA genotyping of Japanese descendants from Curitiba, a city of predominantly European ancestry from Southern Brazil

We hereby report the KIR gene frequencies and the frequencies of HLA ligands of KIR for Brazilians of Japanese ancestry. A total of 51 individuals were genotyped for presence/absence of KIR (killer cell immunoglobulin-like receptors) genes and presence of HLA (human leukocyte antigens) ligands. KIR was genotyped using two pairs of sequence-specific primers (PCR-SSP) and HLA ligands were typed by LABType® SSO reagent kits (One Lambda, USA). These data are fully available in Allele Frequencies Net Database, under the population name "Brazil Curitiba Japanese KIR".

Differential gene expression levels might explain association of LAIR2 polymorphisms with pemphigus

The leukocyte-associated immunoglobulin-like receptor 1 (LAIR-1) is a collagen-binding inhibitory receptor important for the regulation of immune responses, expressed on the majority of peripheral blood mononuclear cells (PBMC). LAIR-2 is a soluble homolog that antagonizes LAIR-1 inhibitory function by binding the same ligands. We sought to investigate whether LAIR1 and LAIR2 single nucleotide polymorphisms (SNP) are associated with differential mRNA expression levels. We analyzed 14 SNPs of LAIR1 (6) and LAIR2 (8) by mass spectrometry-based genotyping and extracted mRNA from PBMC of 177 healthy subjects, followed by quantitative assays. Four SNPs of LAIR1 and two SNPs of LAIR2 mark differential mRNA levels in healthy individuals. To verify the biological relevance of these findings, we analyzed additional 515 individuals (282 patients and 233 controls) to check if LAIR1 and LAIR2 differential mRNA expression could be related to susceptibility to pemphigus foliaceus (PF), an autoimmune blistering skin disease endemic in Brazil. Two LAIR1 variants (rs56802430 G, OR = 1.52, p = 0.0329; rs11084332 C, OR = 0.57, p = 0.0022) and one LAIR2 (rs2287828 T+, OR = 1.9, p = 0.0097) contribute to differential susceptibility to PF. Furthermore, we demonstrate interactions among four LAIR2 SNPs (rs2042287, rs2287828, rs2277974 and rs114834145). A haplotype harboring these SNPs is strongly associated with higher LAIR2 mRNA levels (4.5-fold, p = 0.0069) and with higher susceptibility to PF (OR = 4.02, p = 0.0008). We suggest that LAIR1 and LAIR2 genetic variants are associated with regulation of gene expression and variable PF susceptibility, and show indirect association of LAIR2 differential mRNA expression with PF pathogenesis. Our data demonstrate how this relatively unknown disease can add invaluable knowledge regarding the role of LAIR1 and LAIR2 in immune responses.

Variability of innate immune system genes in Native American populations-relationship with history and epidemiology

OBJECTIVES

The immune system of a host, defending him/her against invading pathogens, has two main subsystems: innate immunity and acquired immunity. There are several evidences showing that Native American populations are immunologically different from non-Native populations. Our aim was to describe the variability of innate immune system genes in Native American populations.

MATERIALS AND METHODS

We investigated heterozygozities and patterns of population differentiation (FST ) of 14 polymorphisms related to the innate immune response in five Native American populations (Aché, Guarani-Kaiowá, Guarani-Ñandeva, Kaingang, and Xavante) and the results were compared with the three major world population data (YRI, CEU, and CHB) available at the 1,000 genomes database.

RESULTS

Mean heterozygosities ranged between 0.241 ± 0.057 (Aché) and 0.343 ± 0.033 (Kaingang), but no significant differences were observed (Friedman test, P = 0.197). Mean heterozygosities were also not significantly different when Amerindians were pooled and compared with the 1000 genomes populations (Friedman test, P = 0.506). When the Native American populations were grouped as Amerindians, a significantly higher FST value (0.194) was observed between the Amerindian and African populations. The Ewens-Watterson neutrality test showed that these markers are not under strong selective pressure.

DISCUSSION

Native American populations present similar levels of heterozygosity as those of other continents, but are different from Africans in the frequency of polymorphisms of innate immune genes. This higher differentiation is probably due to demographic processes that occurred during the out-of-Africa event.

Pemphigus is associated with KIR3DL2 expression levels and provides evidence that KIR3DL2 may bind HLA-A3 and A11 in vivo

Although HLA-A3 and A11 have been reported to be ligands for KIR3DL2, evidence for any in vivo relevance of this interaction is still missing. To explore the functional importance of KIR3DL2 allelic variation, we analyzed the autoimmune disease pemphigus foliaceus, previously associated (lower risk) with activating KIR genes. KIR3DL2*001 was increased in patients (odds ratio (OR) = 2.04; p = 0.007). The risk was higher for the presence of both KIR3DL2*001 and HLA-A3 or A11 (OR = 3.76, p = 0.013), providing the first evidence that HLA-A3 and A11 may interact with KIR3DL2 in vivo. The nonsynonymous single nucleotide polymorphism 1190T (rs3745902) was associated with protection (OR = 0.52, p = 0.018). This SNP results in a threonine-to-methionine substitution. Individuals who have methionine in this position exhibit a lower percentage of KIR3DL2-positive natural killer (NK) cells and also lower intensity of KIR3DL2 on expressing natural killer cells; additionally, we show that the expression of KIR3DL2 is independent of other killer cell immunoglobulin-like receptors. Pemphigus foliaceus is a very unique complex disease strongly associated with immune-related genes. It is the only autoimmune disease known to be endemic, showing a strong correlation with environmental factors. Our data demonstrate that this relatively unknown autoimmune disease may facilitate understanding of the molecular mechanisms of KIR3DL2 ligand recognition.

KIR and HLA under pressure: evidences of coevolution across worldwide populations

KIR (killer cell immunoglobulin-like receptors) and HLA (human leukocyte antigens) are two distinct gene families with remarkable importance for human immune responses. The recognition of HLA molecules by activating and inhibitory KIR promotes a balance of signals that regulates NK cell function and is especially important for the innate defense against pathogens and early placentation. There is no documented gametic association between these two gene families and no evidence of common regulation. However, due to the critical role of KIR recognition for immunity and reproduction, the possibility of KIR-HLA combinations being under selective pressure is not surprising. In this manuscript, we first summarize the HLA-KIR system, the HLA molecules that are the putative ligands for KIR, and then we review the evidences that suggest these two gene families are coevolving as an integrated system.