Allele-Level Haplotype Frequencies and Pairwise Linkage Disequilibrium for 14 KIR Loci in 506 European-American Individuals

The immune responses of natural killer cells are regulated, in part, by killer cell immunoglobulin-like receptors (KIR). The 16 closely-related genes in the KIR gene system have been diversified by gene duplication and unequal crossing over, thereby generating haplotypes with variation in gene copy number. Allelic variation also contributes to diversity within the complex. In this study, we estimated allele-level haplotype frequencies and pairwise linkage disequilibrium statistics for 14 KIR loci. The typing utilized multiple methodologies by four laboratories to provide at least 2x coverage for each allele. The computational methods generated maximum-likelihood estimates of allele-level haplotypes. Our results indicate the most extensive allele diversity was observed for the KIR framework genes and for the genes localized to the telomeric region of the KIR A haplotype. Particular alleles of the stimulatory loci appear to be nearly fixed on specific, common haplotypes while many of the less frequent alleles of the inhibitory loci appeared on multiple haplotypes, some with common haplotype structures. Haplotype structures cA01 and/or tA01 predominate in this cohort, as has been observed in most populations worldwide. Linkage disequilibrium is high within the centromeric and telomeric haplotype regions but not between them and is particularly strong between centromeric gene pairs KIR2DL5∼KIR2DS3S5 and KIR2DS3S5∼KIR2DL1, and telomeric KIR3DL1∼KIR2DS4. Although 93% of the individuals have unique pairs of full-length allelic haplotypes, large genomic blocks sharing specific sets of alleles are seen in the most frequent haplotypes. These high-resolution, high-quality haplotypes extend our basic knowledge of the KIR gene system and may be used to support clinical studies beyond single gene analysis.     

Killer Ig-like receptor haplotype analysis by gene content: evidence for genomic diversity with a minimum of six basic framework haplotypes,each with multiple subsets

Killer Ig-like receptor (KIR) genes constitute a multigene family whose genomic diversity is achieved through differences in gene content and allelic polymorphism. KIR haplotypes containing a single activating KIR gene (A-haplotypes), and KIR haplotypes with multiple activating receptor genes (B-haplotypes) have been described. We report the evaluation of KIR gene content in extended families, sibling pairs, and an unrelated Caucasian panel through identification of the presence or absence of 14 KIR genes and 2 pseudogenes. Haplotype definition included subtyping for the expressed and nonexpressed KIR2DL5 variants, for two alleles of pseudogene 3DP1, and for two alleles of 2DS4, including a novel 2DS4 allele, KIR1D. KIR1D appears functionally homologous to the rhesus monkey KIR1D and likely arose as a consequence of a 22 nucleotide deletion in the coding sequence of 2DS4, leading to disruption of Ig-domain 2D and a premature termination codon following the first amino acid in the putative transmembrane domain. Our investigations identified 11 haplotypes within 12 families. From 49 sibling pairs and 17 consanguineous DNA samples, an additional 12 haplotypes were predicted. Our studies support a model for KIR haplotype diversity based on six basic gene compositions. We suggest that the centromeric half of the KIR genomic region is comprised of three major combinations, while the telomeric half can assume a short form with either 2DS4 or KIR1D or a long form with multiple combinations of several stimulatory KIR genes. Additional rare haplotypes can be identified, and may have arisen by gene duplication, intergenic recombination, or deletions.     

<Emphasis Type="Italic">KIR</Emphasis> haplotypes defined by segregation analysis in 59 Centre d’Etude Polymorphisme Humain (CEPH) families

The killer cell immunoglobulin-like receptor (KIR) gene cluster exhibits extensive allelic and haplotypic diversity. Variation at the locus is associated with an increasing number of human diseases, reminiscent of the HLA loci. Characterization of diversity at the KIR locus has progressed over the past several years, particularly since the sequence of entire KIR haplotypes have become available. To determine the extent of KIR haplotypic variability among individuals of northern European descent, we genotyped 59 CEPH families for presence/absence of all KIR genes and performed limited allelic subtyping at several KIR loci. A total of 20 unique haplotypes differing in gene content were identified, the most common of which was the previously defined A haplotype (f = 0.52). Several unusual haplotypes that probably arose as a consequence of unequal crossing over events were also identified. Linkage disequilibrium (LD) analysis indicated strong negative and positive LD between several pairs of genes, values that may be useful in determining haplotypic structure when family data are not available. These data provide a resource to aid in the interpretation of disease association data involving individuals of European descent. An erratum to this article can be found at     

