Allelic polymorphism in introns 1 and 2 of the HLA-DQA1 gene

Human leukocyte antigen (HLA) class II antigens are highly polymorphic membrane glycoproteins, encoded by the A and B genes of DR, DQ, and DP. The polymorphism is mainly located in exon 2, with the exception of DQA1. Of the 27 DQA1 alleles presently known, 18 cannot be identified on the basis of exon 2 alone, but need additional information from the other exons. DQA1 has been reported to be the most ancient class II gene. For evolutionary comparison and to assess the degree of polymorphism outside the exons, the sequences of introns 1 and 2 were determined from 30 different cell lines, encompassing 15 different DQA1 alleles. The sequences revealed major nucleotide differences between the different lineages, whereas within each lineage few differences were present. Phylogenetic analysis of intron and exon sequences confirmed this lineage specificity. Altogether, the present data indicate that the HLA-DQA1 lineages represent ancient entities. The observed variation of the introns in alleles with identical exon sequences implicates conservative selection of the exons within a given lineage. Intron sequences may provide the means to set up an accurate typing system.     

HLA-DQA1 introns 2 and 3 sequencing: DQA1 sequencing-based typing and characterization of a highly polymorphic microsatellite at intron 3 of DQA1*0505

DQA1 class II gene encodes the alpha-chain of the human leukocyte antigen (HLA)-DQ heterodimer. Sequencing-based typing (SBT) for HLA genes is the most powerful methodology described. However, most of the SBT procedures reported for HLA class II genes are not able to define complete exon 2 region. For that purpose, we have characterized introns 2 and 3 from most DQA1 alleles to design amplification procedures that were able to obtain complete exon 2 and 3 sequences from DQA1 genes. This coding information allowed us to reduce the number of ambiguities for DQA1 typing. DQA1 intron 2 and 3 characterization demonstrated the presence of two polymorphisms for alleles with the same exons 2 and 3 sequence from DQA1*05 group. Different samples including the DQA1*050101 alleles showed a single nucleotide polymorphism at position 53 of intron 2 (G53T). Additional haplotypic analysis showed the possible association of T53 allele with the Ax-Cw5-B18-DR17-DQ2 extended haplotype. On the other hand, DQA1*0505 sequencing from different control samples noticed the existence of a microsatellite (TTTC/AAAG)n located at position 126 of intron 3. Fragment length analysis demonstrated a high polymorphism for this short tandem repeat system (0505STR), defining alleles that ranged from 8 to 20 repetitions in our population.     

The nature of polymorphism of the HLA-DRB intron sequences is lineage specific

Abstract: The sequence database of HLA-DRB genes is mainly derived from mRNA analysis or has focused exclusively on the polymorphism of the 2nd exon. Little is known about the non-coding sequences of the different DRB alleles which represent about 94% of the genes. In this study we have determined the sequence of the 3' 500 bp intron 1 fragment adjacent to exon 2 in all serologically defined HLA-DRB genes and their most frequent allelic subtypes. The intron sequences turned out to be highly polymorphic. Similar to the class I introns, this variability was not characterized by random point mutations but by a highly systematic diversity reflecting the lineage-specific relationship of the HLA-DR alleles. With a few exceptions in DRB1*15, 13 and 08 as well as DRB 4 and 5, the variability mirrors the serological diversity. As well as delivering insight into the genetic relationship between the different DRB alleles, these sequences will provide an extremely valuable basis for developing advanced DRB sequencing strategies for clinical purposes.     

Nucleotide sequence and polymorphism of the caprine major histocompatibility complex class II DQA1 (Cahi-DQA1) gene

