Significant variation between SNP-based HLA imputations in diverse populations: the last mile is the hardest.

Four single nucleotide polymorphism (SNP)-based human leukocyte antigen (HLA) imputation methods (e-HLA, HIBAG, HLA*IMP:02 and MAGPrediction) were trained using 1000 Genomes SNP and HLA genotypes and assessed for their ability to accurately impute molecular HLA-A, -B, -C and -DRB1 genotypes in the Human Genome Diversity Project cell panel. Imputation concordance was high (>89%) across all methods for both HLA-A and HLA-C, but HLA-B and HLA-DRB1 proved generally difficult to impute. Overall, <27.8% of subjects were correctly imputed for all HLA loci by any method. Concordance across all loci was not enhanced via the application of confidence thresholds; reliance on confidence scores across methods only led to noticeable improvement (+3.2%) for HLA-DRB1. As the HLA complex is highly relevant to the study of human health and disease, a standardized assessment of SNP-based HLA imputation methods is crucial for advancing genomic research. Considerable room remains for the improvement of HLA-B and especially HLA-DRB1 imputation methods, and no imputation method is as accurate as molecular genotyping. The application of large, ancestrally diverse HLA and SNP reference data sets and multiple imputation methods has the potential to make SNP-based HLA imputation methods a tractable option for determining HLA genotypes.

Elevations in serum muscle enzyme activities in racehorses due to unaccustomed exercise and training.

Hereditary muscular disease is well described in racehorses, yet little is known about traumatic muscle disease associated with unaccustomed exercise or training. The objective of the study was to compare sedentary horses, racehorses undergoing training for the first time (unaccustomed exercise), and experienced racehorses during a training season (accustomed exercise) to investigate the effect of exercise and training on serum muscle enzyme activities and other variables. Horses were sampled prior to exercise for serum activities of aspartate amino transferase (AST), creatine kinase and other variables γ glutamyl transferase (GGT) and serum amyloid A (SAA) in a three-part study. Serum activities of AST and GGT were higher in fit racehorses (n=47) compared with sedentary horses (n=57) at a single time point (P<0.05). The monthly serum activity of AST in two-year-old racehorses (n=10) increased from month 1 to 4 of unaccustomed training (P<0.05). The serum activities of AST and GGT in three-year-old racehorses (n=12) previously accustomed to exercise sampled fortnightly to monthly showed a linear increase (P<0.05) with cumulative training days over the seven months of training, but showed minimal increase during the first four months of accustomed training. SAA was weakly correlated to cumulative training days. In conclusion, AST activity was increased by unaccustomed exercise and cumulative training stress in the racehorse. GGT appeared to be correlated to cumulative training load. Mild to moderate elevations in serum AST in racehorses may be associated with cumulative muscle damage from training or trauma associated with unaccustomed exercise.

Genotype List String: a grammar for describing HLA and KIR genotyping results in a text string.

Knowledge of an individual's human leukocyte antigen (HLA) genotype is essential for modern medical genetics, and is crucial for hematopoietic stem cell and solid-organ transplantation. However, the high levels of polymorphism known for the HLA genes make it difficult to generate an HLA genotype that unambiguously identifies the alleles that are present at a given HLA locus in an individual. For the last 20 years, the histocompatibility and immunogenetics community has recorded this HLA genotyping ambiguity using allele codes developed by the National Marrow Donor Program (NMDP). While these allele codes may have been effective for recording an HLA genotyping result when initially developed, their use today results in increased ambiguity in an HLA genotype, and they are no longer suitable in the era of rapid allele discovery and ultra-high allele polymorphism. Here, we present a text string format capable of fully representing HLA genotyping results. This Genotype List (GL) String format is an extension of a proposed standard for reporting killer-cell immunoglobulin-like receptor (KIR) genotype data that can be applied to any genetic data that use a standard nomenclature for identifying variants. The GL String format uses a hierarchical set of operators to describe the relationships between alleles, lists of possible alleles, phased alleles, genotypes, lists of possible genotypes, and multilocus unphased genotypes, without losing typing information or increasing typing ambiguity. When used in concert with appropriate tools to create, exchange, and parse these strings, we anticipate that GL Strings will replace NMDP allele codes for reporting HLA genotypes.

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.

