Multiscale networks in multiple sclerosis.

Complex diseases such as Multiple Sclerosis (MS) cover a wide range of biological scales, from genes and proteins to cells and tissues, up to the full organism. In fact, any phenotype for an organism is dictated by the interplay among these scales. We conducted a multilayer network analysis and deep phenotyping with multi-omics data (genomics, phosphoproteomics and cytomics), brain and retinal imaging, and clinical data, obtained from a multicenter prospective cohort of 328 patients and 90 healthy controls. Multilayer networks were constructed using mutual information for topological analysis, and Boolean simulations were constructed using Pearson correlation to identified paths within and among all layers. The path more commonly found from the Boolean simulations connects protein MK03, with total T cells, the thickness of the retinal nerve fiber layer (RNFL), and the walking speed. This path contains nodes involved in protein phosphorylation, glial cell differentiation, and regulation of stress-activated MAPK cascade, among others. Specific paths identified were subsequently analyzed by flow cytometry at the single-cell level. Combinations of several proteins (GSK3AB, HSBP1 or RS6) and immune cells (Th17, Th1 non-classic, CD8, CD8 Treg, CD56 neg, and B memory) were part of the paths explaining the clinical phenotype. The advantage of the path identified from the Boolean simulations is that it connects information about these known biological pathways with the layers at higher scales (retina damage and disability). Overall, the identified paths provide a means to connect the molecular aspects of MS with the overall phenotype.

Quantitative proteomics reveals protein dysregulation during T cell activation in multiple sclerosis patients compared to healthy controls.

BACKGROUND: Multiple sclerosis (MS) is an autoimmune, neurodegenerative disorder with a strong genetic component that acts in a complex interaction with environmental factors for disease development. CD4+ T cells are pivotal players in MS pathogenesis, where peripherally activated T cells migrate to the central nervous system leading to demyelination and axonal degeneration. Through a proteomic approach, we aim at identifying dysregulated pathways in activated T cells from MS patients as compared to healthy controls. METHODS: CD4+ T cells were purified from peripheral blood from MS patients and healthy controls by magnetic separation. Cells were left unstimulated or stimulated in vitro through the TCR and costimulatory CD28 receptor for 24 h prior to sampling. Electrospray liquid chromatography-tandem mass spectrometry was used to measure protein abundances. RESULTS: Upon T cell activation the abundance of 1801 proteins was changed. Among these proteins, we observed an enrichment of proteins expressed by MS-susceptibility genes. When comparing protein abundances in T cell samples from healthy controls and MS patients, 18 and 33 proteins were differentially expressed in unstimulated and stimulated CD4+ T cells, respectively. Moreover, 353 and 304 proteins were identified as proteins exclusively induced upon T cell activation in healthy controls and MS patients, respectively and dysregulation of the Nur77 pathway was observed only in samples from MS patients. CONCLUSIONS: Our study highlights the importance of CD4+ T cell activation for MS, as proteins that change in abundance upon T cell activation are enriched for proteins encoded by MS susceptibility genes. The results provide evidence for proteomic disturbances in T cell activation in MS, and pinpoint to dysregulation of the Nur77 pathway, a biological pathway known to limit aberrant effector T cell responses.

Serum neurofilament light chain concentration predicts disease worsening in multiple sclerosis.

BACKGROUND: Serum neurofilament light (sNfL) chain is a promising biomarker reflecting neuro-axonal injury in multiple sclerosis (MS). However, the ability of sNfL to predict outcomes in real-world MS cohorts requires further validation. OBJECTIVE: The aim of the study is to investigate the associations of sNfL concentration, magnetic resonance imaging (MRI) and retinal optical coherence tomography (OCT) markers with disease worsening in a longitudinal European multicentre MS cohort. METHODS: MS patients (n = 309) were prospectively enrolled at four centres and re-examined after 2 years (n = 226). NfL concentration was measured by single molecule array assay in serum. The patients' phenotypes were thoroughly characterized with clinical examination, retinal OCT and MRI brain scans. The primary outcome was disease worsening at median 2-year follow-up. RESULTS: Patients with high sNfL concentrations (⩾8 pg/mL) at baseline had increased risk of disease worsening at median 2-year follow-up (odds ratio (95% confidence interval) = 2.8 (1.5-5.3), p = 0.001). We found no significant associations of MRI or OCT measures at baseline with risk of disease worsening. CONCLUSION: Serum NfL concentration was the only factor associated with disease worsening, indicating that sNfL is a useful biomarker in MS that might be relevant in a clinical setting.

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.

