Regulatory SVA retrotransposons and classical HLA genotyped-transcripts associated with Parkinson's disease

Introduction

Parkinson's disease (PD) is a neurodegenerative and polygenic disorder characterised by the progressive loss of neural dopamine and onset of movement disorders. We previously described eight SINE-VNTR-Alu (SVA) retrotransposon-insertion-polymorphisms (RIPs) located and expressed within the Human Leucocyte Antigen (HLA) genomic region of chromosome 6 that modulate the differential co-expression of 71 different genes including the HLA classical class I and class II genes in a Parkinson's Progression Markers Initiative (PPMI) cohort.

Aims and methods

In the present study, we (1) reanalysed the PPMI genomic and transcriptomic sequencing data obtained from whole blood of 1521 individuals (867 cases and 654 controls) to infer the genotypes of the transcripts expressed by eight classical HLA class I and class II genes as well as DRA and the DRB3/4/5 haplotypes, and (2) examined the statistical differences between three different PD subgroups (cases) and healthy controls (HC) for the HLA and SVA transcribed genotypes and inferred haplotypes.

Results

Significant differences for 57 expressed HLA alleles (21 HLA class I and 36 HLA class II alleles) up to the three-field resolution and four of eight expressed SVA were detected at p<0.05 by the Fisher's exact test within one or other of three different PD subgroups (750 individuals with PD, 57 prodromes, 60 individuals who had scans without evidence of dopamine deficits [SWEDD]), when compared against a group of 654 HCs within the PPMI cohort and when not corrected by the Bonferroni test for multiple comparisons. Fourteen of 20 significant alleles were unique to the PD-HC comparison, whereas 31 of the 57 alleles overlapped between two or more different subgroup comparisons. Only the expressed HLA-DRA*01:01:01 and -DQA1*03:01:01 protective alleles (PD v HC), the -DQA1*03:03:01 risk (HC v Prodrome) or protective allele (PD v Prodrome), the -DRA*01:01:02 and -DRB4*01:03:02 risk alleles (SWEDD v HC), and the NR_SVA_381 present genotype (PD v HC) at a 5% homozygous insertion frequency near HLA-DPA1, were significant (Pc<0.1) after Bonferroni corrections. The homologous NR_SVA_381 insertion significantly decreased the transcription levels of HLA-DPA1 and HLA-DPB1 in the PPMI cohort and its presence as a homozygous genotype is a risk factor (Pc=0.012) for PD. The most frequent NR_SVA_381 insertion haplotype in the PPMI cohort was NR_SVA_381/DPA1*02/DPB1*01 (3.7%). Although HLA C*07/B*07/DRB5*01/DRB1*15/DQB1*06 was the most frequent HLA 5-loci phased-haplotype (n, 76) in the PPMI cohort, the NR_SVA_381 insertion was present in only six of them (8%).

Conclusions

These data suggest that expressed SVA and HLA gene alleles in circulating white blood cells are coordinated differentially in the regulation of immune responses and the long-term onset and progression of PD, the mechanisms of which have yet to be elucidated.

Nucleotide alterations in the HLA-C class I gene can cause aberrant splicing and marked changes in RNA levels in a polymorphic context-dependent manner

Polymorphisms of HLA genes, which play a crucial role in presenting peptides with diverse sequences in their peptide-binding pockets, are also thought to affect HLA gene expression, as many studies have reported associations between HLA gene polymorphisms and their expression levels. In this study, we devised an ectopic expression assay for the HLA class I genes in the context of the entire gene, and used the assay to show that the HLA-C*03:03:01 and C*04:01:01 polymorphic differences observed in association studies indeed cause different levels of RNA expression. Subsequently, we investigated the C*03:23N null allele, which was previously noted for its reduced expression, attributed to an alternate exon 3 3' splice site generated by G/A polymorphism at position 781 within the exon 3. We conducted a thorough analysis of the splicing patterns of C*03:23N, and revealed multiple aberrant splicing, including the exon 3 alternative splicing, which overshadowed its canonical counterpart. After confirming a significant reduction in RNA levels caused by the G781A alteration in our ectopic assay, we probed the function of the G-rich sequence preceding the canonical exon 3 3' splice site. Substituting the G-rich sequence with a typical pyrimidine-rich 3' splice site sequence on C*03:23N resulted in a marked elevation in RNA levels, likely due to the enhanced preference for the canonical exon 3 3' splice site over the alternate site. However, the same substitution led to a reduction in RNA levels for C*03:03:01. These findings suggested the dual roles of the G-rich sequence in RNA expression, and furthermore, underscore the importance of studying polymorphism effects within the framework of the entire gene, extending beyond conventional mini-gene reporter assays.

HLA Genetics for the Human Diseases

Highly polymorphic human leukocyte antigen (HLA) molecules (alleles) expressed by different classical HLA class I and class II genes have crucial roles in the regulation of innate and adaptive immune responses, transplant rejection and in the pathogenesis of numerous infectious and autoimmune diseases. To date, over 35,000 HLA alleles have been published from the IPD-IMGT/HLA database, and specific HLA alleles and HLA haplotypes have been reported to be associated with more than 100 different diseases and phenotypes. Next generation sequencing (NGS) technology developed in recent years has provided breakthroughs in various HLA genomic/gene studies and transplant medicine. In this chapter, we review the current information on the HLA genomic structure and polymorphisms, as well as the genetic context in which numerous disease associations have been identified in this region.

MHC-DRB alleles with amino acids Val11, Phe13, and the shared epitopes promote collagen-induced arthritis and a rapid IgG1 response in Filipino cynomolgus macaques

Macaques are useful animal models for studying the pathogenesis of rheumatoid arthritis (RA) and the development of anti-rheumatic drugs. The purpose of this study was to identify the major histocompatibility complex (MHC) polymorphisms associated with the pathology of collagen-induced arthritis (CIA) and anti-collagen IgG induction in a cynomolgus macaque model, as MHC polymorphisms affect the onset of CIA in other animal models. Nine female Filipino cynomolgus macaques were immunized with bovine type II collagen (b-CII) to induce CIA, which was diagnosed clinically by scoring the symptoms of joint swelling over 9 weeks. MHC polymorphisms and anti-b-CII antibody titers were compared between symptomatic and asymptomatic macaques. Four of 9 (44%) macaques were defined as the CIA-affected group. Anti-b-CII IgG in the affected group increased in titer approximately 3 weeks earlier compared with the asymptomatic group. The mean plasma IgG1 titer in the CIA-affected group was significantly higher (p < 0.05) than that of the asymptomatic group. Furthermore, the cynomolgus macaque MHC (Mafa)-DRB1*10:05 or Mafa-DRB1*10:07 alleles, which contain the well-documented RA-susceptibility five amino acid sequence known as the shared epitope (SE) in positions 70 to 74, with valine at position 11 (Val11, V11) and phenylalanine at position 13 (Phe13, F13), were detected in the affected group. In contrast, no MHC polymorphisms specific to the asymptomatic group were identified. In conclusion, the presence of V11 and F13 along with SE in the MHC-DRB1 alleles seems essential for the production of IgG1 and the rapid induction of severe CIA in female Filipino cynomolgus macaques.

Regulation of expression quantitative trait loci by SVA retrotransposons within the major histocompatibility complex

Genomic and transcriptomic studies of expression quantitative trait loci (eQTL) revealed that SINE-VNTR-Alu (SVA) retrotransposon insertion polymorphisms (RIPs) within human genomes markedly affect the co-expression of many coding and noncoding genes by coordinated regulatory processes. This study examined the polymorphic SVA modulation of gene co-expression within the major histocompatibility complex (MHC) genomic region where more than 160 coding genes are involved in innate and adaptive immunity. We characterized the modulation of SVA RIPs utilizing the genomic and transcriptomic sequencing data obtained from whole blood of 1266 individuals in the Parkinson's Progression Markers Initiative (PPMI) cohort that included an analysis of human leukocyte antigen (HLA)-A regulation in a subpopulation of the cohort. The regulatory properties of eight SVAs located within the class I and class II MHC regions were associated with differential co-expression of 71 different genes within and 75 genes outside the MHC region. Some of the same genes were affected by two or more different SVA. Five SVA are annotated in the human genomic reference sequence GRCh38.p14/hg38, whereas the other three were novel insertions within individuals. We also examined and found distinct structural effects (long and short variants and the CT internal variants) for one of the SVA (R_SVA_24) insertions on the differential expression of the HLA-A gene within a subpopulation (550 individuals) of the PPMI cohort. This is the first time that many HLA and non-HLA genes (multilocus expression units) and splicing mechanisms have been shown to be regulated by eight structurally polymorphic SVA within the MHC genomic region by applying precise statistical analysis of RNA data derived from the blood samples of a human cohort population. This study shows that SVA within the MHC region are important regulators or rheostats of gene co-expression that might have potential roles in diversity, health, and disease.