16(th) IHIW: immunogenomic data-management methods. report from the immunogenomic data analysis working group (IDAWG).

SUMMARY: The goal of the immunogenomic data analysis working group (IDAWG) is to facilitate the consistent analysis of HLA and KIR data, and the sharing of those data among the immunogenomic and larger genomic communities. However, the data management approaches currently applied by immunogenomic researchers are not widely discussed or reported in the literature, and the effect of different approaches on data analyses is not known. With ASHI's support, the IDAWG developed a 45 question survey on HLA and KIR data generation, data management and data analysis practices. Survey questions detailed the loci genotyped, typing systems used, nomenclature versions reported, computer operating systems and software used to manage and transmit data, the approaches applied to resolve HLA ambiguity and the methods used for basic population-level analyses. Respondents were invited to demonstrate their HLA ambiguity resolution approaches in simulated data sets. By May 2012, 156 respondents from 35 nations had completed the survey. These survey respondents represent a broad sampling of the Immunogenomic community; 52% were European, 30% North American, 10% Asian, 4% South American and 4% from the Pacific. The project will continue in conjunction with the 17th Workshop, with the aim of developing community data sharing standards, ambiguity resolution documentation formats, single-task data Management tools and novel data analysis methods and applications. While additional project details and plans for the 17th IHIW will be forthcoming, we welcome the input and participation in these projects from the histocompatibility and immunogenetics community.

16(th) IHIW: extending the number of resources and bioinformatics analysis for the investigation of HLA rare alleles.

Continuing a project presented at the 15th International HLA and Immunogenetics Workshop (IHIWS) on the rarity of HLA alleles, we sought to expand the number of data sources and bioinformatics tools available in the Allele Frequencies Net Database website (AFND, www.allelefrequencies.net). In this 16th IHIWS Rare Alleles project, HLA alleles described in the latest IMGT/HLA Database (release 3.8.0) were queried against different sources including data from registries (stem cell) and from 74 different laboratories around the world. We demonstrated that approximately 40% of the alleles officially named in the IMGT/HLA Database have been reported only once across all different sources. To facilitate the large-scale analysis of rare alleles, we have produced an online tool called the Rare Allele Detector that simplifies the detection of alleles that are considered to be 'very rare', 'rare' or 'frequent'. Tools and associated data can be accessed via the www.allelefrequencies.net website.

A community standard for immunogenomic data reporting and analysis: proposal for a STrengthening the REporting of Immunogenomic Studies statement.

Modern high-throughput HLA and KIR typing technologies are generating a wealth of immunogenomic data with the potential to revolutionize the fields of histocompatibility and immune-related disease association and population genetic research, much as SNP-based approaches have revolutionized association research. The STrengthening the REporting of Genetic Association studies (STREGA) statement provides community-based data reporting and analysis standards for genomic disease-association studies, identifying specific areas in which adoption of reporting guidelines can improve the consistent interpretation of genetic studies. While aspects of STREGA can be applied to immunogenomic studies, HLA and KIR research requires additional consideration, as the high levels of polymorphism associated with immunogenomic data pose unique methodological and computational challenges to the synthesis of information across datasets. Here, we outline the principle challenges to consistency in immunogenomic studies, and propose that an immunogenomic-specific analog to the STREGA statement, a STrengthening the REporting of Immunogenomic Studies (STREIS) statement, be developed as part of the 16th International HLA and Immunogenetics Workshop. We propose that STREIS extends at least four of the 22 elements of the STREGA statement to specifically address issues pertinent to immunogenomic data: HLA and KIR nomenclature, data-validation, ambiguity resolution, and the analysis of highly polymorphic genetic systems. As with the STREGA guidelines, the intent behind STREIS is not to dictate the design of immunogenomic studies, but to ensure consistent and transparent reporting of research, facilitating the synthesis of HLA and KIR data across studies.

Allele Name Translation Tool and Update NomenCLature: software tools for the automated translation of HLA allele names between successive nomenclatures.

