HLA Association among Thai Patients with Diffuse and Limited Cutaneous Systemic Sclerosis.

This study aimed to clarify the association of HLA Class I and II with dcSSc and lcSSc in Thais. HLA typing for 11 gene loci (Class I: HLA-A, B and C, and Class II [HLA-DR, DP and DQ]) was carried out using the Next Generation DNA Sequencing method (three fields) in 92 Thai patients with systemic sclerosis (55 dcSSc, 37 lcSSc) and 135 healthy controls (HCs). The distribution of HLA alleles in patients with dcSSc and lcSSc was compared. When compared with HCs, the AF of A*24:02:01, A*24:07:01, B*27:04:01 and B*27:06 showed an increasing trend in lcSSc patients without statistical significance. DRB1*15:02:01, DRB5*01:02:01, DQA1*01:01:01, DQB1*05:01:24, DPA1*02:01:01 and DPB1*13:01:01 increased significantly in dcSSc patients. DQB1*05:01:24 and DPB1*13:01:01 also increased significantly in lcSSc patients, but less significantly than in dcSSc patients. The association of DPB1*05:01:01 with lcSSc was significantly protective. HLA-A*24:02:01, B*27:06 and C*03:04:01 formed a three-locus haplotype that also constituted an eight-locus haplotype with DRB1*15:02:01, DQA1*01:01:01, DQB1*05:01:24, DPA1*02:01:01 and DPB1*13:01:01. There was a possibility that HLA Class I would play a role in the pathogenesis of lcSSc, while Class II played more of a role in the dcSSc in Thai patients.

Amino acid polymorphisms in human histocompatibility leukocyte antigen class II and proinsulin epitope have impacts on type 1 diabetes mellitus induced by immune-checkpoint inhibitors.

Introduction: Immune-checkpoint inhibitors are effective in various advanced cancers. Type 1 diabetes mellitus induced by them (ICI-T1DM) is a serious complication requiring prompt insulin treatment, but the immunological mechanism behind it is unclear. Methods: We examined amino acid polymorphisms in human histocompatibility leukocyte antigen (HLA) molecules and investigated proinsulin epitope binding affinities to HLA molecules. Results and Discussion: Twelve patients with ICI-T1DM and 35 patients in a control group without ICI-T1DM were enrolled in the study. Allele and haplotype frequencies of HLA-DRB1*04:05, DQB1*04:01, and most importantly DPB1*05:01 were significantly increased in patients with ICI-T1DM. In addition, novel amino acid polymorphisms in HLA-DR (4 polymorphisms), in DQ (12 polymorphisms), and in DP molecules (9 polymorphisms) were identified. These amino acid polymorphisms might be associated with the development of ICI-T1DM. Moreover, novel human proinsulin epitope clusters in insulin A and B chains were discovered in silico and in vitro peptide binding assays to HLA-DP5. In conclusion, significant amino acid polymorphisms in HLA-class II molecules, and conformational alterations in the peptide-binding groove of the HLA-DP molecules were considered likely to influence the immunogenicity of proinsulin epitopes in ICI-T1DM. These amino acid polymorphisms and HLA-DP5 may be predictive genetic factors for ICI-T1DM.

Association of HLA-DRB1*15:02:01, DQB1*05:01:24 and DPB1*13:01:01 in Thai patients with systemic sclerosis.

HLA studies in patients with systemic sclerosis (SSc) have shown variable results. This study aimed to examine the association of HLA class I and II risk alleles in Thai SSc patients, and clarify the contribution of risk HLA alleles to the pathogenesis and clinical manifestations. Blood samples from 92 SSc patients and 135 healthy controls (HCs) were collected. Eleven loci of the HLA class I (HLA-A, B, and C) and class II (HLA-DR, DP, and DQ) genes were determined by a 3-field (6-digit) analysis using the Next Generation DNA Sequencing (NGS) method. Anti-topoisomerase-I antibodies (ATA) and anti-centromere antibodies (ACA) were identified by ELISA methods. Allele frequencies (AFs) of HLA-DRB1*15:02:01, DRB5*01:02:01, DQB1*05:01:24, DPB1*13:01:01, and DQA1*01:01:01 were increased significantly in the whole SSc and SSc patients with positive ATA, but with negative ACA (SSc/ATA+/ACA-). Of these, DPB1*13:01:01 was the most susceptible allele. The DRB1*15:02:01, DQB1:05:01:24, and DPB1*13:01:01 alleles were estimated to locate on the unique haplotype, and haplotype frequency was estimated to be significantly higher than those in the HCs (p = 0.002). The linkage analysis of DRB1*15/16 revealed that most of the DRB1*15:02:01 alleles were linked to DRB5*01:02:01 or DRB5*01:08:01N. The linkage of DRB1*16:02:01 to DRB5*01:01:01 was observed frequently. The associations of risk alleles with several SSc clinical features were observed. HLA-DRB1*15:02:01, DRB5*01:02:01, DQB1*05:01:24, and DPB1*13:01:01 on the unique haplotype were associated with the pathogenesis and clinical features of SSc in Thai patients. The linkage of DRB1*15:02:01 to DRB5*01:08:01N was observed commonly in northern Thai patients.

Identification of HLA-A*02:06:01 as the primary disease susceptibility HLA allele in cold medicine-related Stevens-Johnson syndrome with severe ocular complications by high-resolution NGS-based HLA typing.