Large-Scale Polymorphism Analysis of Dog Leukocyte Antigen Class I and Class II Genes (DLA-88, DLA-12/88L and DLA-DRB1) and Comparison of the Haplotype Diversity between Breeds in Japan

Polymorphisms of canine leukocyte antigen (DLA) class I (DLA-88 and DLA-12/88L) and class II (DLA-DRB1) genes are important for disease susceptibility studies, but information on the genetic diversity among dog breeds is still lacking. To better elucidate the polymorphism and genetic diversity between breeds, we genotyped DLA-88, DLA-12/88L, and DLA-DRB1 loci using 829 dogs of 59 breeds in Japan. Genotyping by Sanger sequencing identified 89, 43, and 61 alleles in DLA-88, DLA-12/88L, and DLA-DRB1 loci, respectively, and a total of 131 DLA-88-DLA-12/88L-DLA-DRB1 haplotypes (88-12/88L-DRB1) were detected more than once. Of the 829 dogs, 198 were homozygotes for one of the 52 different 88-12/88L-DRB1 haplotypes (homozygosity rate: 23.8%). Statistical modeling suggests that 90% of the DLA homozygotes or heterozygotes with one or other of the 52 different 88-12/88L-DRB1 haplotypes within somatic stem cell lines would benefit graft outcome after 88-12/88L-DRB1-matched transplantation. As previously reported for DLA class II haplotypes, the diversity of 88-12/88L-DRB1 haplotypes varied remarkably between breeds but was relatively conserved within most breeds. Therefore, the genetic characteristics of high DLA homozygosity rate and poor DLA diversity within a breed are useful for transplantation therapy, but they may affect biological fitness as homozygosity progresses.

Association of immune-mediated necrotizing myopathy with HLA polymorphisms

Immune-mediated necrotizing myopathy (IMNM) is a type of autoimmune myositis typically characterized clinically by proximal muscle weakness with elevated creatine kinase levels, pathologically by myofiber necrosis and regeneration with paucity of lymphocytic cell infiltration, and serologically by the presence of either of two myositis-specific autoantibodies, anti-SRP, and anti-HMGCR antibodies. However, the HLA loci and alleles associated with IMNM are still not fully understood at least partly because IMNM was a relatively recently established condition. In this study, we genotyped the six HLA loci (HLA-A, -B, -C, -DRB1, -DQB1 and -DPB1) in 250 patients (237 patients over age 18 years and 13 juvenile patients) diagnosed with IMNM based on clinicopathological features and autoantibody information and performed a case control study with Japanese healthy subjects. In the adult patients, specific HLA alleles associated with IMNM were identified at all HLA loci, with DRB1*08:03 showing the strongest association (OR = 2.5; p = 0.00000017). Furthermore, subgroup analysis with various clinical information showed that C*03:04 (OR = 3.7; p = 0.00012) was a higher risk allele for collagen disease in adult patients, and B*13:01 (OR = 23.2; p = 0.021) and C*03:04 (OR = 5.8; p = 0.0074) were higher risk for juvenile patients with anti-HMGCR antibody-positive IMNM. These findings will help to better understand the HLA genetic background and features of IMNM in designing future studies.

Genetic Association between Farrowing Rates and Swine Leukocyte Antigen Alleles or Haplotypes in Microminipigs

We have previously reported specific swine leukocyte antigen (SLA) haplotype associations with significant effects on several reproduction performance traits in a highly inbred miniature pig population of Microminipigs (MMPs). In this study, to clarify the effects on farrowing rates of SLA similarity between mating partners in the MMP population, we compared the farrowing rates as a measure of reproductive success after 1063-cumulative matings among the following three groups of mating partners: (1) completely sharing SLA class I or class II haplotypes or alleles between partners (CS), (2) only one sharing the haplotypes or alleles (OS), and (3) non-sharing the haplotypes or alleles (NS). Average farrowing rates in CS groups consisting of completely sharing SLA class II haplotypes or DRBI and DQB1 alleles were lowest in the three groups. Moreover, lower farrowing rates were indicated in mating pairs with smaller amino acid pairwise genetic distances of SLA-1, SLA-3, DRB1 and DQB1 alleles between the pairs. These results suggested that the dissimilarity of SLA class I and class II alleles between mating partners markedly improved reproductive performance; therefore, SLA alleles or haplotypes are potentially useful genetic markers for the selection of mating pairs in breeding programs and epistatic studies of reproductive traits of MMPs.

Sequence Variations Within HLA-G and HLA-F Genomic Segments at the Human Leukocyte Antigen Telomeric End Associated With Acute Graft-Versus-Host Disease in Unrelated Bone Marrow Transplantation

Acute graft-versus-host disease (aGVHD) is defined as a syndrome of an immunological response of graft to the host that occurs early after allogeneic hematopoietic stem cell transplantation (HCT). This disease is frequently observed even in HCT matched for human leukocyte antigen (HLA) alleles at multiple gene loci. Although the HLA region represents complex and diverse genomic characteristics, detailed association analysis is required for the identification of uncharacterized variants that are strongly associated with aGVHD. We genotyped three loci, OR2H2, HLA-F-AS1, and HLA-G, that are located in the 460 kb of HLA telomeric region and statistically analyzed the genotypes including HLA-DPB1 with clinical and transplantation outcomes using 338 unrelated bone marrow transplantation (UR-BMT) patient–donor pairs who were matched for HLA-A, HLA-B, HLA-C, HLA-DRB1, and HLA-DQB1 (HLA-10/10). Multivariate analyses demonstrated that HLA-F-AS1 and HLA-DPB1 mismatches were associated with grade II–IV aGVHD (hazard ratio (HR), 1.76; 95% CI, 1.07–2.88; p = 0.026; and HR, 1.59; CI, 1.02–2.49; p = 0.042, respectively). There was no confounding between HLA-F-AS1 and HLA-DPB1 (p = 0.512), suggesting that the HLA-F-AS1 mismatch has a strong effect on aGVHD independently of HLA-DPB1. Moreover, a stratified analysis suggested possible associations of HLA-F-AS1, HLA-DPB1, and/or HLA-G mismatches with grade II–IV aGVHD and the more severe grade III–IV aGVHD. These findings provide new insights into understanding the molecular mechanism of aGVHD caused by HLA-matched UR-BMT.

A novel swab storage gel is superior to dry swab DNA collection, and enables long-range high resolution NGS HLA typing from buccal cell samples

HLA sequence-based DNA typing (SBT) by long-range PCR amplification (LR PCR) and next-generation sequencing (NGS) is a high-throughput DNA sequencing method (LR-NGS-SBT) for the efficient and sensitive detection of novel and null HLA alleles to the field-4 level of allelic resolution without phase ambiguity. However, the accuracy and reliability of the HLA typing results using buccal cells (BCs) and saliva as genetic source materials for the LR-NGS-SBT method are dependent largely on the quality of the extracted genomic DNA (gDNA) because a large degree of gDNA fragmentation can result in insufficient PCR amplification with the incorrect assignment of HLA alleles due to allele dropouts. In this study, we developed a new cost-efficient swab storage gel (SSG) for wet swab collection of BCs (BC-SSG) and evaluated its usefulness by performing different DNA analytical parameters including LR-NGS-SBT to compare the quality and quantity of gDNA extracted from BCs (in SSG or air dried), blood and saliva of 30 subjects. The BC-SSG samples after 5 days of storage revealed qualitative and quantitative DNA values equivalent to that of blood and/or saliva and better than swabs that were only air-dried (BC-nSSG). Moreover, all the gDNA extracted from blood, saliva and BC-SSG samples were HLA-typed successfully to an equivalent total of 408 alleles for each sample type. Therefore, the application of BC-SSG collection media for LR-NGS-SBT has benefits over BC dried samples (dry swabs) such as reducing retesting and the number of untestable BC samples due to insufficient DNA amplification. This article is protected by copyright. All rights reserved.

