KIR gene content in amerindians indicates influence of demographic factors

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.

Killer-Cell Immunoglobulin-like Receptor Diversity in an Admixed South American Population

Natural Killer (NK) cells are innate immune cells that mediate antiviral and antitumor responses. NK cell activation and induction of effector functions are tightly regulated by the integration of activating and inhibitory receptors such as killer immunoglobulin-like receptors (KIR). KIR genes are characterized by a high degree of diversity due to presence or absence, gene copy number and allelic polymorphism. The aim of this study was to establish the distribution of KIR genes and genotypes, to infer the most common haplotypes in an admixed Colombian population and to compare these KIR gene frequencies with some Central and South American populations and worldwide. A total of 161 individuals from Medellin, Colombia were included in the study. Genomic DNA was used for KIR and HLA genotyping. We analyzed only KIR gene-content (presence or absence) based on PCR-SSO. The KIR genotype, most common haplotypes and combinations of KIR and HLA ligands frequencies were estimated according to the presence or absence of KIR and HLA genes. Dendrograms, principal component (PC) analysis and Heatmap analysis based on genetic distance were constructed to compare KIR gene frequencies among Central and South American, worldwide and Amerindian populations. The 16 KIR genes analyzed were distributed in 37 different genotypes and the 7 most frequent KIR inferred haplotypes. Importantly, we found three new genotypes not previously reported in any other ethnic group. Our genetic distance, PC and Heatmap analysis revealed marked differences in the distribution of KIR gene frequencies in the Medellin population compared to worldwide populations. These differences occurred mainly in the activating KIR isoforms, which are more frequent in our population, particularly KIR3DS1. Finally, we observed unique structural patterns of genotypes, which evidences the potential diversity and variability of this gene family in our population, and the need for exhaustive genetic studies to expand our understanding of the KIR gene complex in Colombian populations.

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.

Understanding the genetic history of Malay populations in Peninsular Malaysia via KIR genes diversity

OBJECTIVES

Killer cell immunoglobulin-like receptor (KIR) genes with high polymorphism at genotypic levels are important in providing immune defense and have been expanded towards human population genetics. The aim of this study is to provide supporting information from this new biomarker to strengthen the comprehension of genetic history of the complex Malay population.

METHODS

KIR genotyping for 213 unadmixed Malay individuals from six subethnic groups (Acheh, Bugis, Champa, Mandailing, Minang and Kedah) was carried out using PCR-SSP (sequence specific primers) method in 16 independent reactions.

RESULTS

The most frequent KIR genotype observed is AA1, followed by AB4 and AB5. Five genotypes; AA1, AB4, AB5, AB7 and AB8 were shared among all Malay subethnic groups. The highest frequency of KIR haplotype A was observed in Minang Malays, whereas Acheh and Kedah Malays carry a balanced distribution of A and B KIR haplotypes. PCA for the KIR genes clearly illustrated six ethnogeographical population clusters; Africans, Amerindian, Northeast Asian, South Asian, Oceania and Southeast Asian populations. All six Malay subethnic groups fell within the Southeast Asian cluster.

CONCLUSIONS

The complex array of KIR genotypes observed in the Malays indicates their historical interactions with various populations, especially with the Chinese, Indians and Orang Asli. This study has demonstrated the potential of KIR genes as a genetic marker for deducing population structure and genetic relationship between populations.

Single Nucleotide Polymorphism in KIR2DL1 Is Associated With HLA-C Expression in Global Populations

Regulation of NK cell activity is mediated through killer-cell immunoglobulin-like receptors (KIR) ability to recognize human leukocyte antigen (HLA) class I molecules as ligands. Interaction of KIR and HLA is implicated in viral infections, autoimmunity, and reproduction and there is growing evidence of the coevolution of these two independently segregating gene families. By leveraging KIR and HLA-C data from 1000 Genomes consortium we observed that the KIR2DL1 variant rs2304224 ^* T is associated with lower expression of HLA-C in individuals carrying the ligand HLA-C2 (p = 0.0059). Using flow cytometry, we demonstrated that this variant is also associated with higher expression of KIR2DL1 on the NK cell surface (p = 0.0002). Next, we applied next generation sequencing to analyze KIR2DL1 sequence variation in 109 Euro and 75 Japanese descendants. Analyzing the extended haplotype homozygosity, we show signals of positive selection for rs4806553 ^* G and rs687000 ^* G, which are in linkage disequilibrium with rs2304224 ^* T. Our results suggest that lower expression of HLA-C2 ligands might be compensated for higher expression of the receptor KIR2DL1 and bring new insights into the coevolution of KIR and HLA.