Stevens-Johnson syndrome (SJS) and toxic epidermal necrolysis (TEN) are life-threatening acute inflammatory vesiculobullous reactions of the skin and mucous membranes. These severe cutaneous drug reactions are known to be caused by inciting drugs and infectious agents. Previously, we have reported the association of HLA-A*02:06 and HLA-B*44:03 with cold medicine (CM)-related SJS/TEN with severe ocular complications (SOCs) in the Japanese population. However, the conventional HLA typing method (PCR-SSOP) sometimes has ambiguity in the final HLA allele determination. In this study, we performed HLA-disease association studies in CM-SJS/TEN with SOCs at 3- or 4-field level. 120 CM-SJS/TEN patients with SOCs and 817 Japanese healthy controls are HLA genotyped using the high-resolution next-generation sequencing (NGS)-based HLA typing of HLA class I genes, including HLA-A, HLA-B, and HLA-C. Among the alleles of HLA class I genes, HLA-A*02:06:01 was strongly associated with susceptibility to CM-SJS/TEN (p = 1.15 × 10-18, odds ratio = 5.46). Four other alleles (HLA-A*24:02:01, HLA-B*52:01:01, HLA-B*46:01:01, and HLA-C*12:02:02) also demonstrated significant associations. HLA haplotype analyses indicated that HLA-A*02:06:01 is primarily associated with susceptibility to CM-SJS/TEN with SOCs. Notably, there were no specific disease-causing rare variants among the high-risk HLA alleles. This study highlights the importance of higher resolution HLA typing in the study of disease susceptibility, which may help to elucidate the pathogenesis of CM-SJS/TEN with SOCs.

Aggregation of rare/low-frequency variants of the mitochondria respiratory chain-related proteins in rheumatoid arthritis patients.

Exome sequencings were conducted using 59 patients having rheumatoid arthritis (RA) and 93 controls. After stepwise filtering, 107 genes showed less than 0.05 of P-values by gene-burden tests. Among 107 genes, NDUFA7 which is a subunit of the complex I in the mitochondrial respiratory chain was selected for further analysis based on previous reports. A case-control study was performed on the three single-nucleotide variants (SNVs) of NDUFA7 with 432 cases and 432 controls. An association was observed between NDUFA7 and RA with severe erosive arthritis. These results together with previous reports suggested the involvement of reactive oxygen species (ROS) in the pathogenesis of RA. In the next step, four SNVs from three genes related to the mitochondrial respiratory chain were selected, which is a major source of ROS, and conducted a case-control study. An association was observed based on a pathway-burden test comprising NDUFA7, SDHAF2, SCO1 and ATP5O: P=1.56E-04, odds ratio=2.16, 95% confidence interval=1.43-3.28. Previous reports suggested the involvement of ROS in the pathogenesis of RA. The aggregation of SNVs in the mitochondria respiratory chain suggests the pivotal role of those SNVs in the pathogenesis of RA with severe erosive arthritis.

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.

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.

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.

Exact break point of a 50 kb deletion 8 kb centromeric of the HLA-A locus with HLA-A*24:02: the same deletion observed in other A*24 alleles and A*23:01 allele.

In a structural aberration analysis of patients with arthritis mutilans, a 50 kb deletion near the HLA-A locus with HLA-A*24:02 allele was detected. It was previously reported that HLA-A*24:02 haplotype harbored a large-scale deletion telomeric of the HLA-A gene in healthy individuals. In order to confirm that the deletion are the same in patients with arthritis mutilans and in healthy individuals, and to identify the break point of this deletion, the boundary sequences across the deletion in A*24:02 was amplified by polymerase chain reaction (PCR) as a 3.7 kb genomic fragment and subjected to nucleotide sequence determination. A comparison of these genomic sequences with those of the non-A*24:02 haplotype revealed that the deleted genomic region spanning 50 kb was flanked by 3.7 kb repetitive element-rich segments homologous to each other on both sides in non-A*24. The nucleotide sequences of the PCR products were identical in patients with arthritis mutilans and in healthy individuals, revealing that the deletion linked to A*24:02 is irrelevant to the onset of arthritis mutilans. The deletion was detected in all other A*24 alleles so far examined but not in other HLA-A alleles, except A*23:01. This finding, along with the phylogenic tree of HLA-A alleles and the presence of the 3.7 kb highly homologous segments at the boundary of the deleted genomic region in A*03 and A*32, may suggest that this HLA-A*24:02-linked deletion was generated by homologous recombination within two 3.7 kb homologous segments situated 50 kb apart in the ancestral A*24 haplotype after divergence from the A*03 and A*32 haplotypes.

Particular human leukocyte antigen alleles are associated with biochemical traits in the Japanese population.

We analyzed genetic associations among 7 biochemical traits (fasting plasma glucose, HbA1c, total cholesterol, low-density lipoprotein [LDL] cholesterol, high-density lipoprotein cholesterol, triglyceride, and uric acid) and 6 HLA loci using 1,616 individuals who visited the Health Evaluation and Promotion Center at Tokai University Hospital. Significant differences between the individuals carrying particular HLA alleles and those not carrying the alleles in certain biochemical traits were observed by Mann-Whitney U test. In female subjects, DPB1*03:01 was significantly associated with HbA1c (p = 0.0000665), and DRB1*14:03 was associated with total cholesterol concentration (p = 0.0015). In male subjects, C*14:02 demonstrated significant associations with fasting plasma glucose with p values of 0.0041. By contrast, Fisher's exact test indicated that female DRB1*14:03 was associated with a high concentration of total cholesterol (p = 000323, odds ratio [OR] = 4.32, 95% confidence interval [95% CI] = 1.83-10.36), whereas female DPB1*02:01 had a protective effect against a high concentration of LDL cholesterol (p =0.0043, OR = 0.41, 95% CI = 0.19-0.79). These associations have a statistical power of more than 0.8 and still retain significance after Bonferroni correction.