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.

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.

High KIR diversity in Amerindians is maintained using few gene-content haplotypes.

Interaction between killer cell immunoglobulin-like receptors (KIR) and cognate HLA class I ligands influences the innate and adaptive immune response to infection. The KIR family varies in gene content and allelic polymorphism, thereby, distinguishing individuals and populations. KIR gene content was determined for 230 individuals from three Amerindian tribes from Venezuela: the Yucpa, Bari and Warao. Gene-content haplotypes could be assigned to 212 individuals (92%) because only five different haplotypes were present-group A and four group B. Six different haplotype combinations accounted for >80% of individuals. Each tribe has distinctive genotype frequencies. Despite few haplotypes, all 14 KIR genes are at high frequency in the three tribes, with the exception of 2DS3. Each population has an even frequency of group A and B haplotypes. Allele-level analysis of 3DL1/S1 distinguished five group A haplotypes and six group B haplotypes. The high frequency and divergence of the KIR haplotypes in the Amerindian tribes provide greater KIR diversity than is present in many larger populations. An extreme case being the Yucpa, for whom two gene-content haplotypes account for >90% of the population. These comprise the group A haplotype and a group B haplotype containing all the KIR genes, except 2DS3, that typify the group B haplotypes. Here is clear evidence for balancing selection on the KIR system and the biological importance of both A and B haplotypes for the survival of human populations.

SNP haplotypes and allele frequencies show evidence for disruptive and balancing selection in the human leukocyte receptor complex.

The human leukocyte receptor complex (LRC) of Chromosome 19q13.4 encodes polymorphic and highly homologous genes that are expressed by cells of the immune system and regulate their function. There is an enormous diversity at the LRC, most particularly the variable number of killer cell immunoglobulin-like receptor (KIR) genes. KIR have been associated with several disease processes due to their interaction with polymorphic human leukocyte antigen class I molecules. We have assessed haplotype compositions, linkage disequilibrium patterns and allele frequencies in two Caucasoid population samples (n=54, n=100), using a composite of single-nucleotide polymorphism (SNP) markers and high-resolution, allele-specific molecular genotyping. Particular KIR loci segregated with SNP and other markers, forming two blocks that were separated by a region with a greater history of recombination. The KIR haplotype composition and allele frequency distributions were consistent with KIR having been subject to balancing selection (Watterson's F: P=0.001). In contrast, there was a high inter - population heterogeneity measure for the LRC-encoded leukocyte immunoglobulin-like receptor A3 (LILRA3), indicating pathogen-driven disruptive selection (Wright's FST=0.32). An assessment of seven populations representative of African, Asian and Caucasoid ethnic groups (total n=593) provided little evidence for long-range LRC haplotypes. The different natural selection pressures acting on each locus may have contributed to a lack of linkage disequilibrium between them.

A comparison of HLA-DR and -DQ allele and haplotype frequencies in Trinidadian populations of African, South Asian, and mixed ancestry.

Using polymerase chain reaction-sequence-specific primer (PCR-SSP) typing, this study determined the frequency of human leukocyte antigen (HLA) DR- and -DQ alleles and haplotypes in individuals of African (n = 75), South Asian (n = 98), and mixed (n = 102) ancestry from the Caribbean island of Trinidad. We detected 19 different haplotypes containing DRB3, 8 containing DRB4, 6 containing DRB5, and 6 different haplotypes without DRB3/4/5 genes. Twenty-nine haplotypes were identified in Africans, 24 in the South Asians, and 32 in the mixed group. We detected significant differences between the populations, principally at the DQA1 and DQB1 loci, although the allele frequency for DRB1*0901 was highest in the Africans (p(c) < 0.05). Trinidad African and mixed groups were generally more diverse than the South Asians and displayed a wider range of DRB1-DQB1 associations; DQB1*02 and DQB1*0301 each associated with five to six different DRB1 alleles in the Africans and mixed group but only two in South Asians. In the Africans and the mixed group, DQB1*04 was found with DRB1*0302 and DRB1*04, but only with DRB1*08 in the South Asians. Trinidad Africans revealed consistencies with populations in Western, Central, and Northern Africa, but differed considerably from individual populations on the African continent. Trinidad South Asians displayed similar allele frequencies and associations to other populations from Northern India.