Evaluation of Cell-Specific Epigenetic Age Acceleration in People With Multiple Sclerosis.

BACKGROUND AND OBJECTIVES: In multiple sclerosis (MS), accelerated aging of the immune system (immunosenescence) may be associated with disease onset or drive progression. DNA methylation (DNAm) is an epigenetic factor that varies among lymphocyte subtypes, and cell-specific DNAm is associated with MS. DNAm varies across the life span and can be used to accurately estimate biological age acceleration, which has been linked to a range of morbidities. The objective of this study was to test for cell-specific epigenetic age acceleration (EAA) in people with MS. METHODS: This was a case-control study of EAA using existing DNAm data from several independent previously published studies. Data were included if .idat files from Illumina 450K or EPIC arrays were available for both a case with MS and an age-matched and sex-matched control, from the same study. Multifactor statistical modeling was performed to assess the primary outcome of EAA. We explored the relationship of EAA and MS, including interaction terms to identify immune cell-specific effects. Cell-sorted DNA methylation data from 3 independent datasets were used to validate findings. RESULTS: We used whole blood DNA methylation data from 583 cases with MS and 643 non-MS controls to calculate EAA using the GrimAge algorithm. The MS group exhibited an increased EAA compared with controls (approximately 9 mths, 95% CI 3.6-14.4), p = 0.001). Statistical deconvolution showed that EAA is associated with MS in a B cell-dependent manner (ß int = 1.7, 95% CI 0.3-2.8), p = 0.002), irrespective of B-cell proportions. Validation analysis using 3 independent datasets enriched for B cells showed an EAA increase of 5.1 years in cases with MS compared with that in controls (95% CI 2.8-7.4, p = 5.5 × 10-5). By comparison, there was no EAA difference in MS in a T cell-enriched dataset. We found that EAA was attributed to the DNAm surrogates for Beta-2-microglobulin (difference = 47,546, 95% CI 10,067-85,026; p = 7.2 × 10-5), and smoking pack-years (difference = 8.1, 95% CI 1.9-14.2, p = 0.002). DISCUSSION: This study provides compelling evidence that B cells exhibit marked EAA in MS and supports the hypothesis that premature B-cell immune senescence plays a role in MS. Future MS studies should focus on age-related molecular mechanisms in B cells.

Exploring the role of the multiple sclerosis susceptibility gene CLEC16A in T cells.

C-type lectin-like domain family 16 member A (CLEC16A) is associated with autoimmune disorders, including multiple sclerosis (MS), but its functional relevance is not completely understood. CLEC16A is expressed in several immune cells, where it affects autophagic processes and receptor expression. Recently, we reported that the risk genotype of an MS-associated single nucleotide polymorphism in CLEC16A intron 19 is associated with higher expression of CLEC16A in CD4+ T cells. Here, we show that CLEC16A expression is induced in CD4+ T cells upon T cell activation. By the use of imaging flow cytometry and confocal microscopy, we demonstrate that CLEC16A is located in Rab4a-positive recycling endosomes in Jurkat TAg T cells. CLEC16A knock-down in Jurkat cells resulted in lower cell surface expression of the T cell receptor, however, this did not have a major impact on T cell activation response in vitro in Jurkat nor in human, primary CD4+ T cells.

Genetic variants are major determinants of CSF antibody levels in multiple sclerosis.

Immunological hallmarks of multiple sclerosis include the production of antibodies in the central nervous system, expressed as presence of oligoclonal bands and/or an increased immunoglobulin G index-the level of immunoglobulin G in the cerebrospinal fluid compared to serum. However, the underlying differences between oligoclonal band-positive and -negative patients with multiple sclerosis and reasons for variability in immunoglobulin G index are not known. To identify genetic factors influencing the variation in the antibody levels in the cerebrospinal fluid in multiple sclerosis, we have performed a genome-wide association screen in patients collected from nine countries for two traits, presence or absence of oligoclonal bands (n = 3026) and immunoglobulin G index levels (n = 938), followed by a replication in 3891 additional patients. We replicate previously suggested association signals for oligoclonal band status in the major histocompatibility complex region for the rs9271640*A-rs6457617*G haplotype, correlated with HLA-DRB1*1501, and rs34083746*G, correlated with HLA-DQA1*0301 (P comparing two haplotypes = 8.88 × 10(-16)). Furthermore, we identify a novel association signal of rs9807334, near the ELAC1/SMAD4 genes, for oligoclonal band status (P = 8.45 × 10(-7)). The previously reported association of the immunoglobulin heavy chain locus with immunoglobulin G index reaches strong evidence for association in this data set (P = 3.79 × 10(-37)). We identify two novel associations in the major histocompatibility complex region with immunoglobulin G index: the rs9271640*A-rs6457617*G haplotype (P = 1.59 × 10(-22)), shared with oligoclonal band status, and an additional independent effect of rs6457617*G (P = 3.68 × 10(-6)). Variants identified in this study account for up to 2-fold differences in the odds of being oligoclonal band positive and 7.75% of the variation in immunoglobulin G index. Both traits are associated with clinical features of disease such as female gender, age at onset and severity. This is the largest study population so far investigated for the genetic influence on antibody levels in the cerebrospinal fluid in multiple sclerosis, including 6950 patients. We confirm that genetic factors underlie these antibody levels and identify both the major histocompatibility complex and immunoglobulin heavy chain region as major determinants.

