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.

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 5-year Expanded Disability Status Scale (EDSS) progression in patients with relapsing-remitting MS (RRMS). We hypothesized that network measures, derived from 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 5 years (mean follow-up = 5.0 ± 0.6 years). 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 1 year after baseline. Grey matter atrophy over 1 year and white matter lesion load were determined. A single-subject brain network was reconstructed from T1-weighted scans based on grey matter atrophy measures derived from a statistical parameter mapping-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 grey matter atrophy and white matter 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 5 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 grey matter atrophy. Network degree and global efficiency measurements demonstrated superiority of network measures in the ROC analyses over grey matter atrophy and white matter lesion load in predicting EDSS worsening (all P-values < 0.05). Our findings provide evidence that grey matter network reorganization over 1 year discloses relevant information about subsequent clinical worsening in RRMS. Early grey matter restructuring towards lower network efficiency predicts disability accumulation and outperforms conventional MRI predictors.

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.

Immune cell subpopulations and serum neurofilament light chain are associated with increased risk of disease worsening in multiple sclerosis.

Changes is lymphocyte subpopulations in peripheral blood have been proposed as biomarkers for evaluation of disease activity in multiple sclerosis (MS). Serum neurofilament light chain (sNfL) is a biomarker reflecting neuro-axonal injury in MS that could be used to monitor disease activity, response to drugs and to prognosticate disease course. Here we show a moderate correlation between sNfL and lymphocyte cell subpopulations, and our data furthermore suggest that sNfL and specific immune cell subpopulations together could predict future disease worsening in MS.

Using The Virtual Brain to study the relationship between structural and functional connectivity in patients with multiple sclerosis: a multicenter study.

The relationship between structural connectivity (SC) and functional connectivity (FC) captured from magnetic resonance imaging, as well as its interaction with disability and cognitive impairment, is not well understood in people with multiple sclerosis (pwMS). The Virtual Brain (TVB) is an open-source brain simulator for creating personalized brain models using SC and FC. The aim of this study was to explore SC-FC relationship in MS using TVB. Two different model regimes have been studied: stable and oscillatory, with the latter including conduction delays in the brain. The models were applied to 513 pwMS and 208 healthy controls (HC) from 7 different centers. Models were analyzed using structural damage, global diffusion properties, clinical disability, cognitive scores, and graph-derived metrics from both simulated and empirical FC. For the stable model, higher SC-FC coupling was associated with pwMS with low Single Digit Modalities Test (SDMT) score (F=3.48, P$\lt$0.05), suggesting that cognitive impairment in pwMS is associated with a higher SC-FC coupling. Differences in entropy of the simulated FC between HC, high and low SDMT groups (F=31.57, P$\lt$1e-5), show that the model captures subtle differences not detected in the empirical FC, suggesting the existence of compensatory and maladaptive mechanisms between SC and FC in MS.

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.

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 Retinal Blood Vessel Diameters as Biomarkers in Multiple Sclerosis.

We aimed to determine whether retinal vessel diameters and retinal oxygen saturation in newly diagnosed patients with multiple sclerosis (pwMS) are different from those of a healthy population. Retinal blood vessel diameters were measured using imaging with a spectrophotometric non-invasive retinal oximeter. Twenty-three newly diagnosed untreated relapsing-remitting MS (RRMS) patients (mean age: 32.2 ± 7.5 years, age range = 18-50 years, 56.5% female) were measured and compared to 23 age- and sex-matched healthy controls (HCs) (mean age: 34.8 ± 8.1 years). Patients with Optic Neuritis were excluded. Retinal venular diameter (143.8 µm versus 157.8 µm: mean; p = 0.0013) and retinal arteriolar diameter (112.6 µm versus 120.6 µm: mean; p = 0.0089) were smaller in pwMS when compared with HCs, respectively. There was no significant difference in the oxygen saturation in retinal venules and arterioles in pwMS (mean: 60.0% and 93.7%; p = 0.5980) compared to HCs (mean: 59.3% and 91.5%; p = 0.8934), respectively. There was a significant difference in the median low contrast visual acuity (2.5% contrast) between the pwMS and the HC groups (p = 0.0143) Retinal arteriolar and venular diameter may have potential as objective biomarkers for MS.

Deep neural networks learn general and clinically relevant representations of the ageing brain.

The discrepancy between chronological age and the apparent age of the brain based on neuroimaging data - the brain age delta - has emerged as a reliable marker of brain health. With an increasing wealth of data, approaches to tackle heterogeneity in data acquisition are vital. To this end, we compiled raw structural magnetic resonance images into one of the largest and most diverse datasets assembled (n=53542), and trained convolutional neural networks (CNNs) to predict age. We achieved state-of-the-art performance on unseen data from unknown scanners (n=2553), and showed that higher brain age delta is associated with diabetes, alcohol intake and smoking. Using transfer learning, the intermediate representations learned by our model complemented and partly outperformed brain age delta in predicting common brain disorders. Our work shows we can achieve generalizable and biologically plausible brain age predictions using CNNs trained on heterogeneous datasets, and transfer them to clinical use cases.

Cardiometabolic risk factors associated with brain age and accelerate brain ageing.