Haplotype structures and polymorphisms of dog leukocyte antigen (DLA) class I loci shaped by intralocus and interlocus recombination events

The dog leukocyte antigen (DLA) class I genomic region is located on chromosome 12, and the class I genomic region is composed of at least two distinct haplotypic gene structures, DLA-88-DLA-12 and DLA-88-DLA-88L. However, detailed information of the genomic differences among DLA-88, DLA-12, and DLA-88L are still lacking at the full-length gene level, and therefore, DLA allelic sequences classified for each of these loci are limited in number so far. In this study, we determined the DNA sequence of a 95-kb DLA class I genomic region including DLA-88, DLA-12/88L, and DLA-64 with three DLA homozygous dogs and of 37 full-length allelic gene sequences for DLA-88 and DLA-12/88L loci in 26 DLA class I homozygous dogs. Nucleotide diversity profiles of the 95-kb regions and sequence identity scores of the allelic sequences suggested that DLA-88L is a hybrid gene generated by interlocus and/or intralocus gene conversion between DLA-88 and DLA-12. The putative minimum conversion tract was estimated to be at least an 850-bp segment in length located from the 5´flanking untranslated region to the end of intron 2. In addition, at least one DLA-12 allele (DLA-12*004:01) was newly generated by interlocus gene conversion. In conclusion, the analysis for the occurrence of gene conversion within the dog DLA class I region revealed intralocus gene conversion tracts in 17 of 27 DLA-88 alleles and two of 10 DLA-12 alleles, suggesting that intralocus gene conversion has played an important role in expanding DLA allelic variations.

Identification of three novel HLA-DRA alleles by next-generation sequencing

Three novel HLA-DRA alleles, DRA*01:03, DRA*01:04, and DRA*01:05 alleles with unique amino acid sequences.

Haplotypic Associations and Differentiation of MHC Class II Polymorphic Alu Insertions at Five Loci With HLA-DRB1 Alleles in 12 Minority Ethnic Populations in China

The analysis of polymorphic variations in the human major histocompatibility complex (MHC) class II genomic region on the short-arm of chromosome 6 is a scientific enquiry to better understand the diversity in population structure and the effects of evolutionary processes such as recombination, mutation, genetic drift, demographic history, and natural selection. In order to investigate associations between the polymorphisms of HLA-DRB1 gene and recent Alu insertions (POALINs) in the HLA class II region, we genotyped HLA-DRB1 and five Alu loci (AluDPB2, AluDQA2, AluDQA1, AluDRB1, AluORF10), and determined their allele frequencies and haplotypic associations in 12 minority ethnic populations in China. There were 42 different HLA-DRB1 alleles for ethnic Chinese ranging from 12 alleles in the Jinuo to 28 in the Yugur with only DRB1^∗08:03, DRB1^∗09:01, DRB1^∗12:02, DRB1^∗14:01, DRB1^∗15:01, and DRB1^∗15:02 present in all ethnic groups. The POALINs varied in frequency between 0.279 and 0.514 for AluDPB2, 0 and 0.127 for AluDQA2, 0.777 and 0.995 for AluDQA1, 0.1 and 0.455 for AluDRB1 and 0.084 and 0.368 for AluORF10. By comparing the data of the five-loci POALIN in 13 Chinese ethnic populations (including Han-Yunnan published data) against Japanese and Caucasian published data, marked differences were observed between the populations at the allelic or haplotypic levels. Five POALIN loci were in significant linkage disequilibrium with HLA-DRB1 in different populations and AluDQA1 had the highest percentage association with most of the HLA-DRB1 alleles, whereas the nearby AluDRB1 indel was strongly haplotypic for only DRB1^∗01, DRB1^∗10, DRB1^∗15 and DRB1^∗16. There were 30 five-locus POALIN haplotypes inferred in all populations with H5 (no Alu insertions except for AluDQA1) and H21 (only AluDPB2 and AluDQA1 insertions) as the two predominant haplotypes. Neighbor joining trees and principal component analyses of the Alu and HLA-DRB1 polymorphisms showed that genetic diversity of these genomic markers is associated strongly with the population characteristics of language family, migration and sociality. This comparative study of HLA-DRB1 alleles and multilocus, lineage POALIN frequencies of Chinese ethnic populations confirmed that POALINs whether investigated alone or together with the HLA class II alleles are informative genetic and evolutionary markers for the identification of allele and haplotype lineages and genetic variations within the same and/or different populations.

Stillbirth rates and their association with swine leucocyte antigen class II haplotypes in Microminipigs

Objective

Microminipig (MMP) is a miniature pig with an extra small body size for experimental use. In the present study, the occurrence of stillbirths and their genetic association with swine leukocyte antigen (SLA) class II haplotypes were evaluated in a population of MMPs.

Methods

The occurrences of stillbirth and genetic association with SLA class II haplotypes using 483 stillborn and 2,246 live piglets, and their parents were compared among the three groups of newborn piglet litters; an all stillborn (AS) group consisting of only stillborn piglet litters, a partial stillborn (PS) group consisting of stillborn and live piglet litters, and an all alive (AA) group consisting of only live piglet litters.

Results

The incidence of stillborn piglets was 483/2,729 (17.7%). Distributions of litter sizes, numbers of stillborn piglets in a litter, parities, and gestation periods were distinct among the three groups. The frequencies of low resolution haplotype (Lr)-0.7 or Lr-0.23 were higher in the AS group than in the PS or AA groups. In sires, the frequency of Lr-0.7 associated with the AS group was significantly higher in the AS group than with the AA group. In dams, the frequency of Lr-0.23 was significantly higher in the AS group than in the PS or AA groups, whereas the frequency of Lr-0.7 was not significantly different.

Conclusion

The incidence of stillborn piglets in MMPs appears to be higher than those in other pig breeds. Several traits related with stillbirths such as the number of stillborn piglets and parities of the AS group were different from those of the PS and AA groups. Specific SLA class II haplotypes were associated significantly with a high incidence of stillbirths and could be used as genetic markers to adopt breeding strategies to lower the rate of stillbirth in MMPs.

Haplotype Shuffling and Dimorphic Transposable Elements in the Human Extended Major Histocompatibility Complex Class II Region

The major histocompatibility complex (MHC) on chromosome 6p21 is one of the most single-nucleotide polymorphism (SNP)-dense regions of the human genome and a prime model for the study and understanding of conserved sequence polymorphisms and structural diversity of ancestral haplotypes/conserved extended haplotypes. This study aimed to follow up on a previous analysis of the MHC class I region by using the same set of 95 MHC haplotype sequences downloaded from a publicly available BioProject database at the National Center for Biotechnology Information to identify and characterize the polymorphic human leukocyte antigen (HLA)-class II genes, the MTCO3P1 pseudogene alleles, the indels of transposable elements as haplotypic lineage markers, and SNP-density crossover (XO) loci at haplotype junctions in DNA sequence alignments of different haplotypes across the extended class II region (∼1 Mb) from the telomeric PRRT1 gene in class III to the COL11A2 gene at the centromeric end of class II. We identified 42 haplotypic indels (20 Alu, 7 SVA, 13 LTR or MERs, and 2 indels composed of a mosaic of different transposable elements) linked to particular HLA-class II alleles. Comparative sequence analyses of 136 haplotype pairs revealed 98 unique XO sites between SNP-poor and SNP-rich genomic segments with considerable haplotype shuffling located in the proximity of putative recombination hotspots. The majority of XO sites occurred across various regions including in the vicinity of MTCO3P1 between HLA-DQB1 and HLA-DQB3, between HLA-DQB2 and HLA-DOB, between DOB and TAP2, and between HLA-DOA and HLA-DPA1, where most XOs were within a HERVK22 sequence. We also determined the genomic positions of the PRDM9-recombination suppression sequence motif ATCCATG/CATGGAT and the PRDM9 recombination activation partial binding motif CCTCCCCT/AGGGGAG in the class II region of the human reference genome (NC_ 000006) relative to published meiotic recombination positions. Both the recombination and anti-recombination PRDM9 binding motifs were widely distributed throughout the class II genomic regions with 50% or more found within repeat elements; the anti-recombination motifs were found mostly in L1 fragmented repeats. This study shows substantial haplotype shuffling between different polymorphic blocks and confirms the presence of numerous putative ancestral recombination sites across the class II region between various HLA class II genes.

