Transcriptomic analyses of treatment-naïve pediatric ulcerative colitis patients and exploration of underlying disease pathogenesis.

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.

Epigenome-wide methylation haplotype association analysis identified HLA-DRB1, HLA-DRB5 and HLA-DQB1 as risk factors for rheumatoid arthritis.

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.

Co-Expression Network Analysis Identifies Molecular Determinants of Loneliness Associated with Neuropsychiatric and Neurodegenerative Diseases.

Loneliness and social isolation are detrimental to mental health and may lead to cognitive impairment and neurodegeneration. Although several molecular signatures of loneliness have been identified, the molecular mechanisms by which loneliness impacts the brain remain elusive. Here, we performed a bioinformatics approach to untangle the molecular underpinnings associated with loneliness. Co-expression network analysis identified molecular 'switches' responsible for dramatic transcriptional changes in the nucleus accumbens of individuals with known loneliness. Loneliness-related switch genes were enriched in cell cycle, cancer, TGF-ß, FOXO, and PI3K-AKT signaling pathways. Analysis stratified by sex identified switch genes in males with chronic loneliness. Male-specific switch genes were enriched in infection, innate immunity, and cancer-related pathways. Correlation analysis revealed that loneliness-related switch genes significantly overlapped with 82% and 68% of human studies on Alzheimer's (AD) and Parkinson's diseases (PD), respectively, in gene expression databases. Loneliness-related switch genes, BCAM, NECTIN2, NPAS3, RBM38, PELI1, DPP10, and ASGR2, have been identified as genetic risk factors for AD. Likewise, switch genes HLA-DRB5, ALDOA, and GPNMB are known genetic loci in PD. Similarly, loneliness-related switch genes overlapped in 70% and 64% of human studies on major depressive disorder and schizophrenia, respectively. Nine switch genes, HLA-DRB5, ARHGAP15, COL4A1, RBM38, DMD, LGALS3BP, WSCD2, CYTH4, and CNTRL, overlapped with known genetic variants in depression. Seven switch genes, NPAS3, ARHGAP15, LGALS3BP, DPP10, SMYD3, CPXCR1, and HLA-DRB5 were associated with known risk factors for schizophrenia. Collectively, we identified molecular determinants of loneliness and dysregulated pathways in the brain of non-demented adults. The association of switch genes with known risk factors for neuropsychiatric and neurodegenerative diseases provides a molecular explanation for the observed prevalence of these diseases among lonely individuals.

High-resolution HLA class II sequencing of Swedish multiple sclerosis patients.

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.

Construction and benchmarking of a multi-ethnic reference panel for the imputation of HLA class I and II alleles.

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.

Targeted immune epitope prediction to HHLA2 and MAGEB5 protein variants as therapeutic approach to related viral diseases.

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_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.

Complementary Effects of Carbamylated and Citrullinated LL37 in Autoimmunity and Inflammation in Systemic Lupus Erythematosus.

LL37 acts as T-cell/B-cell autoantigen in Systemic lupus erythematosus (SLE) and psoriatic disease. Moreover, when bound to "self" nucleic acids, LL37 acts as "danger signal," leading to type I interferon (IFN-I)/pro-inflammatory factors production. T-cell epitopes derived from citrullinated-LL37 act as better antigens than unmodified LL37 epitopes in SLE, at least in selected HLA-backgrounds, included the SLE-associated HLA-DRB1*1501/HLA-DRB5*0101 backgrounds. Remarkably, while "fully-citrullinated" LL37 acts as better T-cell-stimulator, it loses DNA-binding ability and the associated "adjuvant-like" properties. Since LL37 undergoes a further irreversible post-translational modification, carbamylation and antibodies to carbamylated self-proteins other than LL37 are present in SLE, here we addressed the involvement of carbamylated-LL37 in autoimmunity and inflammation in SLE. We detected carbamylated-LL37 in SLE-affected tissues. Most importantly, carbamylated-LL37-specific antibodies and CD4 T-cells circulate in SLE and both correlate with disease activity. In contrast to "fully citrullinated-LL37," "fully carbamylated-LL37" maintains both innate and adaptive immune-cells' stimulatory abilities: in complex with DNA, carbamylated-LL37 stimulates plasmacytoid dendritic cell IFN-α production and B-cell maturation into plasma cells. Thus, we report a further example of how different post-translational modifications of a self-antigen exert complementary effects that sustain autoimmunity and inflammation, respectively. These data also show that T/B-cell responses to carbamylated-LL37 represent novel SLE disease biomarkers.