The major histocompatibility class II DQ molecules are dimeric glycoproteins involved in antigen presentation to CD4(+) T cells. In the current work, we have performed the molecular analysis of the goat Cahi-DQA1 gene. Sequencing of the Cahi-DQA1 cDNA revealed a single 768bp open reading frame. The alignment of this sequence with its bovine and ovine DQA1 counterparts revealed a remarkable degree of nucleotide identity (92-93% for the most similar bovine and ovine sequences). Moreover, we amplified a region including the 3'-end of intron 1, exon 2 and the 5'-end of intron 2. We identified seven Cahi-DQA1 alleles that likely correspond to four different allelic lineages. The alignment of these seven Cahi-DQA1 alleles revealed the existence of 23 amino acid polymorphic sites, seven of which (alpha(10), alpha(55), alpha(56), alpha(68), alpha(69), alpha(71) and alpha(76)) are highly polymorphic with at least three amino acid substitutions. Ten of the 23 polymorphic amino acid sites were included in the peptide binding region and consequently they might play a crucial role in immunological processes modulating disease pathogenesis.     

Molecular analysis of HLA-B35 alleles and their relationship to HLA-B15 alleles

The HLA-B35 serotype is one of the largest allelic groups of HLA class I molecules and includes four isotypes. Of the four, the B35 variant isoform is relatively rare and is the most acidic form. DNA sequencing of die rare isoforms revealed three alleles, B*1522, B*3511, and B*3517. A phylogenetic tree of HLA-B15- and HLA-B35-related alleles for the exon 2 and 3 nucleotide sequences showed that exon 2 of B* 1522 clusters with B35 alleles whereas exon 3 clusters with B15 alleles. Branches of the tree suggest that the serodeterminants of B35, B62, B63, and B70 may reside in the αl domain, encoded by exon 2. The B*1520 and B*1522 genes, which type as B62 and B35, respectively, are hybrid molecules alternatively using exon 2 and exon 3 sequences of B*3501 and B*1501. A comparison of intron 2 sequences for B*3501, B*1501 and B*1522 suggests that the recombination site may have been in the region at the 3' end of intron 2. Despite being flanked by two highly polymorphic exons (exons 2 and 3), intron 2 is relatively well conserved in the B-locus, and it is characterized by seven to eight tandem repeats of the CGGGG pentanucleotide. A high degree of sequence homology and repetitive sequences are essential for a significant frequency of recombination. In this report, we reveal more about the complex evolutionary history of the HLA-B alleles     

Mhc-E polymorphism in Pongidae primates: the same allele is found in two different species

Mhc-E intron 1, exon 2, intron 2, and exon 3 from pygmy chimpanzee (Pan paniscus) , chimpanzee (Pan troglodytes) , gorilla (Gorilla gorilla) and orangutan (Pongo pygmaeus) have been sequenced; six new Mhc-E alleles have been obtained but sequence changes are only placed either in introns or in synonymous exonic bases. One pygmy chimpanzee Mhc-E DNA sequence is identical to another sequence from chimpanzee; the fact that no variation is found also at the intronic level suggests that these two species of chimpanzee may have recently separated and/or that both of them might only represent subspecies. Mhc-E phylogenetic trees separate two evolutionary groups: Pongidae , including humans, and Cercopithecinae ; this is also found by studying another non-classical class I gene, Mhc-G. The Mhc-E alleles' invariance at the protein level supports that strong selective forces are operating at the Mhc-E locus, as has also been found in both Cercopithecinae and humans. These allelic and evolutionary data suggest an altogether different functionality for HLA-E (and also HLA-G ) compared with classical class I proteins: i.e., sending negative (tolerogenic) signals to NK and T cells.     

Evolutionary relationships between HLA-B alleles as indicated by an analysis of intron sequences

