HLA-C alleles confer risk for anti-citrullinated peptide antibody-positive rheumatoid arthritis independent of HLA-DRB1 alleles.

OBJECTIVE: The MHC exerts the greatest contribution to RA susceptibility, where certain HLA-DRB1 alleles confer the greatest risk. Interestingly, there is evidence for more risk factors in the MHC with regions surrounding the HLA class I loci, but whether these antigen-presenting loci could be causal risk variants has not been directly investigated. In this study we investigate the HLA association by direct genotyping of the HLA loci. METHODS: Nine hundred and fifty RA patients and 933 healthy controls were genotyped for HLA-A, -B and -C. Eleven single-nucleotide polymorphisms (SNPs) and one insertion/deletion in the MHC were also included. Conditional logistic regression analyses were performed separately in ACPA-positive and -negative RA to identify the strongest susceptibility locus and additional risk loci. RESULTS: In ACPA-positive RA, the most significantly associated locus was HLA-DRB1 (P = 1.58 × 10(-54)), with SE alleles being predisposing. After controlling for HLA-DRB1, the HLA-C locus was found to confer susceptibility (P = 2.32 × 10(-9)), particularly, the HLA-C*03 allele. Also, in ACPA-negative RA, HLA-DRB1 was the most significant locus (P = 7.22 × 10(-9)), but with other risk alleles (particularly DRB1*03). A possible independent involvement of HLA-C was also observed for ACPA-negative RA (P = 0.02). CONCLUSION: HLA-DRB1 was the major MHC risk locus in both ACPA-positive and ACPA-negative RA, but with allelic risk heterogeneity. Joint analyses of the HLA class I loci together with previously proposed SNP associations pointed at HLA-C as a second susceptibility locus in ACPA-positive RA.

Small extracellular vesicles have distinct CD81 and CD9 tetraspanin expression profiles in plasma from rheumatoid arthritis patients.

Extracellular vesicles (EVs) are implicated in the pathogenesis of rheumatoid arthritis (RA) but little is known about the composition of specific small EV (sEV) subpopulations. This study aimed to characterize the CD63, CD81 and CD9 tetraspanin profile in the membrane of single EVs in plasma from treatment naïve RA patients and assess potential discrepancies between methotrexate (MTX) responder groups. EVs isolated from plasma were characterized using transmission electron microscopy, and detection of surface markers (CD63, CD81 and CD9) on single EVs was performed on the ExoView platform. All RA patients (N = 8) were newly diagnosed, treatment naïve, females, ACPA positive and former smokers. The controls (N = 5) were matched for age and gender. After three months of MTX treatment, responders (N = 4) were defined as those with ΔDAS28 > 1.2 and DAS28 ≤ 3.2 post-treatment. The isolated EVs were 50-200 nm in size. The RA patients had a higher proportion of both CD9 and CD81 single positive sEVs compared to healthy controls, while there was a decrease in CD81/CD9 double positive sEVs in patients. Stratification of RA patients into MTX responders and non-responders revealed a distinctly higher proportion of CD81 single positive sEVs in the responder group. The proportion of CD81/CD9 double positive sEVs (anti-CD9 captured) was lower in the non-responders, but increased upon 3 months of MTX treatment. Our exploratory study revealed distinct tetraspanin profiles in RA patients suggesting their implication in RA pathophysiology and MTX treatment response.

Methotrexate Treatment of Newly Diagnosed RA Patients Is Associated With DNA Methylation Differences at Genes Relevant for Disease Pathogenesis and Pharmacological Action.