Retinoic acid enhances the levels of IL-10 in TLR-stimulated B cells from patients with relapsing-remitting multiple sclerosis.

We have explored the beneficial effects of retinoic acid (RA) on B cells from multiple sclerosis (MS) patients. When co-stimulated via the toll-like receptors (TLRs) TLR9 and RP105, MS B cells secreted less of the anti-inflammatory cytokine interleukin 10 (IL-10) compared to B cells from healthy controls. Importantly, RA enhanced the secretion of IL-10 by MS-derived B cells without affecting the levels of the pro-inflammatory cytokine TNF-α. RA revealed the same ability to induce IL-10 as did interferon-ß-1b (IFN-ß-1b), and B-cells from patients treated with glatiramer acetate or IFN-ß-1b still displayed the beneficial effects of RA on the IL-10/TNF-α ratio.

A gain of function mutation in TNFRSF11B encoding osteoprotegerin causes osteoarthritis with chondrocalcinosis.

OBJECTIVE: To identify pathogenic mutations that reveal underlying biological mechanisms driving osteoarthritis (OA). METHODS: Exome sequencing was applied to two distant family members with dominantly inherited early onset primary OA at multiple joint sites with chondrocalcinosis (familial generalised osteoarthritis, FOA). Confirmation of mutations occurred by genotyping and linkage analyses across the extended family. The functional effect of the mutation was investigated by means of a cell-based assay. To explore generalisability, mRNA expression analysis of the relevant genes in the discovered pathway was explored in preserved and osteoarthritic articular cartilage of independent patients undergoing joint replacement surgery. RESULTS: We identified a heterozygous, probably damaging, read-through mutation (c.1205A=>T; p.Stop402Leu) in TNFRSF11B encoding osteoprotegerin that is likely causal to the OA phenotype in the extended family. In a bone resorption assay, the mutant form of osteoprotegerin showed enhanced capacity to inhibit osteoclastogenesis and bone resorption. Expression analyses in preserved and affected articular cartilage of independent OA patients showed that upregulation of TNFRSF11B is a general phenomenon in the pathophysiological process. CONCLUSIONS: Albeit that the role of the molecular pathway of osteoprotegerin has been studied in OA, we are the first to demonstrate that enhanced osteoprotegerin function could be a directly underlying cause. We advocate that agents counteracting the function of osteoprotegerin could comply with new therapeutic interventions of OA.

A meta-analysis of genome-wide association studies identifies novel variants associated with osteoarthritis of the hip.

OBJECTIVES: Osteoarthritis (OA) is the most common form of arthritis with a clear genetic component. To identify novel loci associated with hip OA we performed a meta-analysis of genome-wide association studies (GWAS) on European subjects. METHODS: We performed a two-stage meta-analysis on more than 78,000 participants. In stage 1, we synthesised data from eight GWAS whereas data from 10 centres were used for 'in silico' or 'de novo' replication. Besides the main analysis, a stratified by sex analysis was performed to detect possible sex-specific signals. Meta-analysis was performed using inverse-variance fixed effects models. A random effects approach was also used. RESULTS: We accumulated 11,277 cases of radiographic and symptomatic hip OA. We prioritised eight single nucleotide polymorphism (SNPs) for follow-up in the discovery stage (4349 OA cases); five from the combined analysis, two male specific and one female specific. One locus, at 20q13, represented by rs6094710 (minor allele frequency (MAF) 4%) near the NCOA3 (nuclear receptor coactivator 3) gene, reached genome-wide significance level with p=7.9×10(-9) and OR=1.28 (95% CI 1.18 to 1.39) in the combined analysis of discovery (p=5.6×10(-8)) and follow-up studies (p=7.3×10(-4)). We showed that this gene is expressed in articular cartilage and its expression was significantly reduced in OA-affected cartilage. Moreover, two loci remained suggestive associated; rs5009270 at 7q31 (MAF 30%, p=9.9×10(-7), OR=1.10) and rs3757837 at 7p13 (MAF 6%, p=2.2×10(-6), OR=1.27 in male specific analysis). CONCLUSIONS: Novel genetic loci for hip OA were found in this meta-analysis of GWAS.

GREM1, FRZB and DKK1 mRNA levels correlate with osteoarthritis and are regulated by osteoarthritis-associated factors.

