HLA alleles and haplotypes in Burmese (Myanmarese) and Karen in Thailand.

This is the first report on human leukocyte antigen (HLA) allele and haplotype frequencies at three class I loci and two class II loci in unrelated healthy individuals from two ethnic groups, 170 Burmese and 200 Karen, originally from Burma (Myanmar), but sampled while residing in Thailand. Overall, the HLA allele and haplotype frequencies detected by polymerase chain reaction sequence-specific primer (PCR-SSP) at five loci (A, B, C, DRB1 and DRQB1) at low resolution showed distinct differences between the Burmese and Karen. In Burmese, five HLA-B*15 haplotypes with different HLA-A and HLA-DR/DQ combinations were detected with three of these not previously reported in other Asian populations. The data are important in the fields of anthropology, transplantation and disease-association studies.

Characterization of swine leukocyte antigen alleles and haplotypes on a novel miniature pig line, Microminipig.

Microminipigs are extremely small-sized, novel miniature pigs that were recently developed for medical research. The inbred Microminipigs with defined swine leukocyte antigen (SLA) haplotypes are expected to be useful for allo- and xenotransplantation studies and also for association analyses between SLA haplotypes and immunological traits. To establish SLA-defined Microminipig lines, we characterized the polymorphic SLA alleles for three class I (SLA-1, SLA-2 and SLA-3) and two class II (SLA-DRB1 and SLA-DQB1) genes of 14 parental Microminipigs using a high-resolution nucleotide sequence-based typing method. Eleven class I and II haplotypes, including three recombinant haplotypes, were found in the offspring of the parental Microminipigs. Two class I and class II haplotypes, Hp-31.0 (SLA-1*1502-SLA-3*070102-SLA-2*1601) and Hp-0.37 (SLA-DRB1*0701-SLA-DQB1*0502), are novel and have not so far been reported in other pig breeds. Crossover regions were defined by the analysis of 22 microsatellite markers within the SLA class III region of three recombinant haplotypes. The SLA allele and haplotype information of Microminipigs in this study will be useful to establish SLA homozygous lines including three recombinants for transplantation and immunological studies.

A genome-wide association study of anorexia nervosa.

Anorexia nervosa (AN) is a complex and heritable eating disorder characterized by dangerously low body weight. Neither candidate gene studies nor an initial genome-wide association study (GWAS) have yielded significant and replicated results. We performed a GWAS in 2907 cases with AN from 14 countries (15 sites) and 14 860 ancestrally matched controls as part of the Genetic Consortium for AN (GCAN) and the Wellcome Trust Case Control Consortium 3 (WTCCC3). Individual association analyses were conducted in each stratum and meta-analyzed across all 15 discovery data sets. Seventy-six (72 independent) single nucleotide polymorphisms were taken forward for in silico (two data sets) or de novo (13 data sets) replication genotyping in 2677 independent AN cases and 8629 European ancestry controls along with 458 AN cases and 421 controls from Japan. The final global meta-analysis across discovery and replication data sets comprised 5551 AN cases and 21 080 controls. AN subtype analyses (1606 AN restricting; 1445 AN binge-purge) were performed. No findings reached genome-wide significance. Two intronic variants were suggestively associated: rs9839776 (P=3.01 × 10(-7)) in SOX2OT and rs17030795 (P=5.84 × 10(-6)) in PPP3CA. Two additional signals were specific to Europeans: rs1523921 (P=5.76 × 10(-)(6)) between CUL3 and FAM124B and rs1886797 (P=8.05 × 10(-)(6)) near SPATA13. Comparing discovery with replication results, 76% of the effects were in the same direction, an observation highly unlikely to be due to chance (P=4 × 10(-6)), strongly suggesting that true findings exist but our sample, the largest yet reported, was underpowered for their detection. The accrual of large genotyped AN case-control samples should be an immediate priority for the field.

Variation and linkage disequilibrium between a structurally polymorphic Alu located near the OR12D2 gene of the extended major histocompatibility complex class I region and HLA-A alleles.

