Polymorphism, recombination, and linkage disequilibrium within the HLA class II region.

Thirty-nine CEPH (Centre d'Etude du Polymorphisme Humain) families, comprised of 502 individuals, have been typed for the HLA class II genes DRB1, DQA1, DQB1, and DPB1 using nonradioactive sequence-specific oligonucleotide probes to analyze polymerase chain reaction amplified DNA. This population, which consists of 266 independent chromosomes, contains 27 DRB1, 7 DQA1, 12 DQB1, and 17 DPB1 alleles. Analysis of the distribution of allele frequencies using the homozygosity statistic, which gives an indication of past selection pressures, suggests that balancing selection has acted on the DRB1, DQA1, and DQB1 loci. The distribution of DPB1 alleles, however, suggests a different evolutionary past. Family data permits the estimation of recombination rates and the unambiguous assignment of haplotypes. No recombinants were found between DRB1, DQA1, and DQB1; however, recombinants were detected between DQB1 and DPB1, resulting in an estimated recombination fraction of greater than or equal to 0.008 +/- 0.004. Only 33 distinct DRB1-DQA1-DQB1 haplotypes were found in this population which illustrates the extreme nonrandom haplotypic association of alleles at these loci. A few of these haplotypes are unusual (previously unreported) for a Caucasian population and most likely result from past recombination events between the DR and DQ subregions. Examination of disequilibrium across the HLA region using these data and the available serologic HLA-A and HLA-B types of these samples shows that global disequilibrium between these loci declines with the recombination fraction, approaching statistic nonsignificance at the most distant interval, HLA-A to HLA-DP.DR-DQ haplotypes in linkage disequilibrium with DPB1 and B are noted and, finally, the evolutionary origin of certain class II haplotypes is addressed.     

Analysis of HLA class II haplotypes in the Cayapa Indians of Ecuador: a novel DRB1 allele reveals evidence for convergent evolution and balancing selection at position 86.

PCR amplification, oligonucleotide probe typing, and sequencing were used to analyze the HLA class II loci (DRB1, DQA1, DQB1, and DPB1) of an isolated South Amerindian tribe. Here we report HLA class II variation, including the identification of a new DRB1 allele, several novel DR/DQ haplotypes, and an unusual distribution of DPB1 alleles, among the Cayapa Indians (N = 100) of Ecuador. A general reduction of HLA class II allelic variation in the Cayapa is consistent with a population bottle-neck during the colonization of the Americas. The new Cayapa DRB1 allele, DRB1*08042, which arose by a G-->T point mutation in the parental DRB1*0802, contains a novel Val codon (GTT) at position 86. The generation of DRB1*08042 (Val-86) from DRB1*0802 (Gly-86) in the Cayapa, by a different mechanism than the (GT-->TG) change in the creation of DRB1*08041 (Val-86) from DRB1*0802 in Africa, implicates selection in the convergent evolution of position 86 DR beta variants. The DRB1*08042 allele has not been found in > 1,800 Amerindian haplotypes and thus presumably arose after the Cayapa separated from other South American Amerindians. Selection pressure for increased haplotype diversity can be inferred in the generation and maintenance of three new DRB1*08042 haplotypes and several novel DR/DQ haplotypes in this population. The DPB1 allelic distribution in the Cayapa is also extraordinary, with two alleles, DPB1*1401, a very rare allele in North American Amerindian populations, and DPB1*0402, the most common Amerindian DPB1 allele, constituting 89% of the Cayapa DPB1.(ABSTRACT TRUNCATED AT 250 WORDS)     

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.     

PCR/oligonucleotide probe typing of HLA class II alleles in a Filipino population reveals an unusual distribution of HLA haplotypes.

We have analyzed the distribution of HLA class II alleles and haplotypes in a Filipino population by PCR amplification of the DRB1, DQB1, and DPB1 second-exon sequences from buccal swabs obtained from 124 family members and 53 unrelated individuals. The amplified DNA was typed by using nonradioactive sequence-specific oligonucleotide probes. Twenty-two different DRB1 alleles, including the novel Filipino *1105, and 46 different DRB1/DQB1 haplotypes, including the unusual DRB1*0405-DQB1*0503, were identified. An unusually high frequency (f = .383) of DPB1*0101, a rare allele in other Asian populations, was also observed. In addition, an unusual distribution of DRB1 alleles and haplotypes was seen in this population, with DR2 (f = .415) and DRB1*1502-DQB1*0502 (f = .233) present at high frequencies. This distribution of DRB1 alleles differs from the typical HLA population distribution, in which the allele frequencies are more evenly balanced. The distribution of HLA class II alleles and haplotypes in this Filipino population is different from that of other Asian and Pacific groups: of those populations studied to date; the Indonesian population is the most similar. DRB1*1502-DQB1*0502 was in strong linkage disequilibrium (D'' = .41) with DPB1*0101 (f = .126, for the extended haplotype), which is consistent with selection for this DR, DQ, DP haplotype being responsible for the high frequency of these three class II alleles in this population.     

Discordant patterns of linkage disequilibrium of the peptide-transporter loci within the HLA class II region.

Disequilibrium between genetic markers is expected to decline monotonically with recombinational map distance. We present evidence from the HLA class II region that seems to violate this principle. Pairwise disequilibrium values were calculated from six loci ranging in physical separation from 15 kb to 550 kb. The histocompatibility loci DRB1, DQA1, and DQB1, located on the distal end of the class II region, behave as a single evolutionary unit within which extremely high linkage disequilibrium exists. Lower but still significant levels of disequilibrium are present between these loci and DPB1, located at the proximal edge of the HLA complex. The peptide-transporter loci TAP1 and TAP2, located in the intervening region, reveal no disequilibrium with each other and low or negligible disequilibrium with the flanking loci. The action of two genetic process is required to account for this phenomenon: a recombinational hotspot operating between TAP1 and TAP2, to eliminate disequilibrium between these loci, and at the same time selection operating on particular combinations of alleles across the DR-DP region, to create disequilibrium in the favored haplotypes. The forces producing the patterns of disequilibrium observed here have implications for the mapping of train loci and disease genes: markers of TAP1, for example, would give a false impression as to the influence of DPB1 on a trait known to be associated with DQB1.     

Polymorphism at the ovine major histocompatibility complex class II loci

Southern hybridization analysis of the ovine major histocompatibility complex (MHC) ( MhcOvar ) class II region, using sheep-specific probes for the DQA1, DQA2, DQB and DRA loci, has revealed extensive polymorphism. DQA1 and DQAP had eight and 16 alleles respectively, DQB had six and DRA had three alleles. Little information was derived from the DRB locus owing to extensive cross-hybridization between the DRB probe and the DQB locus. Differences in allele frequency between breeds were revealed. At the DQA1 locus a null allele (DQA1-N) was observed with a frequency of between 27% and 45%, making this the most common DQA1 allele in all breeds examined. The frequency of DQA1-N homozygotes was between 11% and 18%, raising questions as to the functional significance of the DQA1 gene. Linkage analysis between the DQA1, DQA2, DQB and DRA loci did not reveal any recombination.     

The severe phenotype of females with tiny ring X chromosomes is associated with inability of these chromosomes to undergo X inactivation.

Molecular typing of HLA class II loci has been performed on a sample of 196 patients with Hodgkin lymphoma. Division of patients into two histological categories--nodular sclerosing Hodgkin disease versus all other types--shows significant overall association of the nodular sclerosing group with the HLA class II region. Haplotypes and alleles defined for the four loci typed--DRB1, DQA1, DQB1, and DPB1--were present in both excess and deficit in the nodular sclerosing sample. Some of the effects are attributable to particular DRB1 and DQB1 alleles, while other effects are best explained by haplotypes marking the entire class II region. The latter effects might be due to variation in additional, as-yet-unexamined loci in the class II region or to particular combinations of alleles from two or more loci. These data also explain why earlier studies showed HLA linkage but not association, and they substantiate the specific involvement of the immune system in certain neoplastic diseases.     

Major histocompatibility complex variation at three class II loci in the northern elephant seal

Northern elephant seals were hunted to near extinction in the 19th century, yet have recovered remarkably and now number around 175,000. We surveyed 110 seals for single-strand conformation polymorphism (SSCP) and sequence variation at three major histocompatibility (MHC) class II loci (DQA, DQB and DRB) to evaluate the genetic consequences of the population bottleneck at these loci vs. other well-studied genes. We found very few alleles at each MHC locus, significant variation among breeding sites for the DQA locus, and linkage disequilibrium between the DQB and DRB loci. Northern elephant seals are evidently inbred, although there is as yet no evidence of correlative reductions in fitness.     

The role of HLA class II genes in insulin-dependent diabetes mellitus: molecular analysis of 180 Caucasian,multiplex families.

We report here our analysis of HLA class II alleles in 180 Caucasian nuclear families with at least two children with insulin-dependent diabetes mellitus (IDDM). DRB1, DQA1, DQB1, and DPB1 genotypes were determined with PCR/sequence-specific oligonucleotide probe typing methods. The data allowed unambiguous determination of four-locus haplotypes in all but three of the families. Consistent with other studies, our data indicate an increase in DR3/DR4, DR3/DR3, and DR4/DR4 genotypes in patients compared to controls. In addition, we found an increase in DR1/DR4, DR1/DR3, and DR4/DR8 genotypes. While the frequency of DQB1*0302 on DR4 haplotypes is dramatically increased in DR3/DR4 patients, DR4 haplotypes in DR1/DR4 patients exhibit frequencies of DQB1*0302 and DQB1*0301 more closely resembling those in control populations. The protective effect of DR2 is evident in this data set and is limited to the common DRB1*1501-DQB1*0602 haplotype. Most DR2+ patients carry the less common DR2 haplotype DRB1*1601-DQB1*0502, which is not decreased in patients relative to controls. DPB1 also appears to play a role in disease susceptibility. DPB1*0301 is increased in patients (P < .001) and may contribute to the disease risk of a number of different DR-DQ haplotypes. DPB1*0101, found almost exclusively on DR3 haplotypes in patients, is slightly increased, and maternal transmissions of DRB1*0301-DPB1*0101 haplotypes to affected children occur twice as frequently as do paternal transmissions. Transmissions of DR3 haplotypes carrying other DPB1 alleles occur at approximately equal maternal and paternal frequencies. The complex, multigenic nature of HLA class II-associated IDDM susceptibility is evident from these data.     

Polymorphism of the HLA class II loci in Siberian populations

The populations that colonized Siberia diverged from one another in the Paleolithic and evolved in isolation until today. These populations are therefore a rich source of information about the conditions under which the initial divergence of modern humans occurred. In the present study we used the HLA system, first, to investigate the evolution of the human major histocompatibility complex (MHC) itself, and second, to reveal the relationships among Siberian populations. We determined allelic frequencies at five HLA class II loci (DRB1, DQA1, DQB1, DPA1, and DPB1) in seven Siberian populations (Ket, Evenk, Koryak, Chukchi, Nivkh, Udege, and Siberian Eskimo) by the combination of single-stranded conformational polymorphism and DNA sequencing analysis. We then used the gene frequency data to deduce the HLA class II haplotypes and their frequencies. Despite high polymorphism at four of the five loci, no new alleles could be detected. This finding is consistent with a conserved evolution of human class II MHC genes. We found a high number of HLA class II haplotypes in Siberian populations. More haplotypes have been found in Siberia than in any other population. Some of the haplotypes are shared with non-Siberian populations, but most of them are new, and some represent “forbidden” combinations of DQA1 and DQB1 alleles. We suggest that a set of “public” haplotypes was brought to Siberia with the colonizers but that most of the new haplotypes were generated in Siberia by recombination and are part of a haplotype pool that is turning over rapidly. The allelic frequencies at the DRB1 locus divide the Siberian populations into eastern and central Siberian branches; only the former shows a clear genealogical relationship to Amerinds. Received: 18 August 1997 / Accepted: 6 October 1997     

Visualizing human leukocyte antigen class II risk haplotypes in human systemic lupus erythematosus

Human leukocyte antigen (HLA) class I and class II alleles are implicated as genetic risk factors for many autoimmune diseases. However, the role of the HLA loci in human systemic lupus erythematosus (SLE) remains unclear. Using a dense map of polymorphic microsatellites across the HLA region in a large collection of families with SLE, we identified three distinct haplotypes that encompassed the class II region and exhibited transmission distortion. DRB1 and DQB1 typing of founders showed that the three haplotypes contained DRB1*1501/ DQB1*0602, DRB1*0801/ DQB1*0402, and DRB1*0301/DQB1*0201 alleles, respectively. By visualizing ancestral recombinants, we narrowed the disease-associated haplotypes containing DRB1*1501 and DRB1*0801 to an approximately 500-kb region. We conclude that HLA class II haplotypes containing DRB1 and DQB1 alleles are strong risk factors for human SLE.     

Trans-species polymorphism and evidence of selection on class II MHC loci in tuco-tucos (Rodentia: Ctenomyidae)

Balancing selection acting over the evolutionary history of a lineage can result in the retention of alleles among species for longer than expected under neutral evolution. The associated pattern of trans-species polymorphism, in which similar or even identical alleles are shared among species, is often used to infer that balancing selection has occurred. The genes of the major histocompatibility complex (MHC) are thought to be subject to balancing selection that maintains alleles associated with response to specific pathogens. To explore the role of balancing selection in shaping MHC diversity in ctenomyid rodents, we examined allelic variability at the class II DRB and DQA loci in 18 species in the genus Ctenomys. Previous studies of four of these species had revealed significant within-population evidence of positive selection on MHC loci. The current study expands upon these analyses to (1) evaluate among-species evidence of positive selection and (2) explore the potential for balancing selection on MHC genes. Interspecific nucleotide sequence variation revealed significant evidence of positive selection on the DRB and DQA loci. At the same time, comparisons of phylogenetic trees for these MHC loci with a putative species tree based on mitochondrial sequence data revealed multiple examples of trans-specific polymorphism, including sharing of identical DRB and DQA alleles among distantly related species of Ctenomys. These findings suggest that MHC genes in these animals have historically been subject to balancing selection and yield new insights into the complex suite of forces shaping MHC diversity in free-living vertebrates.     

Genetic drift vs. natural selection in a long-term small isolated population: major histocompatibility complex class II variation in the Gulf of California endemic porpoise (Phocoena sinus)

Although many studies confirm long-term small isolated populations (e.g. island endemics) commonly sustain low neutral genetic variation as a result of genetic drift, it is less clear how selection on adaptive or detrimental genes interplay with random forces. We investigated sequence variation at two major histocompatibility complex (Mhc) class II loci on a porpoise endemic to the upper Gulf of California, México (Phocoena sinus, or vaquita). Its unique declining population is estimated around 500 individuals. Single-strand conformation polymorphism analysis revealed one putative functional allele fixed at the locus DQB (n = 25). At the DRB locus, we found two presumed functional alleles (n = 29), differing by a single nonsynonymous nucleotide substitution that could increase the stability at the dimer interface of alphabeta-heterodimers on heterozygous individuals. Identical trans-specific DQB1 and DRB1 alleles were identified between P. sinus and its closest relative, the Burmeister's porpoise (Phocoena spinipinnis). Comparison with studies on four island endemic mammals suggests fixation of one allele, due to genetic drift, commonly occurs at the DQA or DQB loci (effectively neutral). Similarly, deleterious alleles of small effect are also effectively neutral and can become fixed; a high frequency of anatomical malformations on vaquita gave empirical support to this prediction. In contrast, retention of low but functional polymorphism at the DRB locus was consistent with higher selection intensity. These observations indicated natural selection could maintain (and likely also purge) some crucial alleles even in the face of strong and prolonged genetic drift and inbreeding, suggesting long-term small populations should display low inbreeding depression. Low levels of Mhc variation warn about a high susceptibility to novel pathogens and diseases in vaquita.     

Human leukocyte antigen and interleukin 2, 10 and 12p40 cytokine responses to measles: is there evidence of the HLA effect?

HLA class I and class II associations were examined in relation to measles virus-specific cytokine responses in 339 healthy children who had received two doses of live attenuated measles vaccine. Multivariate linear regression modeling analysis revealed suggestions of associations between the expression of DPA1*0201 (p=0.03) and DPA1*0202 (p=0.09) alleles and interleukin-2 (IL-2) cytokine production (global p-value 0.06). Importantly, cytokine production and DQB1 allele associations (global p-value 0.04) revealed that the alleles with the strongest association with IL-10 secretion were DQB1*0302 (p=0.02), DQB1*0303 (p=0.07) and DQB1*0502 (p=0.06). Measles-specific IL-10 secretion associations approached significance with DRB1 and DQA1 loci (both global p-values 0.08). Specifically, suggestive associations were found between DRB1*0701 (p=0.07), DRB1*1103 (p=0.06), DRB1*1302 (p=0.08), DRB1*1303 (p=0.06), DQA1*0101 (p=0.08), and DQA1*0201 (p=0.04) alleles and measles-induced IL-10 secretion. Further, suggestive association was observed between specific DQA1*0505 (p=0.002) alleles and measles-specific IL-12p40 secretion (global p-value 0.09) indicating that cytokine responses to measles antigens are predominantly influenced by HLA class II genes. We found no associations between any of the alleles of HLA A, B, and Cw loci and cytokine secretion. These novel findings suggest that HLA class II genes may influence the level of cytokine production in the adaptive immune responses to measles vaccine.     

Frequencies of the <Emphasis Type="Italic">DRB1</Emphasis>, <Emphasis Type="Italic">DQA1</Emphasis>, <Emphasis Type="Italic">DQB1</Emphasis> and <Emphasis Type="Italic">TNFA</Emphasis> alleles in immigrant population of west Siberia

Based on population analysis of the DRB1, DQA1, DQB1 and TNFA allele frequency distribution patterns, regional features of immunogenetic structure of the population of West Siberia were investigated. Statistically significant linkage disequilibrium within the HLA class II region, as well as between the TNFA and DRB1, DQA1, and DQB1 was demonstrated. Population frequency distribution patterns of two- and multilocus haplotypes were examined.     

MHC class II genes in European wolves: a comparison with dogs

The genome of the grey wolf, one of the most widely distributed land mammal species, has been subjected to both stochastic factors, including biogeographical subdivision and population fragmentation, and strong selection during the domestication of the dog. To explore the effects of drift and selection on the partitioning of MHC variation in the diversification of species, we present nine DQA, 10 DQB, and 17 DRB1 sequences of the second exon for European wolves and compare them with sequences of North American wolves and dogs. The relatively large number of class II alleles present in both European and North American wolves attests to their large historical population sizes, yet there are few alleles shared between these regions at DQB and DRB1. Similarly, the dog has an extensive array of class II MHC alleles, a consequence of a genetically diverse origin, but allelic overlap with wolves only at DQA. Although we might expect a progression from shared alleles to shared allelic lineages during differentiation, the partitioning of diversity between wolves and dogs at DQB and DRB1 differs from that at DQA. Furthermore, an extensive region of nucleotide sequence shared between DRB1 and DQB alleles and a shared motif suggests intergenic recombination may have contributed to MHC diversity in the Canidae.     

Origin and affinities of indigenous Siberian populations as revealed by HLA class II gene frequencies

Gene frequencies of eight Siberian populations (Mansi, Tuva, Todja, Tofalar, Buryat, Okhotsk Evenki, Ulchi, and Negidal) were determined for the three most polymorphic HLA class II loci ( DRB1, DQA1, and DQB1) by a combination of single-stranded conformational polymorphism typing and DNA sequencing. The number of alleles per population ranged from 16 to 25, from seven to eight, and from nine to 14 for the DRB1, DQA1, and DQB1 loci, respectively. The alleles at the three loci occurred in 66 different combinations (haplotypes), most of which appeared to be of ancient origin, but some may have arisen within the Siberian populations. Phylogenetic analysis of the frequency data suggests that the HLA genes of Asian and indigenous American populations stem from a single pool distinct from the gene pools of European and African populations. The Asian populations separate into two clusters, one of which encompasses nearly all the Siberian populations and all the indigenous American populations tested, while the other consists of Central, Eastern, and Southeastern Asian populations. The position of the Tuva people appears to be near the node from which the two clusters diverge. The divergence time of the two clusters is estimated to be 21,000-24,000 years BP. Three different branches of the native Siberian peoples seem to have contributed founders for the indigenous American ethnic groups.     

Two divergent routes of evolution gave rise to the DRw13 haplotypes

The HLA class II genes and haplotypes have evolved over a long period of evolutionary time by mechanisms such as gene conversion, reciprocal recombination and point mutation. The extent of the diversity generated is most clearly evident in an analysis of the HLA class II alleles present within DRw13 haplotypes. This study uses cDNA sequencing to examine the first domains of DRB1, DRB3, DQA1, and DQB1 alleles from several American black individuals expressing seven different DRw13 haplotypes, five with undefined HLA-D specificities (i.e., not Dw18 or Dw19). Two new DRw13 alleles described in this study are the first examples of convergent evolution of DR alleles in which gene conversion has apparently combined segments of DRB1 alleles encoding DRw11 and DRw8 to generate two new DRB1 alleles, DRB1*1303 and DRB1*1304, that encode molecules bearing serologic determinants of a third allele, DRw13. These new DRw13 alleles are found embedded in haplotypes of DRw11 origin distinct from haplotypes encoding previously identified DRw13 alleles, DRB1*1301 and DRB1*1302. These data suggest that two evolutionary pathways may have given rise to two subgroups of alleles encoding molecules that share DRw13 serologic determinants yet which possess different structural and, likely, functional motifs. Reciprocal gene recombination events resulting in different DR, DRw52 and DQ allele combinations also appear to have played a crucial role in augmenting the level of diversity found in DRw13 haplotypes. Recombination has resulted in the association of one of the new DRw13 alleles with a DQw2 allele normally found associated with DR7 and the association of the DRw52c-associated DRw13 allele (DRB1*1302) with three different DQw1 alleles. The seven DRw13 haplotypes that have resulted from the effect of recombination on haplotypes formed by the two pathways of DRw13 allelic diversification have resulted in different repertoires of class II molecules and, most likely, different immune response profiles in individuals with these haplotypes.