Genomic structure and evolution of the human pepsinogen A multigene family

Summary A human cosmid library was screened with a pepsinogen A (PGA) cDNA probe, yielding 18 clones with (parts of) one, two or three PGA genes. By aligning these cosmids a restriction map of a PGA gene quadruplet was obtained in which the four genes are arranged in a highly ordered fashion in a head-to-tail orientation. Using the length in kilobases of the large polymorphic EcoRI fragment of the PGA genes, this quadruplet can be described as 15.0-12.0-12.0-16.6. An AvaII polymorphism allowed us to identify the two PGA haplotypes of the individual whose DNA had been cloned in the cosmid library to be a gene triplet and a gene quadruplet. By comparing the restriction maps of the central 12.0 genes in these multiplets to those of the flanking 15.0 and 16.6 genes, we postulate that these central genes arose from unequal but homologous crossing over between two 15.0–16.6 gene pairs. This hypothesis provides for the creation of a variety of haplotypes by additional cross overs and mutations. Southern blots of family and population material supports the existance of at least five common PGA haplotypes, including a single-gene haplotype, giving rise to a large number of different EcoRI patterns. The single PGA gene is probably the reciprocal crossing over product. Comparison between the DNA and protein polymorphisms suggests further micro-heterogeneity in the different PGA haplotypes.     

Definition of gene content for nine common group B haplotypes of the Caucasoid population: KIR haplotypes contain between seven and eleven KIR genes

The segregation of killer cell immunoglobulin-like receptor ( KIR) genes was determined for a panel of 21 Caucasoid families: 23 different KIR gene patterns were found and could be assigned to combinations of 16 different haplotypes. Four loci were held in common by all haplotypes: KIR2DL4, KIR3DL2, the putative pseudogene KIR3DL3 and KIR2DL2/KIR2DL3, the latter likely being alleles of one gene. Group A haplotypes, which have a unique combination of seven KIR genes, were found at 80% frequency in the family panel, the polygenic group B haplotypes at 65% frequency. KIR gene segregation was fully determined for the nine group B haplotypes, which occurred at highest frequencies in both the family panel and a panel of unrelated individuals. The group B haplotypes carried between seven and 11 KIR genes and encoded inhibitory KIR for one, two, or all three major HLA class I epitopes. Analysis of human leucocyte antigen (HLA) class I genotypes revealed that most, but not all, individuals possess an inhibitory KIR for a self HLA class I epitope. The number of stimulatory KIR genes in group B haplotypes varied considerably between one and five. The data show that group B haplotypes possess a broad spectrum of KIR gene patterns, which is largely complementary to the KIR gene set of group A haplotypes. The results suggest that rapid diversification of group B haplotypes is the result of pathogen-mediated selection for KIR genotypes that have more than the set of KIR genes provided by the group A haplotype.     

Homology Requirements for Unequal Crossing over in Humans

To gain insight into mechanisms of unequal homologous recombination in vivo, genes generated by homologous unequal crossovers in the human beta-globin gene cluster were examined by nucleotide sequencing and hybridization experiments. The naturally occurring genes studied included one delta-beta Lepore-Baltimore fusion gene, one delta-beta Lepore-Hollandia fusion gene, 12 delta-beta Lepore-Boston genes, one A gamma-beta fusion Kenya gene, one A gamma-G gamma fusion (the central gene of a triplication) and one G gamma-A gamma fusion. A comparison of the nucleotide sequences of three Lepore-Boston genes indicates that they were derived from at least two independent homologous but unequal crossover events, although the crossovers occurred within the same 58-bp region. Nine additional Lepore-Boston genes from individuals of various ethnic origins were shown, by hybridization to specific oligonucleotide probes, to have been generated by a crossover in the same region as the sequenced genes. Evidence for gene conversion accompanying a homologous unequal crossover event was found in only one case (although some of the single nucleotide differences observed in other genes in this study may be related to the crossover events in ways that we do not presently understand). Thus, as judged by this limited sample, concurrent gene conversions are not commonly associated with homologous but unequal exchange in humans in vivo. Classification of the recombinant chromosomes by their polymorphic restriction sites in the beta-globin gene cluster indicated that the Lepore-Boston genes are found in at least six different haplotype backgrounds. Therefore the total number of independent examples in this study is at least 6, and at most 12. We have shown that in at least six cases of genes that have arisen by homologous but unequal crossing over in vivo, each event occurred in a relatively extensive region of uninterrupted identity between the parental genes. This preference cannot be explained by a mechanism whereby crossovers occur at random within misaligned related but not identical genes. In general, crossovers occur in regions that are among the largest available stretches of identity for a particular pair of mismatched genes. Our data are in agreement with those of other types of studies of homologous recombination, and support the idea that sequence identity, rather than general homology, is a critical factor in homologous recombination.     