INTRODUCTION: Osteoarthritis is, at least in a subset of patients, associated with hypertrophic differentiation of articular chondrocytes. Recently, we identified the bone morphogenetic protein (BMP) and wingless-type MMTV integration site (WNT) signaling antagonists Gremlin 1 (GREM1), frizzled-related protein (FRZB) and dickkopf 1 homolog (Xenopus laevis) (DKK1) as articular cartilage's natural brakes of hypertrophic differentiation. In this study, we investigated whether factors implicated in osteoarthritis or regulation of chondrocyte hypertrophy influence GREM1, FRZB and DKK1 expression levels. METHODS: GREM1, FRZB and DKK1 mRNA levels were studied in articular cartilage from healthy preadolescents and healthy adults as well as in preserved and degrading osteoarthritic cartilage from the same osteoarthritic joint by quantitative PCR. Subsequently, we exposed human articular chondrocytes to WNT, BMP, IL-1ß, Indian hedgehog, parathyroid hormone-related peptide, mechanical loading, different medium tonicities or distinct oxygen levels and investigated GREM1, FRZB and DKK1 expression levels using a time-course analysis. RESULTS: GREM1, FRZB and DKK1 mRNA expression were strongly decreased in osteoarthritis. Moreover, this downregulation is stronger in degrading cartilage compared with macroscopically preserved cartilage from the same osteoarthritic joint. WNT, BMP, IL-1ß signaling and mechanical loading regulated GREM1, FRZB and DKK1 mRNA levels. Indian hedgehog, parathyroid hormone-related peptide and tonicity influenced the mRNA levels of at least one antagonist, while oxygen levels did not demonstrate any statistically significant effect. Interestingly, BMP and WNT signaling upregulated the expression of each other's antagonists. CONCLUSIONS: Together, the current study demonstrates an inverse correlation between osteoarthritis and GREM1, FRZB and DKK1 gene expression in cartilage and provides insight into the underlying transcriptional regulation. Furthermore, we show that BMP and WNT signaling are linked in a negative feedback loop, which might prove essential in articular cartilage homeostasis by balancing BMP and WNT activity.

A genome-wide linkage scan reveals CD53 as an important regulator of innate TNF-alpha levels.

Cytokines are major immune system regulators. Previously, innate cytokine profiles determined by lipopolysaccharide stimulation were shown to be highly heritable. To identify regulating genes in innate immunity, we analyzed data from a genome-wide linkage scan using microsatellites in osteoarthritis (OA) patients (The GARP study) and their innate cytokine data on interleukin (IL)-1beta, IL-1Ra, IL-10 and tumor necrosis factor (TNF)alpha. A confirmation cohort consisted of the Leiden 85-Plus study. In this study, a linkage analysis was followed by manual selection of candidate genes in linkage regions showing LOD scores over 2.5. An single-nucleotide polymorphism (SNP) gene tagging method was applied to select SNPs on the basis of the highest level of gene tagging and possible functional effects. QTDT was used to identify the SNPs associated with innate cytokine production. Initial association signals were modeled by a linear mixed model. Through these analyses, we identified 10 putative genes involved in the regulation of TNFalpha. SNP rs6679497 in gene CD53 showed significant association with TNFalpha levels (P=0.001). No association of this SNP was observed with OA. A novel gene involved in the innate immune response of TNFalpha is identified. Genetic variation in this gene may have a role in diseases and disorders in which TNFalpha is closely involved.

Incidence of gastrointestinal stromal tumours is underestimated: results of a nation-wide study.

Gastrointestinal stromal tumours (GIST) are identified by their specific morphology added with immunohistochemical staining with anti-CD117, the phenotypic marker for GIST in the majority of cases. In this study the reported incidence of GIST and GIST-like tumours before and after the routine availability of the specific diagnostic marker CD117 antigen was investigated. All patients with GIST or GIST-like tumours were selected from PALGA, a nation-wide network and registry of histo- and cytopathology in the Netherlands, to calculate the incidence in 1995 and longitudinally over time between 1998 and 2003. Pathological reports were retrieved to assess the type of immunostaining used and to assess the risk category for malignant behaviour according to the recently published consensus criteria. The annual incidence of GIST as distilled from the national pathology registry increased from 2.1 per million inhabitants in 1995 to 12.7 per million inhabitants in 2003. The incidence of GIST-like tumours decreased from 17.6 per million inhabitants in 1995 to 12.7 per million inhabitants in 2003. The incidences were stable from 2000 onwards. Additional analysis in 2003 indicated that more than 90% of the GIST tested was CD117 positive, compared to only 4% of the GIST-like tumours. Almost 50% of the GIST was considered to be at high-risk for malignant behaviour, according to the consensus criteria. The increased incidence of GIST 1995-2003 is related to increased understanding of GIST pathobiology and the routine availability of the diagnostic immunohistochemical antibody directed against the CD117 antigen.