Jill A. Hollenbach

Analysis of the frequencies of HLA-A, B, and C alleles and haplotypes in the five major ethnic groups of the United States reveals high levels of diversity in these loci and contrasting distribution patterns in these populations

The HLA system is the most polymorphic of all human genetic systems. The frequency of HLA class I alleles and their linkage disequilibrium patterns differ significantly among human populations as shown in studies using serologic methods. Many DNA-defined alleles with identical serotypes may have variable frequencies in different populations. We typed HLA-A, B, and C loci at the allele level by PCR-based methods in 1,296 unrelated subjects from five major outbred groups living in the U.S.A (African, AFAM; Caucasians, CAU; Asian, ORI; Hispanic, HIS, and North American Natives, NAI). We detected 46, 100 and 32 HLA-A, B, and C alleles, respectively. ORI and HIS presented more alleles at each of these loci. There was lack of correlation between the levels of heterozygosity and the number of alleles detected in each population. In AFAM, heterozygosity (>90%) is maximized at all class I loci. HLA-A had the lowest heterozygosity in all populations but CAU. Tight LD was observed between HLA-B and C alleles. AFAM had weaker or nonexistent associations between alleles of HLA-A and B than other populations. Analysis of the genetic distances between these and other populations showed a close relationship between specific US populations and a population from their original continents. ORI exhibited the largest genetic distance with all the other U.S. groups and were closer to NAI. Evidence of admixture with CAU was observed for AFAM and HIS. HIS also had significant frequencies of AFAM and Mexican Indian alleles. Differences in both LD and heterozygosity levels suggest distinct evolutionary histories of the HLA loci in the geographical regions from where the U.S. populations originated.

HLA diversity, differentiation, and haplotype evolution in Mesoamerican Natives.

Genetic variation of the Human Leukocyte Antigen region (HLA) in three Amerindian populations from the Southern Mexican state of Oaxaca, the Zapotec, Mixtec and the Mixe is examined. Individuals were typed using PCR-SSOP for four class II loci (DRB1, DQA1, DQB1, DPB1) and three class I loci (HLA-A, -B, and -C). Based on known HLA distributions, European admixture ranged from 1% to 10%. Individuals with European alleles were excluded from subsequent analysis. New alleles were revealed at each of the class I loci. In general, genotype frequencies were in Hardy-Weinberg equilibrium, although some deviations were detected. Allele frequency distributions at the DRB1, DQA1, DQB1 and HLA-A loci in all populations were more even than expected under neutrality, supporting a model of balancing selection at these loci. A history of directional selection for DPB1 in all three populations was indicated, as homozygosity values were significantly above expected values. Allele frequency distributions at HLA-B and HLA-C did not differ significantly from neutrality expectations. The data also provide evidence from linkage disequilibrium that strong haplotypic associations are present across the entire HLA region in each of the populations. Significant overall linkage disequilibrium exists between all pairs of loci typed in these populations, except those which include the DPB1 locus. These associations exist despite the fact that the recombination fraction between HLA-A, in the class I region, and DQB1, in the class II region, may exceed 0.02. One explanation is that selective pressures are maintaining the relationships between particular alleles at these loci in these populations. These relationships are maintained in general across the entire HLA region in the Oaxacan Amerindians, with the exception of DPB1.

Juvenile idiopathic arthritis and HLA Class I and Class II interactions and age-at-onset effects

OBJECTIVE The aim of this study was to quantitate risk and to examine heterogeneity for HLA at high resolution in patients with the most common subtypes of juvenile idiopathic arthritis (JIA), IgM rheumatoid factor-negative polyarticular JIA and oligoarticular JIA. Use of 4-digit comprehensive HLA typing enabled great precision, and a large cohort allowed for consideration of both age at disease onset and disease subtype. METHODS Polymerase chain reaction-based high-resolution HLA typing for class I and class II loci was accomplished for 820 patients with JIA and 273 control subjects. Specific HLA epitopes, potential interactions of alleles at specific loci and between loci (accounting for linkage disequilibrium and haplotypic associations), and an assessment of the current International League of Associations for Rheumatology classification criteria were considered. RESULTS An HLA-DRB1/DQB1 effect was shown to be exclusively attributable to DRB1 and was similar between patients with oligoarticular JIA and a younger subgroup of patients with polyarticular JIA. Furthermore, patients with polyarticular JIA showed age-specific related effects, with disease susceptibility in the group older than age 6 years limited to an effect of the HLA-DRB1*08 haplotype, which is markedly different from the additional susceptibility haplotypes, HLA-DRB1*1103/1104, found in the group with oligoarticular JIA and the group of younger patients with polyarticular JIA. Also in contrast to findings for oligoarticular JIA, patients with polyarticular arthritis had no evidence of an HLA class I effect. Markers associated with a reduced risk of disease included DRB1*1501, DRB1*0401, and DRB1*0701. DRB1*1501 was shown to reduce risk across the whole cohort, whereas DRB1*0401 and DRB1*0701 were protective for selected JIA subtypes. Surprisingly, the disease predisposition mediated by DPB1*0201 in individuals without any disease-predisposing DRB1 alleles was great enough to overcome even the very strong protective effect observed for DRB1*1501. CONCLUSION Inherited HLA factors in JIA show similarities overall as well as differences between JIA subtypes.

