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.

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.

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.

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.

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.

The genomic sequence and analysis of the swine major histocompatibility complex.

We describe the generation and analysis of an integrated sequence map of a 2.4-Mb region of pig chromosome 7, comprising the classical class I region, the extended and classical class II regions, and the class III region of the major histocompatibility complex (MHC), also known as swine leukocyte antigen (SLA) complex. We have identified and manually annotated 151 loci, of which 121 are known genes (predicted to be functional), 18 are pseudogenes, 8 are novel CDS loci, 3 are novel transcripts, and 1 is a putative gene. Nearly all of these loci have homologues in other mammalian genomes but orthologues could be identified with confidence for only 123 genes. The 28 genes (including all the SLA class I genes) for which unambiguous orthology to genes within the human reference MHC could not be established are of particular interest with respect to porcine-specific MHC function and evolution. We have compared the porcine MHC to other mammalian MHC regions and identified the differences between them. In comparison to the human MHC, the main differences include the absence of HLA-A and other class I-like loci, the absence of HLA-DP-like loci, and the separation of the extended and classical class II regions from the rest of the MHC by insertion of the centromere. We show that the centromere insertion has occurred within a cluster of BTNL genes located at the boundary of the class II and III regions, which might have resulted in the loss of an orthologue to human C6orf10 from this region.

cDNA cloning and genetic polymorphism of the swine major histocompatibility complex (SLA) class II DMA gene.

cDNA clones corresponding to the swine histocompatibility complex (SLA: swine leucocyte antigen)-DM alpha chain were isolated using the polymerase chain reaction (PCR) products from the third exon in the human HLA-DMA gene as a probe. Amino acid comparative analysis revealed that these clones were more closely related to the bovine and human DMA genes than to the other swine class II genes alpha chain genes, DRA, DQA and DOA. These results suggest that the SLA-DMA gene is expressed and may function, like HLA-DM, as an important modulator in class II restricted antigen processing in swine. Furthermore, based on the sequences and PCR-restriction fragment length polymorphism (PCR-RFLP) patterns in the SLA-DMA gene, no allelic variation was recognized in the second exon, but five allelic variations were recognized in the third exon in five different breeds of swine. These DMA alleles were defined by variation at four nucleotide positions. Two of these alleles resulted in an amino acid substitution. These results suggest that SLA-DMA has little polymorphism as observed in HLA-DMA and mouse H2-Ma.

Genomic anatomy of a premier major histocompatibility complex paralogous region on chromosome 1q21-q22.

Human chromosomes 1q21-q25, 6p21.3-22.2, 9q33-q34, and 19p13.1-p13.4 carry clusters of paralogous loci, to date best defined by the flagship 6p MHC region. They have presumably been created by two rounds of large-scale genomic duplications around the time of vertebrate emergence. Phylogenetically, the 1q21-25 region seems most closely related to the 6p21.3 MHC region, as it is only the MHC paralogous region that includes bona fide MHC class I genes, the CD1 and MR1 loci. Here, to clarify the genomic structure of this model MHC paralogous region as well as to gain insight into the evolutionary dynamics of the entire quadriplication process, a detailed analysis of a critical 1.7 megabase (Mb) region was performed. To this end, a composite, deep, YAC, BAC, and PAC contig encompassing all five CD1 genes and linking the centromeric +P5 locus to the telomeric KRTC7 locus was constructed. Within this contig a 1.1-Mb BAC and PAC core segment joining CD1D to FCER1A was fully sequenced and thoroughly analyzed. This led to the mapping of a total of 41 genes (12 expressed genes, 12 possibly expressed genes, and 17 pseudogenes), among which 31 were novel. The latter include 20 olfactory receptor (OR) genes, 9 of which are potentially expressed. Importantly, CD1, SPTA1, OR, and FCERIA belong to multigene families, which have paralogues in the other three regions. Furthermore, it is noteworthy that 12 of the 13 expressed genes in the 1q21-q22 region around the CD1 loci are immunologically relevant. In addition to CD1A-E, these include SPTA1, MNDA, IFI-16, AIM2, BL1A, FY and FCERIA. This functional convergence of structurally unrelated genes is reminiscent of the 6p MHC region, and perhaps represents the emergence of yet another antigen presentation gene cluster, in this case dedicated to lipid/glycolipid antigens rather than antigen-derived peptides.

