Common and well-documented HLA alleles: 2012 update to the CWD catalogue.

We have updated the catalogue of common and well-documented (CWD) human leukocyte antigen (HLA) alleles to reflect current understanding of the prevalence of specific allele sequences. The original CWD catalogue designated 721 alleles at the HLA-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, and -DPB1 loci in IMGT (IMmunoGeneTics)/HLA Database release 2.15.0 as being CWD. The updated CWD catalogue designates 1122 alleles at the HLA-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, -DPA1 and -DPB1 loci as being CWD, and represents 14.3% of the HLA alleles in IMGT/HLA Database release 3.9.0. In particular, we identified 415 of these alleles as being 'common' (having known frequencies) and 707 as being 'well-documented' on the basis of ~140,000 sequence-based typing observations and available HLA haplotype data. Using these allele prevalence data, we have also assigned CWD status to specific G and P designations. We identified 147/151 G groups and 290/415 P groups as being CWD. The CWD catalogue will be updated on a regular basis moving forward, and will incorporate changes to the IMGT/HLA Database as well as empirical data from the histocompatibility and immunogenetics community. This version 2.0.0 of the CWD catalogue is available online at cwd.immunogenomics.org, and will be integrated into the Allele Frequencies Net Database, the IMGT/HLA Database and National Marrow Donor Program's bioinformatics web pages.

HLA class II alleles and susceptibility and resistance to insulin dependent diabetes mellitus in Mexican-American families.

The role of HLA class II alleles in genetic predisposition to insulin dependent diabetes mellitus (IDDM) was examined by PCR/oligonucleotide probe typing of 42 Mexican-American IDDM families derived from Hispanic Caucasians and Native Americans. All high risk haplotypes (HLA-DR3 and DR4) were of European origin while the most strongly protective haplotype (DRB1*1402) was Native American. Of the 16 DR-DQ DR4 haplotypes identified, only those bearing DQB1*0302 conferred risk; the DRB1 allele, however, also markedly influenced IDDM risk. The general pattern of neutral and protective haplotypes indicates that the presence of Asp-57 in the HLA-DQ beta chain does not confer IDDM protection per se and indicates that both DRB1 and DQB1 influence IDDM susceptibility as well as protection.

HLA DR-DQ associations with cervical carcinoma show papillomavirus-type specificity.

Cervical carcinoma is now known to be associated with human papillomaviruses (HPV), but the evidence for a link with specific HLA loci is controversial. The role of genetic variation at the HLA class II loci and among HPV types in cervical carcinoma was investigated by PCR DNA amplification and oligonucleotide probe typing of paraffin-embedded invasive cervical cancer tissue from Hispanic patients and of cervical swabs from Hispanic controls. Certain HLA class II haplotypes (such as DRB1*1501-DQB1*0602) were associated significantly, while DR13 haplotypes were negatively associated with cervical carcinoma. These associations are HPV16-type specific. These results suggest that specific HLA class II haplotypes may influence the immune response to specific HPV-encoded epitopes and affect the risk of cervical neoplasia.

Genetic differentiation of Jewish populations.

The Jewish diaspora can be viewed as a natural process in population dispersion and differentiation. We extend genetic studies on the Jewish diaspora to an analysis of human leukocyte antigen (HLA) haplotype distributions in the Jewish peoples, and show the value of this information for the design of Jewish marrow donor registries. HLA data from the Hadassah Bone Marrow Registry having parental country-of-origin information comprise samples of geographically discrete regions. We analyzed the HLA allele and haplotype frequencies for each national sample using population genetic and clustering methods. Population differentiation among diaspora populations was shown on the basis of HLA haplotype frequencies, including differences within the more recently diverged European groups. A method of haplotype and population clustering showed patterns of unique haplotype affinities associated with specific Jewish populations. The evidence showed that diaspora Jewish populations can be sorted into distinct clades of which the Ashkenazi are but one. Relationships among Jewish populations are interpretable in light of the historical record. We suggest that a major contributing factor to the genetic divergence between Jewish groups may have been admixture with local host populations, while, at the same time, threads of Eastern Mediterranean ancestry remain evident.

Four-locus high-resolution HLA typing in a sample of Mexican Americans.

Mexicans are the most common minority population of the United States. From a sample of 553 bone marrow donor registrants of self-described Mexican ancestry, human leukocyte antigen (HLA) loci A, C, B and DRB1 were typed by high resolution sequence based typing (SBT) methods. A total of 47, 34, 76 and 46 distinct alleles at A, C, B and DRB1 respectively were identified, including 3 new alleles. The four-locus haplotype frequency distribution was extremely skewed with only 53.9% of 1106 chromosomes present with more than one estimated copy. Haplotypes of Native American origin were identified. These data form an initial basis for determining the requirements for an adequate donor pool for stem cell transplantation in this population.