SNP-Density Crossover Maps of Polymorphic Transposable Elements and HLA Genes Within MHC Class I Haplotype Blocks and Junction

The genomic region (~4 Mb) of the human major histocompatibility complex (MHC) on chromosome 6p21 is a prime model for the study and understanding of conserved polymorphic sequences (CPSs) and structural diversity of ancestral haplotypes (AHs)/conserved extended haplotypes (CEHs). The aim of this study was to use a set of 95 MHC genomic sequences downloaded from a publicly available BioProject database at NCBI to identify and characterise polymorphic human leukocyte antigen (HLA) class I genes and pseudogenes, MICA and MICB, and retroelement indels as haplotypic lineage markers, and single-nucleotide polymorphism (SNP) crossover loci in DNA sequence alignments of different haplotypes across the Olfactory Receptor (OR) gene region (~1.2 Mb) and the MHC class I region (~1.8 Mb) from the GPX5 to the MICB gene. Our comparative sequence analyses confirmed the identity of 12 haplotypic retroelement markers and revealed that they partitioned the HLA-A/B/C haplotypes into distinct evolutionary lineages. Crossovers between SNP-poor and SNP-rich regions defined the sequence range of haplotype blocks, and many of these crossover junctions occurred within particular transposable elements, lncRNA, OR12D2, MUC21, MUC22, PSORS1A3, HLA-C, HLA-B, and MICA. In a comparison of more than 250 paired sequence alignments, at least 38 SNP-density crossover sites were mapped across various regions from GPX5 to MICB. In a homology comparison of 16 different haplotypes, seven CEH/AH (7.1, 8.1, 18.2, 51.x, 57.1, 62.x, and 62.1) had no detectable SNP-density crossover junctions and were SNP poor across the entire ~2.8 Mb of sequence alignments. Of the analyses between different recombinant haplotypes, more than half of them had SNP crossovers within 10 kb of LTR16B/ERV3-16A3_I, MLT1, Charlie, and/or THE1 sequences and were in close vicinity to structurally polymorphic Alu and SVA insertion sites. These studies demonstrate that (1) SNP-density crossovers are associated with putative ancestral recombination sites that are widely spread across the MHC class I genomic region from at least the telomeric OR12D2 gene to the centromeric MICB gene and (2) the genomic sequences of MHC homozygous cell lines are useful for analysing haplotype blocks, ancestral haplotypic landscapes and markers, CPSs, and SNP-density crossover junctions.

Preparation and characterization of monoclonal antibodies recognizing two CD4 isotypes of Microminipigs

Cluster of differentiation 4 (CD4) molecule expressed on the leukocytes is known to function as a co-receptor for class II major histocompatibility complex (MHC) binding to T cell receptor (TCR) on helper T cells. We previously identified two CD4 alleles (CD4.A and CD4.B) in a Microminipig population based on nucleotide sequencing and PCR detection of their gene sequences. However, CD4.B protein expression was not examined because of the unavailability of a reactive antibody to a CD4.B epitope. In this study, we have produced two swine-specific monoclonal antibodies (mAbs) against CD4.B molecules, one that recognizes only CD4.B (b1D7) and the other that recognizes both the CD4.A and CD4.B alleles (x1E10) and that can be used to distinguish CD4 T cell subsets by flow cytometry and immunohistochemistry. Using these two mAbs, we identified CD4.A and CD4.B allele-specific proteins on the surface of CD4.A (+/+) and CD4.B (+/+) T cells at a similar level of expression. Moreover, stimulation of peripheral blood mononuclear cells (PBMCs) derived from CD4.A (+/+) and CD4.B (+/+) swine with toxic shock syndrome toxin-1 (TSST-1) in vitro similarly activated both groups of cells that exhibited a slight increase in the CD4/CD8 double positive (DP) cell ratio. A large portion of the DP cells from the allelic CD4.A (+/+) and CD4.B (+/+) groups enhanced the total CD4 and class I swine leukocyte antigen (SLA) expression. The x1E10 mAb delayed and reduced the TSST-1-induced activation of CD4 T cells. Thus, CD4.B appears to be a functional protein whose expression on activated T cells is analogous to CD4.A.

Identification of Novel Alleles and Structural Haplotypes of Major Histocompatibility Complex Class I and DRB Genes in Domestic Cat (Felis catus) by a Newly Developed NGS-Based Genotyping Method

The major histocompatibility complex (MHC) is a highly polymorphic and duplicated genomic region that encodes transplantation and immune regulatory molecules. Although it is well-known that particular MHC allelic polymorphisms and haplotypes are genetically relate to immune-mediated diseases detailed information of the cat MHC (Feline Leukocyte Antigen; FLA) genetic and haplotypic structure and diversity is limited in comparison to humans and many other species. In this study, to better understand the degree and types of allele and allelic haplotype diversity of FLA-class I (FLA-I) and FLA-DRB loci in domestic cats, we identified six expressible FLA-I loci in peripheral white blood cells by in silico estimation of the coding exons and NGS-based amplicon sequencing using five unrelated cats. We then used a newly developed NGS-based genotyping method to genotype and annotate 32 FLA-I and 16 FLA-DRB sequences in two families of 20 domestic cats. A total of 14 FLA-I and seven FLA-DRB were identified as novel polymorphic sequences. Phylogenetic analyses grouped the sequences into six FLA-I (FLA-E/H/K, FLA-A, FLA-J, FLA-L, FLA-O and a tentatively named FLA-E/H/K_Rec) and four FLA-DRB (FLA-DRB1, FLA-DRB3, FLA-DRB4, and FLA-DRB5) lineages. Pedigree analysis of two cat families revealed eight distinct FLA structural haplotypes (Class I - DRB) with five to eight FLA-I and two to three FLA-DRB transcribed loci per haplotype. It is evident that the eight FLA haplotypes were generated by gene duplications and deletions, and rearrangements by genetic recombination with the accumulation and/or inheritance of novel polymorphisms. These findings are useful for further genetic diversity analysis and disease association studies among cat breeds and in veterinary medicine.

Capturing Differential Allele-Level Expression and Genotypes of All Classical HLA Loci and Haplotypes by a New Capture RNA-Seq Method

The highly polymorphic human major histocompatibility complex (MHC) also known as the human leukocyte antigen (HLA) encodes class I and II genes that are the cornerstone of the adaptive immune system. Their unique diversity (>25,000 alleles) might affect the outcome of any transplant, infection, and susceptibility to autoimmune diseases. The recent rapid development of new next-generation sequencing (NGS) methods provides the opportunity to study the influence/correlation of this high level of HLA diversity on allele expression levels in health and disease. Here, we describe the NGS capture RNA-Seq method that we developed for genotyping all 12 classical HLA loci (HLA-A, HLA-B, HLA-C, HLA-DPA1, HLA-DPB1, HLA-DQA1, HLA-DQB1, HLA-DRA, HLA-DRB1, HLA-DRB3, HLA-DRB4, and HLA-DRB5) and assessing their allelic imbalance by quantifying their allele RNA levels. This is a target enrichment method where total RNA is converted to a sequencing-ready complementary DNA (cDNA) library and hybridized to a complex pool of RNA-specific HLA biotinylated oligonucleotide capture probes, prior to NGS. This method was applied to 161 peripheral blood mononuclear cells and 48 umbilical cord blood cells of healthy donors. The differential allelic expression of 10 HLA loci (except for HLA-DRA and HLA-DPA1) showed strong significant differences (P < 2.1 × 10-15). The results were corroborated by independent methods. This newly developed NGS method could be applied to a wide range of biological and medical questions including graft rejections and HLA-related diseases.