Allelic imbalance of multiple sclerosis susceptibility genes IKZF3 and IQGAP1 in human peripheral blood.

BACKGROUND: Multiple sclerosis is a chronic inflammatory, demyelinating disease of the central nervous system. Recent genome-wide studies have revealed more than 110 single nucleotide polymorphisms as associated with susceptibility to multiple sclerosis, but their functional contribution to disease development is mostly unknown. RESULTS: Consistent allelic imbalance was observed for rs907091 in IKZF3 and rs11609 in IQGAP1, which are in strong linkage disequilibrium with the multiple sclerosis associated single nucleotide polymorphisms rs12946510 and rs8042861, respectively. Using multiple sclerosis patients and healthy controls heterozygous for rs907091 and rs11609, we showed that the multiple sclerosis risk alleles at IKZF3 and IQGAP1 are expressed at higher levels as compared to the protective allele. Furthermore, individuals homozygous for the multiple sclerosis risk allele at IQGAP1 had a significantly higher total expression of IQGAP1 compared to individuals homozygous for the protective allele. CONCLUSIONS: Our data indicate a possible regulatory role for the multiple sclerosis-associated IKZF3 and IQGAP1 variants. We suggest that such cis-acting mechanisms may contribute to the multiple sclerosis association of single nucleotide polymorphisms at IKZF3 and IQGAP1.

Genetic overlap between multiple sclerosis and several cardiovascular disease risk factors.

BACKGROUND: Epidemiological findings suggest a relationship between multiple sclerosis (MS) and cardiovascular disease (CVD) risk factors, although the nature of this relationship is not well understood. OBJECTIVE: We used genome-wide association study (GWAS) data to identify shared genetic factors (pleiotropy) between MS and CVD risk factors. METHODS: Using summary statistics from a large, recent GWAS (total n > 250,000 individuals), we investigated overlap in single nucleotide polymorphisms (SNPs) associated with MS and a number of CVD risk factors including triglycerides (TG), low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol, body mass index, waist-to-hip ratio, type 2 diabetes, systolic blood pressure, and C-reactive protein level. RESULTS AND CONCLUSION: Using conditional enrichment plots, we found 30-fold enrichment of MS SNPs for different levels of association with LDL and TG SNPs, with a corresponding reduction in conditional false discovery rate (FDR). We identified 133 pleiotropic loci outside the extended major histocompatibility complex with conditional FDR < 0.01, of which 65 are novel. These pleiotropic loci were located on 21 different chromosomes. Our findings point to overlapping pathobiology between clinically diagnosed MS and cardiovascular risk factors and identify novel common variants associated with increased MS risk.

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.

A meta-analysis of thyroid-related traits reveals novel loci and gender-specific differences in the regulation of thyroid function.

Thyroid hormone is essential for normal metabolism and development, and overt abnormalities in thyroid function lead to common endocrine disorders affecting approximately 10% of individuals over their life span. In addition, even mild alterations in thyroid function are associated with weight changes, atrial fibrillation, osteoporosis, and psychiatric disorders. To identify novel variants underlying thyroid function, we performed a large meta-analysis of genome-wide association studies for serum levels of the highly heritable thyroid function markers TSH and FT4, in up to 26,420 and 17,520 euthyroid subjects, respectively. Here we report 26 independent associations, including several novel loci for TSH (PDE10A, VEGFA, IGFBP5, NFIA, SOX9, PRDM11, FGF7, INSR, ABO, MIR1179, NRG1, MBIP, ITPK1, SASH1, GLIS3) and FT4 (LHX3, FOXE1, AADAT, NETO1/FBXO15, LPCAT2/CAPNS2). Notably, only limited overlap was detected between TSH and FT4 associated signals, in spite of the feedback regulation of their circulating levels by the hypothalamic-pituitary-thyroid axis. Five of the reported loci (PDE8B, PDE10A, MAF/LOC440389, NETO1/FBXO15, and LPCAT2/CAPNS2) show strong gender-specific differences, which offer clues for the known sexual dimorphism in thyroid function and related pathologies. Importantly, the TSH-associated loci contribute not only to variation within the normal range, but also to TSH values outside the reference range, suggesting that they may be involved in thyroid dysfunction. Overall, our findings explain, respectively, 5.64% and 2.30% of total TSH and FT4 trait variance, and they improve the current knowledge of the regulation of hypothalamic-pituitary-thyroid axis function and the consequences of genetic variation for hypo- or hyperthyroidism.

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.