Association of HLA-DPB1*0301 with IDDM in Mexican-Americans.

Susceptibility to IDDM has been associated with specific alleles at the HLA class II loci in a variety of human populations. Previous studies among Mexican-Americans, a group ancestrally derived from Native Americans and Hispanic whites, showed that the DR4 haplotypes (DRB1*0405-DQB1*0302 and DRB1*0402-DQB1*0302) and the DR3 haplotype (DRB1*0301-DQB1*0201) were increased among patients and suggested a role for both DR and DQ alleles in susceptibility and resistance. Based on the analysis of 42 Mexican-American IDDM families and ethnically matched control subjects by polymerase chain reaction/sequence-specific oligonucleotide probe typing, we report an association of IDDM with the DPB1 allele, *0301 (relative risk = 6.6; P = 0.0012) in this population. The analysis of linkage disequilibrium patterns in this population indicates that the observed increased frequency in DPB1*0301 among patients cannot be attributed simply to linkage disequilibrium with high-risk DR-DQ haplotypes. These data suggest that in addition to alleles at the DRB1 and DQB1 loci, polymorphism at the DPB1 locus may also influence IDDM risk.

Disequilibrium pattern analysis. I. Theory.

We have developed a method, disequilibrium pattern analysis, for examining the disequilibrium distribution of the entire array of two locus multiallelic haplotypes in a population. It is shown that a selected haplotype will produce a distinct pattern of linkage disequilibrium values for all generations while the selection is acting. This pattern will also presumably be maintained for many generations after the selection event, until the disequilibrium pattern is eventually broken down by genetic drift and recombination. Related haplotypes, sharing an allele with a selected haplotype, assume a value of linkage disequilibrium proportional to the frequency of the unshared allele and have a single negative value of the normalized linkage disequilibrium. The analysis assumes zero linkage disequilibrium for all allelic combinations initially. The same basic results continue to apply if the selection involves a new mutant, the occurrence of which creates linkage disequilibrium for some haplotypes. The disequilibrium pattern predicted under selection is robust with respect to the influence of migration and random genetic drift. This method is applicable to population data having linked polymorphic loci including that determined from protein or DNA sequencing.

Disequilibrium pattern analysis. II. Application to Danish HLA A and B locus data.

Disequilibrium pattern analysis, a general method for analyzing evolutionary events acting on pairs of tightly linked polymorphic loci, is applied to a large sample of Danish individuals typed for A and B loci of the HLA (human leukocyte antigen) system. Cases of selection on particular haplotypes are revealed from patterns of linkage disequilibrium among the HLA haplotypes. These patterns cannot be explained by either population admixture or random genetic drift. Six haplotypes out of the total array of 273 haplotypes have been identified which show in varying extents the patterns indicating selection.

Constrained disequilibrium values and hitchhiking in a three-locus system.

Positive selection on a new mutant allele can increase the frequencies of closely linked alleles (through hitchhiking), as well as create linkage disequilibrium between them. Because this disequilibrium is induced by the selected allele, one may be able to identify loci under selection by measuring the influence of a candidate locus on pairwise disequilibrium values at nearby loci. The constrained disequilibrium values (CDV) method approaches this problem by examining differences in pairwise disequilibrium values, which have been normalized for two- and three-locus systems, respectively. We have investigated in detail the reliability of inferences based on CDV, using simulation and analytical methods. Our main results are (i) in some circumstances, CDV may not distinguish well between a selected locus and a neighboring neutral locus, but (ii) CDV seldom indicates "selection" in neutral haplotypes with moderate to large 4Nc. We conclude that, although the CDV method does not appear to precisely locate selected alleles, it can be used to screen for regions in which hitchhiking is a plausible hypothesis. We present a microsatellite data set from human chromosome 6, in which constrained disequilibrium values suggest the action of selection in a region containing the human leukocyte antigen (HLA)-A and myelin oligodendrocyte glycoprotein (MOG) loci. The connection between hitchhiking and disequilibrium has received relatively little attention, so our investigation presents opportunities to address more general issues.

Gene conversion in the evolution of both the H-2 and Qa class I genes of the murine major histocompatibility complex.

In order to better understand the role of gene conversion in the evolution of the class I gene family of the major histocompatibility complex (MHC), we have used a computer algorithm to detect clustered sequence similarities among 24 class I DNA sequences from the H-2, Qa, and Tla regions of the murine MHC. Thirty-four statistically significant clusters were detected; individual analysis of the clusters suggested at least 25 past gene conversion or recombination events. These clusters are comparable in size to the conversions observed in the spontaneously occurring H-2Kbm and H-2Kkm2 mutations, and are distributed throughout all exons of the class I gene. Thus, gene conversion does not appear to be restricted to the regions of the class I gene encoding their antigen-presentation function. Moreover, both the highly polymorphic H-2 loci and the relatively monomorphic Qa and Tla loci appear to have participated as donors and recipients in conversion events. If gene conversion is not limited to the highly polymorphic loci of the MHC, then another factor, presumably natural selection, must be responsible for maintaining the observed differences in level of variation.