Unveiling the Diversity of Immunoglobulin Heavy Constant Gamma (IGHG) Gene Segments in Brazilian Populations Reveals 28 Novel Alleles and Evidence of Gene Conversion and Natural Selection

Even though immunoglobulins are critical for immune responses and human survival, the diversity of the immunoglobulin heavy chain gene (IGH) is poorly known and mostly characterized only by serological methods. Moreover, this genomic region is not well-covered in genomic databases and genome-wide association studies due to particularities that impose technical difficulties for its analysis. Therefore, the IGH gene has never been systematically sequenced across populations. Here, we deliver an unprecedented and comprehensive characterization of the diversity of the IGHG1, IGHG2, and IGHG3 gene segments, which encode the constant region of the most abundant circulating immunoglobulins: IgG1, IgG2, and IgG3, respectively. We used Sanger sequencing to analyze 357 individuals from seven different Brazilian populations, including five Amerindian, one Japanese-descendant and one Euro-descendant population samples. We discovered 28 novel IGHG alleles and provided evidence that some of them may have been originated by gene conversion between common alleles of different gene segments. The rate of synonymous substitutions was significantly higher than the rate of the non-synonymous substitutions for IGHG1 and IGHG2 (p = 0.01 and 0.03, respectively), consistent with purifying selection. Fay and Wu's test showed significant negative values for most populations (p < 0.001), which indicates that positive selection in an adjacent position may be shaping IGHG variation by hitchhiking of variants in the vicinity, possibly the regions that encode the Ig variable regions. This study shows that the variation in the IGH gene is largely underestimated. Therefore, exploring its nucleotide diversity in populations may provide valuable information for comprehension of its evolution, its impact on diseases and vaccine research.

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.

Cost-effective and fast KIR gene-content genotyping by multiplex melting curve analysis

Killer cell immunoglobulin-like receptor (KIR) genes encode cell surface molecules that recognize HLA molecules and modulate the activity of natural killer (NK) cells. KIR genes exhibit presence and absence polymorphism, which generates a variety of gene-content haplotypes in worldwide populations. KIR gene-content variation is implicated in many diseases and is also important for placentation and transplantation. Because of the complexity of KIR polymorphism, variation in this family is still mostly studied at the gene-content level, even with the advent of next-generation sequencing (NGS) methods. Gene-content determination is generally expensive and/or time-consuming. To overcome these difficulties, we developed a method based on multiplex polymerase chain reaction with specific sequence primers (PCR-SSP) followed by melting curve analysis that allows cost-effective, precise and fast generation of results. Our method was 100% concordant with a gel-based method and 99.9% concordant with presence and absence determination by NGS. The limit of detection for accurate typing was 30 ng of DNA (0.42 μM) with 260/230 and 260/280 ratios as low as 0.19 and of 0.44. In addition, we developed a user-friendly Java-based computational application called killerPeak that interprets the raw data generated by Viia7 or QuantStudio 7 quantitative PCR machines and reliably exports the final genotyping results in spreadsheet file format. The combination of a reliable method that requires low amount of DNA with an automated interpretation of results allows scaling the KIR genotyping in large cohorts with reduced turnaround time.

Report from the Killer-cell Immunoglobulin-like Receptors (KIR) component of the 17th International HLA and Immunogenetics Workshop

The goals of the KIR component of the 17th International HLA and Immunogenetics Workshop (IHIW) were to encourage and educate researchers to begin analyzing KIR at allelic resolution, and to survey the nature and extent of KIR allelic diversity across human populations. To represent worldwide diversity, we analyzed 1269 individuals from ten populations, focusing on the most polymorphic KIR genes, which express receptors having three immunoglobulin (Ig)-like domains (KIR3DL1/S1, KIR3DL2 and KIR3DL3). We identified 13 novel alleles of KIR3DL1/S1, 13 of KIR3DL2 and 18 of KIR3DL3. Previously identified alleles, corresponding to 33 alleles of KIR3DL1/S1, 38 of KIR3DL2, and 43 of KIR3DL3, represented over 90% of the observed allele frequencies for these genes. In total we observed 37 KIR3DL1/S1 allotypes, 40 for KIR3DL2 and 44 for KIR3DL3. As KIR allotype diversity can affect NK cell function, this demonstrates potential for high functional diversity worldwide. Allelic variation further diversifies KIR haplotypes. We determined KIR3DL3 ∼ KIR3DL1/S1 ∼ KIR3DL2 haplotypes from five of the studied populations, and observed multiple population-specific haplotypes in each. This included 234 distinct haplotypes in European Americans, 191 in Ugandans, 35 in Papuans, 95 in Egyptians and 86 in Spanish populations. For another 35 populations, encompassing 642,105 individuals we focused on KIR3DL2 and identified another 375 novel alleles, with approximately half of them observed in more than one individual. The KIR allelic level data gathered from this project represents the most comprehensive summary of global KIR allelic diversity to date, and continued analysis will improve understanding of KIR allelic polymorphism in global populations. Further, the wealth of new data gathered in the course of this workshop component highlights the value of collaborative, community-based efforts in immunogenetics research, exemplified by the IHIW.