The HLA-B locus is the most polymorphic of the class I genes encoded within the human major histocompatibility complex. This polymorphism is mainly located in exons 2 and 3, which code for the molecule's alpha1 and alpha2 domains and includes the antigenic peptide binding site. However, information about adjacent non-coding regions (introns 1 and 2) has not been extensively reported but could be very important in establishing an understanding of the evolutionary mechanisms involved in the polymorphism generation of HLA-B and the Mhc loci. In the present work, introns 1 and 2 of 14 HLA-B alleles are studied and their significance is discussed; 10 have been sequenced in our own laboratory and the other 4 have been previously reported by others. Different serological families share the complete intron 1 sequence; at this region, 12 out of 14 HLA-B alleles could be included in four groups with the same intron 1 sequence: a) B*0702, B*4201, B*4801; b) B*27052, B*4002, B*4011; c) B*40012, B*4101, including B*4501, B*5001 (these latter two alleles have specific characteristics in both introns 1 and 2, which may reflect a common evolutionary pathway); and d) B*44031, B*44032. The other alleles, B*1402, and B*1801, do not have identical intron 1 sequences compared to any of the described groups, but share many similarities with them. The B*1801 evolutionary pathway seems to be very specific since it branches separately from other alleles both in intron 1 and intron 2 dendrograms. On the other hand, HLA-B allelic group distribution and similarities according to intron 1 sequences were not confirmed when using intron 2, especially in the cases of B*4002, B*4101 and B*4801. This would suggest that both point mutations fixed by genetic drift and gene conversion events are involved in HLA-B diversification. The latter events could be supported by the strong homology between intron 1 and, to a lesser extent, intron 2, and also the CG content within them. Finally, the precise knowledge of these non-coding regions could be important for developing DNA base typing strategies for the HLA-B alleles.     

Phylogenetic history of hominoid DRB loci and alleles inferred from intron sequences

Summary: The evolutionary relationships among the MHC class II DRB4, DRB5 and DRB6 loci as well as the allelic lineages and alleles of the DRB1 locus were studied based on intron 1 and Intron 2 sequences from humans, chimpanzee (Pan troglodytes). bonobo (Pan paniscus) and gorilla (Gorilla gorilla). The phylogenetic trees for these sequences indicate that most of the DRB 1 allelic lineages predate the separation of the hominoid species studied, consistent with previous analysis of the coding sequences of these lineages. However, the intron sequence variation among alleles within DRB1 allelic lineages is very limited, consistent with the notion that the majority of the contemporary alleles have been generated within the last 250,000 years. The clustering of the DRB1 allelic lineages *08 and *12 with *03 supports a common ancestry for the DR8 and DR52 haplotypes. Similarly, the clustering of DRB1 allelic lineages *15 and *01 with the DRB3 locus is consistent with a common ancestry for the DR1 and DR51 haplotypes. Two cases of recombination around the second exon were observed: 1) the HLA-DRB6 locus appears to have been generated through a recombination between a DRB5 allele and an ancestral DRB6 allele, and 2) the gorilla sequence Gogo-DRB1*03 appears to have been generated through a recombination between the DRB3 locus and an allele from the DRB1*03 allelic lineage. The nucleotide substitution rate of DRB introns was estimated to 0.85–1.63 × 10^-9 per site per year, based on comparisons between the most closely related sequences from different hominoid species. This estimate is similar to the substitution rate for other intronic regions of the primate genome.     

Four novel human leukocyte antigen-DQA1 alleles identified in the Korean population

Four novel human leukocyte antigen (HLA)-DQA1 alleles have been characterized by direct DNA sequencing of coding exons 1-4. All the novel alleles exhibited a single nucleotide substitution either in exon 3 or in exon 4 when compared with previously defined alleles. Thus, it is likely that alleles were generated by point mutation from pre-existing alleles in the population. Substitutions resulted in either a silent (DQA1*010203) or an amino acid change (DQA1*0506, DQA1*0507, and DQA1*0508). The substituted sites were both previously known polymorphic and conserved positions. Putative haplotypes associated with the novel alleles were deduced based on the HLA types shared by the individuals carrying a novel allele or from previously reported population data.     

Identification of two novel DQA1 alleles, DQA1*0107 and DQA1*0602, by sequence-based typing in the GoKinD population

Two novel DQA1 alleles, DQA1*0107 and DQA1*0602, were discovered using DQA1 sequence-based typing (SBT) in participants in the Genetics of Kidneys in Diabetes (GoKinD) Study. The DQA1*0107 allele, found in three unrelated Caucasian participants, contains a novel polymorphism at codon 79 of exon 2 (CGC→TGC), which results in an amino acid change from an arginine to a cysteine. The participants containing this novel polymorphism also had a 1-bp insertion in intron 2 that is common to the *01 alleles. The DQA1*0602 allele, found in one Caucasian participant, contains a novel polymorphism at codon 139 of exon 3 (AGC→CGC), which results in an amino acid change from a serine to an arginine. Additionally, the *0602 allele has a base change in intron 1 that is common to the *06 alleles. Both new alleles were isolated using single-allele amplification SBT and confirmed using sequence-specific primer amplification.     