Selection, hitchhiking and disequilibrium analysis at three linked loci with application to HLA data.

The HLA system has been extensively studied from an evolutionary perspective. Although it is clear that selection has acted on the genes in the HLA complex, the nature of this selection has yet to be fully clarified. A study of constrained disequilibrium values is presented that is applicable to HLA and other less polymorphic systems with three or more linked loci, with the purpose of identifying selection events. The method uses the fact that three locus systems impose additional constraints on the range of possible disequilibrium values for any pair of loci. We have thus examined the behavior of the normalized pairwise disequilibrium measures using two locus (D'), and also three locus (D"), constraints on pairwise disequilibria in a three locus system when one of the three loci is under positive selection. The difference between these measures, delta = magnitude of D' - magnitude of D", has a distribution for the two unselected loci differing from that for the selected locus with either of the unselected loci (the hallmark is a high positive value of delta for the two unselected loci). An examination of genetic drift indicates that positive delta values are unlikely to be found in human populations in the absence of selection when recombination is greater than about 0.1%. This measure can thus provide insight into which allele of several linked loci might have been subject to selection. Application of this method to HLA haplotypes from a large French population study (Provinces Francaise) identifies selected alleles on particular haplotypes. Application of a complementary method, disequilibrium pattern analysis also confirms the action of selection on these haplotypes.

Evolution of HLA class II molecules: Allelic and amino acid site variability across populations.

Analysis of the highly polymorphic beta1 domains of the HLA class II molecules encoded by the DRB1, DQB1, and DPB1 loci reveals contrasting levels of diversity at the allele and amino acid site levels. Statistics of allele frequency distributions, based on Watterson's homozygosity statistic F, reveal distinct evolutionary patterns for these loci in ethnically diverse samples (26 populations for DQB1 and DRB1 and 14 for DPB1). When examined over all populations, the DQB1 locus allelic variation exhibits striking balanced polymorphism (P < 10(-4)), DRB1 shows some evidence of balancing selection (P < 0.06), and while there is overall very little evidence for selection of DPB1 allele frequencies, there is a trend in the direction of balancing selection (P < 0.08). In contrast, at the amino acid level all three loci show strong evidence of balancing selection at some sites. Averaged over polymorphic amino acid sites, DQB1 and DPB1 show similar deviation from neutrality expectations, and both exhibit more balanced polymorphic amino acid sites than DRB1. Across ethnic groups, polymorphisms at many codons show evidence for balancing selection, yet data consistent with directional selection were observed at other codons. Both antigen-binding pocket- and non-pocket-forming amino acid sites show overall deviation from neutrality for all three loci. Only in the case of DRB1 was there a significant difference between pocket- and non-pocket-forming amino acid sites. Our findings indicate that balancing selection at the MHC occurs at the level of polymorphic amino acid residues, and that in many cases this selection is consistent across populations.

Genetic variability and linkage disequilibrium within the DP region in the CEPH families.

A PCR-based SSO-assay has been developed to characterize the allelic polymorphism at the HLA-DPA1 locus. To validate the performance of this assay, 77 samples were typed side by side in a blinded fashion by the SSO assay and sequencing-based typing (SBT); 100% concordance was seen between the two methods. To address questions of genetic variability and linkage disequilibrium within the class II region, 478 members of the 37 original Caucasian Centre d'Etude du Polymorphisme Humain (CEPH) families were typed for DPA1 using the SSO assay providing information on 247 independent chromosomes. Six of the eight known DPA1 alleles were detected in this population; DPA1*0103 was the most frequent allele. Analysis of the distribution of allele and haplotype frequencies using the homozygosity statistic suggests that balancing selection does not appear to be acting on the DPA1 locus nor on the functional DP heterodimer in this population. Family data permits the unambiguous assignment of haplotypes. Of the 247 independent chromosomes analyzed, 24 distinct DPA1-DPB1 haplotypes were identified with DPA1*0103-DPB1*0401 being the most common. Twelve of the 18 DPB1 alleles identified in this population have an exclusive association with one DPA1 allele. Of the remaining six DPB1 alleles, four are present at a frequency of >3% and show preferential association with just one DPA1 allele. Calculation of the normalized disequilibrium parameter (D') shows 13 DP haplotypes to be in significant positive disequilibrium. These data suggest there is strong linkage disequilibrium between the DPA1 and DPB1 loci in this Caucasoid population and provide a basis with which to study linkage disequilibrium in other ethnic groups as well as analyze the evolutionary forces which govern allelic and haplotypic variation.