Background: Methotrexate (MTX) is the first line treatment of rheumatoid arthritis (RA), and methylation changes in bulk T cells have been reported after treatment with MTX. We have investigated cell-type specific DNA methylation changes across the genome in naïve and memory CD4+ T cells before and after MTX treatment of RA patients. DNA methylation profiles of newly diagnosed RA patients (N=9) were assessed by reduced representation bisulfite sequencing. Results: We found that MTX treatment significantly influenced DNA methylation levels at multiple CpG sites in both cell populations. Interestingly, we identified differentially methylated sites annotated to two genes; TRIM15 and SORC2, previously reported to predict treatment outcome in RA patients when measured in bulk T cells. Furthermore, several of the genes, including STAT3, annotated to the significant CpG sites are relevant for RA susceptibility or the action of MTX. Conclusion: We detected CpG sites that were associated with MTX treatment in CD4+ naïve and memory T cells isolated from RA patients. Several of these sites overlap genetic regions previously associated with RA risk and MTX treatment outcome.

MicroRNA Expression Differences in Blood-Derived CD19+ B Cells of Methotrexate Treated Rheumatoid Arthritis Patients.

Rheumatoid arthritis (RA) is a complex disease with a wide range of underlying susceptibility factors. Recently, dysregulation of microRNAs (miRNAs) in RA have been reported in several immune cell types from blood. However, B cells have not been studied in detail yet. Given the autoimmune nature of RA with the presence of autoantibodies, CD19+ B cells are a key cell type in RA pathogenesis and alterations in CD19+ B cell subpopulations have been observed in patient blood. Therefore, we aimed to reveal the global miRNA repertoire and to analyze miRNA expression profile differences in homogenous RA patient phenotypes in blood-derived CD19+ B cells. Small RNA sequencing was performed on CD19+ B cells of newly diagnosed untreated RA patients (n=10), successfully methotrexate (MTX) treated RA patients in remission (MTX treated RA patients, n=18) and healthy controls (n=9). The majority of miRNAs was detected across all phenotypes. However, significant expression differences between MTX treated RA patients and controls were observed for 27 miRNAs, while no significant differences were seen between the newly diagnosed patients and controls. Several of the differentially expressed miRNAs were previously found to be dysregulated in RA including miR-223-3p, miR-486-3p and miR-23a-3p. MiRNA target enrichment analysis, using the differentially expressed miRNAs and miRNA-target interactions from miRTarBase as input, revealed enriched target genes known to play important roles in B cell activation, differentiation and B cell receptor signaling, such as STAT3, PRDM1 and PTEN. Interestingly, many of those genes showed a high degree of correlated expression in CD19+ B cells in contrast to other immune cell types. Our results suggest important regulatory functions of miRNAs in blood-derived CD19+ B cells of MTX treated RA patients and motivate for future studies investigating the interactive mechanisms between miRNA and gene targets, as well as the possible predictive power of miRNAs for RA treatment response.

Human Leukocyte Antigen alleles associated with Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS).

The etiology and pathogenesis of Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) are unknown, and autoimmunity is one of many proposed underlying mechanisms. Human Leukocyte Antigen (HLA) associations are hallmarks of autoimmune disease, and have not been thoroughly investigated in a large ME/CFS patient cohort. We performed high resolution HLA -A, -B, -C, -DRB1, -DQB1 and -DPB1 genotyping by next generation sequencing in 426 adult, Norwegian ME/CFS patients, diagnosed according to the Canadian Consensus Criteria. HLA associations were assessed by comparing to 4511 healthy and ethnically matched controls. Clinical information was collected through questionnaires completed by patients or relatives. We discovered two independent HLA associations, tagged by the alleles HLA-C*07:04 (OR 2.1 [95% CI 1.4-3.1]) and HLA-DQB1*03:03 (OR 1.5 [95% CI 1.1-2.0]). These alleles were carried by 7.7% and 12.7% of ME/CFS patients, respectively. The proportion of individuals carrying one or both of these alleles was 19.2% in the patient group and 12.2% in the control group (OR 1.7 [95% CI 1.3-2.2], pnc = 0.00003). ME/CFS is a complex disease, potentially with a substantial heterogeneity. We report novel HLA associations pointing toward the involvement of the immune system in ME/CFS pathogenesis.

Rheumatoid Arthritis Patients, Both Newly Diagnosed and Methotrexate Treated, Show More DNA Methylation Differences in CD4+ Memory Than in CD4+ Naïve T Cells.