Different patterns of evolution in the centromeric and telomeric regions of group A and B haplotypes of the human killer cell Ig-like receptor locus

The fast evolving human KIR gene family encodes variable lymphocyte receptors specific for polymorphic HLA class I determinants. Nucleotide sequences for 24 representative human KIR haplotypes were determined. With three previously defined haplotypes, this gave a set of 12 group A and 15 group B haplotypes for assessment of KIR variation. The seven gene-content haplotypes are all combinations of four centromeric and two telomeric motifs. 2DL5, 2DS5 and 2DS3 can be present in centromeric and telomeric locations. With one exception, haplotypes having identical gene content differed in their combinations of KIR alleles. Sequence diversity varied between haplotype groups and between centromeric and telomeric halves of the KIR locus. The most variable A haplotype genes are in the telomeric half, whereas the most variable genes characterizing B haplotypes are in the centromeric half. Of the highly polymorphic genes, only the 3DL3 framework gene exhibits a similar diversity when carried by A and B haplotypes. Phylogenetic analysis and divergence time estimates, point to the centromeric gene-content motifs that distinguish A and B haplotypes having emerged ~6 million years ago, contemporaneously with the separation of human and chimpanzee ancestors. In contrast, the telomeric motifs that distinguish A and B haplotypes emerged more recently, ~1.7 million years ago, before the emergence of Homo sapiens. Thus the centromeric and telomeric motifs that typify A and B haplotypes have likely been present throughout human evolution. The results suggest the common ancestor of A and B haplotypes combined a B-like centromeric region with an A-like telomeric region.     

Familial clustering of IGHC deletions and duplications: functional and molecular analysis

The human immunoglobulin heavy chain constant region locus (IGHC) comprises nine genes and two pseudogenes clustered in a 350 kilobase (kb) region on chromosome 14q32. Several IGHC haplotypes with single or multiple gene deletions and duplications have been characterized. The most likely mechanism accounting for these unusual haplotypes is the unequal crossing-over between homologous regions within the locus. Here we report the analysis of an unusual case of familial clustering of deletions/duplications. In the two branches of the BON family, three duplicated and two deleted haplotypes, all probably independent in origin, have been characterized. The structure of the haplotypes, one of which is described here for the first time, supports the hypothesis of homologous unequal crossing-over as the origin of recombinant haplotypes. The analysis of serological markers in a subject carrying one deleted and one duplicated haplotype allowed us the first direct inferences concerning the functions of the duplicated IGHC haplotypes.     

<Emphasis Type="Italic">KIR</Emphasis> haplotype content at the allele level in 77 Northern Irish families