The immunogenetics of multiple sclerosis: A comprehensive review.

Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system and common cause of non-traumatic neurological disability in young adults. The likelihood for an individual to develop MS is strongly influenced by her or his ethnic background and family history of disease, suggesting that genetic susceptibility is a key determinant of risk. Over 100 loci have been firmly associated with susceptibility, whereas the main signal genome-wide maps to the class II region of the human leukocyte antigen (HLA) gene cluster and explains up to 10.5% of the genetic variance underlying risk. HLA-DRB1*15:01 has the strongest effect with an average odds ratio of 3.08. However, complex allelic hierarchical lineages, cis/trans haplotypic effects, and independent protective signals in the class I region of the locus have been described as well. Despite the remarkable molecular dissection of the HLA region in MS, further studies are needed to generate unifying models to account for the role of the MHC in disease pathogenesis. Driven by the discovery of combinatorial associations of Killer-cell Immunoglobulin-like Receptor (KIR) and HLA alleles with infectious, autoimmune diseases, transplantation outcome and pregnancy, multi-locus immunogenomic research is now thriving. Central to immunity and critically important for human health, KIR molecules and their HLA ligands are encoded by complex genetic systems with extraordinarily high levels of sequence and structural variation and complex expression patterns. However, studies to-date of KIR in MS have been few and limited to very low resolution genotyping. Application of modern sequencing methodologies coupled with state of the art bioinformatics and analytical approaches will permit us to fully appreciate the impact of HLA and KIR variation in MS.

Immunogenetics as a tool in anthropological studies

The genes coding for the main molecules involved in the human immune system – immunoglobulins, human leucocyte antigen (HLA) molecules and killer‐cell immunoglobulin‐like receptors (KIR) – exhibit a very high level of polymorphism that reveals remarkable frequency variation in human populations. ‘Genetic marker’ (GM) allotypes located in the constant domains of IgG antibodies have been studied for over 40 years through serological typing, leading to the identification of a variety of GM haplotypes whose frequencies vary sharply from one geographic region to another. An impressive diversity of HLA alleles, which results in amino acid substitutions located in the antigen‐binding region of HLA molecules, also varies greatly among populations. The KIR differ between individuals according to both gene content and allelic variation, and also display considerable population diversity. Whereas the molecular evolution of these polymorphisms has most likely been subject to natural selection, principally driven by host–pathogen interactions, their patterns of genetic variation worldwide show significant signals of human geographic expansion, demographic history and cultural diversification. As current developments in population genetic analysis and computer simulation improve our ability to discriminate among different – either stochastic or deterministic – forces acting on the genetic evolution of human populations, the study of these systems shows great promise for investigating both the peopling history of modern humans in the time since their common origin and human adaptation to past environmental (e.g. pathogenic) changes. Therefore, in addition to mitochondrial DNA, Y‐chromosome, microsatellites, single nucleotide polymorphisms and other markers, immunogenetic polymorphisms represent essential and complementary tools for anthropological studies.

Co-evolution of Human Leukocyte Antigen (HLA) Class I Ligands with Killer-Cell Immunoglobulin-Like Receptors (KIR) in a Genetically Diverse Population of Sub-Saharan Africans

Interactions between HLA class I molecules and killer-cell immunoglobulin-like receptors (KIR) control natural killer cell (NK) functions in immunity and reproduction. Encoded by genes on different chromosomes, these polymorphic ligands and receptors correlate highly with disease resistance and susceptibility. Although studied at low-resolution in many populations, high-resolution analysis of combinatorial diversity of HLA class I and KIR is limited to Asian and Amerindian populations with low genetic diversity. At the other end of the spectrum is the West African population investigated here: we studied 235 individuals, including 104 mother-child pairs, from the Ga-Adangbe of Ghana. This population has a rich diversity of 175 KIR variants forming 208 KIR haplotypes, and 81 HLA-A, -B and -C variants forming 190 HLA class I haplotypes. Each individual we studied has a unique compound genotype of HLA class I and KIR, forming 1–14 functional ligand-receptor interactions. Maintaining this exceptionally high polymorphism is balancing selection. The centromeric region of the KIR locus, encoding HLA-C receptors, is highly diverse whereas the telomeric region encoding Bw4-specific KIR3DL1, lacks diversity in Africans. Present in the Ga-Adangbe are high frequencies of Bw4-bearing HLA-B*53:01 and Bw4-lacking HLA-B*35:01, which otherwise are identical. Balancing selection at key residues maintains numerous HLA-B allotypes having and lacking Bw4, and also those of stronger and weaker interaction with LILRB1, a KIR-related receptor. Correspondingly, there is a balance at key residues of KIR3DL1 that modulate its level of cell-surface expression. Thus, capacity to interact with NK cells synergizes with peptide binding diversity to drive HLA-B allele frequency distribution. These features of KIR and HLA are consistent with ongoing co-evolution and selection imposed by a pathogen endemic to West Africa. Because of the prevalence of malaria in the Ga-Adangbe and previous associations of cerebral malaria with HLA-B*53:01 and KIR, Plasmodium falciparum is a candidate pathogen.