Underlying molecular mechanisms of DIO2 susceptibility in symptomatic osteoarthritis.

OBJECTIVES: To investigate how the genetic susceptibility gene DIO2 confers risk to osteoarthritis (OA) onset in humans and to explore whether counteracting the deleterious effect could contribute to novel therapeutic approaches. METHODS: Epigenetically regulated expression of DIO2 was explored by assessing methylation of positional CpG-dinucleotides and the respective DIO2 expression in OA-affected and macroscopically preserved articular cartilage from end-stage OA patients. In a human in vitro chondrogenesis model, we measured the effects when thyroid signalling during culturing was either enhanced (excess T3 or lentiviral induced DIO2 overexpression) or decreased (iopanoic acid). RESULTS: OA-related changes in methylation at a specific CpG dinucleotide upstream of DIO2 caused significant upregulation of its expression (ß=4.96; p=0.0016). This effect was enhanced and appeared driven specifically by DIO2 rs225014 risk allele carriers (ß=5.58, p=0.0006). During in vitro chondrogenesis, DIO2 overexpression resulted in a significant reduced capacity of chondrocytes to deposit extracellular matrix (ECM) components, concurrent with significant induction of ECM degrading enzymes (ADAMTS5, MMP13) and markers of mineralisation (ALPL, COL1A1). Given their concurrent and significant upregulation of expression, this process is likely mediated via HIF-2α/RUNX2 signalling. In contrast, we showed that inhibiting deiodinases during in vitro chondrogenesis contributed to prolonged cartilage homeostasis as reflected by significant increased deposition of ECM components and attenuated upregulation of matrix degrading enzymes. CONCLUSIONS: Our findings show how genetic variation at DIO2 could confer risk to OA and raised the possibility that counteracting thyroid signalling may be a novel therapeutic approach.

Genes expressed in blood link osteoarthritis with apoptotic pathways.

OBJECTIVE: To identify novel gene expression networks in blood of osteoarthritis patients compared to controls. METHODS: A comprehensive exploration of gene expression in peripheral blood was performed by microarray analysis for a subset of osteoarthritis patients from the Genetics osteoARthritis and Progression (GARP) study in comparison with sex and age-matched healthy controls. To identify pathways, we performed gene enrichment analyses (database for annotation, visualisation and integrated discovery and search tool for the retrieval of interacting genes). Quantitative PCR analysis in overlapping and in additional osteoarthritis samples was performed for prioritised genes to validate and replicate findings. Classification of cases and controls was explored by applying statistical models. RESULTS: 741 probes representing 694 unique genes were differentially expressed between cases and controls, including 86 genes expressed with at least a 1.5-fold difference. ATF4, GPR18 and H3F3B were among the top genes identified (p<4.5 × 10(-8)). We found that in the blood of osteoarthritis patients the apoptosis pathway, including the well-known gene CASP3, was significantly enriched among the differentially expressed genes. Our findings were validated in independent samples and when using a small subset of the validated genes, we could accurately distinguish patients from controls (area under the curve 98%). CONCLUSIONS: In the current study, we have identified specific gene expression networks, in the easily accessible tissue blood, which associated consistently with osteoarthritis among GARP study cases. Our data further hint at the relevance of apoptosis as an aetiological factor in osteoarthritis onset, thereby qualifying expression profiling of blood as a useful tool to understand the underlying molecular mechanisms of osteoarthritis.

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.

Oligoclonal bands and age at onset correlate with genetic risk score in multiple sclerosis.