Background: Differences in DNA methylation have been reported in B and T lymphocyte populations, including CD4+ T cells, isolated from rheumatoid arthritis (RA) patients when compared to healthy controls. CD4+ T cells are a heterogeneous cell type with subpopulations displaying distinct DNA methylation patterns. In this study, we investigated DNA methylation using reduced representation bisulfite sequencing in two CD4+ T cell populations (CD4+ memory and naïve cells) in three groups: newly diagnosed, disease modifying antirheumatic drugs (DMARD) naïve RA patients (N = 11), methotrexate (MTX) treated RA patients (N = 18), and healthy controls (N = 9) matched for age, gender and smoking status. Results: Analyses of these data revealed significantly more differentially methylated positions (DMPs) in CD4+ memory than in CD4+ naïve T cells (904 vs. 19 DMPs) in RA patients compared to controls. The majority of DMPs (72%) identified in newly diagnosed and DMARD naïve RA patients with active disease showed increased DNA methylation (39 DMPs), whereas most DMPs (80%) identified in the MTX treated RA patients in remission displayed decreased DNA methylation (694 DMPs). Interestingly, we also found that about one third of the 101 known RA risk loci overlapped (±500 kb) with the DMPs. Notably, introns of the UBASH3A gene harbor both the lead RA risk SNP and two DMPs in CD4+ memory T cells. Conclusion: Our results suggest that RA associated DNA methylation differences vary between the two T cell subsets, but are also influenced by RA characteristics such as disease activity, disease duration and/or MTX treatment.

Transcriptomes of antigen presenting cells in human thymus.

Antigen presenting cells (APCs) in the thymus play an essential role in the establishment of central tolerance, i.e. the generation of a repertoire of functional and self-tolerant T cells to prevent autoimmunity. In this study, we have compared the transcriptomes of four primary APCs from human thymus (mTECs, CD19+ B cells, CD141+ and CD123+ DCs). We investigated a set of genes including the HLA genes, genes encoding transcriptional regulators and finally, tissue-enriched genes, i.e, genes with a five-fold higher expression in a particular human tissue. We show that thymic CD141+ DCs express the highest levels of all classical HLA genes and 67% (14/21) of the HLA class I and II pathway genes investigated in this study. CD141+ DCs also expressed the highest levels of the transcriptional regulator DEAF1, whereas AIRE and FEZF2 expression were mainly found in primary human mTECs. We found expression of "tissue enriched genes" from the Human Protein Atlas (HPA) in all four APC types, but the mTECs were clearly dominating in the number of uniquely expressed tissue enriched genes (20% in mTECs, 7% in CD19+ B cells, 4% in CD123+ DCs and 2% in CD141+ DCs). The tissue enriched genes also overlapped with reported human autoantigens. This is, to our knowledge, the first study that performs RNA sequencing of mTECs, CD19+ B cells, CD141+ and CD123+ DCs isolated from the same individuals and provides insight into the transcriptomes of these human thymic APCs.

Association analysis in type 1 diabetes of the PRSS16 gene encoding a thymus-specific serine protease.

We have previously mapped a separate type 1 diabetes (T1D) association in the extended MHC class I region, marked by D6S2223, on the DRB1*03-DQA1*0501-DQB1*0201 haplotype. The associated region encompasses a gene encoding a thymus-specific serine protease (PRSS16), presumably involved in positive selection of T cells or in T-cell regulation. Fourteen PRSS16 polymorphisms were genotyped in two steps using a total of six T1D family data sets, as well as case-control materials for both T1D and celiac disease (CD). An association with a 15 base-pair deletion in exon 12 of PRSS16 was found on the DRB1*03-DQA1*0501-DQB1*0201 haplotype for both T1D and CD, but it could not explain the more pronounced disease associations observed at marker D6S2223. We compared the performance of the 14 tested PRSS16 polymorphisms, selected after our previous comprehensive screen, against HapMap selected tag SNPs. Use of a HapMap based SNP selection strategy would result in loss of a large proportion of the genetic variation in PRSS16. Our data suggest that it is unlikely that polymorphisms within the PRSS16 gene are involved in the predisposition to T1D. However, we cannot rule out that regulatory polymorphisms located some distance away from the gene may be involved.