In this brief communication, we describe the Allele Name Translation Tool (antt) and Update NomenCLature (uncl), free programs developed to facilitate the translation of human leukocyte antigen (HLA) allele names recorded using the December 2002 version of the HLA allele nomenclature (e.g. A*01010101) to those recorded using the colon-delimited version of the HLA allele nomenclature (e.g. A*01:01:01:01) that was adopted in April 2010. In addition, the antt and uncl translate specific HLA allele-name changes (e.g. DPB1*0502 is translated to DPB1*104:01), as well as changes to the locus prefix for HLA-C (i.e. Cw* is translated to C*). The antt and uncl will also translate allele names that have been truncated to two, four, or six digits, as well as ambiguous allele strings. The antt is a locally installed and run application, while uncl is a web-based tool that requires only an Internet connection and a modern browser. The antt accepts a variety of HLA data-presentation and allele-name formats. In addition, the antt can translate using user-defined conversion settings (e.g. the names of alleles that encode identical peptide binding domains can be translated to a common 'P-code'), and can serve as a preliminary data-sanity tool. The antt is available for download, and uncl for use, at www.igdawg.org/software.

Association analysis of SNPs in the IL4R locus with type I diabetes.

The Type I Diabetes Genetics Consortium (T1DGC) has collected thousands of multiplex and simplex families with type I diabetes (T1D) with the goal of identifying genes involved in T1D susceptibility. These families have all been genotyped for the HLA class I and class II loci and a subset of samples has been typed for an major histocompatibility complex (MHC) single-nucleotide polymorphism (SNP) panel. In addition, the T1DGC has genotyped SNPs in candidate genes to evaluate earlier reported T1D associations. Individual SNPs and SNP haplotypes in IL4R, which encodes the alpha-chain of the IL4 and IL13 receptors, have been associated with T1D in some reports, but not in others. In this study, 38 SNPs in IL4R were genotyped using the Sequenom iPLEX Gold MassARRAY technology in 2042 multiplex families from nine cohorts. Association analyses (transmission-disequilibrium test and parental-disequilibrium test) were performed on individual SNPs and on three-SNP haplotypes. Analyses were also stratified on the high-risk HLA DR3/DR4-DQB1*0302 genotype. A modest T1D association in HBDI families (n=282) was confirmed in this larger collection of HBDI families (n=424). The variant alleles at the non-synonymous SNPs (rs1805011 (E400A), rs1805012 (C431R), and rs1801275 (Q576R)), which are in strong linkage disequilibrium, were negatively associated with T1D risk. These SNPs were more associated with T1D among non-DR3/DR4-DQB1*0302 genotypes than DR3/DR4-DQB1*0302 genotypes. This association was stronger, both in terms of odds ratio and P-values, than the initial report of the smaller collection of HBDI families. However, the IL4R SNPs and the three-SNP haplotype containing the variant alleles were not associated with T1D in the total data. Thus, in the overall families, these results do not show evidence for an association of SNPs in IL4R with T1D.

HLA-A, -B, -C, and -DRB1 allele and haplotype frequencies distinguish Eastern European Americans from the general European American population.