BACKGROUND: Many genetic risk variants are now well established in multiple sclerosis (MS), but the impact on clinical phenotypes is unclear. OBJECTIVE: To investigate the impact of established MS genetic risk variants on MS phenotypes, in well-characterized MS cohorts. METHODS: Norwegian MS patients (n = 639) and healthy controls (n = 530) were successfully genotyped for 61 established MS-associated single nucleotide polymorphisms (SNPs). Data including and excluding Major Histocompatibility Complex (MHC) markers were summed to a MS Genetic Burden (MSGB) score. Study replication was performed in a cohort of white American MS patients (n = 1997) and controls (n = 708). RESULTS: The total human leukocyte antigen (HLA) and the non-HLA MSGB scores were significantly higher in MS patients than in controls, in both cohorts (P << 10(-22)). MS patients, with and without cerebrospinal fluid (CSF) oligoclonal bands (OCBs), had a higher MSGB score than the controls; the OCB-positive patients had a slightly higher MSGB than the OCB-negative patients. An early age at symptom onset (AAO) also correlated with a higher MSGB score, in both cohorts. CONCLUSION: The MSGB score was associated with specific clinical MS characteristics, such as OCBs and AAO. This study underlines the need for well-characterized, large cohorts of MS patients, and the usefulness of summarizing multiple genetic risk factors of modest effect size in genotype-phenotype analyses.

Predicting disease severity in multiple sclerosis using multimodal data and machine learning.

BACKGROUND: Multiple sclerosis patients would benefit from machine learning algorithms that integrates clinical, imaging and multimodal biomarkers to define the risk of disease activity. METHODS: We have analysed a prospective multi-centric cohort of 322 MS patients and 98 healthy controls from four MS centres, collecting disability scales at baseline and 2 years later. Imaging data included brain MRI and optical coherence tomography, and omics included genotyping, cytomics and phosphoproteomic data from peripheral blood mononuclear cells. Predictors of clinical outcomes were searched using Random Forest algorithms. Assessment of the algorithm performance was conducted in an independent prospective cohort of 271 MS patients from a single centre. RESULTS: We found algorithms for predicting confirmed disability accumulation for the different scales, no evidence of disease activity (NEDA), onset of immunotherapy and the escalation from low- to high-efficacy therapy with intermediate to high-accuracy. This accuracy was achieved for most of the predictors using clinical data alone or in combination with imaging data. Still, in some cases, the addition of omics data slightly increased algorithm performance. Accuracies were comparable in both cohorts. CONCLUSION: Combining clinical, imaging and omics data with machine learning helps identify MS patients at risk of disability worsening.

Identification of new susceptibility loci for osteoarthritis (arcOGEN): a genome-wide association study.

BACKGROUND: Osteoarthritis is the most common form of arthritis worldwide and is a major cause of pain and disability in elderly people. The health economic burden of osteoarthritis is increasing commensurate with obesity prevalence and longevity. Osteoarthritis has a strong genetic component but the success of previous genetic studies has been restricted due to insufficient sample sizes and phenotype heterogeneity. METHODS: We undertook a large genome-wide association study (GWAS) in 7410 unrelated and retrospectively and prospectively selected patients with severe osteoarthritis in the arcOGEN study, 80% of whom had undergone total joint replacement, and 11,009 unrelated controls from the UK. We replicated the most promising signals in an independent set of up to 7473 cases and 42,938 controls, from studies in Iceland, Estonia, the Netherlands, and the UK. All patients and controls were of European descent. FINDINGS: We identified five genome-wide significant loci (binomial test p≤5·0×10(-8)) for association with osteoarthritis and three loci just below this threshold. The strongest association was on chromosome 3 with rs6976 (odds ratio 1·12 [95% CI 1·08-1·16]; p=7·24×10(-11)), which is in perfect linkage disequilibrium with rs11177. This SNP encodes a missense polymorphism within the nucleostemin-encoding gene GNL3. Levels of nucleostemin were raised in chondrocytes from patients with osteoarthritis in functional studies. Other significant loci were on chromosome 9 close to ASTN2, chromosome 6 between FILIP1 and SENP6, chromosome 12 close to KLHDC5 and PTHLH, and in another region of chromosome 12 close to CHST11. One of the signals close to genome-wide significance was within the FTO gene, which is involved in regulation of bodyweight-a strong risk factor for osteoarthritis. All risk variants were common in frequency and exerted small effects. INTERPRETATION: Our findings provide insight into the genetics of arthritis and identify new pathways that might be amenable to future therapeutic intervention. FUNDING: arcOGEN was funded by a special purpose grant from Arthritis Research UK.

