Prognostic value of single-subject grey matter networks in early multiple sclerosis

The identification of prognostic markers in early multiple sclerosis (MS) is challenging and requires reliable measures that robustly predict future disease trajectories. Ideally, such measures should make inferences at the individual level to inform clinical decisions. This study investigated the prognostic value of longitudinal structural networks to predict five-year EDSS progression in patients with relapsing-remitting MS (RRMS). We hypothesized that network measures, derived from magnetic resonance imaging (MRI), outperform conventional MRI measurements at identifying patients at risk of developing disability progression. This longitudinal, multicentre study within the Magnetic Resonance Imaging in MS (MAGNIMS) network included 406 patients with RRMS (mean age = 35.7 ± 9.1 years) followed up for five years (mean follow-up = 5.0 ± 0.6 years). Expanded Disability Status Scale (EDSS) was determined to track disability accumulation. A group of 153 healthy subjects (mean age = 35.0 ± 10.1 years) with longitudinal MRI served as controls. All subjects underwent MRI at baseline and again one year after baseline. Grey matter (GM) atrophy over one year and white matter (WM) lesion load were determined. A single-subject brain network was reconstructed from T1-weighted scans based on GM atrophy measures derived from a statistical parameter mapping (SPM)-based segmentation pipeline. Key topological measures, including network degree, global efficiency and transitivity, were calculated at single-subject level to quantify network properties related to EDSS progression. Areas under receiver operator characteristic (ROC) curves were constructed for GM atrophy, WM lesion load and the network measures, and comparisons between ROC curves were conducted. The applied network analyses differentiated patients with RRMS who experience EDSS progression over five years through lower values for network degree [H(2)=30.0, p<0.001] and global efficiency [H(2)=31.3, p<0.001] from healthy controls but also from patients without progression. For transitivity, the comparisons showed no difference between the groups (H(2)= 1.5, p=0.474). Most notably, changes in network degree and global efficiency were detected independent of disease activity in the first year. The described network reorganization in patients experiencing EDSS progression was evident in the absence of GM atrophy. Network degree and global efficiency measurements demonstrated superiority of network measures in the ROC analyses over GM atrophy and WM lesion load in predicting EDSS worsening (all p-values < 0.05). Our findings provide evidence that GM network reorganization over one year discloses relevant information about subsequent clinical worsening in RRMS. Early GM restructuring towards lower network efficiency predicts disability accumulation and outperforms conventional MRI predictors

Cross-Sectional and Longitudinal MRI Brain Scans Reveal Accelerated Brain Aging in Multiple Sclerosis

Multiple sclerosis (MS) is an inflammatory disorder of the central nervous system. By combining longitudinal MRI-based brain morphometry and brain age estimation using machine learning, we tested the hypothesis that MS patients have higher brain age relative to chronological age than healthy controls (HC) and that longitudinal rate of brain aging in MS patients is associated with clinical course and severity. Seventy-six MS patients [71% females, mean age 34.8 years (range 21–49) at inclusion] were examined with brain MRI at three time points with a mean total follow up period of 4.4 years (±0.4 years). We used additional cross-sectional MRI data from 235 HC for case-control comparison. We applied a machine learning model trained on an independent set of 3,208 HC to estimate individual brain age and to calculate the difference between estimated and chronological age, termed brain age gap (BAG). We also assessed the longitudinal change rate in BAG in individuals with MS. MS patients showed significantly higher BAG (4.4 ± 6.6 years) compared to HC (Cohen's D = 0.69, p = 4.0 × 10−6). Longitudinal estimates of BAG in MS patients showed high reliability and suggested an accelerated rate of brain aging corresponding to an annual increase of 0.41 (SE = 0.15) years compared to chronological aging (p = 0.008). Multiple regression analyses revealed higher rate of brain aging in patients with more brain atrophy (Cohen's D = 0.86, p = 4.3 × 10−15) and increased white matter lesion load (WMLL) (Cohen's D = 0.55, p = 0.015). On average, patients with MS had significantly higher BAG compared to HC. Progressive brain aging in patients with MS was related to brain atrophy and increased WMLL. No significant clinical associations were found in our sample, future studies are warranted on this matter. Brain age estimation is a promising method for evaluation of subtle brain changes in MS, which is important for predicting clinical outcome and guide choice of intervention

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

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

Pathways validation for canonical pathway.

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.</div

Sys4MS cohort: Clinical and imaging variables of MS patients and healthy controls.

Disability scales are shown as the mean ± SD, except for the EDSS which is displayed as the median (range).</p

Venn diagram describing the overlap between the paths identified in the single-cell analysis and the paths identified in the UNIPROT database.

(1) CD56 Neg > INL—mRNFL > EDSS—T25WT; (2) Total CD8 > NGMV—T2LV > EDSS - 9HPT–SDMT; (3) MK03 > Total T Cells > mRNFL > T25WT; (4) HSPB1 > B Memory > NBV > T25WT; (5) STAT6 > Th17 > NGMV Change > Years with Disease; (6) KS6B1—LCK > Total T Cells—Th1 Non Classic > NGMV—T2LV> LCVA Change—MSSS—Years since Relapse; (7) MP2K1—STAT6 > Th17 > mRNFL > T25WT—ARMSS (8) MP2K1—STAT6 > Th17 > INL > EDSS Change; (9) MP2K1 > CD8 Treg > GCIPL > EDSS Change; (10) Atypical B Memory–B Memory–Th1 Classic > mRNFL–T2LV > EDSS–T25WT.</p

Linear regression for HSPB1 node.

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.</div

Difference in Pearson correlation between healthy and infected cases.

The networks shown contain paths that were identified from the Boolean simulations in the infected network. Furthermore, each path contains at least two nodes from two different layers that are present in the acute phase response signaling biological pathway. The same paths do not necessarily appear in the healthy network, so edges with Pearson correlation are shown. There is a notable increase in the strength of the connections, both positive and negative, in the infected case.</p