The structure and integrity of the ageing brain is interchangeably linked to physical health, and cardiometabolic risk factors (CMRs) are associated with dementia and other brain disorders. In this mixed cross-sectional and longitudinal study (interval mean = 19.7 months), including 790 healthy individuals (mean age = 46.7 years, 53% women), we investigated CMRs and health indicators including anthropometric measures, lifestyle factors, and blood biomarkers in relation to brain structure using MRI-based morphometry and diffusion tensor imaging (DTI). We performed tissue specific brain age prediction using machine learning and performed Bayesian multilevel modeling to assess changes in each CMR over time, their respective association with brain age gap (BAG), and their interaction effects with time and age on the tissue-specific BAGs. The results showed credible associations between DTI-based BAG and blood levels of phosphate and mean cell volume (MCV), and between T1-based BAG and systolic blood pressure, smoking, pulse, and C-reactive protein (CRP), indicating older-appearing brains in people with higher cardiometabolic risk (smoking, higher blood pressure and pulse, low-grade inflammation). Longitudinal evidence supported interactions between both BAGs and waist-to-hip ratio (WHR), and between DTI-based BAG and systolic blood pressure and smoking, indicating accelerated ageing in people with higher cardiometabolic risk (smoking, higher blood pressure, and WHR). The results demonstrate that cardiometabolic risk factors are associated with brain ageing. While randomized controlled trials are needed to establish causality, our results indicate that public health initiatives and treatment strategies targeting modifiable cardiometabolic risk factors may also improve risk trajectories and delay brain ageing.

LesionQuant for Assessment of MRI in Multiple Sclerosis-A Promising Supplement to the Visual Scan Inspection.

Background and Goals: Multiple sclerosis (MS) is a central nervous system inflammatory disease where magnetic resonance imaging (MRI) is an important tool for diagnosis and disease monitoring. Quantitative measurements of lesion volume, lesion count, distribution of lesions, and brain atrophy have a potentially significant value for evaluating disease progression. We hypothesize that utilizing software designed for evaluating MRI data in MS will provide more accurate and detailed analyses compared to the visual neuro-radiological evaluation. Methods: A group of 56 MS patients (mean age 35 years, 70% females and 96% relapsing-remitting MS) was examined with brain MRI one and 5 years after diagnosis. The T1 and FLAIR brain MRI sequences for all patients were analyzed using the LesionQuant (LQ) software. These data were compared with data from structured visual evaluations of the MRI scans performed by neuro-radiologists, including assessments of atrophy, and lesion count. The data from LQ were also compared with data from other validated research methods for brain segmentation, including assessments of whole brain volume and lesion volume. Correlations with clinical tests like the timed 25-foot walk test (T25FT) were performed to explore additional value of LQ analyses. Results: Lesion count assessments by LQ and by the neuro-radiologist were significantly correlated one year (cor = 0.92, p = 2.2 × 10-16) and 5 years (cor = 0.84, p = 2.7 × 10-16) after diagnosis. Analyzes of the intra- and interrater variability also correlated significantly (cor = 0.96, p < 0.001, cor = 0.97, p < 0.001). Significant positive correlation was found between lesion volume measured by LQ and by the software Cascade (cor = 0.7, p < 0.001. LQ detected a reduction in whole brain percentile >10 in 10 patients across the time-points, whereas the neuro-radiologist assessment identified six of these. The neuro-radiologist additionally identified five patients with increased atrophy in the follow-up period, all of them displayed decreasing low whole brain percentiles (median 11, range 8-28) in the LQ analysis. Significant positive correlation was identified between lesion volume measured by LQ and test performance on the T25FT both at 1 and 5 years after diagnosis. Conclusion: For the number of MS lesions at both time-points, we demonstrated strong correlations between the assessments done by LQ and the neuro-radiologist. Lesion volume evaluated with LQ correlated with T25FT performance. LQ-analyses classified more patients to have brain atrophy than the visual neuro-radiological evaluation. In conclusion, LQ seems like a promising supplement to the evaluation performed by neuro-radiologists, providing an automated tool for evaluating lesions in MS patients and also detecting early signs of atrophy in both a longitudinal and cross-sectional setting.

The genetic architecture of human brainstem structures and their involvement in common brain disorders.

Brainstem regions support vital bodily functions, yet their genetic architectures and involvement in common brain disorders remain understudied. Here, using imaging-genetics data from a discovery sample of 27,034 individuals, we identify 45 brainstem-associated genetic loci, including the first linked to midbrain, pons, and medulla oblongata volumes, and map them to 305 genes. In a replication sample of 7432 participants most of the loci show the same effect direction and are significant at a nominal threshold. We detect genetic overlap between brainstem volumes and eight psychiatric and neurological disorders. In additional clinical data from 5062 individuals with common brain disorders and 11,257 healthy controls, we observe differential volume alterations in schizophrenia, bipolar disorder, multiple sclerosis, mild cognitive impairment, dementia, and Parkinson's disease, supporting the relevance of brainstem regions and their genetic architectures in common brain disorders.

Publisher Correction: Common brain disorders are associated with heritable patterns of apparent aging of the brain.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Common brain disorders are associated with heritable patterns of apparent aging of the brain.

Common risk factors for psychiatric and other brain disorders are likely to converge on biological pathways influencing the development and maintenance of brain structure and function across life. Using structural MRI data from 45,615 individuals aged 3-96 years, we demonstrate distinct patterns of apparent brain aging in several brain disorders and reveal genetic pleiotropy between apparent brain aging in healthy individuals and common brain disorders.

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.