Southern Bahia, Brazil: KIR gene-content distribution in the highly admixed population from Ilhéus

We here describe the KIR gene-content genotypes in the highly admixed population from Ilhéus, a city located in Bahia State, Northeast Brazil. Bahia is the State with the largest contribution of African ethnicity in comparison to other Brazilian areas and this is the first report of a population from this region. A total of 32 KIR gene-content genotypes have been found, from which 18 were observed in single individuals. The distinct KIR gene-content distribution observed in Ilhéus, along with the elevated diversity of genotypes point to importance of describing KIR polymorphism in unique populations from Brazil.

How strong was the bottleneck associated to the peopling of the Americas? New insights from multilocus sequence data

In spite of many genetic studies that contributed for a deep knowledge about the peopling of the Americas, no consensus has emerged about important parameters such as the effective size of the Native Americans founder population. Previous estimates based on genomic datasets may have been biased by the use of admixed individuals from Latino populations, while other recent studies using samples from Native American individuals relied on approximated analytical approaches. In this study we use resequencing data for nine independent regions in a set of Native American and Siberian individuals and a full-likelihood approach based on isolation-with-migration scenarios accounting for recent flow between Asian and Native American populations. Our results suggest that, in agreement with previous studies, the effective size of the Native American population was small, most likely in the order of a few hundred individuals, with point estimates close to 250 individuals, even though credible intervals include a number as large as ~4,000 individuals. Recognizing the size of the genetic bottleneck during the peopling of the Americas is important for determining the extent of genetic markers needed to characterize Native American populations in genome-wide studies and to evaluate the adaptive potential of genetic variants in this population.

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.

Association between KIR genes and dust mite sensitization in a Brazilian population

BACKGROUND

Killer cell immunoglobulin-like receptors (KIRs), found on the surface of natural killer (NK) cells, play a key role in controlling the innate response. Such response depends on a series of cellular interactions between these receptors and HLA activating/inhibiting ligands. Atopic diseases have been associated with genes that regulate cytokine production and HLA genes, which may either protect or predispose to hypersensitivity.

OBJECTIVE

To verify an association study of KIR genes with sensitization to the following mites: Dermatophagoides farinae, Dermatophagoides pteronyssinus, and Blomia tropicalis.

METHODS

A total of 341 children aged up to 14 years, were classified as mite-sensitive or mite-insensitive after undergoing a skin prick test for immediate allergic reactions. The presence/absence of KIR genes and their human leukocyte antigen (HLA) ligands was determined by polymerase chain reaction-sequence specific oligonucleotide (PCR-SSO) with the commercial kit LabType™ using Luminex™.

RESULTS

The frequencies of KIR genes and their respective class I HLA ligands and the frequency of haplotypes were performed in sensitive and insensitive individuals, and no significant differences were found.

CONCLUSION

Our results suggest no influence of KIR genes on resistance/susceptibility to sensitization to dust mites.

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.

KIR genes and HLA class I ligands in a Caucasian Brazilian population with colorectal cancer

Colorectal cancer (CRC) can occur anywhere in the colon or rectum and represents the third most common cancer in the world in both sexes. Natural killer cells (NK) are part of the innate immune system recognizing class I HLA molecules on target cells through their membrane receptors, called killer cell immunoglobulin-like receptors (KIR). The aim of our study was to evaluate the association between the KIR genes and HLA ligands in patients with colorectal cancer and healthy controls. We examined the polymorphism of 16 KIR genes and their HLA ligands in 154 caucasoid CRC patients and 216 controls. When both groups were compared, no significant differences were found for HLA ligands and KIR genes after Bonferroni correction. However, the Bx haplotypes (heterozygous and homozygous for the haplotype B) were more frequent in controls, when compared with patients. These findings suggest that individuals with Bx haplotypes could have some protection to colorectal cancer. The hypothesis is not related with the presence of a special KIR gene and HLA ligand related to the disease, but to the presence of several activating genes in the individuals with no better action of one in relation to other. Further studies to confirm this observation are warranted.

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.