A comprehensive search for segregation distortion in HLA.

Segregation distortion, the non-Mendelian segregation of gametes, has been well documented among diverse groups of organisms. These cases are characterized by extreme segregation ratios found only in males. Previous reports have suggested the existence of segregation distortion operating in the HLA system of humans, a tightly linked complex of genes which regulates the immune system. In mice, some alleles of the T/t complex, which is linked to H-2 (the HLA homologue of mice), cause extreme segregation distortion in wild mice populations. Here we report on the examination of a large body of pedigree data on non-diseased families, scored for the alleles of five HLA region loci. We searched for segregation distortion on the basis of five different models of inheritance: allelic, haplotypic, genotypic, diffuse occurrence in families, and autosomal effects on the sex ratio. There was no clear evidence for segregation distortion. In particular, the possibility of extreme levels of segregation distortion was firmly rejected in the populations examined, thus reducing the likelihood of common distortion-causing HLA associated haplotypes in our species.

Molecular analysis of the role of the HLA class II genes DRB1, DQA1, DQB1, and DPB1 in susceptibility to Lyme arthritis.

Using the polymerase chain reaction and sequence-specific oligonucleotide probes, we have determined the distribution of the alleles at the polymorphic class II loci, DRB1, DQA1, DQB1, and DPB1 in 25 patients with Lyme arthritis. Comparison of the frequencies of the DRB1 and DPB1 alleles in Lyme arthritis patients to the Centre d'Etude du Polymorphisme Humaine control panel revealed statistically significant differences between the two groups. An increase in DRB1*1301 is noted along with a unique distribution of the DR4 subtypes within the 9 DR4+ patients. In addition, an increase in the rate DPB1*1001 and the more common DPB1*0201 alleles is reported.

Evolution of the CCR5 Delta32 mutation based on haplotype variation in Jewish and Northern European population samples.

The chemokine receptor 5 (CCR5) serves as a fusion cofactor for macrophage-tropic strains of HIV-1. In addition, CCR5 has been shown to mediate the entry of poxviruses into target cells. Individuals homozygous for the Delta32 deletion-mutation have no surface expression of CCR5 and are highly protected against HIV-1 infection. To gain insights into the evolution of the mutation in modern populations, the relatively high frequency of the Delta32-ccr5 allele in some European and Jewish populations is explored here by examining haplotypes of 3p21.3 constructed of five polymorphic marker loci surrounding CCR5. By sampling Ashkenazi, non-Ashkenazi and non-Jewish populations, we utilize the natural experiment that occurred as a consequence of the Jewish Diaspora, and demonstrate that a single mutation was responsible for all copies of Delta32. This mutation must have moved from Northern European populations to the Ashkenazi Jews where evidence suggests that Delta32 carriers of both groups were favored by repeated occurrence of epidemic small pox beginning in the 8th century AD.

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.

Analysis of genotype frequencies and interlocus association for the PM, DQA1, and D1S80 loci in four populations.

Allele frequencies of the LDLR, HBGG, GYPA, D7S8, GC, DQA1, and D1S80 loci are presented and genotypes are analyzed for each of four ethnic groups: African Americans (n = 200), US Caucasians (n = 200), US Hispanics (n = 200), and Japanese (n = 89). Hardy-Weinberg genotypic proportions were observed in all but two of the 28 population-locus tests undertaken. Those two instances are attributable to type I statistical error. Gametic equilibrium among loci is an assumption invoked for application of the product rule to utilize the discriminatory power from two or more loci simultaneously. Two statistical methods, a genotype matching statistic and log-linear modeling, were used to evaluate gametic disequilibrium. The match statistic, comparing observed to expected likelihood of genotypic identity for seven loci among pairs of individuals within the database, revealed only one statistically significant deviation among 20 tests. As expected, the probability of match was generally lowest in the test on all ethnic groups combined, indicating that allele frequencies differ among ethnic groups for some of the loci. This was confirmed with the statistic theta to measure ethnic stratification, in which about 0.10 of the genetic variation is apportioned among the four ethnic groups for four of the structural loci (LDLR, HBGG, GC, and DQA1), while for GYPA, D7S8, and D1S80, variation is more uniformly distributed among ethnic groups. Log-linear modeling was also applied to the five PM loci. The most parsimonious log-linear model included only three higher order terms: the two-way interactions of three of the PM loci with ethnic group. These three instances (LDLR, HBGG, and GC) indicated differences in allele frequencies between ethnic groups. No two or higher way interaction (disequilibrium) was observed among loci. In summary, the assumptions of Hardy-Weinberg and gametic equilibrium that facilitate the use of the five PM loci, DQA1 and D1S80 in forensic applications are consistent with the allele and genotype frequencies observed in these populations.