Sequence-based typing was used to identify human leukocyte antigen (HLA)-A, -B, -C, and -DRB1 alleles from 558 consecutively recruited US volunteers with Eastern European ancestry for an unrelated hematopoietic stem cell registry. Four of 31 HLA-A alleles, 29 HLA-C alleles, 59 HLA-B alleles, and 42 HLA-DRB1 alleles identified (A*0325, B*440204, Cw*0332, and *0732N) are novel. The HLA-A*02010101g allele was observed at a frequency of 0.28. Two-, three-, and four-locus haplotypes were estimated using the expectation-maximization algorithm. The highest frequency extended haplotypes (A*010101g-Cw*070101g-B*0801g-DRB1*0301 and A*03010101g-Cw*0702-B*0702-DRB1*1501) were observed at frequencies of 0.04 and 0.03, respectively. Linkage disequilibrium values (Dij') of the constituent two-locus haplotypes were highly significant for both extended haplotypes (P values were less than 8 x 10(-10)) but were consistently higher for the more frequent haplotype. Balancing selection was inferred to be acting on all the four loci, with the strongest evidence of balancing selection observed for the HLA-C locus. Comparisons of the A-C-B haplotypes and DRB1 frequencies in this population with those for African, European, and western Asian populations showed high degrees of identity with Czech, Polish, and Slovenian populations and significant differences from the general European American population.

14th International HLA and Immunogenetics Workshop: report of progress in methodology, data collection, and analyses.

The Biostatistics Component of the 13th International Histocompatibility Workshop (IHWS) developed the PyPop (Python for Population Genomics) software framework for high-throughput analysis and quality control (QC) assessments of highly polymorphic genotype data. Since its initial release, the software has had several new analysis modules added to it. These additions, combined with improved data filtering and QC modules, facilitate analyses of data at different levels (allele, haplotype, amino acid sequence, and nucleotide sequence). Since the 13th IHWS, much of the human leukocyte antigen (HLA) data from the workshop, QCed via PyPop and other methods, have been made publicly available through the Major Histocompatibility Complex database web site at the National Center for Biotechnology Information (http://ncbi.nih.gov/mhc/). The Anthropology/Human Genetic Diversity component (AHGDC) data have been used in a variety of studies. Prugnolle et al. used this data to corroborate a model of pathogen-driven selection as a factor related to high levels of diversity at HLA loci. Using a comparative genomics approach contrasting results for HLA and non-HLA markers, Meyer et al. analyzed a subset of the 13th IHWS AHGDC data and showed that HLA loci show detectable signs of both natural selection and the demographic history of populations.

Population samples and genotyping technology.

The 14th International HLA (human leukocyte antigen) Immunogenetics Workshop (14th-IHIWS) Biostatistics and Anthropology/Human Genetic Diversity project continues the population sampling, genotype data generation, and biostatistic analyses of the 13th International Histocompatibility Workshop Anthropology/Human Genetic Diversity Component, with the overall goal of further characterizing global HLA allele and haplotype diversity and better describing the relationships between major histocompatibility complex diversity, geography, linguistics, and population history. Since the 13th Workshop, new investigators have and continue to be recruited to the project and new high-resolution class I and class II genotype data are being generated for 112 population samples from around the world.

HLA-A, -B, -C, -DRB1 allele and haplotype frequencies in an African American population.

Sequence-based typing was used to identify human leukocyte antigen (HLA)-A, -B, -C, and -DRB1 alleles from 564 consecutively recruited African American volunteers for an unrelated hematopoietic stem cell registry. The number of known alleles identified at each locus was 42 for HLA-A, HLA-B 67, HLA-C 33, and HLA-DRB1 44. Six novel alleles (A*260104, A*7411, Cw*0813, Cw*1608, Cw*1704, and DRB1*130502) not observed in the initial sequence-specific oligonucleotide probe testing were characterized. The action of balancing selection, shaping more 'even' than expected allele frequency distributions, was inferred for all four loci and significantly so for the HLA-A and DRB1 loci. Two-, three-, and four-locus haplotypes were estimated using the expectation maximization algorithm. Comparisons with other populations from Africa and Europe suggest that the degree of European admixture in the African American population described here is lower than that in other African American populations previously reported, although HLA-A:B haplotype frequencies similar to those in previous studies of African American individuals were also noted.

Genetic variation among the Golla pastoral caste subdivisions of Andhra Pradesh, India, according to the HLA system.

The HLA allele frequency distributions have been characterized for the HLA class I and class II loci of the Golla pastoral caste, from Southeast India, subdivided into the subcastes (Puja, Punugu, Kurava, Pokanati, Karnam, and Doddi). Genetic distances, neighbor-joining, correspondence, and haplotype analyses all indicate that the subcastes exhibit a high haplotype variability and that their genetic substratum may be the result of European-Middle East/Asian admixture with the autochthonous populations. The Karnam subcaste seems to be the one that has undergone a higher degree of admixture, when compared with the other subcastes. The Golla speak an old Indian Dravidian language and should theoretically represent the basic Indian substratum that existed before the postulated "Aryan" invasion.

Genetic variability and linkage disequilibrium within the HLA-DP region: analysis of 15 different populations.

In order to understand the forces governing the evolution of the genetic diversity in the HLA-DP molecule, polymerase chain reaction (PCR)-based methods were used to characterize genetic variation at the DPA1 and DPB1 loci encoding this heterodimer on 2,807 chromosomes from 15 different populations including individuals of African, Asian, Amerindian, Indian and European origin. These ethnically diverse samples represent a variety of population substructures and include small, isolated populations as well as larger, presumably admixed populations. Ten DPA1 and 39 DPB1 alleles were identified and observed on 87 distinct DP haplotypes, 34 of which were found to be in significant positive linkage disequilibrium in at least one population. Some haplotypes were found in all ethnic groups while others were confined to a single ethnic group or population. Strong positive global linkage disequilibrium (Wn) between DPA1 and DPB1 was present in all 15 populations. The African populations displayed the lowest values of Wn whereas the Amerindian populations displayed near absolute disequilibrium. Analysis of the distribution of haplotypes using the normalized deviate of the Ewens-Watterson homozygosity statistic, F, suggests that DP haplotypes encoding the functional heterodimer are subject to much lower degrees of balancing selection than other loci within the HLA region. Finally, neighbor joining tree analyses demonstrate the power of haplotype diversity for inferring the relationships between the different populations.

Diversity is demonstrated in class I HLA-A and HLA-B alleles in Cameroon, Africa: description of HLA-A*03012, *2612, *3006 and HLA-B*1403, *4016, *4703.

To examine the genetic diversity in west Africa, class I HLA-A and HLA-B alleles of 92 unrelated individuals from two areas in the Cameroon, the capital Yaoundé and the village of Etoa, were identified by direct automated DNA sequencing of exons 2 and 3 of the HLA-B locus alleles and sequence-specific oligonucleotide probe (SSOP) and/or sequencing of the HLA-A locus alleles. HLA-A*2301 (18.7%), A*2902 (10.4%), B*5301 (10.9%), and B*5802 (10.9%) were the most frequently detected alleles, present in at least 10% of the population. A total of 30 HLA-A locus and 33 HLA-B locus alleles, including six novel alleles, were detected. The novel alleles were HLA-A*03012, A*2612, A*3006 and HLA-B*1403, B*4016, and B*4703. HLA-B*4703 contains a novel amino acid sequence that is a combination of the first 5 amino acids of the Bw6 epitope and the last 2 residues of the Bw4 epitope. The addition of 6 alleles to the ever-expanding number of known class I HLA alleles supports our hypothesis that extensive genetic diversity, including previously undescribed alleles, would be observed in this African population. In the Yaoundé population, the allele frequency distribution at the HLA-A locus is consistent with distributions indicative of balancing selection. Extensive HLA-A-B haplotypes were observed in this population suggesting that only a fraction of the Cameroon HLA-A-B haplotype diversity has been observed.

Evolution of Pacific/Asian populations inferred from HLA class II allele frequency distributions.

The allele frequency distributions for the HLA class II loci, DRB1, DQB1 and DPB1, in eight Pacific/Asian populations: Hawaiian, Samoan, Malay, Papua New Guinea (PNG) Highlands, and two Indonesian and PNG Lowland groups, were determined using high-resolution polymerase chain reaction/sequence-specific oligonucleotide probe (PCR/SSOP) typing methods. The allele frequency distributions for the HLA-DRB1 locus were determined for a third Indonesian population as well as for an additional Filipino population. DRB1 alleles in the DR2 serogroup (or allelic lineage) are very common in this region; in some populations, more than 50% of the alleles belong to this serogroup. The DRB1*1502 allele is frequent in nine of the ten populations studied, reaching a frequency of 0.48 in one Indonesian population and among Filipinos. Extensive DR-DQ haplotype diversity was detected in these populations. Seven different DR2-DQB1 haplotypes were observed in the Indonesian and PNG Lowland populations, eight in the PNG Highlands and ten in Malays and Filipinos. The DRB1*0410 allele, commonly observed in Australia, is observed in the PNG Highlands at a low frequency (f=0.03) and is absent in the other populations. Two additional DRB1 alleles commonly observed in Australia, DRB1*0405 and *1407, are also observed in the PNG Highlands at high frequencies (f=0.132 and 0.126), while they are rare in the PNG Lowlands (f=0.039 and 0.013). These alleles are generally rare or absent in the other populations. The DPB1*0501 allele, common in Chinese and Japanese populations, is most frequent in the Samoan, Hawaiian, Indonesian, and Malay populations, and the *0401 allele is the most frequent DPB1 allele in the PNG Lowlands. Both of these alleles have the same very high frequency (f=0.34) in the PNG Highlands. Analyses of homozygosity (the Ewens-Watterson F statistic) in these and other populations indicate that, while most allele frequency distributions are consistent with balancing selection, values of F for the Indonesian and Javan populations may reflect positive directional selection. Phylogenetic trees constructed using the allele frequencies at the DRB1 locus of the populations reported here, as well as those for additional Pacific, Asian, and Australian populations, indicate that the PNG Highland population is more closely related to Australian populations than to PNG Lowland populations, while the PNG Lowlands are more closely related to other Melanesian populations.

HLA class I allele distributions in six Pacific/Asian populations: evidence of selection at the HLA-A locus.

The distributions of HLA-A alleles in six Pacific/Asian populations, Malay, Papua New Guinea (PNG) Highlands, two Indonesian groups, and two PNG Lowland groups, as well as the distribution of the HLA-B alleles in the PNG Highlands population, were determined using polymerase chain reaction (PCR) immobilized sequence-specific oligonucleotide (SSO) probe typing methods. The allele frequency distributions at the HLA class II loci, DRB1, DQB1 and DPB1 were also determined by PCR/SSO methods in an additional study of the same populations. In most of these populations, the HLA-A*2402 allele was the most frequent, attaining a frequency of 0.78 in the PNG Highlands. A*1101 was the next most frequent allele, followed in frequency by the *3401 allele. The HLA-B*1506, *4001, *5601 and *5602 alleles comprised 73% of the allele diversity at the B-locus in the PNG Highlands. Two previously unreported HLA-A alleles were identified in Indonesians and Malays, based on novel probe reactivity patterns. Cloning and sequencing identified these as A*1104 and *2410. Sequence comparisons show that these new alleles differ at codon 187 from their putative parental alleles (*1101 and *2403) by dinucleotide changes in the first two codon positions. These changes involve a Thr to Arg (CG to AC) and an Arg to Thr substitution (AC to CG) at position 187; residues at this position participate in pocket A of the peptide binding groove. Comparison of the HLA-A allele frequency distributions indicate that Malays are the most diverse (heterozygosity (h)=0.88) and the PNG Highlanders are, by far, the least diverse (h=0.37) of the groups studied. However, the diversity of B-locus alleles in the PNG Highlanders (h=0.91) was greater than that observed at the A-locus of any of the populations reported here. The remarkably high allele frequency of A*2402 in the PNG Highlands could reflect founder effects and population bottlenecks, genetic drift, or positive directional selection. The distribution of the HLA-B locus alleles and class II alleles, as well as mtDNA sequence data in the PNG Highlands indicates a reasonably high level of diversity at other loci, arguing that the high frequency of A*2402 cannot be attributed entirely to founder effects, bottlenecks, or drift and suggests the operation of positive selection for the A*2402 allele in this population.

Tracing the origin of HLA-DRB1 alleles by microsatellite polymorphism.

We analyzed the origin of allelic diversity at the class II HLA-DRB1 locus, using a complex microsatellite located in intron 2, close to the polymorphic second exon. A phylogenetic analysis of human, gorilla, and chimpanzee DRB1 sequences indicated that the structure of the microsatellite has evolved, primarily by point mutations, from a putative ancestral (GT)x(GA)y-complex-dinucleotide repeat. In all contemporary DRB1 allelic lineages, with the exception of the human *04 and the gorilla *08 lineages, the (GA)y repeat is interrupted, often by a G-->C substitution. In general, the length of the 3' (GA)y repeat correlates with the allelic lineage and thus evolves more slowly than a middle (GA)z repeat, whose length correlates with specific alleles within the lineage. Comparison of the microsatellite sequence from 30 human DRB1 alleles showed the longer 5' (GT)x to be more variable than the shorter middle (GA)z and 3' (GA)y repeats. Analysis of multiple samples with the same exon sequence, derived from different continents, showed that the 5' (GT)x repeat evolves more rapidly than the middle (GA)z and the 3' (GA)y repeats, which is consistent with findings of a higher mutation rate for longer tracts. The microsatellite-repeat-length variation was used to trace the origin of new DRB1 alleles, such as the new *08 alleles found in the Cayapa people of Ecuador and the Ticuna people of Brazil.