There has been an explosion in population studies determining the frequency of KIR genes. However, there is still limited knowledge of allele and haplotype frequencies in different populations. The present study aims to determine the haplotype frequencies using allele information on ten genes and presence/absence of the other seven genes in the parents of 77 families. There were 26 of 154 different genotypes without using allele information and 143 of 154 different genotypes using allele information. These genotypes came from 96 of 308 different haplotypes. Of these, 41 were A and 55 were B. Forty-nine haplotypes occurred only once. In total, 181 (58.8%) of haplotypes were A and 127 (41.2%) were B. Three different haplotypes carried two copies of KIR2DL4, two different haplotypes were truncated with both KIR2DL4 and KIR3DL1/S1 missing, and three different haplotypes were negative for both KIR2DL2 and KIR2DL3; two of these haplotypes carried KIR2DS2. A further haplotype, present in two individuals, appeared to have two alleles of KIR2DL5A present. The percentages of individuals who were homozygous for the A haplotype, heterozygous for the A and B haplotype and homozygous for the B haplotype were 35.1%, 47.4% and 17.5% respectively. The genes KIR3DL1, KIR2DS4 and KIR2DL3 were present on 31, 32 and 15 different B haplotypes, respectively, and 64, 65 and 40 of the total B haplotypes, respectively. Sixty B haplotypes had both KIR3DL1 and KIR2DS4, and four haplotypes had KIR2DS4 and KIR2DL3. However, in 40 of 41 different and 180 of 181 total A haplotypes, KIR3DL1, KIR2DS4 and KIR2DL3 were all present (we did not allele-type for KIR2DL1 and therefore could not determine presence/absence on those haplotypes). At the allele level, homozygosity was found in 22.1%, 9.7% and 12.6% for KIR2DL4, KIR3DL2 and KIR3DL1 genes, respectively, but 62.6% and 53% for KIR2DL3 and KIR2DS4 genes, respectively, despite the fact that no one allele dominated the frequency in any of these genes.     

Predominance of group A KIR haplotypes in Japanese associated with diverse NK cell repertoires of KIR expression

Genomic DNA from a panel of 41 healthy unrelated Japanese individuals was typed for the presence or absence of 16 KIR genes and pseudogenes. Only eight different KIR genotypes were found. Group A haplotypes outnumbered group B haplotypes in frequency by approximately 3:1, with individuals having two group A haplotypes accounting for 56% of the panel. The frequency of A haplotypes in the Japanese is higher than that observed in other populations. Flow cytometric comparison of KIR expression in 19 panel members showed considerable diversity in NK cell repertoire, which was also seen within the group of individuals having two A haplotypes. This diversity is likely due to allelic polymorphism in expressed genes of the A haplotype. In comparison to other populations, the Japanese appear less heterogeneous in KIR genotype as assessed by gene content.     

Haplotype inference for present-absent genotype data using previously identified haplotypes and haplotype patterns

MOTIVATION: Killer immunoglobulin-like receptor (KIR) genes vary considerably in their presence or absence on a specific regional haplotype. Because presence or absence of these genes is largely detected using locus-specific genotyping technology, the distinction between homozygosity and hemizygosity is often ambiguous. The performance of methods for haplotype inference (e.g. PL-EM, PHASE) for KIR genes may be compromised due to the large portion of ambiguous data. At the same time, many haplotypes or partial haplotype patterns have been previously identified and can be incorporated to facilitate haplotype inference for unphased genotype data. To accommodate the increased ambiguity of present-absent genotyping of KIR genes, we developed a hybrid approach combining a greedy algorithm with the Expectation-Maximization (EM) method for haplotype inference based on previously identified haplotypes and haplotype patterns. RESULTS: We implemented this algorithm in a software package named HAPLO-IHP (Haplotype inference using identified haplotype patterns) and compared its performance with that of HAPLORE and PHASE on simulated KIR genotypes. We compared five measures in order to evaluate the reliability of haplotype assignments and the accuracy in estimating haplotype frequency. Our method outperformed the two existing techniques by all five measures when either 60% or 25% of previously identified haplotypes were incorporated into the analyses. AVAILABILITY: The HAPLO-IHP is available at http://www.soph.uab.edu/Statgenetics/People/KZhang/HAPLO-IHP/index.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Evolution of killer cell Ig-like receptor (KIR) genes: definition of an orangutan KIR haplotype reveals expansion of lineage III KIR associated with the emergence of MHC-C