Transethnic, Genome-Wide Analysis Reveals Immune-Related Risk Alleles and Phenotypic Correlates in Pediatric Steroid-Sensitive Nephrotic Syndrome.

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.

Evaluation of the Common Molecular Basis in Alzheimer's and Parkinson's Diseases.

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.

Association of genetic variants in RAB23 and ANXA11 with uveitis in sarcoidosis.

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.

A novel approach for multi-SNP GWAS and its application in Alzheimer's disease.

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.

Association of Parkinson disease with structural and regulatory variants in the HLA region.

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_C(∗)07:02_DRB5(∗)01_DRB1(∗)15:01_DQA1(∗)01:02_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.

Genome-wide epistatic expression quantitative trait loci discovery in four human tissues reveals the importance of local chromosomal interactions governing gene expression.

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.

Multiple mismatches at the low expression HLA loci DP, DQ, and DRB3/4/5 associate with adverse outcomes in hematopoietic stem cell transplantation.

A single mismatch in highly expressed HLA-A, -B, -C, and -DRB1 loci (HEL) is associated with worse outcomes in hematopoietic stem cell transplantation, while less is known about the cumulative impact of mismatches in the lesser expressed HLA loci DRB3/4/5, DQ, and DP (LEL). We studied whether accumulation of LEL mismatches is associated with deleterious effects in 3853 unrelated donor transplants stratified according to number of matches in the HEL. In the 8/8 matched HEL group, LEL mismatches were not associated with any adverse outcome. Mismatches at HLA-DRB1 were associated with occurrence of multiple LEL mismatches. In the 7/8 HEL group, patients with 3 or more LEL mismatches scored in the graft-versus-host vector had a significantly higher risk of mortality (1.45 and 1.43) and transplant-related mortality (1.68 and 1.54) than the subgroups with 0 or 1 LEL mismatches. No single LEL locus had a more pronounced effect on clinical outcome. Three or more LEL mismatches are associated with lower survival after 7/8 HEL matched transplantation. Prospective evaluation of matching for HLA-DRB3/4/5, -DQ, and -DP loci is warranted to reduce posttransplant risks in donor-recipient pairs matched for 7/8 HEL.

HLA-DRB1, -DRB3, -DRB4 and -DRB5 genotyping at a super-high resolution level by long range PCR and high-throughput sequencing.

Super high-resolution single molecule sequence-based typing (SS-SBT) is a human leukocyte antigen (HLA) DNA typing method to the field 4 level of allelic resolution (formerly known as eight-digit typing) to efficiently detect new and null alleles without phase ambiguity by combination of long ranged polymerase chain reaction (PCR) amplification and next-generation sequencing (NGS) technologies. We previously reported the development and application of the SS-SBT method for the eight classical HLA loci, A, B, C, DRB1, DQA1, DQB1, DPA1 and DPB1. In this article, we describe the development of the SS-SBT method for three DRB1 linked loci, DRB3, DRB4 and DRB5 (DRB3/4/5) and characterization of DRB1-DRB3/4/5 haplotype structures to the field 4 level. Locus specific PCR primers for DRB3/4/5 were designed to amplify the gene regions from intron 1 to exon 6 [3' untranslated region (3'UTR)]. In total 20 DRB1 and 13 DRB3/4/5 allele sequences were determined by the SS-SBT to the field 4 level without phase ambiguity using 19 DR51, DR52 and DR53 positive genomic DNA samples obtained from Japanese. Moreover, 18 DRB1-DRB3/4/5 haplotypes were estimated to the field 4 level by the SS-SBT method in contrast to 10 haplotypes estimated by conventional methods to the field 1 level (formerly known as two digit typing). Therefore, DRB1-DRB3/4/5 haplotyping by SS-SBT is expected to provide informative data for improved HLA matching in medical research, transplantation procedures, HLA-related disease studies and human population diversity studies.

Myelin basic protein-specific TCR/HLA-DRB5*01:01 transgenic mice support the etiologic role of DRB5*01:01 in multiple sclerosis.

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.

Copy number variations of HLA-DRB5 is associated with systemic lupus erythematosus risk in Chinese Han population.

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.

Discovery of potential recombinant alleles HLA-DRB3*02:171 and HLA-DRB5*01:140.

Two novel alleles that appear to be generated from recombination between exons 2 and 3.

Detection of the novel HLA-DRB5*01:01:01:07 allele by Pacific Biosciences HiFi sequencing in a Chinese individual.

HLA-DRB5*01:01:01:07 differs from HLA-DRB5*01:01:01:01 by two nucleotide changes in intron 1 and intron 2.