We investigated the genetic structure and population frequency of an Alu repeat dimorphism (absence or presence) located near the OR12D2 gene within the olfactory receptor gene region telomeric of the alpha HLA class I region (HLA-J, -A, -G, -F). The structurally polymorphic Alu insertion (POALIN) locus rs33972478 that we designated as AluOR and its allele and haplotype frequencies and association with HLA-A and six other structurally polymorphic retroelements (3 Alu, 2 SVA and an HERVK9) were determined in 100 Japanese, 174 Caucasians and 100 African American DNA samples. The AluOR insertion varied in population frequency between 14.4% and 31.5% with significant differences between the Japanese and Caucasians, but not between the Caucasian and African Americans. Although AluOR is located 600 kb from the HLA-A gene, there was a significant linkage disequilibrium between the two loci and a high percentage association of the AluOR insertion with HLA-A29 (79%) in Caucasians and HLA-A31 (69.4%) in Japanese. Inferred haplotypes among three-locus to eight-locus haplotype structures showed maximum differences between the populations with the eight-locus haplotypes. The most frequent multilocus haplotype shared between the populations was the HLA-A2 allele in combination with the AluHG insertion. The AluOR whether investigated alone or together with the HLA class I alleles and other dimorphic retroelements is an informative ancestral marker for the identification of lineages and variations within the same and/or different populations and for examining the linkage and crossing-over between the HLA and OR genomic regions in the extended MHC.

HLA-DRB1, -DRB3, -DRB4 and -DRB5 genotyping at a super-high resolution level by long range PCR and high-throughput sequencing.

Super high-resolution single molecule sequence-based typing (SS-SBT) is a human leukocyte antigen (HLA) DNA typing method to the field 4 level of allelic resolution (formerly known as eight-digit typing) to efficiently detect new and null alleles without phase ambiguity by combination of long ranged polymerase chain reaction (PCR) amplification and next-generation sequencing (NGS) technologies. We previously reported the development and application of the SS-SBT method for the eight classical HLA loci, A, B, C, DRB1, DQA1, DQB1, DPA1 and DPB1. In this article, we describe the development of the SS-SBT method for three DRB1 linked loci, DRB3, DRB4 and DRB5 (DRB3/4/5) and characterization of DRB1-DRB3/4/5 haplotype structures to the field 4 level. Locus specific PCR primers for DRB3/4/5 were designed to amplify the gene regions from intron 1 to exon 6 [3' untranslated region (3'UTR)]. In total 20 DRB1 and 13 DRB3/4/5 allele sequences were determined by the SS-SBT to the field 4 level without phase ambiguity using 19 DR51, DR52 and DR53 positive genomic DNA samples obtained from Japanese. Moreover, 18 DRB1-DRB3/4/5 haplotypes were estimated to the field 4 level by the SS-SBT method in contrast to 10 haplotypes estimated by conventional methods to the field 1 level (formerly known as two digit typing). Therefore, DRB1-DRB3/4/5 haplotyping by SS-SBT is expected to provide informative data for improved HLA matching in medical research, transplantation procedures, HLA-related disease studies and human population diversity studies.

Association study of genetic variants on chromosome 7q31 with susceptibility to normal tension glaucoma in a Japanese population.

The caveolin 1 to caveolin 2 (CAV1-CAV2) gene region on chromosome 7q31 has been reported to be associated with susceptibility to primary open angle glaucoma (POAG) and normal tension glaucoma (NTG) in previous studies. We investigated whether genetic variants in the CAV1-CAV2 region are associated with NTG in Japanese patients. Two hundred and ninety-two Japanese patients with NTG and 352 Japanese healthy controls were recruited. We genotyped three single-nucleotide polymorphisms; that is, rs1052990, rs4236601, and rs7795356, in the CAV1-CAV2 gene region and assessed the allelic diversity among cases and controls. The frequency of the minor allele (G) of rs1052990 was significantly decreased in NTG cases compared with controls (P=0.014, OR=0.71), whereas NTG or POAG cases had a significantly higher frequency of the allele than controls in previous studies. Conversely, rs7795356 did not show any significant association with NTG cases, and rs4236601 was monomorphic in the Japanese study population. Our findings did not correspond with previous positive results, suggesting that CAV1-CAV2 variants studied in the present study are not important risk factors for NTG susceptibility in all populations. Further studies are needed to elucidate the possible contribution of the CAV1-CAV2 region to the development of glaucoma.

Super high resolution for single molecule-sequence-based typing of classical HLA loci at the 8-digit level using next generation sequencers.