MIC gene polymorphism and haplotype diversity in rhesus macaques

Rhesus macaques (Macaca mulatta) mainly originating from India were analysed for their major histocompatibility complex class I-related (MIC) gene repertoire. Thus far, three distinct genes, designated MIC1, MIC2 and MIC3, have been identified in the rhesus macaque. In addition, an MICD pseudogene has been described mapping apart from the other loci in a telomeric direction. Genomic comparisons and the presence of a characteristic microsatellite in exon 5 suggest that the MIC1 gene is the equivalent of the human MICA gene. Hence, the MIC2 gene, lacking the microsatellite - as do humans -, is considered to be the equivalent of human MICB. The MIC3 gene, a hybrid of MICA and MICB, seems to be generated by a crossing-over event with one breakpoint in intron 3 and accordingly is named MICA/B. Apart from their human counterparts, MICA, MICB and MICA/B cluster in separate branches in the phylogenetic tree, confirming the hybrid character of the MICA/B gene. Population analyses have shown that the various genes display polymorphism, and six MICA, five MICB and three MICA/B alleles have been identified. In the panel of homozygous typing cells, two distinct haplotype configurations have been defined by segregation analyses. Each haplotype comprises an MICB gene in conjunction with either an MICA or an MICA/B gene. Furthermore, the presence of a polymorphic microsatellite in the MICA and MICA/B alleles facilitates speedy and accurate haplotyping.     

Comprehensive sequence analysis of HLA-DQA1 and -DQB1 alleles and characterization of DRB1-QAP-DQA1-DQB1 haplotypes

It is well known that both chain and β chain of HLA-DQ are highly polymorphic. However the polymorphisms outside the hypervariable region were not fully examined so far. To further clarify the polymorphisms in DQ genes, we determined the nucleotide sequences of full length cDNA, spanning from the leader sequence to the stop codon, from 15 DQA1 alleles and 15 DQB1 alleles. We identified several new DQ alleles which had identical exon 2 sequence and were different in other exons. On the basis of the sequence analyses, a comprehensive PCR-based oligotyping system for DQA1 gene was established. We then characterized DRB1-QAP(DQA1 promoter)-DQA1-DQB1 haplotypes of B-lymphoblastoid cell lines homozygous for HLA and healthy unrelated Japanese and Norwegian populations. It was revealed that DQA1 alleles, which were identical in exon 2 but different in other exons, showed close linkage disequilibrium with diferent characteristic DRB1, QAP and DQB1 alleles. These results suggest that DR-DQ haplotypes have been generated in the early stage of molecular evolution.     

HLA-DQB1 sequencing-based typing using newly identified conserved nucleotide sequences in introns 1 and 2

Sequencing-based typing (SBT) human leukocyte antigen (HLA) class I and II genes should examine entire exon sequences where polymorphisms lie. Primers for the amplification of complete exons therefore anneal in introns and their design relies on accurate intron sequences being available. We decided to develop a SBT method for HLA-DQB1 using amplification primers which anneal in introns 1 and 2, yet the amount of intron sequence data previously available in databases was sparse. Therefore, we undertook a systematic sequencing of introns 1 and 2 using DNA from cell lines homozygous for DQB1. This study confirmed an earlier report that the non-coding regions of this gene are the most polymorphic seen in the human genome. Intron sequences within an allele group were largely identical, the exceptions being DQB1*0301 differing from other DQB1*03 allele groups and DQB1*0601 differing from all other DQB1*06 alleles. A retroviral Alu element, related to the AluYa5a2 subfamily, was identified uniquely inserted in intron 2 of DQB1*02 alleles. For the typing approach, six amplification primers were designed based on conserved allele group sequences covering all of the HLA DQB antigens, and two sequencing primers were also designed which anneal in intron 2. This method has proved to be very robust and has been used as part of a referral DNA sequencing service for a number of years.     