Orangutan (Pongo pygmaeus) MHC-C appears less evolved than human HLA-C: Popy-C is not fixed and its alleles encode only one (C1) of the two motifs for killer cell Ig-like receptor (KIR) ligands. To assess the structure and complexity of the orangutan KIR locus, the complete nucleotide sequence of an orangutan KIR haplotype was determined. The PopyKIR locus is flanked by LILR and FCAR and consists of seven genes and pseudogenes, two novel and five corresponding to known cDNA. Distinguishing all KIRs in this rapidly evolving KIR locus from the KIR3DX1 gene is an LTR33A/MLT1D element in intron 3. These two forms of KIR represent lineages that originated by duplication of a common ancestor. The conserved, framework regions of primate KIR loci comprise the 5' part of a lineage V KIR, the 3' part of a pseudogene, the complete 2DL4 gene, and the 3' part of a lineage II KIR. Although previously defined PopyKIR2DL4 alleles contain premature termination codons, the sequenced haplotype's PopyKIR2DL4 allele encodes a full-length protein. A model for KIR evolution is proposed. Distinguishing the orangutan KIR haplotype from the proposed common ancestor of primate KIR haplotypes is an increased number to give three lineage III KIR genes in the centromeric part of the locus, the site for most human lineage III genes encoding HLA-C specific KIR. Thus, expansion of lineage III KIR is associated with emergence of MHC-C.     

Allelic polymorphism synergizes with variable gene content to individualize human KIR genotype

Killer Ig-like receptor (KIR) genes are a multigene family on human chromosome 19. KIR genes occur in various combinations on different haplotypes. Additionally, KIR genes are polymorphic. To examine how allelic polymorphism diversifies KIR haplotypes with similar or identical combinations of KIR genes, we devised methods for discriminating alleles of KIR2DL1, -2DL3, -3DL1, and -3DL2. These methods were applied to 143 individuals from 34 families to define 98 independent KIR haplotypes at the allele level. Three novel 3DL2 alleles and a chimeric 3DL1/3DL2 sequence were also identified. Among the A group haplotypes were 22 different combinations of 2DL1, 2DL3, 3DL1, and 3DL2 alleles. Among the B group haplotypes that were unambiguously determined were 15 distinct haplotypes involving 9 different combinations of KIR genes. A and B haplotypes both exhibit strong linkage disequilibrium (LD) between 2DL1 and 2DL3 alleles, and between 3DL1 and 3DL2 alleles. In contrast, there was little LD between the 2DL1/2DL3 and 3DL1/3DL2 pairs that define the two halves of the KIR gene complex. The synergistic combination of allelic polymorphism and variable gene content individualize KIR genotype to an extent where unrelated individuals almost always have different KIR types. This level of diversity likely reflects strong pressure from pathogens on the human NK cell response.     

Origin of rat beta-globin haplotypes containing three and five genes

We have reported in rat three adult beta-gene haplotypes containing either five or three genes. Detailed sequence analysis reveals that the leftmost gene is the major gene and that at the opposite end downstream lies the minor gene. All of the genes lying between them are minor-major hybrids indicating their origin by unequal crossing-over. In two haplotypes beta-globin genes were found with an L1(1) element inserted directly into IVS2. The described results allow the formulation of a pathway of mutational events leading from the ancient two-beta-gene rodent ancestor through a three-gene haplotype to five-gene haplotypes, one of which is postulated to have arisen in common laboratory strains since their capture in the wild.     

The E beta hot spot of recombination in wild-derived natural recombinant MHC haplotypes. Cross-over site mapping and the identification of a 1.0-kb E beta deletion in the p and w14 haplotypes

Detailed molecular analysis of three wild-derived MHC haplotypes provided evidence for an important role of the E beta recombinational hot spot in the recent evolution of the mouse I region. Examination of RFLP and restriction maps of cloned DNA permitted the mapping of the natural cross-over events in the haplotypes carried by strains B10.GAA37 (w21) and B10.KPB128 (w19) to a fragment of DNA not exceeding 4.1 kb, which lies almost entirely within the intron separating the beta 1 and beta 2 exons of the E beta gene. In the w14 haplotype (strain B10.STC77), which appears to be a natural recombinant between a p-like parental haplotype and another wild-derived haplotype, the site of crossing over can be mapped to a segment between the beta 2 exon of the E beta gene (left border) and the E beta 2 gene (right border). This segment containing the cross-over site in the w14 haplotype includes the E beta hot spot. In addition, the w14 haplotype as well as the standard p haplotype contain a deletion of approximately 1.0 kb in the second intron of the E beta gene, which may represent the product of an unequal cross-over event in a E beta recombinational hot spot.     