Current human leukocyte antigen (HLA) DNA typing methods such as the sequence-based typing (SBT) and sequence-specific oligonucleotide (SSO) methods generally yield ambiguous typing results because of oligonucleotide probe design limitations or phase ambiguity for HLA allele assignment. Here we describe the development and application of the super high-resolution single-molecule sequence-based typing (SS-SBT) of HLA loci at the 8-digit level using next generation sequencing (NGS). NGS which can determine an HLA allele sequence derived from a single DNA molecule is expected to solve the phase ambiguity problem. Eight classical HLA loci-specific polymerase chain reaction (PCR) primers were designed to amplify the entire gene sequences from the enhancer-promoter region to the 3' untranslated region. Phase ambiguities of HLA-A, -B, -C, -DRB1 and -DQB1 were completely resolved and unequivocally assigned without ambiguity to single HLA alleles. Therefore, the SS-SBT method described here is a superior and effective HLA DNA typing method to efficiently detect new HLA alleles and null alleles without ambiguity.

Associations between six classical HLA loci and rheumatoid arthritis: a comprehensive analysis.

Although the HLA region contributes to one-third of the genetic factors affecting rheumatoid arthritis (RA), there are few reports on the association of the disease with any of the HLA loci other than the DRB1. In this study we examined the association between RA and the alleles of the six classical HLA loci including DRB1. Six HLA loci (HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1) of 1659 Japanese subjects (622 cases; 488 anti-cyclic citrullinated peptides (CCP) antibody (Ab) positive (82.6%); 103 anti-CCP Ab negative (17.4%); 31 not known and 1037 controls) were genotyped. Disease types and positivity/negativity for CCP autoantibodies were used to stratify the cases. Statistical and genetic assessments were performed by Fisher's exact tests, odds ratio, trend tests and haplotype estimation. None of the HLA loci were significantly associated with CCP sero-negative cases after Bonferroni correction and we therefore limited further analyses to using only the anti CCP-positive RA cases and both anti-CCP positive and anti-CCP negative controls. Some alleles of the non-DRB1 HLA loci showed significant association with RA, which could be explained by linkage disequilibrium with DRB1 alleles. However, DPB1*02:01, DPB1*04:01 and DPB1*09:01 conferred RA risk/protection independently from DRB1. DPB1*02:01 was significantly associated with the highly erosive disease type. The odds ratio of the four HLA-loci haplotypes with DRB1*04:05 and DQB1*04:01, which were the high-risk HLA alleles in Japanese, varied from 1.01 to 5.58. C*07:04, and B*15:18 showed similar P-values and odds ratios to DRB1*04:01, which was located on the same haplotype. This haplotype analysis showed that the DRB1 gene as well as five other HLA loci is required for a more comprehensive understanding of the genetic association between HLA and RA than analyzing DRB1 alone.

Genetic variation and hitchhiking between structurally polymorphic Alu insertions and HLA-A, -B, and -C alleles and other retroelements within the MHC class I region.

We investigated structurally polymorphic Alu insertions (POALINs) at five loci in the major histocompatibility complex (MHC) class I genomic region to determine their allele and haplotype frequencies and associations with the human leukocyte antigen (HLA)-A, -B, and -C genes in three populations, the Australian Caucasians, Japanese, and African Americans. The POALINs varied in allelic frequency between 0% and 42.3% with significant differences between populations at three of the five loci. The linkage disequilibrium (LD) between Alu insertions and the HLA-A, -B, or -C alleles and previously published polymorphic retroelements (four SVA and human endogenous retrovirus type 9 (HERVK9) loci) within the class I region of the MHC were calculated in pairwise analyses of haplotypes to show strong allelic associations and possible crossing-over events between some loci. Each POALIN was in significant LD with a variety of HLA-A, -B, or -C two-digit alleles probably as a result of hitchhiking. The POALINs helped to further stratify the HLA-A:B:C haplotypes into different POALIN:HLA-A:B:C haplotype frequencies. Of the multilocus haplotype analyses, the seven- and eight-locus haplotypes showed the largest number of differences between the populations, and fewer matched haplotypes between populations that ranged, for example, from 49% for HLA-B:HLA-A haplotypes to 7% for AluMICB:HLA-B:HLA-C:AluTF:AluHJ:HLA-A:AluHG:AluTF haplotypes in the Japanese. This comparative study of multilocus POALINs in the HLA class I region of three ethnic populations shows that POALINs alone or together with the HLA class I alleles and other retroelements are informative ancestral markers for assessing the interrelationship of HLA class I haplotype lineages, LD, and genetic diversity within the same and/or different populations.