Genetic Association between Swine Leukocyte Antigen Class II Haplotypes and Reproduction Traits in Microminipigs

The effects of swine leukocyte antigen (SLA) molecules on numerous production and reproduction performance traits have been mainly reported as associations with specific SLA haplotypes that were assigned using serological typing methods. In this study, we intended to clarify the association between SLA class II genes and reproductive traits in a highly inbred population of 187 Microminipigs (MMP), that have eight different types of SLA class II haplotypes. In doing so, we compared the reproductive performances, such as fertility index, gestation period, litter size, and number of stillbirth among SLA class II low resolution haplotypes (Lrs) that were assigned by a polymerase chain reaction-sequence specific primers (PCR-SSP) typing method. Only low resolution haplotypes were used in this study because the eight SLA class II high-resolution haplotypes had been assigned to the 14 parents or the progenitors of the highly inbred MMP herd in a previous publication. The fertility index of dams with Lr-0.13 was significantly lower than that of dams with Lr-0.16, Lr-0.17, Lr-0.18, or Lr-0.37. Dams with Lr-0.23 had significantly smaller litter size at birth than those with Lr-0.17, Lr-0.18, or Lr-0.37. Furthermore, litter size at weaning of dams with Lr-0.23 was also significantly smaller than those dams with Lr-0.16, Lr-0.17, Lr-0.18, or Lr-0.37. The small litter size of dams with Lr-0.23 correlated with the smaller body sizes of these MMPs. These results suggest that SLA class II haplotypes are useful differential genetic markers for further haplotypic and epistatic studies of reproductive traits, selective breeding programs, and improvements in the production and reproduction performances of MMPs.

Long Noncoding RNA HCP5, a Hybrid HLA Class I Endogenous Retroviral Gene: Structure, Expression, and Disease Associations

The HCP5 RNA gene (NCBI ID: 10866) is located centromeric of the HLA-B gene and between the MICA and MICB genes within the major histocompatibility complex (MHC) class I region. It is a human species-specific gene that codes for a long noncoding RNA (lncRNA), composed mostly of an ancient ancestral endogenous antisense 3′ long terminal repeat (LTR, and part of the internal pol antisense sequence of endogenous retrovirus (ERV) type 16 linked to a human leukocyte antigen (HLA) class I promoter and leader sequence at the 5′-end. Since its discovery in 1993, many disease association and gene expression studies have shown that HCP5 is a regulatory lncRNA involved in adaptive and innate immune responses and associated with the promotion of some autoimmune diseases and cancers. The gene sequence acts as a genomic anchor point for binding transcription factors, enhancers, and chromatin remodeling enzymes in the regulation of transcription and chromatin folding. The HCP5 antisense retroviral transcript also interacts with regulatory microRNA and immune and cellular checkpoints in cancers suggesting its potential as a drug target for novel antitumor therapeutics.

MHC class I polymorphic Alu insertion (POALIN) allele and haplotype frequencies in the Arabs of the United Arab Emirates and other world populations

Polymorphic Alu insertions (POALINs) are found throughout the human genome and have been used in various studies to infer geographic origin of human populations. The main aim of this study was to determine the allele and haplotype frequencies of five POALINs, AluHF, AluHG, AluHJ, AluTF and AluMICB, within the major histocompatibility complex (MHC) class I region of 95 UAE Arabs, and correlate their frequencies to those of the HLA-A, HLA-C and HLA-B class I allele lineages. Evolutionary relationships between the POALINs of the Arabs and those previously studied in populations of African, Asian and European descent were compared. At each of the five Alu loci (AluHF, AluHG, AluHJ, AluTF and AluMICB), Alu insertion was designated as Alu(locus)*02 and absence was Alu(locus)*01. The AluHG insertion (AluHG*02) had the highest frequency (0.332), followed by AluHF*02 (0.300), AluHJ*02 (0.263), AluMICB*02 (0.111) and AluTF*02 (0.058). Of the 270 Alu-HLA haplotypes pairs in the UAE Arabs, 110 had no Alu insertion, and 54 had an Alu insertion at >50% per haplotype. An Alu insertion >75% per haplotype was found between AluMICB*02 and HLA-B*14, HLA-B*22, HLA-B*44, HLA-B*55, HLA-B*57 and HLA-B*73, and with HLA-C*01 and HLA-C*18; AluHJ*02 with HLA-A*01, HLA-A*19, HLA-A*24 and HLA-A*32; AluHG*02 with HLA-A*02 and HLA-B*18; and AluHF*02 with HLA-A*10. The genotyped allele and haplotype frequencies of the MHC POALINs in UAE Arabs were compared with the results of 30 previously published Asian, European, American and African populations. Phylogenetic and multidimensional scaling (MDS) analysis of the relative MHC POALINs allele and haplotype frequencies revealed that the UAE Arabs have a similar lineage to Caucasians and the most distant genetic relationship to the Waorani native American population of Ecuador. The structure of both the phylogenetic tree and the MDS analysis supports the Out of Africa theory of human evolution. The nature of the clusters suggests the Arabian Middle East represents a crossroads from which human populations migrated towards Asia in the east and Europe to the north-west.

HLA class I allele lineages and haplotype frequencies in Arabs of the United Arab Emirates

The high degree of polymorphism of the HLA system provides suitable genetic markers to study the diversity and migration of different world populations and is beneficial for forensic identification, anthropology, transplantation and disease associations. Although the United Arab Emirates (UAE) population of about nine million people is heterogeneous, information is limited for the HLA class I allele and haplotype frequencies of the Bedouin ethnic group. We performed low-resolution PCR-SSP genotyping of three HLA class I loci at HLA-A, -B and -C for 95 unrelated healthy Bedouins from the cities of Al Ain and Abu Dhabi in the UAE. A total of 54 HLA allele lineages were detected; the most frequent low-resolution allele lineages at each HLA locus were A*02 (0.268), B*51 (0.163) and C*07 (0.216). The inferred estimates for the two most frequent HLA-A and HLA-B haplotypes were HLA-A*02 ~ HLA-B*50 (0.070) and HLA-A*02 ~ HLA-B*51 (0.051), and the most frequent 3-locus haplotype was HLA-A*02 ~ HLA-B*50 ~ HLA-C*06 (0.068). The HLA allele lineage frequencies of the UAE Arabs were compared to those previously reported for 70 other world populations, and a strong genetic similarity was detected between the UAE Arabs and the Saudi Arabians from the west with evidence of a limited gene flow between the UAE Arabs and Pakistani across the Gulf from the east, and the UAE Arabs and Omani from the south of the Gulf Peninsula.

Super High Resolution for Single Molecule-Sequence-Based Typing of Classical HLA Loci Using Ion Torrent PGM

Super high resolution-single molecule-sequence-based typing (SS-SBT) is an HLA DNA typing method to the field 4 level of allelic resolution (formerly known as 8-digit typing) to efficiently detect novel and null alleles without phase ambiguity by combination of long ranged PCR amplification and next-generation sequencing (NGS) technologies. In this chapter, we describe three basic steps, long ranged PCR, NGS, and allele assignment.

Reference Grade Characterization of Polymorphisms in Full-Length HLA Class I and II Genes With Short-Read Sequencing on the ION PGM System and Long-Reads Generated by Single Molecule, Real-Time Sequencing on the PacBio Platform

Although NGS technologies fuel advances in high-throughput HLA genotyping methods for identification and classification of HLA genes to assist with precision medicine efforts in disease and transplantation, the efficiency of these methods are impeded by the absence of adequately-characterized high-frequency HLA allele reference sequence databases for the highly polymorphic HLA gene system. Here, we report on producing a comprehensive collection of full-length HLA allele sequences for eight classical HLA loci found in the Japanese population. We augmented the second-generation short read data generated by the Ion Torrent technology with long amplicon spanning consensus reads delivered by the third-generation SMRT sequencing method to create reference grade high-quality sequences of HLA class I and II gene alleles resolved at the genomic coding and non-coding level. Forty-six DNAs were obtained from a reference set used previously to establish the HLA allele frequency data in Japanese subjects. The samples included alleles with a collective allele frequency in the Japanese population of more than 99.2%. The HLA loci were independently amplified by long-range PCR using previously designed HLA-locus specific primers and subsequently sequenced using SMRT and Ion PGM sequencers. The mapped long and short-reads were used to produce a reference library of consensus HLA allelic sequences with the help of the reference-aware software tool LAA for SMRT Sequencing. A total of 253 distinct alleles were determined for 46 healthy subjects. Of them, 137 were novel alleles: 101 SNVs and/or indels and 36 extended alleles at a partial or full-length level. Comparing the HLA sequences from the perspective of nucleotide diversity revealed that HLA-DRB1 was the most divergent among the eight HLA genes, and that the HLA-DPB1 gene sequences diverged into two distinct groups, DP2 and DP5, with evidence of independent polymorphisms generated in exon 2. We also identified two specific intronic variations in HLA-DRB1 that might be involved in rheumatoid arthritis. In conclusion, full-length HLA allele sequencing by third-generation and second-generation technologies has provided polymorphic gene reference sequences at a genomic allelic resolution including allelic variations assigned up to the field-4 level for a stronger foundation in precision medicine and HLA-related disease and transplantation studies.

Production of a Locus- and Allele-Specific Monoclonal Antibody for the Characterization of SLA-1*0401 mRNA and Protein Expression Levels in MHC-Defined Microminipigs

The class I major histocompatibility complex (MHC) presents self-developed peptides to specific T cells to induce cytotoxity against infection. The MHC proteins are encoded by multiple loci that express numerous alleles to preserve the variability of the antigen-presenting ability in each species. The mechanism regulating MHC mRNA and protein expression at each locus is difficult to analyze because of the structural and sequence similarities between alleles. In this study, we examined the correlation between the mRNA and surface protein expression of swine leukocyte antigen (SLA)-1*0401 after the stimulation of peripheral blood mononuclear cells (PBMCs) by Staphylococcus aureus superantigen toxic shock syndrome toxin-1 (TSST-1). We prepared a monoclonal antibody (mAb) against a domain composed of Y102, L103 and L109 in the α2 domain. The Hp-16.0 haplotype swine possess only SLA-1*0401, which has the mAb epitope, while other haplotypes possess 0 to 3 SLA classical class I loci with the mAb epitopes. When PBMCs from SLA-1*0401 homozygous pigs were stimulated, the SLA-1*0401 mRNA expression level increased until 24 hrs and decreased at 48 hrs. The kinetics of the interferon regulatory transcription factor-1 (IRF-1) mRNA level were similar to those of the SLA-1*0401 mRNA. However, the surface protein expression level continued to increase until 72 hrs. Similar results were observed in the Hp-10.0 pigs with three mAb epitopes. These results suggest that TSST-1 stimulation induced both mRNA and surface protein expression of class I SLA in the swine PBMCs differentially and that the surface protein level was sustained independently of mRNA regulation.

Identification and characterization of two CD4 alleles in Microminipigs

Background

We previously identified two phenotypes of CD4+ cells with and without reactions to anti-pig CD4 monoclonal antibodies by flow cytometry in a herd of Microminipigs. In this study, we analyzed the coding sequences of CD4 and certified the expression of CD4 molecules in order to identify the genetic sequence variants responsible for the positive and negative PBMCs reactivity to anti-pig CD4 monoclonal antibodies.

Results

We identified two CD4 alleles, CD4.A and CD4.B, corresponding to antibody positive and negative, respectively, by nucleotide sequencing of PCR products using CD4 specific primer pairs. In comparison with the swine CD4 amino-acid sequence [GenBank: NP_001001908], CD4.A had seven amino-acid substitutions and CD4.B had 15 amino-acid substitutions. The amino-acid sequences within domain 1 of CD4.B were identical to the swine CD4.2 [GenBank: CAA46584] sequence that had been reported previously to be a modified CD4 molecule that had lost reactivity with an anti-pig CD4 antibody in NIH miniature pigs. Homozygous and heterozygous CD4.A and CD4.B alleles in the Microminipigs herd were characterised by using the RFLP technique with the restriction endonuclease, BseRI. The anti-pig CD4 antibody recognized pig PBMCs with CD4.AA and CD4.AB, but did not recognized those with CD4.BB. We transfected HeLa cells with the FLAG-tagged CD4.A or CD4.B vectors, and certified that transfected HeLa cells expressed FLAG in both vectors. The failure of cells to react with anti-CD4 antibodies in CD4.B pigs was associated to ten amino-acid substitutions in domain 1 and/or one amino-acid substitution in joining region 3 of CD4.B. We also found exon 8 was defective in some CD4.A and CD4.B resulting in the loss of the transmembrane domain, which implies that these CD4 proteins are secreted from helper T cells into the circulation.

Conclusions

We identified that amino-acids substitutions of domain 1 in CD4.B gave rise to the failure of some CD4 expressing cells to react with particular anti-pig CD4 monoclonal antibodies. In addition, we developed a PCR-RFLP method that enabled us to simply identify the CD4 sequence variant and the positive and negative PBMCs reactivity to our anti-pig CD4 monoclonal antibodies without the need to use flow cytometric analysis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12917-016-0856-8) contains supplementary material, which is available to authorized users.