Characterization of alkaline phosphatase genes expressed in seminoma by cDNA cloning.

Two members of a placental alkaline phosphatase (PLAP) family, PLAP and PLAP-like or germ cell alkaline phosphatase, are aberrantly expressed in tumors of ecotropic origin. To characterize alkaline phosphatase induced in seminoma, alkaline phosphatase cDNA clones were isolated from a cDNA library constructed from seminoma cells and characterized by nucleotide sequence determination. Thus, isolated cDNA clones were classified into two types, germ cell alkaline phosphatase (PLAP-like) and liver/bone/kidney-type alkaline phosphatase (L/B/K AP). These results suggest that other than the PLAP family members, the expression of L/B/K AP is enhanced in seminoma and can serve as a tumor marker in seminoma.

Extended HLA haplotypes in Japanese homozygous typing cells.

We have defined extended HLA haplotypes including the HLA class II genes, the non-HLA genes such as TAP1, TAP2 and LMP2, and the (CTG)n microsatellite repeats within the NOTCH4 gene between DRA and 21OH in 33 Japanese HLA homozygous typing cells (HTC). These conserved haplotypes characterized by unique linkage might be maintained as a result of functional co-operation among them in the antigen presentation pathway. These HTCs can be served as an original and ethnic-specific standard panel, providing useful genetic markers in haplotypic diversity, disease association, and anthropology studies.

Triplet repeat polymorphism in the NOTCH4 gene with the human major histocompatibility complex in a healthy population and patients with a salivary gland tumor in Japan.

The NOTCH4 gene, the human counterpart of the mouse mammary tumor gene, int-3, has been recently localized near the boundary of the HLA class II and class III regions. This gene is one of candidates for development of salivary gland tumor. Microsatellite polymorphism of (CTG)n repeat in the signal peptide domain of NOTCH4 was analyzed in Japanese including the patients with salivary gland tumor. Four alleles consisting of 6, 9, 10 and 11 repetitions of CTG (Leu) were observed and found to be in linkage disequilibria with HLA class I and class II alleles. No significant association of this microsatellite polymorphism with the disease were observed in 26 samples of salivary gland tumor. In this neoplasia, neither large-scale deletion nor translocation was detected around the NOTCH4 gene using genomic Southern hybridization analysis by the NOTCH4 cDNA as a probe.

Triplet repeat polymorphism within the NOTCH4 gene located near the junction of the HLA class II and class III regions in narcolepsy.

A polymorphic (CTG)n microsatellite repeat was found in the signal peptide domain of the NOTCH4 gene located near the junction of the class II and class III regions of the human major histocompatibility complex. This gene belongs to a multigene family of NOTCH originally identified as a differential factor of neuronal cells. To ascertain whether the NOTCH4 gene is involved in the development of neurogenic disease, narcolepsy, which is known to be tightly associated with HLA-DR15, this microsatellite polymorphism of the (CTG)n repeat was analyzed in Japanese patients with narcolepsy One allele, 9 repetitions of CTG (Leu) was significantly increased in the patient group. However, the significant increase of this allele in the patient group could be explained by a strong linkage disequilibrium with the HLA class II alleles, DRB1*1501, DQA1*0102 and DQB1*0602, which were more strongly associated with the disease. These results suggest that the (CTG)n repeat polymorphism in NOTCH4 does not primarily determine the susceptibility to narcolepsy.

cDNA cloning of the human homologues of the mouse Ke4 and Ke6 genes at the centromeric end of the human MHC region.

cDNA clones corresponding to the HKE4 and HKE6 genes at the centromeric end of the HLA region on human chromosome 6p21.3 were isolated and characterized. The predicted amino acid sequences of HKE4 and HKE6 exhibited 81.5 and 85.6% identity to the mouse homologues, Ke4 and Ke6, respectively. HKE4 may encode a membrane protein with histidine-rich charge clusters. HKE6 possesses remarkable amino acid sequence conservation with several bacterial proteins with oxidoreductase function and also shows significant homology with the two unique functional domains containing the nucleotide cofactor binding site and the consensus motif characteristic of the members of the superfamily of short-chain alcohol dehydrogenases such as human and rat steroid and prostaglandin dehydrogenases.