SLA-DRB1 and -DQB1 genotyping by the PCR-SSOP-Luminex method.

A simple and novel genotyping method was developed to detect alleles at the swine leukocyte antigen (SLA)-DRB1 and -DQB1 class II loci by using polymerase chain reaction (PCR)-fluorescently labeled sequence-specific oligonucleotide probes (SSOPs) and Luminex 100 xMAP detection. The PCR-SSOP-Luminex method exhibited accuracy of 95% for both SLA-DRB1 and -DQB1 in 6 homozygous and 16 heterozygous pig samples as confirmed by sequencing the PCR products of the same samples. In addition, 12 low-resolution SLA class II haplotypes consisting of 7 and 9 DRB1 and DQB1 alleles were identified, respectively, in one population of 283 Landrace pigs. This genotyping method facilitates the rapid and accurate identification of two- or four-digit alleles at the SLA-DRB1 and -DQB1 loci.

NFKBIL1 confers resistance to experimental autoimmune arthritis through the regulation of dendritic cell functions.

We and others have reported that human NF-κB inhibitor-like-1 (NFKBIL1) was a putative susceptible gene for autoimmune diseases such as rheumatoid arthritis (RA). However, its precise role in the pathogenesis of RA is still largely unknown. In this study, we generated transgenic mice expressing human NFKBIL1 (NFKBIL1-Tg) and examined whether NFKBIL1 plays some role(s) in the development of autoimmune arthritis. In both a collagen-induced arthritis model and a collagen antibody-induced arthritis model, NFKBIL1-Tg mice showed resistance to arthritis compared to control mice, indicating that the gene product of NFKBIL1 was involved in the control of thusly induced arthritis. Total spleen cells of NFKBIL1-Tg mouse showed decreased proliferation to mitogenic stimuli, consistent with its resistance to arthritis. Unexpectedly, purified T cells of NFKBIL1-Tg mouse showed increased proliferation and cytokine production. This apparent discrepancy was accounted for by the impaired functions of antigen-presenting cells of NFKBIL1-Tg mouse; both T/B cell-depleted spleen cells and bone marrow-derived dendritic cells of the Tg mouse induced less prominent proliferation and IL-2 production of T cells. Furthermore, dendritic cells (DCs) derived from NFKBIL1-Tg mouse showed lower expression of co-stimulatory molecules and decreased production of inflammatory cytokines when they were activated by lipopolysaccharide. Taken together, these results indicated that NFKBIL1 affected the pathogenesis of RA at least in part through the regulation of DC functions.

Cytochrome P450 genetic polymorphisms influence the serum concentration of calcineurin inhibitors in allogeneic hematopoietic SCT recipients.

Calcineurin inhibitors are necessary as immunosuppressants during hematopoietic SCT (HSCT) to prevent alloreactivity, but have unfortunate toxicities. So, we investigated the association of gene polymorphisms with the initial calcineurin inhibitor concentration and the subsequent drug dose from day 1 to day 28 among patients who underwent HSCT at a single institution. We analyzed 58 serial cases of Japanese patients receiving GVHD prophylaxis with CsA (21 cases) or tacrolimus (37 cases). We investigated eight single-nucleotide polymorphisms: rs4244285 (CYP2C19), rs15524, rs4646450, rs3800959, rs776746 (CYP3A5), rs1128503, rs2032582 and rs1045642 (MDR1). The CsA concentration was significantly higher when the genotype of CYP3A5 rs15524 was T/T (P=0.044) or rs776746 was G/G (P=0.027). The CYP3A5 rs776746 and rs4646450 genotypes were also associated with tacrolimus concentration (P=0.013 and P=0.0058, respectively). The dosage of tacrolimus was remarkably reduced from day -1 to day 28 when the patient had the CYP3A5 rs4646450 C/C and/or rs776746 G/G genotype (P=0.0010 and P=0.0021, respectively). In this study, we show that genetic variation has a predictable effect on the pharmacological responses to calcineurin inhibitors in HSCT patients.

Cytotoxic T-lymphocyte antigen 4 haplotype correlates with relapse and survival after allogeneic hematopoietic SCT.

