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.

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.

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.

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.

Identification of novel polymorphisms and two distinct haplotype structures in dog leukocyte antigen class I genes: DLA-88, DLA-12 and DLA-64.

The current information on the polymorphism variation and haplotype structure of the domestic dog leukocyte antigen (DLA) genes is limited in comparison to other experimental animals. In this paper, to better elucidate the degree and types of polymorphisms and genetic differences for DLA-88, DLA-12 and DLA-64, we genotyped four families of 38 beagles and another 404 unrelated dogs representing 49 breeds by RT-PCR based Sanger sequencing. We also sequenced and analyzed the genomic organization of the DLA-88 and DLA-12 gene segments to better define these two-gene DLA haplotypes more precisely. We identified 45 alleles for DLA-88, 15 for DLA-12 and six for DLA-64, of which 20, 14 and six, respectively, were newly described alleles. Therefore, this study shows that the DLA-12 and DLA-64 loci are far more polymorphic than previously reported. Phylogenetic analysis strongly supported that the DLA-88, DLA-12 and DLA-64 alleles were independently generated after the original divergence of the DLA-79 alleles. Two distinct DLA-88 and DLA-12 haplotype structures, tentatively named DLA-88-DLA-12 and DLA-88-DLA-88L, were identified, and the novel haplotype DLA-88-DLA-88L contributed to 32.7% of the unrelated dogs. Quantitative real-time PCR analysis showed that the gene expression levels of DLA-88L and DLA-88 were similar, and that the gene expression level of DLA-12 was significantly lower. In addition, haplotype frequency estimations using frequently occurring alleles revealed 45 different DLA-class I haplotypes (88-88L/12-64) overall, and 22 different DLA-class I haplotypes in homozygous dogs for 18 breeds and mongrels.

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.

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.

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.

Discovery of novel MHC-class I alleles and haplotypes in Filipino cynomolgus macaques (Macaca fascicularis) by pyrosequencing and Sanger sequencing: Mafa-class I polymorphism.

Although the low polymorphism of the major histocompatibility complex (MHC) transplantation genes in the Filipino cynomolgus macaque (Macaca fascicularis) is expected to have important implications in the selection and breeding of animals for medical research, detailed polymorphism information is still lacking for many of the duplicated class I genes. To better elucidate the degree and types of MHC polymorphisms and haplotypes in the Filipino macaque population, we genotyped 127 unrelated animals by the Sanger sequencing method and high-resolution pyrosequencing and identified 112 different alleles, 28 at cynomolgus macaque MHC (Mafa)-A, 54 at Mafa-B, 12 at Mafa-I, 11 at Mafa-E, and seven at Mafa-F alleles, of which 56 were newly described. Of them, the newly discovered Mafa-A8*01:01 lineage allele had low nucleotide similarities (<86%) with primate MHC class I genes, and it was also conserved in the Vietnamese and Indonesian populations. In addition, haplotype estimations revealed 17 Mafa-A, 23 Mafa-B, and 12 Mafa-E haplotypes integrated with 84 Mafa-class I haplotypes and Mafa-F alleles. Of these, the two Mafa-class I haplotypes, F/A/E/B-Hp1 and F/A/E/B-Hp2, had the highest haplotype frequencies at 10.6 and 10.2%, respectively. This suggests that large scale genetic screening of the Filipino macaque population would identify these and other high-frequency Mafa-class I haplotypes that could be used as MHC control animals for the benefit of biomedical research.

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.

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.

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.

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.

Association analysis of the HLA-C gene in Japanese alopecia areata.

Alopecia areata (AA) is an organ-specific and cell-mediated autoimmune disease involving hair loss, but its pathogenesis remains poorly understood. Many autoimmune diseases are genetically associated with alleles of the human leukocyte antigen (HLA) genes within the major histocompatibility complex. Associations between AA and HLA genes were previously observed in some different ethnic groups. However, the results were inconsistent, and a primary susceptibility HLA gene and/or region has not yet been assigned for AA. The aim of this study was to evaluate whether an allele of the HLA-C locus, HLA-C*07:04, which was strongly associated with AA in Chinese Hans, could be replicated in the Japanese population. The HLA-C locus was genotyped by the SSO method using 156 AA patients and 560 healthy controls. As a consequence, among the 17 alleles detected, only two alleles, C*04:01 (OR = 2.25, CI 95 % = 1.35-3.75, P = 1.84E-03) and C*15:02 (OR = 2.52, CI 95 % = 1.37-4.64, P = 2.90E-03), were significantly associated with AA after Bonferroni correction. Further, the stratification analysis suggested that C*04:01, C*07:02, and C*15:02 represented different AA genetic risk factors in each sub-phenotype.

IL12B and IL23R gene SNPs in Japanese psoriasis.

Psoriasis is a common human skin disease whereby abnormal production of inflammatory mediators is believed to play an important role in its pathogenesis. The IL12B gene, which encodes the shared IL-12p40 subunit in two cytokines, IL-12 and IL-23, and the IL23R gene, which encodes a subunit of the receptor for IL-23, were identified as psoriasis-susceptibility genetic factors in recent candidate gene and genome-wide association studies of Chinese and Europeans. Since there are significant differences in the incidence of psoriasis between Europeans and Japanese suggesting a genetic ethnic effect, we examined the association of IL12B and IL23R gene polymorphisms with psoriasis in a cohort of Japanese. In this study, we genotyped two SNPs (rs3212227 and rs6887695) in the IL12B gene and one SNP (rs11209026) in the IL23R gene using 560 Japanese psoriasis cases and 560 controls and compared our results with those previously published for Europeans and Asians. Our study showed significant associations between psoriasis and both IL12B gene SNPs, rs3212227 (odds ratio (OR) = 1.35, P = 4.94E-04) and rs6887695 (OR = 1.32, P = 2.00E-03), but no significant association between psoriasis and the IL23R SNP, rs11209026. Furthermore, a significant haplotype association was found for the IL12B gene protective haplotype C-C (OR = 0.71, P = 1.84E-04) in Japanese, as previously elucidated in the studies of European ancestry.

Exome sequencing identifies novel rheumatoid arthritis-susceptible variants in the BTNL2.

The butyrophilin-like protein 2 gene (BTNL2) within the class III region of the major histocompatibility complex genomic region was identified as a rheumatoid arthritis (RA) susceptibility gene by exome sequencing (19 RA cases) with stepwise filtering analysis, and then validated by Sanger sequencing and association analysis using 432 cases and 432 controls. Logistic regression of the Sanger-sequenced single-nucleotide variants in an association study of 432 cases and 432 controls showed that 12 non-synonymous single-nucleotide polymorphisms (SNPs) in BTNL2 were significantly associated with RA. The lowest P-values were obtained from three SNPs, rs41521946, rs28362677 and rs28362678, which were in absolute linkage disequilibrium: P=4.55E-09, odds ratio=1.88, 95% confidence interval=1.52-2.33. The BTNL2 locates on chromosome 6 between HLA-DRB1 and NOTCH4, and is 170 kb apart from these two genes. Although DRB1 and NOTCH4 were reported to be RA-susceptible, the three BTNL2 SNPs retained significant association with RA when evaluated by the logistic regression with the adjustment for RA-susceptible HLA-DRB1 alleles in Japanese or rs2071282-T in NOTCH4: P=0.0156 and P=0.00368, respectively. These results suggest that the three non-synonymous SNPs in BTNL2 confer RA risk independently from HLA-DRB1 and NOTCH4.

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.

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.

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.

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.