Comparative genomics of the human, macaque and mouse major histocompatibility complex

The MHC is a highly polymorphic genomic region that encodes the transplantation and immune regulatory molecules. It receives special attention for genetic investigation because of its important role in the regulation of innate and adaptive immune responses and its strong association with numerous infectious and/or autoimmune diseases. The MHC locus was first discovered in the mouse and for the past 50 years it has been studied most intensively in both mice and humans. However, in recent years the macaque species have emerged as some of the more important and advanced experimental animal models for biomedical research into MHC with important human immunodeficiency virus/simian immunodeficiency virus and transplantation studies undertaken in association with precise MHC genotyping and haplotyping methods using Sanger sequencing and next-generation sequencing. Here, in this special issue on 'Macaque Immunology' we provide a short review of the genomic similarities and differences among the human, macaque and mouse MHC class I and class II regions, with an emphasis on the association of the macaque class I region with MHC polymorphism, haplotype structure and function.

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.

The transcript repeat element: the human Alu sequence as a component of gene networks influencing cancer

A small percentage (3%) of the 1.3 million copies of Alu sequences in the human genome is expressed individually or as part of various gene transcripts with potential regulatory and pathophysiological importance. In order to better understand the role of repetitive elements within transcripts, this review focuses on Alu-containing transcripts of normal and cancerous tissue in a transcriptome-wide survey of the H-Invitational human transcript database on 106,825 tissue-derived transcripts expressed at 29,979 loci. The Alu elements in transcripts of cancerous tissues are significantly underrepresented in comparison to those in normal tissues. In this review, we propose a model for Alu-mediated siRNA down-regulation of Alu-containing transcripts in cancer tissues. In cancer or other rapidly dividing tissues, hypomethylation of repeat element regions triggers the expression of transposon elements including Alu, which can potentially form double-stranded RNA molecules for use as templates to generate Alu-derived siRNAs (Alu-siRNAs). The generated Alu-siRNAs target endogenous messenger RNAs harbouring sequence similarity to Alu elements. This model correlates with the observation that there is substantial under-representation of Alu-containing mRNAs in cancer cells. This new perspective of gene regulation in disease conditions can provide a basis for starting to account for changes in complex gene network in cancer.

Differentiation ability of multipotent hematopoietic stem/progenitor cells detected by a porcine specific anti-CD117 monoclonal antibody