CTLA-4 is a negative regulator of activated T cells and the association of CTLA-4 polymorphisms with autoimmune diseases and transplant outcome has been reported. We evaluated the effect of donor CTLA-4 polymorphisms on outcome after allogeneic hematopoietic SCT (HSCT). We analyzed 147 Japanese HLA-matched sibling recipients and their donors who had undergone allogeneic HSCT. Genotyping of three single-nucleotide polymorphisms in CTLA-4 (-318, +49, CT60) was performed using TaqMan-PCR. According to the international HapMap database, only these three CTLA-4 haplotypes, classified as C-G-G, C-A-A and T-A-G, are present in the Japanese population. In this study, percentage expression of the C-G-G, C-A-A and T-A-G haplotypes was 59.5, 30.6 and 9.9%, respectively. Recipients of the C-A-A haplotype donor showed a significantly lower risk of relapse (HR: 0.54, 95% CI: 0.30-0.97, P=0.040) and a trend toward higher OS (HR: 0.61, 95% CI: 0.36-1.0, P=0.054) than did recipients of a donor without the C-A-A haplotype. The presence or absence of the C-A-A haplotype did not affect GVHD or non-relapse mortality. As the presence of the C-A-A haplotype reduced relapse risk and improved survival after allogeneic HSCT, this CTLA-4 haplotype may provide useful information for donor selection.

Polymorphic major histocompatibility complex class II Alu insertions at five loci and their association with HLA-DRB1 and -DQB1 in Japanese and Caucasians.

We investigated polymorphic Alu insertion (POALIN) frequencies at five loci in the major histocompatibility complex (MHC) class II genomic region to determine their allele and haplotype frequencies and associations with the human leukocyte antigen (HLA)-DRB1 and -DQB1 genes for 100 Japanese, 174 Australian Caucasians and 67 HLA reference cell lines obtained from different ethnic groups. The POALINs varied in frequency between 11% and 57% with significant differences between the Japanese and Caucasians at three loci. One POALIN locus deviated significantly from Hardy-Weinberg equilibrium (HWE) and four POALIN loci were in significant linkage disequilibrium and had a high percentage association with a variety of HLA-DRB1 or -DQB1 two-digit alleles. Inferred haplotype analysis among two-locus, five-locus and seven-locus haplotype structures showed maximum differences between the Japanese and Caucasians with the seven-locus haplotypes. The most common multilocus haplotype in Caucasians was DRB1*1501/DQB1*0602/AluDQ1/AluDRB1/AluORF10/AluDPB2 (6.7%), whereas the second most common allele HLA-DRB1*15 (17.5%) in Japanese was associated with three or four Alu insertions. The HLA class II POALINs also differentiated within and between HLA-DRB1 super-haplotypes DR1, DR8, DR51, DR52 and DR53. This is the first comparative population study of multilocus POALINs in the HLA class II region, which shows that POALINs whether investigated alone or together with the HLA class II alleles are informative genetic markers for the identification of allele and haplotype lineages and variations within the same and/or different populations.

Mapping of susceptibility locus for endometriosis within the HLA region using microsatellite markers in Japanese women.

Endometriosis is a female disorder characterized by the presence of uterine endometrial tissue in ectopic loci. Previous studies reported a higher prevalence of particular human leukocyte antigen (HLA) in endometriosis. In order to confirm the association between endometriosis and the HLA region, 15 polymorphic microsatellite markers distributed in the HLA class II to class III region were subjected to association analysis by polymerase chain reaction (PCR)-based DNA typing of 89 patients and 136 healthy controls. Statistical analysis of the allelic frequency at each microsatellite locus showed that there were no statistically significant differences in the allele frequency distributions between the cases and controls. This finding suggests that the etiology of endometriosis does not involve the HLA class II genomic region and a portion of class III genomic region in the Japanese population.

Investigation of the association between the GLC3A locus and normal tension glaucoma in Japanese patients by microsatellite analysis.