The role of variants from the innate immune system genes in tuberculosis and skin test response in a Native American population

Native American populations show higher tuberculosis (TB) mortality and infectivity rates than non-Native populations. Variants in the innate immune system seem to have an important role on TB susceptibility. The role of some innate immune system variants in TB susceptibility and/or skin test response (PPD) were investigated in the Aché, a Native American population. Complement receptor 1 and toll like receptor 9 variants were associated with anergy to PPD and protection to TB, respectively. These findings demonstrate an important role of the innate immune system variants in TB 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.

A General Model of Negative Frequency Dependent Selection Explains Global Patterns of Human ABO Polymorphism

The ABO locus in humans is characterized by elevated heterozygosity and very similar allele frequencies among populations scattered across the globe. Using knowledge of ABO protein function, we generated a simple model of asymmetric negative frequency dependent selection and genetic drift to explain the maintenance of ABO polymorphism and its loss in human populations. In our models, regardless of the strength of selection, models with large effective population sizes result in ABO allele frequencies that closely match those observed in most continental populations. Populations must be moderately small to fall out of equilibrium and lose either the A or B allele (N(e) ≤ 50) and much smaller (N(e) ≤ 25) for the complete loss of diversity, which nearly always involved the fixation of the O allele. A pattern of low heterozygosity at the ABO locus where loss of polymorphism occurs in our model is consistent with small populations, such as Native American populations. This study provides a general evolutionary model to explain the observed global patterns of polymorphism at the ABO locus and the pattern of allele loss in small populations. Moreover, these results inform the range of population sizes associated with the recent human colonization of the Americas.

Differing evolutionary histories of the ACTN3*R577X polymorphism among the major human geographic groups

It has been proposed that the functional ACTN3*R577X polymorphism might have evolved due to selection in Eurasian human populations. To test this possibility we surveyed all available population-based data for this polymorphism and performed a comprehensive evolutionary analysis of its genetic diversity, in order to assess the action of adaptive and random mechanisms on its variation across human geographical distribution. The derived 577X allele increases in frequency with distance from Africa, reaching the highest frequencies on the American continent. Positive selection, detected by an extended haplotype homozygosisty test, was consistent only with the Eurasian data, but simulations with neutral models could not fully explain the results found in the American continent. It is possible that particularities of Native American population structure could be responsible for the observed allele frequencies, which would have resulted from a complex interaction between selective and random factors.

Killer cell immunoglobulin-like receptor (KIR) gene diversity in a population naturally exposed to malaria in Porto Velho, Northern Brazil

Killer cell immunoglobulin-like receptors (KIR) are expressed mainly in natural killer cells and specifically recognize human leukocyte antigen (HLA) class I molecules. The repertoire of KIR genes and KIR-HLA pairs is known to play a key role in the susceptibilities to and the outcomes of several diseases, including malaria. The aim of this study was to investigate the distribution of KIR genes, KIR genotypes and KIR-HLA pair combinations in a population naturally exposed to malaria from Brazilian Amazon. All 16 KIR genes investigated were present in the studied population. Overall, 46 KIR genotypes were defined. The two most common genotypes in the Porto Velho communities, genotypes 1 and 2, were present at similar frequencies as in the Americas. Principal component analysis based on the frequencies of the KIR genes placed the Porto Velho population closer to the Venezuela Mestizos, USA California hispanic and Brazil Paraná Mixed in terms of KIR gene frequencies. This analysis highlights the multi-ethnic profile of the Porto Velho population. Most of the individuals were found to have at least one inhibitory KIR-HLA pair. Seventy-five KIR-HLA pair combinations were identified. The KIR-2DL2/3_HLA-C1, KIR3DL1_HLA-Bw4 and KIR2DL1_HLA-C2 pairs were the most common. There was no association between KIR genes, KIR genotypes or KIR-HLA pair combinations and malaria susceptibility in the studied population. This is the first report on the distribution of KIR and known HLA ligands in the Porto Velho population. Taken together, these results should provide baseline information that will be relevant to population evolutionary history, malaria and other diseases studies in populations of the Brazilian Amazon.