CD117 is a cytokine receptor expressed on the surface of hematopoietic stem cells with a likely role in cell survival, proliferation and differentiation. In order to study the differentiation activity of porcine CD117 hematopoietic cells in vitro and in vivo we prepared an anti-swine CD117 Mab (2A1) with high specificity for flow-cytometrical analysis. The 2A1 Mab did not recognize mouse or human mast cells suggesting that 2A1 is species-specific. Swine bone marrow (BM) CD117+ cells differentiated in vitro mainly into erythroid and monocyte lineages in the methylcellulose-based colony assay. When the swine BM CD117+ cells were transplanted in vivo into immunodeficient NOG (NOD/SCID/IL-2gc-null) mice, a significant amount of swine CD45+ leukocytes, including CD3 positive T cells, were developed in the mice. These results revealed that the swine BM CD117+ cells possess hematopoietic stem/progenitor activity and when monitored in immunodeficient mice or in vitro they can develop into lymphoid, erythroid, and myeloid cells efficiently with the new monoclonal antibody.

Cost-efficient multiplex PCR for routine genotyping of up to nine classical HLA loci in a single analytical run of multiple samples by next generation sequencing

Background

HLA genotyping by next generation sequencing (NGS) requires three basic steps, PCR, NGS, and allele assignment. Compared to the conventional methods, such as PCR-sequence specific oligonucleotide primers (SSOP) and -sequence based typing (SBT), PCR-NGS is extremely labor intensive and time consuming. In order to simplify and accelerate the NGS-based HLA genotyping method for multiple DNA samples, we developed and evaluated four multiplex PCR methods for genotyping up to nine classical HLA loci including HLA-A, HLA-B, HLA-C, HLA-DRB1/3/4/5, HLA-DQB1, and HLA-DPB1.

Results

We developed multiplex PCR methods using newly and previously designed middle ranged PCR primer sets for genotyping different combinations of HLA loci, (1) HLA-DRB1/3/4/5, (2) HLA-DQB1 (3.8 kb to 5.3 kb), (3) HLA-A, HLA-B, HLA-C, and (4) HLA-DPB1 (4.6 kb to 7.2 kb). The primer sets were designed to genotype polymorphic exons to the field 3 level or 6-digit typing. When we evaluated the PCR method for genotyping all nine HLA loci (9LOCI) using 46 Japanese reference subjects who represented a distribution of more than 99.5% of the HLA alleles at each of the nine HLA loci, all of the 276 alleles genotyped, except for HLA-DRB3/4/5 alleles, were consistent with known alleles assigned by the conventional methods together with relevant locus balance and no excessive allelic imbalance. One multiplex PCR method (9LOCI) was able to provide precise genotyping data even when only 1 ng of genomic DNA was used for the PCR as a sample template.

Conclusions

In this study, we have demonstrated that the multiplex PCR approach for NGS-based HLA genotyping could serve as an alternative routine HLA genotyping method, possibly replacing the conventional methods by providing an accelerated yet robust amplification step. The method also could provide significant merits for clinical applications with its ability to amplify lower quantity of samples and the cost-saving factors.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1514-4) contains supplementary material, which is available to authorized users.

Association of sick building syndrome with neuropathy target esterase (NTE) activity in Japanese

Sick building syndrome (SBS) is a set of several clinically recognizable symptoms reported by occupants of a building without a clear cause. Neuropathy target esterase (NTE) is a membrane bound serine esterase and its reaction with organophosphates (OPs) can lead to OP-induced delayed neuropathy (OPIDN) and nerve axon degeneration. The aim of our study was to determine whether there was a difference in NTE activity in the peripheral blood mononuclear cells (PBMCs) of Japanese patients with SBS and healthy controls and whether PNPLA6 (alias NTE) gene polymorphisms were associated with SBS. We found that the enzymatic activity of NTE was significantly higher (P < 0.0005) in SBS patients compared with controls. Moreover, population with an AA genotype of a single nucleotide polymorphism (SNP), rs480208, in intron 21 of the PNPLA6 gene strongly reduced the activity of NTE. Fifty-eight SNP markers within the PNPLA6 gene were tested for association in a case-control study of 188 affected individuals and 401 age-matched controls. Only one SNP, rs480208, was statistically different in genotype distribution (P = 0.005) and allele frequency (P = 0.006) between the cases and controls (uncorrected for testing multiple SNP sites), but these were not significant by multiple corrections. The findings of the association between the enzymatic activity of NTE and SBS in Japanese show for the first time that NTE activity might be involved with SBS.

Multiple deletions in mitochondrial DNA in a patient with progressive external ophthalmoplegia, leukoencephalopathy and hypogonadism

Progressive external ophthalmoplegia (PEO) is one of a number of major types of mitochondrial disorders. Most sporadic PEO patients have a heteroplasmic large deletion of mitochondrial DNA (mtDNA) in the mitochondria in skeletal muscles. We herein analyzed mtDNA deletions using sub-cloning and Sanger sequencing of PCR products in a 31-year-old Japanese man with multiple symptoms, including PEO, muscle weakness, hearing loss, leukoencephalopathy and hypogonadism. A large number of multiple deletions was detected, as well as four kinds of deletion breakpoints identified in different locations, including m.3347_12322, m.5818_13964, m.5829_13964 and m.5837_13503.

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 and differentiation of MHC class I and II polymorphic Alu insertions and HLA-A, -B, -C and -DRB1 alleles in the Chinese Han population

In order to investigate the polymorphism of Alu insertions (POALINs) in the HLA region, we genotyped ten Alu loci (AluMICB, AluTF, AluHJ, AluHG, AluHF in the HLA class I region and AluDPB2, AluDQA2, AluDQA1, AluDRB1, AluORF10 in the HLA class II region) to determine their allele frequencies and associations with the HLA-A, HLA-B, HLA-C and HLA-DRB1 genes in the Chinese Han population. Our results showed the ten-loci POALINs varied in frequency between 0.003 and 0.425. By comparing the data of the ten-loci POALIN in Chinese Han with Japanese and Caucasian data, marked differences were observed between the three ethnic groups at the allelic or haplotypic levels. Each POALIN was in significant linkage disequilibrium with a variety of HLA-A, -B, -C and -DRB1 alleles, and was associated with a variety of HLA-A, -B, -C and -DRB1 allele in Chinese Han. This comparative study of multilocus POALINs in the HLA class I and II regions of the Chinese Han population shows that POALINs alone or as haplotypes together with the HLA class I and II alleles are informative genetic markers for the identification of HLA class I and II allele and variations, such as crossing over events within the same and/or different populations.

Regulation of CD93 Cell Surface Expression by Protein Kinase C Isoenzymes

Human CD93, also known as complement protein 1, q subcomponent, receptor (C1qRp), is selectively expressed by cells with a myeloid lineage, endothelial cells, platelets, and microglia and was originally reported to be involved in the complement protein 1, q subcomponent (C1q)-mediated enhancement of phagocytosis. The intracellular molecular events responsible for the regulation of its expression on the cell surface, however, have not been determined. In this study, the effect of protein kinases in the regulation of CD93 expression on the cell surface of a human monocyte-like cell line (U937), a human NK-like cell line (KHYG-1), and a human umbilical vein endothelial cell line (HUV-EC-C) was investigated using four types of protein kinase inhibitors, the classical protein kinase C (cPKC) inhibitor Go6976, the novel PKC (nPKC) inhibitor Rottlerin, the protein kinase A (PKA) inhibitor H-89 and the protein tyrosine kinase (PTK) inhibitor herbimycin A at their optimum concentrations for 24 hr. CD93 expression was analyzed using flow cytometry and glutaraldehyde-fixed cellular enzyme-linked immunoassay (EIA) techniques utilizing a CD93 monoclonal antibody (mAb), mNI-11, that was originally established in our laboratory as a CD93 detection probe. The nPKC inhibitor Rottlerin strongly down-regulated CD93 expression on the U937 cells in a dose-dependent manner, whereas the other inhibitors had little or no effect. CD93 expression was down-regulated by Go6976, but not by Rottlerin, in the KHYG-1 cells and by both Rottlerin and Go6976 in the HUV-EC-C cells. The PKC stimulator, phorbol myristate acetate (PMA), strongly up-regulated CD93 expression on the cell surface of all three cell-lines and induced interleukin-8 (IL-8) production by the U937 cells and interferon-gamma (IFN-gamma) production by the KHYG-1 cells. In addition, both Go6976 and Rottlerin inhibited the up-regulation of CD93 expression induced by PMA and IL-8 or IFN-gamma production in the respective cell-lines. Whereas recombinant tumor necrosis factor-alpha (rTNF-alpha) slightly up-regulated CD93 expression on the U937 cells, recombinant interleukin-1beta (rIL-1beta), recombinant interleukin-2 (rIL-2), recombinant interferon-gamma (rIFN-gamma) and lipopolysaccharide (LPS) had no effect. Taken together, these findings indicate that the regulation of CD93 expression on these cells involves the PKC isoenzymes.