PURPOSE: To investigate whether the GLC3A locus harboring the CYP1B1 gene is associated with normal tension glaucoma (NTG) in Japanese patients. MATERIALS AND METHODS: One hundred forty-two Japanese patients with NTG and 101 Japanese healthy controls were recruited. Patients exhibiting a comparatively early onset were selected as this suggests that genetic factors may show stronger involvement. Genotyping and assessment of allelic diversity was performed on 13 highly polymorphic microsatellite markers in and around the GLC3A locus. RESULTS: There were decreased frequencies of the 444 allele of D2S0416i and the 258 allele of D2S0425i in cases compared to controls (P = 0.022 and P = 0.034, respectively). However, this statistical significance disappeared when corrected (Pc > 0.05). We did not find any significant association between the remaining 11 microsatellite markers, including D2S177, which may be associated with CYP1B1, and NTG (P > 0.05). CONCLUSIONS: Our study showed no association between the GLCA3 locus and NTG, suggesting that the CYP1B1 gene, which is reportedly involved in a range of glaucoma phenotypes, may not be an associated factor in the pathogenesis of NTG.

HLA-B*15 subtypes in Burmese population by sequence-based typing.

Human leukocyte antigen (HLA)-B*15 encompasses an increasing number of subtypes of more than 150. Frequency studies and a strong genetic association between HLA subtypes and susceptibility to drug hypersensitivity have been reported in different ethnic populations. To identify HLA-B*15 subtypes in Burmese using sequence-based typing (SBT) method, we selected 65 HLA-B*15-positive samples from 170 unrelated healthy Burmese who were genotyped HLA-B* by polymerase chain reaction with the sequence-specific primer method. The frequency of HLA-B*15 in Burmese was found to be 38.2%. By the SBT method, results showed 10 alleles of HLA-B*15 subtypes. Four common alleles, B*1502 (45.2%), B*1532 (16.4%), B*1525 (12.3%), and B*1501 (8.2%), were found in 82.1% of HLA-B*15-positive Burmese. Whereas the B*1501 was the highest in the Caucasians, Koreans, and Japanese, the highest frequency of HLA-B*15 alleles in Burmese was B*1502 (45.2%) that is similar to the frequency found in northeastern Thais and Vietnamese. This study is the first report of HLA-B*15 subtypes in Burmese. These results will provide the basic data in the further study in transplantations, genetic association with diseases, and drug hypersensitivity.

New susceptibility locus for high myopia is linked to the uromodulin-like 1 (UMODL1) gene region on chromosome 21q22.3.

PURPOSE: To ascertain and define the position of a potential disease susceptibility gene around D21S0083i prioritized during our previous whole genome case-control association analysis with 27,158 microsatellite markers, in Japanese high-myopia patients. METHODS: 520 high myopic patients and 520 healthy controls were genotyped using 39 SNPs distributed around D21S0083i on chromosome 21q22.3. RESULTS: Only 1 SNP (rs2839471) of 39 SNPs was significant after correction for multiple testing (allele T: P=0.00027, Pc=0.01, OR=1.684). The SNP (rs2839471) did not reside in haplotype blocks constructed by the pair-wise linkage disequilibrium between the SNPs. CONCLUSIONS: The SNP (rs2839471) is suggested to be located in the frequent recombinant region within UMODL1. Together this region might play a critical role for susceptibility to high myopia, and warrants further confirming studies and investigations as to the mechanisms by which UMODL1 may contribute to myopia.

Re-evaluation of heterogeneity in HLA-B*510101 associated with Behçet's disease.

Behçet's disease (BD) is a chronic inflammatory disease characterized by oral aphthous ulcers, genital ulcers, uveitis and skin lesions. Etiology and pathogenesis of BD are not fully elucidated, but the association with human leukocyte antigen (HLA)-B51 or B*5101 has been repeatedly reported. Previous studies have shown that there are few sequence variations in the protein-coding region of B51, while there is a report on many variations in the 5'-flanking region and intron. In this study, HLA-B*5101 gene from 37 individuals including Japanese, Turkish, Jordanian and Iranian patients and healthy controls were fully sequenced to further clarify the B*5101 gene in association with BD. We found that all the patients and healthy controls carried B*510101 with no variation in the 5'-flanking region, exon and intron. However, seven polymorphisms were found in the 3'-flanking region. These polymorphisms composed of six haplotypes that were shared and stretched over the ethnic groups, suggesting that the susceptibility to BD was conferred by the B*510101 itself and not by any genes in linkage disequilibrium with B*510101. In addition, phylogenetic analyses of B*510101 showed that the 3'-flanking sequences followed an evolutional divergence differently from that of the other regions, implying that a unifying selection might operate to conserve B*510101.