Search

HLA-DRB5*02:02:05 differs from HLA-DRB5*02:02:01 by one nucleotide substitution in codon 134 in exon 3., (© 2024 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Two novel alleles that appear to be generated from recombination between exons 2 and 3., (© 2024 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

HLA-DRB5*01:01:01:07 differs from HLA-DRB5*01:01:01:01 by two nucleotide changes in intron 1 and intron 2., (© 2024 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

The aim of this study was to compare nonrandom associations between physically adjacent single methylation polymorphism loci among rheumatoid arthritis (RA) and normal subjects for investigating RA-risk methylation haplotypes (meplotype). With 354 ACPA-positive RA patients and 335 normal controls selected from a case-control study based on Swedish population, we conducted the first RA epigenome-wide meplotype association study using our software EWAS2.0, mainly including (i) converted the β value to methylation genotype (menotype) data, (ii) identified methylation disequilibrium (MD) block, (iii) calculated frequent of each meplotypes in MD block and performed case-control association test and (iv) screened for RA-risk meplotypes by odd ratio (OR) and p-values. Ultimately, 545 meplotypes on 334 MD blocks were identified significantly associated with RA (p-value < .05). These meplotypes were mapped to 329 candidate genes related to RA. Subsequently, combined with gene optimization, eight RA-risk meplotypes were identified on three risk genes: HLA-DRB1, HLA-DRB5 and HLA-DQB1. Our results reported the relationship between DNA methylation pattern on HLA-DQB1 and the risk of RA for the first time, demonstrating the co-demethylation of 'cg22984282' and 'cg13423887' on HLA-DQB1 gene (meplotype UU, p-value = 2.90E - 6, OR = 1.68, 95% CI = [1.35, 2.10]) may increase the risk of RA. Our results demonstrates the potential of methylation haplotype analysis to identify RA-related genes from a new perspective and its applicability to the study of other disease., (© 2023 John Wiley & Sons Ltd.)

Differences in gene-wide DNA methylation of the Alzheimer's disease (AD)-associated genes BIN1, HLA-DRB5, SORL1, SLC24A4 , and ABCA7 are reported to be associated with AD in post-mortem brain samples. We investigated whether the same associations could be found in leukocytes collected pre-mortem. Using cohort data of 544 Swedish twins (204 dementia diagnoses), we replicated the findings in HLA-DRB5 and SLC24A4 at P < 0.05. However, co-twin control analyses indicated that the associations were partly explained by familial confounding. Thus, DNA methylation differences in HLA-DRB5 and SLC24A4 are present in both neuronal cells and leukocytes, and not fully explained familial factors.

The novel allele HLA-DRB5*02:35 differs from HLA-DRB5*02:02:01 by one non-synonymous nucleotide substitution in exon 2., (© 2023 The Authors. HLA: Immune Response Genetics published by John Wiley & Sons Ltd.)

Background: Ulcerative colitis (UC) is a form of chronic inflammatory bowel disease of nonspecific origin. This study used an RNA-Sequencing (RNA-Seq) approach to evaluate the transcriptomic landscape of a well-stratified treatment-naïve pediatric UC patient population by comparing them with healthy control children. The data were analyzed to evaluate the mechanisms driving UC-related intestinal inflammation and fibrosis., Methods: Intestinal mucosal samples from five pediatric UC patients and five healthy controls were analyzed by RNA-Seq, and results were verified by qPCR. A CRISPR/Cas9 approach was used to knock out the expression of HLA-DRB5, and molecular biology techniques were used for additional mechanistic studies., Results: In these analyses, 2290 genes were found to be differentially expressed between the UC and control samples, of which 1258 and 1032 were upregulated and downregulated, respectively. Gene Ontology analysis showed that these genes were enriched in extracellular matrix (ECM)-related processes and that 7 of 8 differentially expressed genes of interest (PIK3CD, IL1β, IL1α, TIMP1, MMP1, MMP12, COL6A3, and HLADRB5) were upregulated and involved in ECM-receptor interaction and inflammatory bowel disease-related pathways. Increased HLA-DRB5 expression driven by intestinal bacteria was found to promote IL-1α secretion, leading to intestinal inflammation and fibrosis, suggesting a possible target for the treatment of UC., Conclusion: These data suggest that intestinal inflammation is present in pediatric UC patients for extended periods before the onset of symptoms, and intestinal fibrosis begins even during the early stages of UC. Intestinal bacteria were also found to trigger intestinal inflammation and fibrosis, with HLA-DRB5 playing a central role in this process., (© 2023. The Author(s).)

Background: Sarcoidosis is a multisystem granulomatous disorder whose etiology is related to genetic and immunological factors. Familial aggregation and ethnic prevalence suggest a genetic predisposition and inherited susceptibility to sarcoidosis., Objective: This study aimed to identify suspected risk loci for familial sarcoidosis patients., Methods: We conducted whole exome sequencing on two sarcoidosis patients and five healthy family members in a Chinese family for a case-control study. The two sarcoidosis patients were siblings who showed chronic disease., Results: The Gene Ontology results showed single nucleotide polymorphisms in three genes, including human leukocyte antigen ( HLA)-DRB1, HLA-DRB5, and KIR2DL4 , associated with both 'antigen processing and presentation' and 'regulation of immune response.' Sanger sequencing verified two nonsynonymous mutations in HLA-DRB5 (rs696318 and rs115817940) located on 6p21.3 in the major histocompatibility complex (MHC) class II beta 1 region. The structural model simulated on Prot- Param protein analysis by the Expert Protein Analysis System predicted that the hydropathy index changed at two mutation sites (rs696318: p.F96L, -1.844 to -1.656 and rs115817940: p.T106N, -0.322 to -0.633), which indicated the probability of changes in peptide-binding selectivity., Conclusion: Our results indicated that two nonsynonymous mutations of HLA-DRB5 have been identified in two sarcoidosis siblings, while their healthy family members do not have the mutations. The two HLA-DRB5 alleles may influence genetic susceptibility and chronic disease progression through peptide mutations on the MHC class II molecule among the two affected family members., (Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.net.)

Multiple sclerosis (MS) is a chronic neurological disease believed to be caused by autoimmune pathogenesis. The aetiology is likely explained by a complex interplay between inherited and environmental factors. Genetic investigations into MS have been conducted for over 50 years, yielding >100 associations to date. Globally, the strongest linkage is with the human leukocyte antigen (HLA) HLA-DRB5*01:01:01-DRB1*15:01:01-DQA1*01:02:01-DQB1*06:02:01 haplotype. Here, high-resolution sequencing of HLA was used to determine the alleles of DRB3, DRB4, DRB5, DRB1, DQA1, DQB1, DPA1 and DPB1 as well as their extended haplotypes and genotypes in 100 Swedish MS patients. Results were compared to 636 population controls. The heterogeneity in HLA associations with MS was demonstrated; among 100 patients, 69 extended HLA-DR-DQ genotypes were found. Three extended HLA-DR-DQ genotypes were found to be correlated to MS; HLA-DRB5*01:01:01-DRB1*15:01:01-DQA1*01:02:01-DQB1*06:02:01 haplotype together with (A) HLA-DRB4*01:01:01//DRB4*01:01:01:01-DRB1*07:01:01-DQA1*02:01//02:01:01-DQB1*02:02:01, (B) HLA-DRBX*null-DRB1*08:01:01-DQA1*04:01:01-DQB1*04:02:01, and (C) HLA-DRB3*01:01:02-DRB1*03:01:01-DQA1*05:01:01-DQB1*02:01:01. At the allelic level, HLA-DRB3*01:01:02 was considered protective against MS. However, when combined with HLA-DRB3*01:01:02-DRB1*03:01:01-DQA1*05:01:01-DQB1*02:01:01, this extended haplotype was considered a predisposing risk factor. This highlights the limitations as included with investigations of single alleles relative to those of extended haplotypes/genotypes. In conclusion, with 69 genotypes presented among 100 patients, high-resolution sequencing was conducted to underscore the wide polymorphisms present among MS patients. Additional studies in larger cohorts will be of importance to define MS among the patient group not associated with HLA-DRB5*01:01:01-DRB1*15:01:01-DQA1*01:02:01-DQB1*06:02:01., (© 2022 The Authors. International Journal of Immunogenetics published by John Wiley & Sons Ltd.)

Objectives: We aimed to identify differentially methylated genes and related signaling pathways in autism spectrum disorder (ASD)., Methods: First, the DNA methylation profile in the brain samples (GSE131706 and GSE80017) and peripheral blood samples (GSE109905) was downloaded from the Gene Expression Omnibus database (GEO) dataset, followed by identification of differentially methylated genes and functional analysis. Second, the GSE109905 data set was used to further validate the methylation state and test the ability to diagnose disease of identified differentially methylated genes. Third, expression measurement of selected differentially methylated genes was performed in whole blood from an independent sample. Finally, protein-protein interaction (PPI) network of core differentially methylated genes was constructed., Results: Totally, 74 differentially methylated genes were identified in ASD, including 38 hypermethylated genes and 36 hypomethylated genes. 15 differentially methylated genes were further identified after validation in the GSE109905 data set. Among these, major histocompatibility complex (HLA)-DQA1 was involved in the molecular function of myosin heavy chain class II receptor activity; HLA-DRB5 was involved in the signaling pathways of cell adhesion molecules, Epstein-Barr virus infection and antigen processing and presentation. In the PPI analysis, the interaction pairs of HLA-DQA1 and HLA-DRB5, FMN2 and ACTR3, and CALCOCO2 and BAZ2B were identified. Interestingly, FMN2, BAZ2B, HLA-DRB5, CALCOCO2 and DUSP22 had a potential diagnostic value for patients with ASD. The expression result of four differentially methylated genes (HLA-DRB5, NTM, IL16 and COL5A3) in the independent sample was consistent with the integrated analysis., Conclusions: Identified differentially methylated genes and enriched signaling pathway could be associated with ASD., (Copyright © 2022 Wolters Kluwer Health, Inc. All rights reserved.)

Introduction: Indeterminate acute liver failure (IND-ALF) is a rare clinical syndrome with a high mortality rate. Lacking a known etiology makes rapid evaluation and treatment difficult, with liver transplantation often considered as the only therapeutic option. Our aim was to identify genetic variants from whole exome sequencing data that might be associated with IND-ALF clinical outcomes., Methods: Bioinformatics analysis was performed on whole exome sequencing data for 22 patients with IND-ALF. A 2-tier approach was used to identify significant single-nucleotide polymorphisms (SNPs) associated with IND-ALF clinical outcomes. Tier 1 identified the SNPs with a higher relative risk in the IND-ALF population compared with those identified in control populations. Tier 2 determined the SNPs connected to transplant-free survival and associated with model for end-stage liver disease serum sodium and Acute Liver Failure Study Group prognostic scores., Results: Thirty-one SNPs were found associated with a higher relative risk in the IND-ALF population compared with those in controls, of which 11 belong to the human leukocyte antigen (HLA) class II genes but none for the class I. Further analysis showed that 5 SNPs: rs796202376, rs139189937, and rs113473719 of HLA-DRB5; rs9272712 of HLA-DQA1; and rs747397929 of IDO1 were associated with a higher probability of IND-ALF transplant-free survival. Using 3 selected SNPs, a model for the polygenic risk score was developed to predict IND-ALF prognoses, which are comparable with those by model for end-stage liver disease serum sodium and Acute Liver Failure Study Group prognostic scores., Discussion: Certain gene variants in HLA-DRB5, HLA-DQA1, and IDO1 were found associated with IND-ALF transplant-free survival. Once validated, these identified SNPs may help elucidate the mechanism of IND-ALF and assist in its diagnosis and management., (Copyright © 2022 Written work prepared by employees of the Federal Government as part of their official duties is, under the U.S. Copyright Act, a “work of the United States Government” for which copyright protection under Title 17 of the United States Code is not available. As such, copyright does not extend to the contributions of employees of the Federal Government.)

Background and Objectives: Associations between HLA alleles and susceptibility to M-type phospholipase A2 receptor (PLA2R)-related membranous nephropathy have been well defined previously in Chinese patients. However, the relationships between HLA alleles and kidney outcome remain unclear., Design, Setting, Participants, & Measurements: Five HLA genes (DRB1, DQA1, DQB1, DRB3, and DRB5) were genotyped in a prospective cohort of 392 patients with PLA2R-related membranous nephropathy. The associations between HLA alleles and kidney outcomes were studied., Results: A total of 79 HLA alleles were identified in this study. Four HLA alleles, DRB1*13:01 ( n =12; hazard ratio, 3.7; 95% confidence interval, 1.8 to 7.8; P <0.001), DQB1*06:03 ( n =12; hazard ratio, 3.7; 95% confidence interval, 1.8 to 7.8; P <0.001), DRB1*04:05 ( n =12; hazard ratio, 3.8; 95% confidence interval, 1.5 to 9.5; P =0.004), and DQB1*03:02 ( n =21; hazard ratio, 3.1; 95% confidence interval, 1.4 to 6.7; P =0.005), were associated with a ≥40% eGFR decline during follow-up. DRB1*13:01 and DQB1*06:03 were tightly linked with each other. Forty-four of the 392 patients (11%) carried at least one of the four identified risk HLA alleles in this study. Compared with patients who were negative for all risk HLA alleles, those carrying at least one risk HLA allele had a significant risk of a ≥40% eGFR decline during follow-up (hazard ratio, 3.9; 95% confidence interval, 2.3 to 6.7; P <0.001). After adjusting for age, sex, proteinuria, albumin, eGFR, and anti-PLA2R antibody levels, multivariable Cox analysis showed that patients carrying any of the four risk HLA alleles remained associated with a higher risk of a ≥40% decline in eGFR (hazard ratio, 4.1; 95% confidence interval, 2.3 to 7.1; P <0.001)., Conclusions: Carrying any of the HLA alleles, DRB1*13:01/DQB1*06:03, DRB1*04:05, and DQB1*03:02, was independently associated with poor prognosis in Chinese patients with PLA2R-related membranous nephropathy., (Copyright © 2021 by the American Society of Nephrology.)

Background: Targeted immunotherapy is mostly associated with cancer treatment wherein designed molecules engage signaling pathways and mutant proteins critical to the survival of the cell. One of several genetic approaches is the use of in silico methods to develop immune epitopes targeting specific antigenic regions on related mutant proteins. In a recent study we showed a functional association between the gamma retrovirus HERV-H Long Terminal Associating (HHLA1, HHLA2 and HHLA3) proteins and melanoma associated antigen of the B class proteins (MAGEB5), with a resultant decrease in expression of HLA class I and II immune variants. HLA-C and HLA-DRB5 were the main HLA class I and II Immune variants, respectively, that showed expression changes across viral samples of interest. Specific immune variants for HLA-C and HLA-DRB5 were filtered for the top ten based on their relative frequency of counts across the samples., Results: Protein variants for HHLA1, HHLA2, HHLA3 and MAGEB5 were used to predict antigenic epitope peptides to immune peptide-MHC class I and II binding using artificial neural networks. For IC50 peptide scores (PS) ≥ 0.5 with a transformed binding ability between 0 and 1, the top 5 epitopes identified for all targeted genes HHLA1,2 & 3 and MAGEB5 were qualified as strong or weak binders according to the threshold. Domain analysis using NCBI Conserved Domain Database (CDD) identified HHLA2 with immunoglobulin-like domains (Ig&#95;C1-set) and MAGEB5 with the MAGE Homology Domain (MHD). Linear regression showed a statistical correlation (P < 0.001) for HHLA2 and MAGEB5 predicted epitope peptides to HLA-C but not HLA-DRB5. The prediction model identified HLA-C variant 9 (HLA-C9, BAA08825.1 HLA-B*1511) at 1.1% as the most valuable immune target for clinical considerations. Identification of the 9-mer epitope peptide within the domain showed for HHLA2: YANRTSLFY (PS = 0.5837) and VLAYYLSSSQNTIIN (PS = 0.77) for HLA-C and HLA-DRB5, respectively and for MAGEB5, peptides: FVRLTYLEY (PS = 0.5293) and YPAHYQFLWGPRAYT (PS = 0.62) for HLA-C and HLA-DRB5, respectively., Conclusion: Specific immune responses to targeted epitope peptides and their prediction models, suggested co-expression and co-evolution for HHLA2 and MAGEB5 in viral related diseases. HHLA2 and MAGEB5 could be considered markers for virus related tumors and targeted therapy for oncogenic diseases., (© 2021. The Author(s).)

Allele frequencies and haplotype frequencies of HLA-A, -B, -C, -DRB1, -DRB3/4/5, -DQA1, -DQB1, -DPA1, and -DPB1 have been rarely reported in South Koreans using unambiguous, phase-resolved next generation DNA sequencing. In this study, HLA typing of 11 loci in 173 healthy South Koreans were performed using next generation DNA sequencing with long-range PCR, TruSight® HLA v2 kit, Illumina MiSeqDx platform system, and Assign™ for TruSight™ HLA software. Haplotype frequencies were calculated using the PyPop software. Direct counting methods were used to investigate the association with DRB1 for samples with only one copy of a particular secondary DRB locus. We compared these allele types with the ambiguous allele combinations of the IPD-IMGT/HLA database. We identified 20, 40, 26, 31, 19, 16, 4, and 16 alleles of HLA-A, HLA-B, HLA-C, HLA-DRB1, HLA-DQA1, HLA-DQB1, HLA-DPA1, and HLA-DPB1, respectively. The number of HLA-DRB3/4/5 alleles was 4, 5, and 3, respectively. The haplotype frequencies of most common haplotypes were as follows: A*33:03:01-B*44:03:01-C*14:03-DRB1*13:02:01-DQB1*06:04:01-DPB1*04:01:01 (2.89%), A*33:03:01-B*44:03:01-C*14:03 (4.91%), DRB1*08:03:02-DQA1*01:03:01-DQB1*06:01:01-DPA1*02:02:02-DPB1*05:01:01 (5.41%), DRB1*04:05:01-DRB4*01:03:01 (12.72%), DQA1*01:03:01-DQB1*06:01:01 (13.01%), and DPA1*02:02:02-DPB1*05:01:01 (30.83%). In samples with only one copy of a specific secondary DRB locus, we examined its association with DRB1. We, thus, resolved 10 allele ambiguities in HLA-B, -C (each exon 2+3), -DRB1, -DQB1, -DQA1, and -DPB1 (each exon 2) of the IPD-IMGT/HLA database. Korean population was geographically close to Japanese and Han Chinese populations in the genetic distances by multidimensional scaling (MDS) plots. The information obtained by HLA typing of the 11 extended loci by next generation sequencing may be useful for more exact diagnostic tests on various transplantations and the genetic population relationship studies in South Koreans., Competing Interests: The authors have declared that no competing interests exist.

The incidence of endometrial cancer is increasing each year, and treatment effects are poor for patients with advanced and specific subtypes. Exploring immune infiltration-related factors in endometrial cancer can aid in the prognosis of patients and provide new immunotherapy targets. We downloaded immune metagene and functional data of patients with different subtypes of endometrial cancer from The Cancer Genome Atlas database and selected the lymphocyte-specific kinase (LCK) metagene as a representative genetic marker of the immune microenvironment in endometrial cancer. The results showed that LCK metagene expression is related to the prognosis of patients with endometrioid endometrial adenocarcinoma subtypes and highly correlated with the PTEN and PIK3CA mutational status. A search for LCK-related modules returned seven independent genetic predictors of survival in patients with endometrial cancer. The TIMER algorithm showed that the expression of these seven genes was positively correlated with the infiltration levels of six types of immune cells. The diagnostic value of these markers was validated using real-time quantitative PCR and immunohistochemical methods. Our results identified CD74, HLA-DRB5, CD52, HLA-DPB1 and HLA-DRB1 as possible valuable genetic markers for the diagnosis and prognosis of endometrial cancer and provided a theoretical basis for immunotherapy targets for its clinical treatment.

Both heritability and environment contribute to risk for schizophrenia. However, the molecular mechanisms of interactions between genetic and non-genetic factors remain unclear. Epigenetic regulation of neuronal genome may be a presumable mechanism in pathogenesis of schizophrenia. Here, we performed analysis of open chromatin landscape of gene promoters in prefrontal cortical (PFC) neurons from schizophrenic patients. We cataloged cell-type-based epigenetic signals of transcriptional start sites (TSS) marked by histone H3-K4 trimethylation (H3K4me3) across the genome in PFC from multiple schizophrenia subjects and age-matched control individuals. One of the top-ranked chromatin alterations was found in the major histocompatibility (MHC) locus on chromosome 6 highlighting the overlap between genetic and epigenetic risk factors in schizophrenia. The chromosome conformation capture (3C) analysis in human brain cells revealed the architecture of multipoint chromatin interactions between the schizophrenia-associated genetic and epigenetic polymorphic sites and distantly located HLA-DRB5 and BTNL2 genes. In addition, schizophrenia-specific chromatin modifications in neurons were particularly prominent for non-coding RNA genes, including an uncharacterized LINC01115 gene and recently identified BNRNA&#95;052780. Notably, protein-coding genes with altered epigenetic state in schizophrenia are enriched for oxidative stress and cell motility pathways. Our results imply the rare individual epigenetic alterations in brain neurons are involved in the pathogenesis of schizophrenia.

Next-generation targeted sequencing of HLA-DRB1 and HLA-DRB3, -DRB4, and -DRB5 (abbreviated as DRB345) provides high resolution of functional variant positions to investigate their associations with type 1 diabetes risk and with autoantibodies against insulin (IAA), GAD65 (GADA), IA-2 (IA-2A), and ZnT8 (ZnT8A). To overcome exceptional DR sequence complexity as a result of high polymorphisms and extended linkage disequilibrium among the DR loci, we applied a novel recursive organizer (ROR) to discover disease-associated amino acid residues. ROR distills disease-associated DR sequences and identifies 11 residues of DRB1, sequences of which retain all significant associations observed by DR genes. Furthermore, all 11 residues locate under/adjoining the peptide-binding groove of DRB1, suggesting a plausible functional mechanism through peptide binding. The 15 residues of DRB345, located respectively in the β49-55 homodimerization patch and on the face of the molecule shown to interact with and bind to the accessory molecule CD4, retain their significant disease associations. Further ROR analysis of DR associations with autoantibodies finds that DRB1 residues significantly associated with ZnT8A and DRB345 residues with GADA. The strongest association is between four residues (β14, β25, β71, and β73) and IA-2A, in which the sequence ERKA confers a risk association (odds ratio 2.15, P = 10 -18 ), and another sequence, ERKG, confers a protective association (odds ratio 0.59, P = 10 -11 ), despite a difference of only one amino acid. Because motifs of identified residues capture potentially causal DR associations with type 1 diabetes, this list of residuals is expected to include corresponding causal residues in this study population., (© 2019 by the American Diabetes Association.)

Alzheimer's disease (AD) and Parkinson's disease (PD) are the most common neurodegenerative disorders related to aging. Though several risk factors are shared between these two diseases, the exact relationship between them is still unknown. In this paper, we analyzed how these two diseases relate to each other from the genomic, epigenomic, and transcriptomic viewpoints. Using an extensive literature mining, we first accumulated the list of genes from major genome-wide association (GWAS) studies. Based on these GWAS studies, we observed that only one gene ( HLA-DRB5 ) was shared between AD and PD. A subsequent literature search identified a few other genes involved in these two diseases, among which SIRT1 seemed to be the most prominent one. While we listed all the miRNAs that have been previously reported for AD and PD separately, we found only 15 different miRNAs that were reported in both diseases. In order to get better insights, we predicted the gene co-expression network for both AD and PD using network analysis algorithms applied to two GEO datasets. The network analysis revealed six clusters of genes related to AD and four clusters of genes related to PD; however, there was very low functional similarity between these clusters, pointing to insignificant similarity between AD and PD even at the level of affected biological processes. Finally, we postulated the putative epigenetic regulator modules that are common to AD and PD.

Next Generation Sequencing allows for testing and typing of entire genes of the HLA region. A better and comprehensive sequence assessment can be achieved by the inclusion of full gene sequences of all the common alleles at a given locus. The common alleles of DRB5 are under-characterized with the full exon-intron sequence of two alleles available. In the present study the DRB5 genes from 18 subjects alleles were cloned and sequenced; haplotype analysis showed that 17 of them had a single copy of DRB5 and one consanguineous subject was homozygous at all HLA loci. Methodological approaches including robust and efficient long-range PCR amplification, molecular cloning, nucleotide sequencing and de novo sequence assembly were combined to characterize DRB5 alleles. DRB5 sequences covering from 5'UTR to the end of intron 5 were obtained for DRB5*01:01, 01:02 and 02:02; partial coverage including a segment spanning exon 2 to exon 6 was obtained for DRB5*01:03, 01:08N and 02:03. Phylogenetic analysis of the generated sequences showed that the DRB5 alleles group together and have distinctive differences with other DRB loci. Novel intron variants of DRB5*01:01:01, 01:02 and 02:02 were identified. The newly characterized DRB5 intron variants of each DRB5 allele were found in subjects harboring distinct associations with alleles of DRB1, B and/or ethnicity. The new information provided by this study provides reference sequences for HLA typing methodologies. Extending sequence coverage may lead to identify the disease susceptibility factors of DRB5 containing haplotypes while the unexpected intron variations may shed light on understanding of the evolution of the DRB region., (Copyright © 2019 American Society for Histocompatibility and Immunogenetics. Published by Elsevier Inc. All rights reserved.)

Next Generation Sequencing (NGS) has become a major technology in HLA typing. The expectations are that highly accurate and unambiguous typing results will be obtained. However, HLA typing by NGS has some limitations caused by imperfections in the PCR amplification. The accuracy of NGS data is investigated by analyzing the Short Tandem Repeats (STR) regions. For this analysis HLA-DRB5 is used as the model. The repeat length in a sample highly influences the repeat length distribution present in the reads of NGS data. With a repeat length of 20 only 50% of all reads were of the estimated repeat length, seriously hampering distinguishing allelic differences in this region correctly. Our findings are confirmed by doing the same analysis in HLA-DRB1. Despite the uncertainty of determining the repeat lengths, several new HLA-DRB5 alleles have been identified in this paper., (Copyright © 2019. Published by Elsevier Inc.)

Genotype imputation of the human leukocyte antigen (HLA) region is a cost-effective means to infer classical HLA alleles from inexpensive and dense SNP array data. In the research setting, imputation helps avoid costs for wet lab-based HLA typing and thus renders association analyses of the HLA in large cohorts feasible. Yet, most HLA imputation reference panels target Caucasian ethnicities and multi-ethnic panels are scarce. We compiled a high-quality multi-ethnic reference panel based on genotypes measured with Illumina's Immunochip genotyping array and HLA types established using a high-resolution next generation sequencing approach. Our reference panel includes more than 1,300 samples from Germany, Malta, China, India, Iran, Japan and Korea and samples of African American ancestry for all classical HLA class I and II alleles including HLA-DRB3/4/5. Applying extensive cross-validation, we benchmarked the imputation using the HLA imputation tool HIBAG, our multi-ethnic reference and an independent, previously published data set compiled of subpopulations of the 1000 Genomes project. We achieved average imputation accuracies higher than 0.924 for the commonly studied HLA-A, -B, -C, -DQB1 and -DRB1 genes across all ethnicities. We investigated allele-specific imputation challenges in regard to geographic origin of the samples using sensitivity and specificity measurements as well as allele frequencies and identified HLA alleles that are challenging to impute for each of the populations separately. In conclusion, our new multi-ethnic reference data set allows for high resolution HLA imputation of genotypes at all classical HLA class I and II genes including the HLA-DRB3/4/5 loci based on diverse ancestry populations., (© The Author(s) 2018. Published by Oxford University Press.)

HLA-DRB5*02:21 differs from HLA-DRB5*02:02 by a single-nucleotide substitution (A → G) at position 3785., (© 2019 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

The association of systemic sclerosis with anti-Topoisomerase 1 antibody (ATASSc) with specific alleles of human leukocyte antigen (HLA)-DR has been observed among various ethnics. The anti-Topoisomerase 1 antibody is a common autoantibody in SSc with diffuse cutaneous scleroderma, which is one of the clinical subtypes of SSc. On the other hand, an immunodominant peptide of topoisomerase 1 (Top1) self-protein (residues 349-368) was reported to have strong association with ATASSc. In this study, molecular dynamics simulation was performed on the complexes of Top1 peptide with various HLA-DR subtypes divided into ATASSc-associated alleles (HLA-DRB1*08:02, HLA-DRB1*11:01 and HLA-DRB1*11:04), suspected allele (HLA-DRB5*01:02), and non-associated allele (HLA-DRB1*01:01). The unique interaction for each system was compared to the others in terms of dynamical behaviors, binding free energies and solvation effects. Our results showed that three HLA-DR/Top1 complexes of ATASSc association mostly exhibited high protein stability and increased binding efficiency without solvent interruption, in contrast to non-association. The suspected case (HLA-DRB5*01:02) binds Top1 as strongly as the ATASSc association case, which implied a highly possible risk for ATASSc development. This finding might support ATASSc development mechanism leading to a guideline for the treatment and avoidance of pathogens like Top1 self-peptide risk for ATASSc.

HLA-DRB5*02:21 differs from HLA-DRB5*02:02:01 by one nucleotide substitution at codon 203 in exon 4., (© 2018 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

DRB5*01:01:01:02 differs from DRB5*01:01:01 by four nucleotide changes in intron 2., (© 2018 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Background Steroid-sensitive nephrotic syndrome (SSNS) is a childhood disease with unclear pathophysiology and genetic architecture. We investigated the genomic basis of SSNS in children recruited in Europe and the biopsy-based North American NEPTUNE cohort. Methods We performed three ancestry-matched, genome-wide association studies (GWAS) in 273 children with NS (Children Cohort Nephrosis and Virus [NEPHROVIR] cohort: 132 European, 56 African, and 85 Maghrebian) followed by independent replication in 112 European children, transethnic meta-analysis, and conditional analysis. GWAS alleles were used to perform glomerular cis -expression quantitative trait loci studies in 39 children in the NEPTUNE cohort and epidemiologic studies in GWAS and NEPTUNE (97 children) cohorts. Results Transethnic meta-analysis identified one SSNS-associated single-nucleotide polymorphism (SNP) rs1063348 in the 3' untranslated region of HLA-DQB1 ( P =9.3×10 -23 ). Conditional analysis identified two additional independent risk alleles upstream of HLA-DRB1 (rs28366266, P =3.7×10 -11 ) and in the 3' untranslated region of BTNL2 (rs9348883, P =9.4×10 -7 ) within introns of HCG23 and LOC101929163 These three risk alleles were independent of the risk haplotype DRB1*07:01-DQA1*02:01-DQB1*02:02 identified in European patients. Increased burden of risk alleles across independent loci was associated with higher odds of SSNS. Increased burden of risk alleles across independent loci was associated with higher odds of SSNS, with younger age of onset across all cohorts, and with increased odds of complete remission across histologies in NEPTUNE children. rs1063348 associated with decreased glomerular expression of HLA-DRB1, HLA-DRB5, and HLA-DQB1. Conclusions Transethnic GWAS empowered discovery of three independent risk SNPs for pediatric SSNS. Characterization of these SNPs provide an entry for understanding immune dysregulation in NS and introducing a genomically defined classification., (Copyright © 2018 by the American Society of Nephrology.)

Purpose: Uveitis occurs in a subset of patients with sarcoidosis. The purpose of this study was to determine whether genetic variants that have been associated previously with overall sarcoidosis are associated with increased risk of developing uveitis., Methods: Seventy-seven subjects were enrolled, including 45 patients diagnosed with sarcoidosis-related uveitis as cases and 32 patients with systemic sarcoidosis without ocular involvement as controls. Thirty-eight single nucleotide polymorphisms (SNPs) previously associated with sarcoidosis, sarcoidosis severity, or other organ-specific sarcoidosis involvement were identified. Allele frequencies in ocular sarcoidosis cases versus controls were compared using the chi-square test, and p values were corrected for multiple hypotheses testing using permutation. All analyses were conducted with PLINK., Results: SNPs rs1040461 and rs61860052, in ras-related protein RAS23 ( RAB23 ) and annexin A11 ( ANXA11 ) genes, respectively, were associated with sarcoidosis-associated uveitis. The T allele of rs1040461 and the A allele of rs61860052 were found to be more prevalent in ocular sarcoidosis cases. These associations remained after correction for the multiple hypotheses tested (p=0.01 and p=0.02). In a subanalysis of Caucasian Americans only, two additional variants within the major histocompatibility complex (MHC) genes on chromosome 6, in HLA-DRB5 and HLA-DRB1 , were associated with uveitis as well (p=0.009 and p=0.04)., Conclusions: Genetic variants in RAB23 and ANXA11 genes were associated with an increased risk of sarcoidosis-associated uveitis. These loci have previously been associated with overall sarcoidosis risk.

To evaluate the shared genetic etiology of type 2 diabetes (T2D) and coronary heart disease (CHD), we conducted a genome-wide, multi-ancestry study of genetic variation for both diseases in up to 265,678 subjects for T2D and 260,365 subjects for CHD. We identify 16 previously unreported loci for T2D and 1 locus for CHD, including a new T2D association at a missense variant in HLA-DRB5 (odds ratio (OR) = 1.29). We show that genetically mediated increase in T2D risk also confers higher CHD risk. Joint T2D-CHD analysis identified eight variants-two of which are coding-where T2D and CHD associations appear to colocalize, including a new joint T2D-CHD association at the CCDC92 locus that also replicated for T2D. The variants associated with both outcomes implicate new pathways as well as targets of existing drugs, including icosapent ethyl and adipocyte fatty-acid-binding protein.

This paper presents an overview of genetic variation associated with the onset and progression of 14 neurological disorders, focusing primarily on association studies. The 14 disorders are heterogeneous in terms of their frequency, age of onset, etiology and progression. There is substantially less evidence on progression than onset. With regard to onset, the conditions are diverse in terms of their epidemiology and patterns of familial aggregation. While the muscular dystrophies and Huntington's disease are monogenic diseases, for the other 12 conditions only a small proportion of cases is associated with specific genetic syndromes or mutations. Excluding these, some familial aggregation remains for the majority of cases. There is considerable variation in the volume of evidence by condition, and by gene within condition. The volume of evidence is greatest for Alzheimer's disease, Parkinson's disease, multiple sclerosis and amyotrophic lateral sclerosis. As for common complex chronic diseases, genome wide association studies have found that validated genomic regions account for a low proportion of heritability. Apart from multiple sclerosis, which shares several susceptibility loci with other immune-related disorders, variation at HLA-DRB5 being associated both with Parkinson's disease and Alzheimer's disease, and the association of the C9orf72 repeat expansion with ALS and frontotemporal degeneration, there was little evidence of gene loci being consistently associated with more than one neurological condition or with other conditions. With the exception of spina bifida, for which maternal MTHFR genotype is associated with risk in the offspring, and corroborates other evidence of the importance of folate in etiology, there was little evidence that the pathways influenced by genetic variation are related to known lifestyle or environmental exposures., (Copyright © 2016. Published by Elsevier B.V.)

Background: HLA-DRB3, DRB4 and DRB5 (DRB3/4/5) are paralogues of HLA-DRB1. They have important roles in transplantation and have been reported to be related to many diseases. HLA typing methods for DRB3/4/5 based on NGS data have many limitations now, such as need of polymerase chain reaction (PCR) or low accuracy., Materials and Methods: We present a HLA typing method for DRB3/4/5 based on read mapping and haplotype assembly from NGS data. Also, copy number of DRB3/4/5 is determined by a k-means clustering method according to ratio of sequencing depth between DRB3/4/5 and DRB1., Results: We achieved 100%, 100%, 100% accuracy on simulated data and 95.88%, 98.89%, 99.34% accuracy on MHC capture Illumina sequencing data at 4-digit resolution with 30-fold coverage for DRB3/4/5 separately. We also explored the DRB3/4/5 profiles in five continental populations through low coverage WGS data generated by the 1000 Genome Project. We found that frequency of DRB4 in African were significantly lower than that in all other populations., Conclusion: Our method for DRB3/4/5 typing has high accuracy. It is a good supplement to regular HLA typing and could help in disease studies, medical applications and human population diversity studies., (© 2017 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Introduction: Frontotemporal dementia (FTD), amyotrophic lateral sclerosis (ALS) and Parkinson's disease (PD) are neurodegenerative diseases that share common genetic risk factors. A recent genome-wide association study has linked risk of FTD with polymorphisms in the HLA-DRA/HLA-DRB5 gene (rs9268877, rs9268856), BTNL2 gene (rs1980493), and RAB38/CTSC gene (rs302668)., Methods: We used the SNPscan™ Kit to genotype these variants in 400 Chinese patients with sporadic ALS, 554 with sporadic PD and 634 healthy controls., Results: The AA genotype at rs9268856 increased risk of ALS (P=0.005). Mean survival time was significantly shorter in patients with the AA genotype (24.8±16.2months) than in patients with other genotypes (36.9±19.9months; P<0.001). Kaplan-Meier curves and Cox analysis indicated significantly lower survival probability for patients carrying the AA genotype (P<0.001)., Conclusion: Our results suggest that the AA genotype at rs9268856 is an independent risk factor and prognostic factor for ALS in Han Chinese from southwest China., (Copyright © 2016 Elsevier B.V. All rights reserved.)

Background: Genome-wide association studies (GWAS) have effectively identified genetic factors for many diseases. Many diseases, including Alzheimer's disease (AD), have epistatic causes, requiring more sophisticated analyses to identify groups of variants which together affect phenotype., Results: Based on the GWAS statistical model, we developed a multi-SNP GWAS analysis to identify pairs of variants whose common occurrence signaled the Alzheimer's disease phenotype., Conclusions: Despite not having sufficient data to demonstrate significance, our preliminary experimentation identified a high correlation between GRIA3 and HLA-DRB5 (an AD gene). GRIA3 has not been previously reported in association with AD, but is known to play a role in learning and memory.

Importance: Late-onset Alzheimer disease (AD), the most common form of dementia, places a large burden on families and society. Although epidemiological and clinical evidence suggests a relationship between inflammation and AD, their relationship is not well understood and could have implications for treatment and prevention strategies., Objective: To determine whether a subset of genes involved with increased risk of inflammation are also associated with increased risk for AD., Design, Setting, and Participants: In a genetic epidemiology study conducted in July 2015, we systematically investigated genetic overlap between AD (International Genomics of Alzheimer's Project stage 1) and Crohn disease, ulcerative colitis, rheumatoid arthritis, type 1 diabetes, celiac disease, and psoriasis using summary data from genome-wide association studies at multiple academic clinical research centers. P values and odds ratios from genome-wide association studies of more than 100 000 individuals were from previous comparisons of patients vs respective control cohorts. Diagnosis for each disorder was previously established for the parent study using consensus criteria., Main Outcomes and Measures: The primary outcome was the pleiotropic (conjunction) false discovery rate P value. Follow-up for candidate variants included neuritic plaque and neurofibrillary tangle pathology; longitudinal Alzheimer's Disease Assessment Scale cognitive subscale scores as a measure of cognitive dysfunction (Alzheimer's Disease Neuroimaging Initiative); and gene expression in AD vs control brains (Gene Expression Omnibus data)., Results: Eight single-nucleotide polymorphisms (false discovery rate P < .05) were associated with both AD and immune-mediated diseases. Of these, rs2516049 (closest gene HLA-DRB5; conjunction false discovery rate P = .04 for AD and psoriasis, 5.37 × 10-5 for AD, and 6.03 × 10-15 for psoriasis) and rs12570088 (closest gene IPMK; conjunction false discovery rate P = .009 for AD and Crohn disease, P = 5.73 × 10-6 for AD, and 6.57 × 10-5 for Crohn disease) demonstrated the same direction of allelic effect between AD and the immune-mediated diseases. Both rs2516049 and rs12570088 were significantly associated with neurofibrillary tangle pathology (P = .01352 and .03151, respectively); rs2516049 additionally correlated with longitudinal decline on Alzheimer's Disease Assessment Scale cognitive subscale scores (β [SE], 0.405 [0.190]; P = .03). Regarding gene expression, HLA-DRA and IPMK transcript expression was significantly altered in AD brains compared with control brains (HLA-DRA: β [SE], 0.155 [0.024]; P = 1.97 × 10-10; IPMK: β [SE], -0.096 [0.013]; P = 7.57 × 10-13)., Conclusions and Relevance: Our findings demonstrate genetic overlap between AD and immune-mediated diseases and suggest that immune system processes influence AD pathogenesis and progression.

The possible contribution of HLA-DRB3, -DRB4, and -DRB5 alleles to type 1 diabetes risk and to insulin autoantibody (IAA), GAD65 (GAD autoantibody [GADA]), IA-2 antigen (IA-2A), or ZnT8 against either of the three amino acid variants R, W, or Q at position 325 (ZnT8RA, ZnT8WA, and ZnT8QA, respectively) at clinical diagnosis is unclear. Next-generation sequencing (NGS) was used to determine all DRB alleles in consecutively diagnosed patients ages 1-18 years with islet autoantibody-positive type 1 diabetes (n = 970) and control subjects (n = 448). DRB3, DRB4, or DRB5 alleles were tested for an association with the risk of DRB1 for autoantibodies, type 1 diabetes, or both. The association between type 1 diabetes and DRB1*03:01:01 was affected by DRB3*01:01:02 and DRB3*02:02:01. These DRB3 alleles were associated positively with GADA but negatively with ZnT8WA, IA-2A, and IAA. The negative association between type 1 diabetes and DRB1*13:01:01 was affected by DRB3*01:01:02 to increase the risk and by DRB3*02:02:01 to maintain a negative association. DRB4*01:03:01 was strongly associated with type 1 diabetes (P = 10(-36)), yet its association was extensively affected by DRB1 alleles from protective (DRB1*04:03:01) to high (DRB1*04:01:01) risk, but its association with DRB1*04:05:01 decreased the risk. HLA-DRB3, -DRB4, and -DRB5 affect type 1 diabetes risk and islet autoantibodies. HLA typing with NGS should prove useful to select participants for prevention or intervention trials., (© 2016 by the American Diabetes Association. Readers may use this article as long as the work is properly cited, the use is educational and not for profit, and the work is not altered.)

The genetic association of HLA-DRB1 with rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE) is well documented, but association with other HLA-DR beta genes (HLA-DRB3, HLA-DRB4 and HLA-DRB5) has not been thoroughly studied, despite their similar functions and chromosomal positions. We examined variants in all functional HLA-DR beta genes in RA and SLE patients and controls, down to the amino-acid level, to better understand disease association with the HLA-DR locus. To this end, we improved an existing HLA reference panel to impute variants in all protein-coding HLA-DR beta genes. Using the reference panel, HLA variants were inferred from high-density SNP data of 9,271 RA-control subjects and 5,342 SLE-control subjects. Disease association tests were performed by logistic regression and log-likelihood ratio tests. After imputation using the newly constructed HLA reference panel and statistical analysis, we observed that HLA-DRB1 variants better accounted for the association between MHC and susceptibility to RA and SLE than did the other three HLA-DRB variants. Moreover, there were no secondary effects in HLA-DRB3, HLA-DRB4, or HLA-DRB5 in RA or SLE. Of all the HLA-DR beta chain paralogs, those encoded by HLA-DRB1 solely or dominantly influence susceptibility to RA and SLE.

Background: Epistasis (synergistic interaction) among SNPs governing gene expression is likely to arise within transcriptional networks. However, the power to detect it is limited by the large number of combinations to be tested and the modest sample sizes of most datasets. By limiting the interaction search space firstly to cis-trans and then cis-cis SNP pairs where both SNPs had an independent effect on the expression of the most variable transcripts in the liver and brain, we greatly reduced the size of the search space., Results: Within the cis-trans search space we discovered three transcripts with significant epistasis. Surprisingly, all interacting SNP pairs were located nearby each other on the chromosome (within 290 kb-2.16 Mb). Despite their proximity, the interacting SNPs were outside the range of linkage disequilibrium (LD), which was absent between the pairs (r(2) < 0.01). Accordingly, we redefined the search space to detect cis-cis interactions, where a cis-SNP was located within 10 Mb of the target transcript. The results of this show evidence for the epistatic regulation of 50 transcripts across the tissues studied. Three transcripts, namely, HLA-G, PSORS1C1 and HLA-DRB5 share common regulatory SNPs in the pre-frontal cortex and their expression is significantly correlated. This pattern of epistasis is consistent with mediation via long-range chromatin structures rather than the binding of transcription factors in trans. Accordingly, some of the interactions map to regions of the genome known to physically interact in lymphoblastoid cell lines while others map to known promoter and enhancer elements. SNPs involved in interactions appear to be enriched for promoter markers., Conclusions: In the context of gene expression and its regulation, our analysis indicates that the study of cis-cis or local epistatic interactions may have a more important role than interchromosomal interactions.

Importance: Recent large-scale genome-wide association studies have discovered several genetic variants associated with Alzheimer disease (AD); however, the extent to which DNA methylation in these AD loci contributes to the disease susceptibility remains unknown., Objective: To examine the association of brain DNA methylation in 28 reported AD loci with AD pathologies., Design, Setting, and Participants: Ongoing community-based clinical pathological cohort studies of aging and dementia (the Religious Orders Study and the Rush Memory and Aging Project) among 740 autopsied participants 66.0 to 108.3 years old., Exposures: DNA methylation levels at individual CpG sites generated from dorsolateral prefrontal cortex tissue using a bead assay., Main Outcomes and Measures: Pathological diagnosis of AD by National Institute on Aging-Reagan criteria following a standard postmortem examination., Results: Overall, 447 participants (60.4%) met the criteria for pathological diagnosis of AD. Brain DNA methylation in SORL1, ABCA7, HLA-DRB5, SLC24A4, and BIN1 was associated with pathological AD. The association was robustly retained after replacing the binary trait of pathological AD with 2 quantitative and molecular specific hallmarks of AD, namely, Aβ load and paired helical filament tau tangle density. Furthermore, RNA expression of transcripts of SORL1 and ABCA7 was associated with paired helical filament tau tangle density, and the expression of BIN1 was associated with Aβ load., Conclusions and Relevance: Brain DNA methylation in multiple AD loci is associated with AD pathologies. The results provide further evidence that disruption of DNA methylation is involved in the pathological process of AD.

Systemic lupus erythematosus (SLE) is a polygenic, systemic, autoimmune disease. Copy number variants (CNVS) have been discovered to be associated with a number of complex disorders. We undertook the current study to explore the potential associations between genomic CNVS and SLE in Chinese Han population. In the discovery stage, seven SLE patients were examined with the high-density comparative genomic hybridization microarrays in the screening test for SLE associated CNVS. Then, in the validation stage, 135 SLE patients and 219 matched healthy subjects were investigated for the CNVS of gene HLA-DRB5 by AccuCopy™ technology. Quantitative polymerase chain reaction was carried out to determine the copy number (CN) and mRNA level of HLA-DRB5 in SLE patients. Although the mRNA level of HLA-DRB5 between the CN deletion group and the CN normal group in SLE patients was not statistically positive (P = 0.46), our results still showed more CN of HLA-DRB5 in SLE patients than in healthy controls (P = 3.98 × 10(-6)). Odds ratio for CN deletion was 0.38 (95% confidence interval (CI), 0.23-0.61, P = 7.79 × 10(-5)) and for CN duplication was 1.89 (95% CI, 0.56-7.66, P = 0.37), respectively. These findings indicated that CNVS of HLA-DRB5 was associated with the risk of SLE, and CN deletion appeared to be protective for SLE.

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., (© 2013 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Purpose: This study investigates the extent of the human transcriptome that can be quantified from conjunctival impression cytology extracts. The aim is to determine if sufficient RNA can be isolated from a patient's conjunctival surface to identify differences in gene expression between dry eye and normal patients of (a) an array of 96 inflammatory biomarkers and associated receptors, and (b) if this comparison can be expanded to the entire transcriptome., Materials and Methods: CIC was used to collect conjunctival surface cells from 53 qualifying normal and dry eye patients. Based on prior optimization of all assay steps, RNA was isolated from the samples using a Qiagen RNeasy Plus Mini Kit and qRT-PCR was used to determine gene expression of 96 genes using TaqMan Low Density Array cards. Samples from six normal and six dry eye patients were then assayed on an Illumina Human HT-12 BeadChip., Results: Optimization steps yielded an RNA processing procedure that improved yield from an initial 12 genes through 96, then to the entire human transcriptome. For the HT-12 BeadChip, more than 30 genes differed by a factor of >1.5 between the dry eye and normal groups and seven genes were down-regulated by a factor of >2.0 in the dry eye group: HLA-DRB5, PSCA, FOS, lysozyme, TSC22D1, CAPN13 and CXCL6., Conclusions: Conjunctival impression cytology can be used to collect sufficient RNA from conjunctival surface cells that, when processed optimally, allows successful transcriptome-wide expression analysis. While the current transcriptome analysis used a limited patient group, larger studies of patients with various types and severities of dry eye should reveal significant gene expression trends that can then be targeted to improve dry eye treatment options.

Historically, association of disease with the major histocompatibility complex (HLA) genes has been tested with HLA alleles that encode antigen-binding affinity. The association with Parkinson disease (PD), however, was discovered with noncoding SNPs in a genome-wide association study (GWAS). We show here that several HLA-region SNPs that have since been associated with PD form two blocks tagged by rs3129882 (p = 9 × 10(-11)) and by rs9268515 and/or rs2395163 (p = 3 × 10(-11)). We investigated whether these SNP-associations were driven by HLA-alleles at adjacent loci. We imputed class I and class II HLA-alleles for 2000 PD cases and 1986 controls from the NeuroGenetics Research Consortium GWAS and sequenced a subset of 194 cases and 204 controls. We were therefore able to assess accuracy of two imputation algorithms against next-generation-sequencing while taking advantage of the larger imputed data sets for disease study. Additionally, we imputed HLA alleles for 843 cases and 856 controls from another GWAS for replication. PD risk was positively associated with the B(∗)07:02&#95;C(∗)07:02&#95;DRB5(∗)01&#95;DRB1(∗)15:01&#95;DQA1(∗)01:02&#95;DQB1(∗)06:02 haplotype and negatively associated with the C(∗)03:04, DRB1(∗)04:04 and DQA1(∗)03:01 alleles. The risk haplotype and DQA1(∗)03:01 lost significance when conditioned on the SNPs, but C(∗)03:04 (OR = 0.72, p = 8 × 10(-6)) and DRB1(∗)04:04 (OR = 0.65, p = 4 × 10(-5)) remained significant. Similarly, rs3129882 and the closely linked rs9268515 and rs2395163 remained significant irrespective of HLA alleles. rs3129882 and rs2395163 are expression quantitative trait loci (eQTLs) for HLA-DR and HLA-DQ (9 × 10(-5) ≥ PeQTL ≥ 2 × 10(-79)), suggesting that HLA gene expression might influence PD. Our data suggest that PD is associated with both structural and regulatory elements in HLA. Furthermore, our study demonstrates that noncoding SNPs in the HLA region can be associated with disease irrespective of HLA alleles, and that observed associations with HLA alleles can sometimes be secondary to a noncoding variant., (Copyright © 2013 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

A genetic study, particularly in HLA-DRs, has never been performed in Thai patients with systemic sclerosis (SSc). This study was performed to investigate the association between the HLA-DR series in Thai SSc patients. HLA-DR subtypes were determined in 50 Thai SSc patients and 99 healthy controls (HCs). All SSc patients met the ACR classification criteria for SSc. HLA-DR typing was performed using INNO-LiPA HLA-DRB Decoder kits (INNOGENETICS) and reconfirmed using MICRO SSP HLA DNA Typing kits (ONE LAMBDA). The allele frequency (AF) of HLA-DR*15, compared with HC, was significantly higher in all SSc patients (41.0 vs 21.7%, Pc = 0.0083) and SSc patients with anti-Scl70 antibody positive (anti-Scl70+) (47.1%, Pc = 0.0018). Among the HLA-DR*15 alleles, the AF of the DRB1*15:02 was increased significantly in all SSc patients (29.0 vs 12.6%, Pc = 0.0219) and SSc patients with anti-Scl70+ (32.4 vs 12.6%, Pc = 0.0196). The AF of the HLA-DRB5*01:02 allele was also increased in all SSc patients (27.0 vs 12.6%, Pc = 0.0166) and in SSc patients with anti-Scl70+ (29.4%, Pc = 0.0124). The AF of the DR*04 was significantly lower in the SSc patients (1.0 vs 9.6%, Pc = 0.0399). However, the AF of the DRB1*15:02 and DRB5*01:02 was not different among SSc patients with or without clinical manifestations (pulmonary fibrosis, digital pitting scar, sclerodactyly, myositis, and sicca symptoms). In addition, there was no significant association between clinical manifestations among individuals who carried HLA-DRB1*15:02 or DRB5*01:02. HLA-DRB1*15:02 and DRB5*01:02 alleles were significantly elevated in Thai SSc patients, especially in those with anti-Scl70+. The HLA-DRB1*04 was a protective allele against Thai SSc patients.

Human leucocyte antigen (HLA) study in patients with systemic lupus erythematosus (SLE) has been investigated in various countries, but the results are still inconclusive. This study was performed to investigate the association between HLA-DR and SLE in patients in northern Thailand. HLA-DR subtyping was performed in 70 patients with SLE and 99 normal healthy controls living in northern Thailand using the INNO-LiPA HLA-DR Decoder kit (Innogenetics) and MICRO SSP HLA DNA Typing kit (One Lambda) for reconfirmation. The allele frequency (AF) of DRB5*01:01 in SLE was significantly higher than in the controls [25.7% vs. 14.6%, P = 0.012, Pc = 0.048, OR = 2.02 (95%CI = 1.17-3.48)]. The AF of DRB1*15:01 and DRB1*16:02 showed a nonsignificant tendency to be higher in SLE (10.7% vs. 8.1%, and 17.9% vs. 11.1%). Interestingly, the DRB5*01:01 allele linked to DRB1*16:02 in 47.2% of SLE and 37.9% of controls, and the prevalence of the DRB1*16:02-DRB5*01:01 haplotype was higher in the patients with SLE [12.1% vs. 5.6%, P = 0.044, OR = 2.35 (95%CI = 1.06-5.19)]. The DRB1*16:02 linked to DRB5*02:02 and *02:03 in 18.2% and 31.8% of controls, respectively, and linked to DRB5*02:03 in 32.0% of SLE patients. The frequency of DRB1*03:01 and *15:02 alleles was not increased in Thai SLE. There was no significant association between DRB5*01:01 and any auto-antibodies or clinical manifestations of SLE. DRB5*01:01 is associated with Thai SLE, and the association is stronger than that of DRB1*15:01. The genetic contribution of DRB5*01:01 is due partially to the linkage disequilibrium between DRB1*16:02 and DRB5*01:01 in the northern Thai population., (© 2012 Blackwell Publishing Ltd.)

We compared the carriage frequencies of HLA-DRB3 and its major alleles and of HLA-DRB4 and HLA-DRB5 in an Australian sIBM cohort and a population control group who had previously been genotyped for the HLA-DRB1*03:01 risk allele. There was a strong disease association with the carriage of the DRB3*01:01 allele which was accounted for by its linkage disequilibrium with DRB1*03:01. The carriage of HLA-DRB4 was found to be strongly protective and abrogated the risk effect of HLA-DRB1*03:01. The findings indicate that haplotypic combinations of alleles at the HLA-DRB1 and secondary HLA-DRB loci have important risk modifying effects in sIBM., (Copyright © 2012 Elsevier B.V. All rights reserved.)

Objective: Interstitial lung disease (ILD) is a serious complication of SSc. We aimed to identify markers associated with SSc-related ILD., Methods: RNA was prepared from the peripheral blood mononuclear cells of 14 SSc patients, divided into four different RNA pools according to the presence or absence of ILD and to the treatment, and subjected to microarray analysis. Real-time quantitative PCR was used to confirm the microarray results in 43 SSc patients, 42 autoimmune controls and 10 healthy controls. Genomic DNA samples were collected from 149 patients with SSc (70 in Hokkaido and 79 in Tokyo) who underwent a high-resolution CT for the evaluation of ILD and from 230 healthy controls. Genotyping was performed using sequence-specific primers., Results: The microarray analysis revealed HLA-DRB5 to be the only gene commonly up-regulated in patients with ILD compared with those without ILD in both comparison groups. High expression levels of HLA-DRB5 in SSc patients with ILD were confirmed by real-time quantitative PCR. The prevalence of HLA-DRB5 gene carriers increased in the SSc patients with ILD relative to those without ILD or to healthy controls in both cohorts. Among the four detected alleles, the HLA-DRB5*01:05 allele was significantly more frequent in SSc patients with ILD than in SSc patients without ILD or in healthy controls. These associations were confirmed in the second cohort., Conclusion: HLA-DRB5 was highly expressed in PBMCs from patients with SSc-related ILD. The HLA-DRB5*01:05 allele is a risk factor for ILD in patients with SSc.

Genetic susceptibility to multiple sclerosis (MS) has been linked to the HLA-DR15 haplotype consisting of DRB1*15:01(DR2b) and DRB5*01:01(DR2a) alleles. Given almost complete linkage disequilibrium of the two alleles, recent studies suggested differential roles in susceptibility (DR2b) or protection from MS (DR2a). Our objective was to assess the potential contribution of DR2a to disease etiology in MS using a humanized model of autoimmunity. To assess the potential contribution of DR2a to disease etiology, we created DR2a humanized transgenic (Tg) mice and subsequently crossed them to Tg mice expressing TL3A6, an MS patient-derived myelin basic protein 83-99-specific TCR. In TL3A6/DR2a Tg mice, CD4 Tg T cells escape thymic and peripheral deletion and initiate spontaneous experimental autoimmune encephalomyelitis (EAE) at low rates, depending on the level of DR2a expression. The ability to induce active EAE was also increased in animals expressing higher levels of DR2a. Inflammatory infiltrates and neuronal damage were present throughout the spinal cord, consistent with a classical ascending EAE phenotype with minor involvement of the cerebellum, brainstem, and peripheral nerve roots in spontaneous, as well as actively induced, disease. These studies emphasize the pathologic contribution of the DR2a allele to the development of autoimmunity when expressed as the sole MHC class II molecule, as well as strongly argue for DR2a as a contributor to the CNS autoimmunity in MS.

Keloid disease (KD) is a common fibroproliferative disorder that can occur following cutaneous injury in genetically susceptible individuals. Familial predisposition, high prevalence in certain populations and occurrence in twins suggest a strong genetic component to KD. However, to date no single causative gene has been identified. Copy number variations (CNVs) in genes have been associated with several human diseases including common skin disorders. The objective of this study was, therefore to determine if CNVs in the human genome may contribute to the development of KD. Agilent SurePrint G3 one Million microarrays were used to detect DNA copy number differences in keloid scars of four Caucasian females and compared to commercial reference DNA samples. Subsequent validation was performed with quantitative polymerase chain reactions (qPCR) using 15 KD cases and 27 Caucasian controls. Further validation using a second cohort was carried out with an additional 11 Caucasian controls and 10 KD cases developed from minor skin puncture wounds (caused by acne, vaccination or chickenpox). HLA-DRB1*15 was typed using allele-specific primers in PCR. Five CNV regions located at chromosome (chr) 6p21.32, chr11q11, chr17q12, chr8p23.1, chr22q13.1, chr19p13.1 and chr2q14.3 were selected for validation with qPCR. When comparing controls to subgroups of KD, according to their cause of scarring, chr6p21.32 (primer design targetting HLA-DRB5) was significantly (P < 0.005) associated with KD developed from minor skin puncture wounds, which is further validated using a larger sample size (P < 0.001). The presence of HLA-DRB5 was associated with HLA-DRB1*15 status; 18 out of 19 individuals were positive for both HLA-DRB5 and HLA-DRB1*15 allele. In conclusion, these preliminary findings further support the possible contribution of the HLA genes in KD pathogenesis.

Although many studies have analyzed HLA allele frequencies in several ethnic groups in patients with scleroderma (SSc), none has been done in French Caucasian patients and none has evaluated which one of the common amino acid sequences, (67)FLEDR(71), shared by HLA-DRB susceptibility alleles, or (71)TRAELDT(77), shared by HLA-DQB1 susceptibility alleles in SSc, was the most important to develop the disease. HLA-DRB and DQB typing was performed for a total of 468 healthy controls and 282 patients with SSc allowing FLEDR and TRAELDT analyses. Results were stratified according to patient's clinical subtypes and autoantibody status. Moreover, standardized HLA-DRß1 and DRß5 reverse transcriptase Taqman PCR assays were developed to quantify ß1 and ß5 mRNA in 20 subjects with HLA-DRB1*15 and/or DRB1*11 haplotypes. FLEDR motif is highly associated with diffuse SSc (χ(2) = 28.4, p<10-6) and with anti-topoisomerase antibody (ATA) production (χ(2) = 43.9, p<10-9) whereas TRAELDT association is weaker in both subgroups (χ(2) = 7.2, p = 0.027 and χ(2) = 14.6, p = 0.0007 respectively). Moreover, FLEDR motif- association among patients with diffuse SSc remains significant only in ATA subgroup. The risk to develop ATA positive SSc is higher with double dose FLEDR than single dose with respectively, adjusted standardised residuals of 5.1 and 2.6. The increase in FLEDR motif is mostly due to the higher frequency of HLA-DRB1*11 and DRB1*15 haplotypes. Furthermore, FLEDR is always carried by the most abundantly expressed ß chain: ß1 in HLA DRB1*11 haplotypes and ß5 in HLA-DRB1*15 haplotypes.In French Caucasian patients with SSc, FLEDR is the main presenting motif influencing ATA production in dcSSc. These results open a new field of potential therapeutic applications to interact with the FLEDR peptide binding groove and prevent ATA production, a hallmark of severity in SSc.

The human leukocyte antigen (HLA) DRB1*1501 has been consistently associated with multiple sclerosis (MS) in nearly all populations tested. This points to a specific antigen presentation as the pathogenic mechanism though this does not fully explain the disease association. The identification of expression quantitative trait loci (eQTL) for genes in the HLA locus poses the question of the role of gene expression in MS susceptibility. We analyzed the eQTLs in the HLA region with respect to MS-associated HLA-variants obtained from genome-wide association studies (GWAS). We found that the Tag of DRB1*1501, rs3135388 A allele, correlated with high expression of DRB1, DRB5 and DQB1 genes in a Caucasian population. In quantitative terms, the MS-risk AA genotype carriers of rs3135388 were associated with 15.7-, 5.2- and 8.3-fold higher expression of DQB1, DRB5 and DRB1, respectively, than the non-risk GG carriers. The haplotype analysis of expression-associated variants in a Spanish MS cohort revealed that high expression of DRB1 and DQB1 alone did not contribute to the disease. However, in Caucasian, Asian and African American populations, the DRB1*1501 allele was always highly expressed. In other immune related diseases such as type 1 diabetes, inflammatory bowel disease, ulcerative colitis, asthma and IgA deficiency, the best GWAS-associated HLA SNPs were also eQTLs for different HLA Class II genes. Our data suggest that the DR/DQ expression levels, together with specific structural properties of alleles, seem to be the causal effect in MS and in other immunopathologies rather than specific antigen presentation alone.