Distinct HLA allele and haplotype distributions in four ethnic groups of China

Distinct human leukocyte antigen (HLA) allele and haplotype distributions occur in the northern and southern Han populations of China. However, different ethnic groups in China show limited regional distributions for many HLA alleles and haplotypes. Therefore, it is necessary and meaningful to study the differences in HLA allele and haplotype distribution for northern and southern ethnic groups of China. A total of 428 unrelated individuals from the Lisu, Nu, Tu and Yugur ethnic populations were genotyped for HLA-A, -B, -C and -DRB1 alleles using the PCR-Luminex typing method. The frequencies of HLA alleles and statistically inferred haplotypes were calculated. A total of 29 HLA-A, 54 HLA-B, 27 HLA-C and 41 HLA-DRB1 alleles were spread throughout these four populations with distinct allele and deduced haplotype frequencies between populations. Some alleles and deduced haplotypes exhibited significantly different distributions between northern (Tu and Yugur) and southern groups (Lisu and Nu). A phylogenetic tree and principal component analysis were used to compare the HLA polymorphism between our dataset and 19 other eastern and southeastern Asian populations. This analysis showed that Lisu and Nu belong to a cluster of southern ethnic groups, while Tu and Yugur are most closely related to other northern groups. Thus, distinct ethnic population histories were revealed by analyzing HLA allelic polymorphisms with the HLA profiles of the Lisu and Nu southern Chinese ethnic groups clearly different from the Tu and Yugur northern ethnic groups. The results will be useful for future association studies of infectious disease and contribute toward a more efficient search of organ/tissue matches for transplantation.

Improved loop-mediated isothermal amplification for HLA-DRB1 genotyping using RecA and a restriction enzyme for enhanced amplification specificity

Our aim was to test and develop the use of loop-mediated isothermal amplification (LAMP) for HLA-DRB1 genotyping. Initially, we found that the conventional LAMP protocols produced non-specific and variable amplification results depending on the sample DNA conditions. Experiments with different concentrations of DNase in the reaction mixture with and without T4 DNA ligase-treated samples suggested that the strand displacement activity of DNA polymerase in LAMP, at least in part, started from randomly existing nicks because T4 DNA ligase treatment of sample DNA resulted in no amplification. Such non-specific amplification due to the randomly existing nicks was improved specifically by the addition of RecA of Escherichia coli and a restriction enzyme, for example, PvuII, to the reaction mixture. We applied the modified LAMP (mLAMP) (1) to detect specific HLA-DRB1 alleles by using only specific primers for amplification or (2) for genotyping in multiple samples with a multi-probe typing system. In the latter case, HLA-DRB1 genotyping was developed by combining the mLAMP with amplicon capture using polymorphic region-specific probes fixed onto the bottom of the wells of a 96-well plate and the captured amplicons visualized as a black spot at the bottom of the well. The multi-probe human leukocyte antigen (HLA) typing method and the specific HLA allele detection method could be applied for point-of-care testing due to no requirement for specific and expensive instruments.

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.

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.

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.

The association and differentiation of MHC class I polymorphic Alu insertions and HLA-B/Cw alleles in seven Chinese populations

We investigated polymorphic Alu insertion (POALIN) frequencies 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)-B and -Cw genes in seven different Chinese ethnic populations, the Han, Bulang, Wa, Dai, Maonan, Hani and Jinuo. The POALINs varied in frequency between 0% and 42.3% with significant differences between populations at all of the loci. Each POALIN was in significant linkage disequilibrium with a variety of HLA-B or -Cw four-digit alleles. The percentage association between Alu insertions and the HLA-B or -Cw alleles was calculated in pairwise analyses of haplotypes to show possible crossing over events between loci. The POALIN insertions also helped to further stratify the HLA-B:-Cw haplotypes into different POALIN:HLA-B:HLA-Cw haplotype frequencies. Of the two-locus, five-locus and seven-locus haplotype analyses, the seven-locus haplotypes showed the largest number of differences between the populations. The most common multilocus haplotype in Han was MICB*1:B*4601:Cw*0102:TF*1:HJ*1:HG*2:HF*1 (15.6%) associated with the AluHG insertion, whereas the second most common multilocus haplotype in Han was MICB*1:B*1502:Cw*0801:TF*1:HJ*2:HG*1:HF*1 (11.8%) associated with the AluHJ insertion. This comparative study of multilocus POALINs in the HLA class I region of seven Chinese ethnic populations shows that POALINs alone or together with the HLA class I alleles are informative genetic markers for the identification of HLA class I allele and haplotype lineages and variations such as crossing over events 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.

Flow cytometric identification of CD93 expression on naive T lymphocytes (CD4(+)CD45RA (+) cells) in human neonatal umbilical cord blood

Human CD93 has a molecular weight of about 100 kDa and is selectively expressed by myeloid cell lineages in peripheral blood (PB) mononuclear cells. Although CD93 was initially identified as a receptor for complement component 1, subcomponent q phagocytosis (C1qRp) involved in the C1q-mediated enhancement of the phagocytosis of various antigens, several recent studies have reported that CD93 is not a receptor for the C1q-mediated enhancement of phagocytosis. The expression patterns of CD93 have been previously investigated in PB mononuclear cells (lymphocytes, monocytes, and granulocytes) from adult PB and neonatal umbilical cord blood (UCB), and the expression of CD93 was not found on lymphocytes from either normal adult PB or neonatal UCB. However, the detection of CD93 expression in neonatal UCB using CD93 monoclonal antibodies (mAbs) that recognize different antigenic epitopes remains poorly understood. In this study, we examined the expression of CD93 on lymphocytes, monocytes, and granulocytes from neonatal UCB using four different types of CD93 mAb detection probes, mNI-11, R139, R3, and X-2, using flow cytometric and western blot analyses. We found that CD93, as defined using all four mAbs, was expressed on monocytes and granulocytes in PB mononuclear cells from adult PB and neonatal UCB. On the other hand, we observed for the first time that the expression of CD93 on lymphocytes in neonatal UCB can only be detected using the mNI-11 mAb, established in our laboratory, and not with commercially available CD93 mAbs (R139, R3, and X-2). However, CD93 expression on lymphocytes from normal adults was not detected using any of the four CD93 mAbs. Two-color flow cytometric analyses showed that the CD93 recognized by mNI-11 mAb was expressed on CD3(+) T lymphocytes (mainly CD4(+) helper T lymphocytes), but not on CD19(+) B lymphocytes or on CD8(+) suppressor/cytotoxic T lymphocytes from neonatal UCB. In addition, CD93 was expressed on CD45RA(+) (naive antigen) lymphocytes from neonatal UCB, but not on CD45RO(+) (memory antigen) lymphocytes from neonatal UCB or on CD45RA(+) and CD45RO(+) lymphocytes from normal adult PB. Three-color flow cytometric analysis showed that CD93 was co-expressed on naive T lymphocytes (CD4(+)CD45RA(+)) from neonatal UCB. In a western blot analysis, the CD93 mAb (mNI-11) immunoprecipitated at a molecular weight of 98 kDa, identified as a CD93 molecule, in the CD4(+)CD45RA(+) cells from neonatal UCB but not from adult PB, similar to the results in the human monocyte-like cell line U937 (human CD93-positive cells). Taken together, these results provide the first direct evidence of a novel/naive cell population (CD4(+)CD45RA(+)CD93(+)) in neonatal UCB that may have an important role in cell biology, transplantation, and immature/mature immune responses.

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.