Increased type II deiodinase protein in OA-affected cartilage and allelic imbalance of OA risk polymorphism rs225014 at DIO2 in human OA joint tissues.

OBJECTIVE: Genetic variation at the type II deiodinase (D2) gene (DIO2) was previously identified as osteoarthritis (OA) risk factor. To investigate mechanisms possibly underlying this association, we assessed D2 protein in healthy and OA-affected cartilage and investigated allelic balance of the OA risk polymorphism rs225014 at DIO2 in human OA joints. METHODS: Immunohistochemical staining of healthy and OA-affected cartilage was performed for D2. We then assessed allelic balance of DIO2 mRNA within OA-affected cartilage both at and away from the lesion, ligaments and subchondral bone. Allelic balance was measured by the amount of alleles 'C' and 'T' of the intragenic OA risk polymorphism rs225014 in heterozygous carriers. RESULTS: A markedly higher amount of D2 positive cells and staining intensity was observed in OA cartilage. A significant, 1.3-fold higher presence was observed for the OA-associated rs225014 'C' allele relative to the 'T' allele of DIO2, which was significant in 28 of 31 donors. CONCLUSION: In OA cartilage, D2 protein presence is increased. The allelic imbalance of the DIO2 mRNA transcript, with the OA risk allele 'C' of rs225014 more abundant than the wild-type 'T' allele in heterozygote carriers provides a possible mechanism by which genetic variation at DIO2 confers OA risk.

Brain disconnectome mapping derived from white matter lesions and serum neurofilament light levels in multiple sclerosis: A longitudinal multicenter study.

BACKGROUND AND OBJECTIVES: Connectivity-based approaches incorporating the distribution and magnitude of the extended brain network aberrations caused by lesions may offer higher sensitivity for axonal damage in patients with multiple sclerosis (MS) than conventional lesion characteristics. Using individual brain disconnectome mapping, we tested the longitudinal associations between putative imaging-based brain network aberrations and levels of serum neurofilament light chain (NfL) as a neuroaxonal injury biomarker. METHODS: MS patients (n = 312, mean age 42.9 years, 71 % female) and healthy controls (HC) (n = 59, mean age 39.9 years, 78 % female) were prospectively enrolled at four European MS centres, and reassessed after two years (MS, n = 242; HC, n = 30). Post-processing of 3 Tesla (3 T) MRI data was performed at one centre using a harmonized pipeline, and disconnectome maps were calculated using BCBtoolkit based on individual lesion maps. Global disconnectivity (GD) was defined as the average disconnectome probability in each patient's white matter. Serum NfL concentrations were measured by single molecule array (Simoa). Robust linear mixed models (rLMM) with GD or T2-lesion volume (T2LV) as dependent variables, patient as a random factor, serum NfL, age, sex, timepoint for visit, diagnosis, treatment, and center as fixed factors were run. RESULTS: rLMM revealed significant associations between GD and serum NfL (t = 2.94, p = 0.003), age (t = 4.21, p = 2.5 × 10-5), and longitudinal changes in NfL (t = -2.29, p = 0.02), but not for sex (t = 0.63, p = 0.53) or treatments (t = 0.80-0.83, p = 0.41-0.42). Voxel-wise analyses revealed significant associations between dysconnectivity in cerebellar and brainstem regions and serum NfL (t = 7.03, p < 0.001). DISCUSSION: In our prospective multi-site MS cohort, rLMMs demonstrated that the extent of global and regional brain disconnectivity is sensitive to a systemic biomarker of axonal damage, serum NfL, in patients with MS. These findings provide a neuroaxonal correlate of advanced disconnectome mapping and provide a platform for further investigations of the functional and potential clinical relevance of brain disconnectome mapping in patients with brain disorders.

Meta-analyses of genes modulating intracellular T3 bio-availability reveal a possible role for the DIO3 gene in osteoarthritis susceptibility.