Killer Immunoglobulin-Like Receptor (KIR) Centromeric-AA Haplotype Is Associated with Ethnicity and Tuberculosis Disease in a Canadian First Nations Cohort

Killer immunoglobulin-like receptors (KIR) on natural killer (NK) cells interact with other immune cells to monitor the immune system and combat infectious diseases, such as tuberculosis (TB). The balance of activating and inhibiting KIR interactions helps determine the NK cell response. In order to examine the enrichment or depletion of KIRs as well as to explore the association between TB status and inhibitory/stimulatory KIR haplotypes, we performed KIR genotyping on samples from 93 Canadian First Nations (Dene, Cree, and Ojibwa) individuals from Manitoba with active, latent, or no TB infection, and 75 uninfected Caucasian controls. There were significant differences in KIR genes between Caucasians and First Nations samples and also between the First Nations ethnocultural groups (Dene, Cree, and Ojibwa). When analyzing ethnicity and tuberculosis status in the study population, it appears that the KIR profile and centromeric haplotype are more predictive than the presence or absence of individual genes. Specifically, the decreased presence of haplotype B centromeric genes and increased presence of centromeric-AA haplotypes in First Nations may contribute to an inhibitory immune profile, explaining the high rates of TB in this population.     

Complexity in cattle <Emphasis Type="Italic">KIR</Emphasis> genes: transcription and genome analysis

Killer immunoglobulin (Ig)-like receptors (KIRs) are the major functional natural killer (NK) cell receptors in human. The presence of KIR genes has only recently been demonstrated in other (non-primate) species, and their expression, genomic arrangement, and function in these species have yet to be investigated. In this study, we describe the KIR gene family in cattle. KIR sequences were amplified from cDNA derived from four animals. Seventeen new sequences were identified in total. Some are alleles of two previously described genes, and the remainder are representative of at least four additional genes. These cDNA data, together with analysis of the cattle genome sequence, confirm that, as in humans, cattle have multiple inhibitory and activating KIR genes, with variable haplotype composition, and putative framework genes. In contrast to human, the majority of the cattle KIR genes encode three Ig-domain KIRs; most of the inhibitory genes encode only one immunoreceptor tyrosine-based inhibitory motif (ITIM), and the activating genes encode molecules with arginine rather than the more usual lysine in the transmembrane domain. A divergent gene, 2DL1, encodes a two Ig-domain KIR with an unusual D0-D2 structure, and a distinct signaling domain with two ITIMs. Similarity to pig and human two Ig-domain (D0-D2) KIRs suggest these may be more related to an ancestral gene than the other cattle KIR genes. Cattle have multiple NKG2A-related genes and at least one Ly49 gene; thus, the data presented here suggest that they have the potential to express more major histocompatibility complex-binding NK receptors than other species.     

Diversity of killer cell immunoglobulin-like receptor genes in four ethnic groups in China

Killer cell immunoglobulin-like receptors (KIRs) show extensive variation in terms of gene content and allelic polymorphisms among different populations. The aim of this study was to analyze the distribution of KIR genes in the Bulang, Nu, Yugu, and Zhuang ethnic groups, which belong to four different language families in China, and thus to provide basic KIR gene and genotype data for these Chinese ethnic groups. Genotyping of 16 KIR genes was performed in 425 unrelated individuals using the polymerase chain reaction-sequence-specific oligonucleotide probe method with the Luminex MultiAnalyte Profiling System. The four framework KIR genes were detected in all four ethnic groups. The activating KIR genes as well as the inhibitory KIR genes showed extreme diversity among these four populations. A total of 35 distinct KIR genotypes were identified, one of which was previously unknown. The four most common genotypes were identified in all four populations and comprised 66.1~91.1% of all the genotypes. The group A haplotype occurred more frequently than the group B haplotype in the Nu, Yugu, and Zhuang populations, as in other East Asian populations. In contrast, the group A and group B haplotypes occurred equally in the Bulang population. The results of the present study suggested that the KIR genes and genotypes are diverse in these four ethnic groups, and each ethnic group has its own characteristic KIR distribution. The findings with respect to KIR gene diversity in these four populations should provide relevant genomic diversity data for the future study of viral infections, autoimmune diseases, and reproductive fitness.