Tracking human migrations by the analysis of the distribution of HLA alleles, lineages and haplotypes in closed and open populations

The human leucocyte antigen (HLA) system shows extensive variation in the number and function of loci and the number of alleles present at any one locus. Allele distribution has been analysed in many populations through the course of several decades, and the implementation of molecular typing has significantly increased the level of diversity revealing that many serotypes have multiple functional variants. While the degree of diversity in many populations is equivalent and may result from functional polymorphism(s) in peptide presentation, homogeneous and heterogeneous populations present contrasting numbers of alleles and lineages at the loci with high-density expression products. In spite of these differences, the homozygosity levels are comparable in almost all of them. The balanced distribution of HLA alleles is consistent with overdominant selection. The genetic distances between outbred populations correlate with their geographical locations; the formal genetic distance measurements are larger than expected between inbred populations in the same region. The latter present many unique alleles grouped in a few lineages consistent with limited founder polymorphism in which any novel allele may have been positively selected to enlarge the communal peptide-binding repertoire of a given population. On the other hand, it has been observed that some alleles are found in multiple populations with distinctive haplotypic associations suggesting that convergent evolution events may have taken place as well. It appears that the HLA system has been under strong selection, probably owing to its fundamental role in varying immune responses. Therefore, allelic diversity in HLA should be analysed in conjunction with other genetic markers to accurately track the migrations of modern humans.

Report from the killer immunoglobulin-like receptor (KIR) anthropology component of the 15th International Histocompatibility Workshop: worldwide variation in the KIR loci and further evidence for the co-evolution of KIR and HLA

The killer immunoglobulin-like receptor (KIR) anthropology component of the 15th International Histocompatibility Workshop (IHIWS) sought to explore worldwide population variation in the KIR loci, and to examine the relationship between KIR genes and their human leukocyte antigen (HLA) ligands. Fifteen laboratories submitted KIR genotype and HLA ligand data in 27 populations from six broad ethnic groups. Data were analyzed for correlations between the frequencies of KIR and their known HLA ligands. In addition, allelic typing was performed for KIR2DL2 and 3DL1 in a subset of populations. Strong and significant correlations were observed between KIR2DL2, 2DL3 genotype frequencies and the frequency of their ligand, HLA-C1. In contrast, only weak associations were seen for 3DL1, 3DS1 and the HLA-Bw4 ligand. Although some aspects of the correlations observed here differ from those reported in other populations, these data provide additional evidence of linked evolutionary histories for some KIR and HLA loci. Investigation of allele-level variation for the B haplotype locus KIR 2DL2 showed that two alleles, *001 and *003, predominate in all populations in this study. Much more allelic variation was observed for the A haplotype locus 3DL1, with several alleles observed at moderate frequencies and extensive variation observed between populations.

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.

Susceptibility to Crohn’s Disease is Mediated by KIR2DL2/KIR2DL3 Heterozygosity and the HLA-C Ligand

In the present study, we investigated the relationship between the KIR loci and the genes encoding their HLA ligands and genetic susceptibility to Crohn’s disease (CD). Analyses of the interactions between KIR3DL1, KIR2DL1, KIR2DL2, and KIR2DL3 with their respective HLA ligands indicate that there is a protective effect for KIR2DL2 in the absence of its HLA ligand C1. Given that KIR2DL2 and KIR2DL3 segregate as alleles, we compared their genotypic distributions to expectations under Hardy–Weinberg Equilibrium (HWE) with regard to the HLA ligand C1 status. While all the genotypic distributions conform to expectations under HWE in controls, in C2 ligand homozygous cases there is significant deviation from HWE, with a reduction of KIR2DL2, KIR2DL3 heterozygotes. KIR2DL2, KIR2DL3 heterozygosity is the only genotypic combination that confers protection from CD. In addition to the protective effect (OR = 0.44, CI = 0.22–0.87; p = 0.018) observed in C2 ligand homozygotes, the KIR2DL2, KIR2DL3 genotype is predisposing (OR = 1.34, CI = 1.03–4.53; p = 0.031) in the presence of C1 ligand. A test for trend of HLA class I C ligand group genotypes with KIR2DL2, KIR2DL3 heterozygosity in cases and controls indicates that C1, C2 ligand group heterozygotes have an intermediate effect on predisposition. These results show for the first time that disease susceptibility may be related to heterozygosity at a specific KIR locus, and that HLA ligand genotype influences the relative effect of the KIR genotype.