OBJECTIVE: To study whether common genetic variants of the genes involved in the complex regulatory mechanism determining the intracellular bio-availability of T3 influence osteoarthritis onset. METHODS: In total 17 genetic variants within the genes encoding WD40-repeat/SOCS-box protein 1, ubiquitin specific protease 33, thyroid hormone receptor α, deiodinase, iodothyronine, type III (DIO3) and Indian hedgehog were measured and associated with osteoarthritis in a meta-analyses in European populations from the UK, The Netherlands, Greece and Spain containing a total of 3252 osteoarthritis cases and 2132 controls. RESULTS: The minor allele of the DIO3 variant rs945006 showed suggestive evidence for protective association in the overall meta-analyses, which was supported by individual osteoarthritis studies and osteoarthritis subtypes. The association appeared most significant in cases with knee and/or hip with an allelic OR of 0.81 (95% CI 0.70 to 0.930) with a nominal p value of 0.004 and a permutation-based corrected p value for multiple testing of 0.039. CONCLUSION: The findings suggest that the DIO3 gene modulates osteoarthritis disease risk; however, additional studies are necessary to replicate our findings. To elucidate the molecular mechanisms focus should be on the local adaptation to T3 availability either during the endochondral ossification process or during ageing of the articular cartilage.

Meta-analysis of genome-wide association studies confirms a susceptibility locus for knee osteoarthritis on chromosome 7q22.

OBJECTIVES: Osteoarthritis (OA) is the most prevalent form of arthritis and accounts for substantial morbidity and disability, particularly in older people. It is characterised by changes in joint structure, including degeneration of the articular cartilage, and its aetiology is multifactorial with a strong postulated genetic component. METHODS: A meta-analysis was performed of four genome-wide association (GWA) studies of 2371 cases of knee OA and 35 909 controls in Caucasian populations. Replication of the top hits was attempted with data from 10 additional replication datasets. RESULTS: With a cumulative sample size of 6709 cases and 44 439 controls, one genome-wide significant locus was identified on chromosome 7q22 for knee OA (rs4730250, p=9.2 × 10⁻9), thereby confirming its role as a susceptibility locus for OA. CONCLUSION: The associated signal is located within a large (500 kb) linkage disequilibrium block that contains six genes: PRKAR2B (protein kinase, cAMP-dependent, regulatory, type II, ß), HPB1 (HMG-box transcription factor 1), COG5 (component of oligomeric golgi complex 5), GPR22 (G protein-coupled receptor 22), DUS4L (dihydrouridine synthase 4-like) and BCAP29 (B cell receptor-associated protein 29). Gene expression analyses of the (six) genes in primary cells derived from different joint tissues confirmed expression of all the genes in the joint environment.

Interleukin-1 gene cluster variants with innate cytokine production profiles and osteoarthritis in subjects from the Genetics, Osteoarthritis and Progression Study.

OBJECTIVE: To assess whether genetic variation in the interleukin-1 (IL-1) gene cluster contributes to familial osteoarthritis (OA) by influencing innate ex vivo production of IL-1beta or IL-1 receptor antagonist (IL-1Ra). METHODS: Innate ex vivo IL-1beta and IL-1Ra production upon lipopolysaccharide (LPS) stimulation of whole blood cells was measured in subjects from the Genetics, Osteoarthritis and Progression (GARP) Study, which includes sibling pairs in which at least one sibling has symptomatic OA at multiple sites. Radiographic OA (ROA) was assessed by Kellgren/Lawrence score. Subjects from the GARP Study and controls from the Rotterdam Study were genotyped for 7 single-nucleotide polymorphisms (SNPs) encompassing the IL-1 gene cluster on chromosome 2q13. Linkage disequilibrium analysis and genotype and haplotype association analysis were performed to assess the relationship between the IL-1 gene cluster SNPs, innate ex vivo cytokine production, and OA. RESULTS: Among subjects in the GARP Study, the haplotype variable-number tandem repeat in intron 2/T+8006C/T+11100C 2/2/1 of the IL1RN gene was significantly associated with reduced innate ex vivo bioavailability of IL-1beta upon LPS stimulation (P = 0.026) and with ROA at the highest number of joint locations. CONCLUSION: These results show that genetic variation at the IL-1 gene cluster is associated with lower IL-1beta bioavailability and with OA at a large number of joint locations. The data further indicate that, among subjects with OA affecting the highest number of joints, the innate immune system may be activated, thereby obscuring possible underlying mechanisms.