Oligoclonal band status in Scandinavian multiple sclerosis patients is associated with specific genetic risk alleles.

The presence of oligoclonal bands (OCB) in cerebrospinal fluid (CSF) is a typical finding in multiple sclerosis (MS). We applied data from Norwegian, Swedish and Danish (i.e. Scandinavian) MS patients from a genome-wide association study (GWAS) to search for genetic differences in MS relating to OCB status. GWAS data was compared in 1367 OCB positive and 161 OCB negative Scandinavian MS patients, and nine of the most associated SNPs were genotyped for replication in 3403 Scandinavian MS patients. HLA-DRB1 genotypes were analyzed in a subset of the OCB positive (n = 2781) and OCB negative (n = 292) MS patients and compared to 890 healthy controls. Results from the genome-wide analyses showed that single nucleotide polymorphisms (SNPs) from the HLA complex and six other loci were associated to OCB status. In SNPs selected for replication, combined analyses showed genome-wide significant association for two SNPs in the HLA complex; rs3129871 (p = 5.7×10(-15)) and rs3817963 (p = 5.7×10(-10)) correlating with the HLA-DRB1*15 and the HLA-DRB1*04 alleles, respectively. We also found suggestive association to one SNP in the Calsyntenin-2 gene (p = 8.83×10(-7)). In HLA-DRB1 analyses HLA-DRB1*15∶01 was a stronger risk factor for OCB positive than OCB negative MS, whereas HLA-DRB1*04∶04 was associated with increased risk of OCB negative MS and reduced risk of OCB positive MS. Protective effects of HLA-DRB1*01∶01 and HLA-DRB1*07∶01 were detected in both groups. The groups were different with regard to age at onset (AAO), MS outcome measures and gender. This study confirms both shared and distinct genetic risk for MS subtypes in the Scandinavian population defined by OCB status and indicates different clinical characteristics between the groups. This suggests differences in disease mechanisms between OCB negative and OCB positive MS with implications for patient management, which need to be further studied.

Late onset myasthenia gravis is associated with HLA DRB1*15:01 in the Norwegian population.

BACKGROUND: Acquired myasthenia gravis (MG) is a rare antibody-mediated autoimmune disease caused by impaired neuromuscular transmission, leading to abnormal muscle fatigability. The aetiology is complex, including genetic risk factors of the human leukocyte antigen (HLA) complex and unknown environmental factors. Although associations between the HLA complex and MG are well established, not all involved components of the HLA predisposition to this heterogeneous disease have been revealed. Well-powered and comprehensive HLA analyses of subgroups in MG are warranted, especially in late onset MG. METHODOLOGY/PRINCIPAL FINDINGS: This case-control association study is of a large population-based Norwegian cohort of 369 MG patients and 651 healthy controls. We performed comprehensive genotyping of four classical HLA loci (HLA-A, -B, -C and -DRB1) and showed that the DRB1*15:01 allele conferred the strongest risk in late onset MG (LOMG; onset ≥ 60 years) (OR 2.38, p(c)7.4 × 10(-5)). DRB1*13:01 was found to be a protective allele for both early onset MG (EOMG) and LOMG (OR 0.31, p(c) 4.71 × 10(-4)), a finding not previously described. No significant association was found to the DRB1*07:01 allele (p(nc) = 0.18) in a subset of nonthymomatous anti-titin antibody positive LOMG as reported by others. HLA-B*08 was mapped to give the strongest contribution to EOMG, supporting previous studies. CONCLUSION: The results from this study provide important new information concerning the susceptibility of HLA alleles in Caucasian MG, with highlights on DRB1*15:01 as being a major risk allele in LOMG.