DQA103 subtypes have different associations with DRB1 and DQB1 alleles

Polymorphisms outside the hypervariable regions of HLA class II alleles that do not affect the peptide-binding site are probably not under selective pressure and could therefore be useful as markers of the evolutionary pathways of the HLA class II haplotypes. We have analyzed such a polymorphism in the variants of DQA1^*03, which differ at residue 160 encoded in exon 3. Our study included homozygous BCLs of the 10th IHWS and samples of a multiracial panel of 723 unrelated subjects which were also typed for allelic variations in exon 2 by hybridization with SSOP. BCLs having DQA1^*03 and 131 selected DQA1^*03-positive samples were typed for the dimorphism in exon 3 that distinguishes DQA1^*0301 and DQA1^*0302. DQA1^*0301 was found to be exclusively associated with DQB1^*0302, while samples carrying DQB1^*0201, 0301, 0303, and 0401 always had DQA1^*0302. A few haplotypes carrying DQB1^*0302 had DQA1^*0302. The fact that DQA1^*0301 is completely included in DQB1^*0302, and not vice versa, suggests that DQA1^*0301 may have arisen from a mutation in a haplotype containing DQA1^*0302-DQB1^*0302. DQB1^*0302 was found to be associated with all DR4 subtypes, suggesting possibly that the current variants of DRB1-DR4 may be of more recent origin. DRB1^*0405 was the only subtype of DR4 which was not associated with DQA1^*0301 and had multiple associations with the DQB1 alleles, therefore, perhaps representing the oldest allele of this group.     

A simple and rapid method for HLA-DQA1 genotyping by digestion of PCR-amplified DNA with allele specific restriction endonucleases

The second exon of the HLA-DQA1 genes was selectively amplified from genomic DNAs of 72 HLA-homozygous B cell lines by the polymerase chain reaction (PCR). Amplified DNAs were digested with HaeIII, Ddel, ScrFI, FokI and RsaI, which recognize allelic sequence variations in the polymorphic segments of the DQA1 second exon, and then subjected to electrophoresis in polyacrylamide gels. Eight different polymorphic patterns of restriction fragments were obtained, and seven were identical to patterns predicted from the known DNA sequences, correlating with each HLA-DQw type defined by serological typing. The remaining one pattern cannot be explained from the sequence data, suggesting the presence of a novel DQA1 allele at the nucleotide level. This PCR-RFLP method provides a simple and rapid technique for accurate definition of the HLA-DQ types at the nucleotide level, eliminating the need for radioisotope as well as allele specific oligonucleotide probes and can be extended and applied to HLA-DR, -Dw DP typing.     

Genomic characterization of KIR2DL4 in families and unrelated individuals reveals extensive diversity in exon and intron sequences including a common frameshift variation occurring in several alleles

The KIR2DL4 gene including a portion of exon 1 through exon 9 was sequenced from two families and eight cell lines from the International Histocompatibility Workshop (IHWS). Two known alleles and eight variants were detected. Overall, there were five synonymous and three non-synonymous changes when the variants were compared to the coding sequences of the most closely related known alleles plus a common frameshift change in five of the variant alleles. Alignment of the new variants with all known alleles showed that the regions encoding the extracellular region and the cytoplasmic tail were the most polymorphic. Two non-synonymous changes, P146H and L161V, occurred in an extracellular immunoglobulin-like domain. Five of the eight variants had a single adenine deletion in the exon encoding the transmembrane region, potentially resulting in a truncated protein lacking the cytoplasmic tail. The distribution of the deletion variant among many KIR2DL4 alleles may explain the high frequency of this variation in the population. Four of the eight consanguineous IHWS cell lines were found to be heterozygous for KIR2DL4 carrying two alleles that differed from one another by a few nucleotide substitutions. Analysis of intron sequences in the families revealed the nature and distribution of interspersed repeat elements which comprise 46% of the KIR2DL4 nucleotide sequence and consist of 12 elements including six SINEs (13.73% of the total length), one LINE (12.41%), and five LTR elements (19.51%). The results revealed the presence of extensive diversity in the KIR2DL4 gene. This is the first extensive report providing both exon and intron data in related individuals.     

PCR-RFLP-based Mamu-DQB1 typing of rhesus monkeys: characterization of two novel alleles

Up to now 19 allelic sequences of the rhesus monkey DQB1 locus have been published. Referring to these sequences, we have developed a typing protocol for Mamu-DQB1 alleles which was verified by additional cloning, sequence analysis and segregation studies. The protocol is based on the amplification of the second exon with only one specific primer pair followed by the digestion of the PCR products with up to 10 different restriction endonucleases. The alleles can be identified in homozyous and heterozygous combinations since most amplified second exon sequences give unique band patterns after digestion with at least one of the selected restriction endonucleases. By the use of this protocol we analyzed DNA-samples from 182 rhesus monkeys. Among these samples two novel Mamu-DQB1 alleles were detected, subsequently cloned and their nucleic sequence determined. Since we typed four complete breeding groups consisting of two generations we were able to identify several DQ haplotypes by segregation analysis using the previously developed typing protocol for DQA1.     

DQA-DQB linkage in Old World monkeys

Abstract: The MHC DQ region in nonhuman primates, as in humans, consists of a and β chains that are polymorphic with strong linkage disequilibrium between certain DQA-DQB alleles. Not only are contemporary HLA class II allelic variants present in evolutionarily distant species, but we demonstrate that linkages between loci also bear ancient roots. In unrelated baboons ( Papio cynocephalus anubis ) and family segregation analysis of pig-tailed macaques ( Macaca nemestrina ) we found cis -linkages between DQA1*O1 and DQB1*05 or *06, between DQA1*05 and DQB1*03, and between DQA1*03 and DQB1*03 alleles, all of which are also prominent in modern humans. In contrast, one linkage that has not been seen in humans, between DQA1*05 and DQB1*06 alleles, was also found. These patterns of selective linkage disequilibrium imply evolutionary mechanisms following the divergence of species that constrain the diversity of haplotypes which evolve.     

A simple and rapid method for HLA-DP genotyping by digestion of PCR-amplified DNA with allele-specific restriction endonucleases

We previously reported a simple and rapid method for HLA-DQA genotyping by digestion of polymerase chain reaction-amplified DQA genes with allele-specific restriction endonucleases. Here we report the application of this method to DP genotyping. The second exon of the HLA-DPB genes was selectively amplified from genomic DNAs of 72 HLA-D homozygous B-cell lines by the polymerase chain reaction method. Amplified DNAs were digested with ApaI, SacI, BstUI, FokI, and RsaI, which can recognize allelic sequence variations in the polymorphic segments of the DPB second exon and then subjected to electrophoresis in polyacrylamide gels. Sixteen different polymorphic patterns of the restriction fragments were found, and twelve were identical to patterns predicted from the known DNA sequences correlating with each HLA-DPw specificity defined by cellular typing. The other four patterns were distinct from those of the known DPw specificities, suggesting the presence of novel DP alleles. This polymerase chain reaction-restriction fragment length polymorphism method provides a simple and rapid technique for accurate definition of HLA-DP types at the nucleotide level, replacing the technically demanding method of primed lymphocyte typing.     

Nucleotide sequences of MHC class I introns 1, 2, and 3 in humans and intron 2 in nonhuman primates

HLA-class I genes are the most polymorphic genetic system yet known. The polymorphic substitutions are mostly located in exon 2 and 3, encoding α1 and α2 domains, respectively, which are involved in peptide binding and T cell receptor interaction. In this study, we present the sequences of the introns neighboring the polymorphic exons in humans with few examples from non-human primates. In general, intron sequences are found to be less polymorphic than the adjacent exons, displaying numerous locus-specific and group-specific sites. These sequences will provide important information for developing DNA based typing strategies for HLA-class I alleles.