Skip to main content

Combination protein biomarkers predict multiple sclerosis diagnosis and outcomes

Abstract

Establishing biomarkers to predict multiple sclerosis diagnosis and prognosis has been challenging using a single biomarker approach. We hypothesised that a combination of biomarkers would increase the accuracy of prediction models to differentiate multiple sclerosis from other neurological disorders and enhance prognostication for people with multiple sclerosis. We measured 24 fluid biomarkers in the blood and cerebrospinal fluid of 77 people with multiple sclerosis and 80 people with other neurological disorders, using ELISA or Single Molecule Array assays. Primary outcomes were multiple sclerosis versus any other diagnosis, time to first relapse, and time to disability milestone (Expanded Disability Status Scale 6), adjusted for age and sex. Multivariate prediction models were calculated using the area under the curve value for diagnostic prediction, and concordance statistics (the percentage of each pair of events that are correctly ordered in time for each of the Cox regression models) for prognostic predictions. Predictions using combinations of biomarkers were considerably better than single biomarker predictions. The combination of cerebrospinal fluid [chitinase-3-like-1 + TNF-receptor-1 + CD27] and serum [osteopontin + MCP-1] had an area under the curve of 0.97 for diagnosis of multiple sclerosis, compared to the best discriminative single marker in blood (osteopontin: area under the curve 0.84) and in cerebrospinal fluid (chitinase-3-like-1 area under the curve 0.84). Prediction for time to next relapse was optimal with a combination of cerebrospinal fluid[vitamin D binding protein + Factor I + C1inhibitor] + serum[Factor B + Interleukin-4 + C1inhibitor] (concordance 0.80), and time to Expanded Disability Status Scale 6 with cerebrospinal fluid [C9 + Neurofilament-light] + serum[chitinase-3-like-1 + CCL27 + vitamin D binding protein + C1inhibitor] (concordance 0.98). A combination of fluid biomarkers has a higher accuracy to differentiate multiple sclerosis from other neurological disorders and significantly improved the prediction of the development of sustained disability in multiple sclerosis. Serum models rivalled those of cerebrospinal fluid, holding promise for a non-invasive approach. The utility of our biomarker models can only be established by robust validation in different and varied cohorts.

Introduction

Multiple sclerosis is one of the most common causes of neurological disability in young adults and long-term disease outcomes are highly variable and difficult to predict at onset. Since the emergence of disease modifying therapies for multiple sclerosis in the 1990s, more than 15 treatments with variable efficacy are now available. However, uncertainty remains over selection and sequencing of disease modifying therapies; aggressive suppression of early multiple sclerosis activity may improve longer-term disease outcomes, but most high-efficacy disease modifying therapies carry risks of serious adverse events [1, 2]. Therefore, there is an urgent need to identify multiple sclerosis biomarkers than can inform disease prognosis, treatment response and potential for adverse events, to enable a more personalised approach to therapeutic interventions.

The hallmark of multiple sclerosis pathology is inflammatory demyelination within the white and grey matter of the brain and spinal cord. Foci of acute inflammation and demyelination are thought to explain the short-lived flare-ups of symptoms (relapses), seen early in multiple sclerosis, but a more complex interplay between immune cells, glia and neurons is likely to explain accumulation of longer-term and/or progressive disability [3, 4]. Recent studies supporting the use of the neuronal protein neurofilament light as a prognostic multiple sclerosis biomarker have offered some optimism for personalising treatment decisions. The prognostic value of neurofilament light in multiple sclerosis is well-established when measured in CSF using ELISA-based assays [5], but this represents an impractical tissue to sample repeatedly. Serum neurofilament light concentrations are an order of magnitude lower and poorly detected using ELISA but can now be accurately measured using single molecule array (SiMoA) assays [6]. Good correlation is observed between serum and CSF concentrations of neurofilament light, a feature which has also been demonstrated for some other candidate multiple sclerosis biomarkers. As a result serum/plasma analysis of selected biomarkers suggests promise for a less invasive approach to disease monitoring [7, 8]. However, group effects observed in neurofilament light studies appear to translate to relatively modest individual predictions, especially when using a single cross-sectional measurement, raising questions over its clinical utility [9]. Other proteins including glial marker chitinase-3-like-1 [10], and inflammatory cytokine osteopontin [11], have also been identified as potential candidates for predicting multiple sclerosis outcomes, and may offer complementary information to neurofilament light. Candidate multiple sclerosis biomarkers are often studied in isolation, and comparing studies that have used different clinical outcomes is challenging [12]. Given the complexity of multiple sclerosis biology, a combination of markers may improve predictions versus any marker in isolation.

In this study, we investigate 31 candidate biomarkers for their ability to distinguish people with multiple sclerosis from people with other neurological conditions, and predict short- and medium-term clinical outcomes in those with multiple sclerosis. We use paired CSF and serum samples from 157 people to study the optimal combination of markers for diagnostic and prognostic predictions.

Materials and methods

We identified candidate biomarkers during a literature search (conducted in 2018, using PubMed database and search terms “multiple sclerosis” and “biomarker” ± “CSF” or “serum” or “plasma”). We identified candidate tissue biomarkers which had been shown to differentiate between multiple sclerosis and controls, between multiple sclerosis subtypes, or to predict severity of the disease (Table 1) [12]. Candidate biomarkers included pro- and anti-inflammatory cytokines (e.g. IL-6, IL-18, and IL-10), proteins and factors involved in the CNS (such as brain-derived neurotrophic factor (BDNF), osteopontin and glial fibrillary acidic protein (GFAP)), chemokines (such as CXCL12 and CXCL13), complement proteins and proteins that are involved in other pathways relevant to multiple sclerosis (such as Vitamin D binding protein).

Table 1 Description of candidate biomarkers studied in this cohort

Participants

People with multiple sclerosis were recruited as part of a long running observational study in South East Wales, UK (Ethics REC Ref: 05/ WSE03/111, 19/WA/0289), which has been described previously [13]. Briefly, participants are invited to annual clinic visits and data are collected at each clinical encounter, including measurement of expanded disability status scale (EDSS) [14], relapse history, disease subtype and disease modifying therapy review. People with new neurological symptoms lasting > 24 h in the absence of infection are reviewed by clinicians to confirm relapse (in the absence of an alternative explanation). Annual clinic data are complemented by a postal questionnaire, which includes an assessment of current disability using a validated self-reported EDSS tool [15]. Where possible, blood and CSF samples are acquired at the time of diagnostic lumbar puncture and archived in a biorepository. We studied paired serum and CSF samples from 77 people who were being investigated for demyelinating disease, all of whom ultimately reached a diagnosis of multiple sclerosis according to contemporary criteria, and 80 controls who were undergoing CSF examination for reasons other than suspected demyelinating disease; none of whom were diagnosed with multiple sclerosis during follow-up. Control samples were obtained from the Welsh Neuroscience Research Tissue Bank (WNRTB Ethics Rec Ref: 19/WA/0058). The demographic and clinical characteristics of the cohort are summarised in Table 2.

Table 2 Demographic and clinical characteristics of the cohort

Samples

Venous blood samples were collected in a BD vacutainer (Gold SST II Advance) between 04/2007 and 01/2020. Whole blood and paired CSF samples were immediately centrifuged at 4500 rpm for 10min at 4 °C within 2–3 h of collection. The resulting supernatant (serum or CSF) was subsequently split into 300 μl aliquots and stored at – 80 °C until use.

Assays

In-house ELISA assays were used for the measurement of complement biomarkers (terminal complement complex (TCC), iC3b, C3, C5, C9, factor H, factor B, factor I) in serum and CSF. C1inhibitor/C1s complex (C1inh/C1s) was measured using a commercially available assay donated by Hycult Biotech (Cat No HK399). For all other biomarkers commercial assays from R&D were used (Additional file 1: Table S1). Dilutions of serum and CSF for all assays are described in Additional file 1: Table S1. For the complement assays, Nunc Maxisorp (VWR, UK) plates were coated with affinity-purified capture antibody overnight at 4 °C. Wells were subsequently washed 1 × with phosphate-buffered saline 0.1% Tween20 (Sigma Aldrich, Germany) (PBS-T) and blocked for 1 h at 37 °C with 1% bovine serum albumin (BSA) in PBS-T. After 1 wash, purified protein standards or serum samples (in duplicate), optimally diluted in 0.1% BSA in PBS-T, 5mM ethylenediaminetetraacetic acid (EDTA), were added to wells in duplicate and incubated for 1.5 h at 37 °C. Wells were washed 3 × with PBS-T then incubated for 1 h at 37 °C with detection antibody (either unlabelled or labelled in-house using horseradish peroxidase (HRP, where EZ-Link Plus Activated Peroxidase Kit (ThermoFisher Scientific, #31489, UK)) or biotin for the TCC assay (#21327) and washed 3 times with PBS-T. For unlabelled detection antibodies, HRP-labelled secondary antibody (anti-mouse or anti-rabbit IgG as appropriate (Jackson ImmunoResearch #715-035-151, #711-035-152, USA), or streptavidin-HRP (R&D systems, #DY995, UK) for TCC detection was added to wells, incubated and washed as above. Signals were detected using 3,3',5,5'-tetramethylbenzidine (TMB) (1-Step™ Ultra TMB-ELISA Substrate Solution, Thermofisher Scientific, UK) and after the addition of 2M sulfuric acid, absorbance (450nm) was measured. In each plate, protein standards were included, and samples were randomly assigned to eliminate assay bias. A nonlinear regression model was used to fit standard curves generated by ELISA. Biomarkers levels were automatically calculated by reference to the standard curve using GraphPad Prism. To account for any difference in the degree of blood–brain barrier permeability between people with multiple sclerosis and controls, the CSF total protein concentration was measured using micro-BCA Protein Assay Kit #23235, ThermoFisher Scientific, UK) [16], and CSF biomarker concentrations were adjusted accordingly. For the non-complement assays (Table 2 for details), manufacturer’s instructions were followed.

During optimisation, it was found that analytes were undetectable in the majority of CSF and serum samples using 5 assays (IFNγ, TGFβ, TNFα, IL-6, IL-12), so these were not taken forward into the main analysis. Analytes from six assays (GFAP, BDNF, IL-10, IL-18, LIF, CCL27, CXCL13) were undetectable in the CSF, and were thereafter used only in serum samples. One analyte (IL-8) was undetectable in serum, and thereafter only used in CSF. For all ELISA assays, intra-assay CVs were less than 10%, and inter-assay CVs were under 15%. For values outside the standard curve, the samples were re-measured using different dilutions. For neat samples, values below the standard curve were imputed (see Statistical analysis section and Additional file 1: Table S2).

Single molecule array (SiMoA) technology was employed to quantify neurofilament light in serum samples in view of the known insensitivity of ELISA for serum neurofilament light [6]. We used an HD-X analyser, a fully automated multiplex digital immunoassay instrument providing ultra-sensitive measurements (fg/ml) of proteins of biological interest, over a wide dynamic range and with low coefficients of variance. Serum samples were quantified using the human neurofilament light single bead-based advantage assay (Quanterix cat#103186), array discs (cat#103347) and assay buffers (cat#100488), following the manufacturers protocol.

Statistical analysis

The demographic and clinical characteristics of the cohort were presented using descriptive statistics. Due to group differences between the multiple sclerosis and control cohorts, all analyses were corrected for age and sex. The concentration of each biomarker for both serum and CSF was individually compared between multiple sclerosis and controls using the non-parametric Mann–Whitney test. The Benjamini–Hochberg correction was applied to the p-values to control for the probability of type 1 error with repeated tests. Concentration of each biomarker (CSF and serum) were individually converted to Z-scores to allow all their model coefficients to be directly compared. Z-scores were generated by subtracting the mean and dividing by the standard deviation for all the valid data for each biomarker. To avoid zero values, levels for samples which were below the minimum detection limit were imputed using a value equal to half of the minimum detectable value (using our assays and samples). Due to scarce CSF or serum for some patients, it was not possible to measure all analytes in all samples. The frequency of missing and imputed data is shown in Additional file 1: Table S2, but did not differ significantly between MS and control cases.

The Z-scores were used to model the utility of single and combination biomarkers of predicting disease status and clinical outcome(s). For determination of MS versus control status, a cross-validation approach was taken to avoid over-fitting of a model. The cohort was randomly divided into four roughly equally sized and distinct groups (Additional file 1: Table S3 and S4). Three of these four groups were separately used as “training datasets” to produce a model using roughly 75% of the data, which was then tested on the remaining approximately 25% of the data (the “test dataset”). This was repeated 4 times with each of the 4 groups being the test dataset with the other 3 producing the model. In investigating the value of biomarkers to distinguish multiple sclerosis from control categories, receiver operating characteristic curve (ROC) curves were produced and the area under the curve (AUC) was used to calculate the ability of single biomarkers (CSF and serum) and combinations of biomarkers, to distinguish multiple sclerosis from control cases. After looping through all combinations of biomarkers, a mean test AUC was calculated from the AUCs from the 4 test cohorts. Mean AUC values were ordered from lowest to highest, and the optimum model was selected when addition of a further analyte resulted in an AUC increase < 0.01 (worked example provided in Additional file 1). CSF and serum biomarkers were explored independently in models adjusted for age and sex, and finally data on both serum and CSF markers were combined in a single adjusted model. To explore the reliability of the results, we re-ran the models 1000 times on random selections of Train and Test data, to generate data on the range of AUCs produced.

We performed a sensitivity analysis to explore the possible association between biomarker concentrations with age and sex. We built a linear model of biomarkers levels adjusted for age and sex as predictors in the control cohort. We used this to generate predicted values for each biomarker and subtracted those from the measured biomarker concentrations in both control and MS cohorts, to adjust for any physiological effect of age and sex on biomarker concentration, so that only the MS effects of the biomarkers would be reflected by the model. Using the ‘corrected’ biomarker concentrations, we repeated the modelling as above.

To investigate the prognostic value of single and combinations of biomarkers, we used Cox regression models. We used two different dependent variables (outcomes) to reflect different aspects of multiple sclerosis biology. Firstly, we studied time from sampling (biomarker measurement) to next relapse for all people with multiple sclerosis. To mitigate the effect of disease modifying therapy on suppressing relapses, we censored follow-up at the time of disease modifying therapy commencement. Secondly, we studied time to EDSS 6 (need for unilateral assistance to walk 100m), sustained for at least 6 months, adjusted according to whether people had ever versus never (1 versus 0) received disease modifying therapy during follow-up. Only EDSS scores measured > 30 days out of relapse were used, to account for temporary fluctuations in disability that are likely to recover. The concordance statistic, i.e. the percentage of each pair of events that are correctly ordered in time for each of the Cox regression models, was used to rank best fit. After looping through all combinations of biomarkers, a mean concordance statistic was calculated from the concordance statistics from the 4 test cohorts. Mean concordance values were ordered from lowest to highest, and the optimum model was selected when addition of a further analyte resulted in a concordance increase < 0.01.

All statistical analysis was performed using R version 4.2.1 (2022-06-23 ucrt), R packages pROC (1.18.0), survival (3.3-1), survminer (0.4.9), and SPSS 27.

Results

CSF and serum biomarkers in people with multiple sclerosis versus controls

Eight markers were found to be significantly differently expressed in CSF of multiple sclerosis versus control cases: chitinase-3-like-1, soluble CD27 (sCD27), neurofilament light, osteopontin, C5, iC3b, C9 and TCC (Additional file 1: Table S4). Those with the highest fold difference in CSF between multiple sclerosis and controls were sCD27 (4.1-fold increase in multiple sclerosis versus controls), neurofilament light (3.4-fold increase in multiple sclerosis) and iC3b (2.6-fold increase in multiple sclerosis). Seven markers were found to be significantly differently expressed in the serum of multiple sclerosis versus controls: osteopontin, Factor B, vitamin D binding protein, C5, iC3b, CRP and neurofilament light (Fig. 1, Additional file 1: Table S5). Neurofilament-light (1.6-fold increase in multiple sclerosis) and osteopontin (1.54-fold increase in multiple sclerosis) were the serum biomarkers with the highest fold difference.

Fig. 1
figure 1

Concentrations of CSF and serum biomarkers in multiple sclerosis versus controls. A CSF and B serum samples, values expressed as Z-scores

Statistical models of CSF and serum analytes differentiate multiple sclerosis from controls

The AUC for predicting multiple sclerosis versus control status using age and sex alone was 0.77. sCD27 and neurofilament light were the single CSF biomarkers that produced highest AUC (adjusted for age and sex) to distinguish multiple sclerosis from controls cases (sCD27: Test AUC 0.85; neurofilament light: Test AUC 0.85; Table 3 and Additional file 1: Table S4). Addition of further CSF biomarkers to the model incrementally improved the AUC up to a combination of five CSF biomarkers: Ch3L1 + Tumour necrosis factor receptor 1 (TNFR1) + sCD27 + C9 + TCC (Train/Test AUC 0.96/0.94; Table 3 and Additional file 1: Table S6).

Table 3 Optimal models to predict clinical multiple sclerosis versus non-multiple sclerosis status

Osteopontin and neurofilament light were the single serum biomarkers with highest AUC to distinguish multiple sclerosis from control cases (osteopontin: Test AUC 0.83; neurofilament light: Test AUC 0.80). Addition of further serum biomarkers to the model incrementally improved AUC up to a combination of four serum biomarkers: osteopontin + iC3b + monocyte chemoattractant protein 1 (MCP1) + CCL27 (Train/Test AUC 0.93/0.91, Table 3 and Additional file 1: Table S6). When combining serum and CSF biomarkers in a single model, the AUC plateaued at a combination of five markers: CSF[chitinase-3-like-1 + TNFR1 + sCD27] + serum[osteopontin + MCP1] (Train/Test AUC 0.97/0.95; Table 3 and Additional file 1: Table S6 and S7).

Figure 2 illustrates example violin plots and ROC curves for the optimum combinations of CSF, serum and mixed CSF + serum for their ability to distinguish MS from non-MS cases. There were no obvious patterns of clinical or demographic characteristics of non-MS cases who were misclassified as controls (data not shown). Re-running modelling for the optimum combinations of biomarkers in CSF, serum and CSF and serum 1000 times in random selections of Train and Test data demonstrated mean AUCs > 0.9 for all models (Additional file 1: Fig. S1).

Fig. 2
figure 2

Logistic regression model output of multiple sclerosis versus non-multiple sclerosis. Illustration of the model expressed according to best combination of A CSF, B serum, and C combined CSF and serum biomarkers. Top left: violin plot illustrating the predicted probability of multiple sclerosis (y-axis) by cohort (x-axis), according to the optimum biomarker model adjusted for sex and age. The width of the violin plots indicates the relative number of probabilities at the probability. Horizontal dashed line illustrates an arbitrary cut-off value that could be applied to distinguish MS from non-MS cases, where red dots above the line indicate false positives and blue dots below the line false negatives. Top right: receiver operator curve illustrating performance of the optimum model. Bottom right: regression table indicating the magnitude and direction of effect of each variable. Note: a cross-validation approach was taken whereby the cohort was randomly divided into four roughly equally sized and distinct groups (see Methods), which were separately used as “training datasets” to produce a model using roughly 75% of the data, which was then tested on the remaining approximately 25% of the data (the “test dataset”). This was repeated 4 times to generate a mean test AUC to determine the optimum combination of biomarkers. The ROC curve, regression table and AUCs shown here were produced using one of the train/test sub-cohorts

Given the difference in age and sex between control and MS cohorts, we performed a sensitivity analysis repeating the modelling of AUC using biomarker values that were corrected for the potential effects of age and sex (measured in the controls). This did not lead to any substantial difference in outcomes (Additional file 1: Table S8).

Cox regression modelling predicts time to next relapse

We used univariate Cox regression, adjusted for age and sex and censored for commencement of disease modifying therapy (since disease-modifying therapy is expected to mask new relapses). Thirty-six (47%) people with multiple sclerosis experienced a relapse during follow-up and prior to disease modifying therapy commencement. The CSF biomarkers with highest concordance for predicting time to next relapse were vitamin D binding protein (concordance 0.65, p = 0.083) and neurofilament light (concordance 0.61, P = 0.42; Additional file 1: Table S9). The addition of further biomarkers in a multivariate model improved concordance up to a combination of six biomarkers: vitamin D binding protein + neurofilament light + CXCL12 + Factor B + sCD27 + MCP1 (concordance 0.72; P = 0.075; Table 4 and Additional file 1: Table S9).

Table 4 Optimal models to predict clinical outcomes

In univariate analysis adjusted for age and sex, and censored for commencement of disease modifying therapy, the serum biomarkers with highest concordance for predicting time to next relapse were C5 (concordance 0.61, P = 0.04) and C1inh/C1s (concordance 0.60, P = 0.05; Additional file 1: Table S9). The addition of further biomarkers to the model improved concordance up to a combination of six serum markers: C1inh/C1s + C3 + Factor H + CXCL12 + sCD27 + vitamin D binding protein (concordance 0.72, P = 0.02, Table 4 and Additional file 1: Table S10). Combining CSF and serum markers together in a single model showed concordance for time to next relapse prediction that improved up to a combination of six markers: CSF[vitamin D binding protein + Factor I + C1inh/C1s] + serum[Factor B + IL4 + C1inh/C1s] (concordance 0.80, P < 0.001).

Cox regression modelling predicts time to disability

Four people with multiple sclerosis had already reached EDSS 6 by the time of sampling so were excluded from the analysis. Nineteen of the remaining 73 (26%) people with multiple sclerosis reached EDSS 6 during follow-up. In univariate Cox regression, adjusted for age and sex, the CSF biomarkers with highest concordance for predicting time to EDSS 6 were neurofilament light (concordance 0.92, P = 0.047) and C1inh/C1s (concordance 0.91, P = 0.79; Additional file 1: Table S11). The addition of biomarkers in a multivariate model improved concordance up to a combination of four biomarkers: neurofilament light + TNFR1 + iC3b + CXCL12 (concordance 0.93; p < 0.001; Table 4 and Additional file 1: Table S10).

In univariate analysis adjusted for age and sex, the serum biomarkers with highest signal for concordance for predicting time to EDSS 6 were C1inh/C1s (concordance 0.90, P = 0.91) and sCD27 (concordance 0.86, P = 0.86). The addition of further biomarkers to the model improved concordance up to a combination of four serum markers: CCL27 + Factor I + osteopontin + C1inh/C1s (concordance 0.94, P < 0.001; Table 4 and Additional file 1: Table S10). Combining CSF and serum markers together in a single model showed concordance for EDSS 6 prediction that improved up to a combination of six markers: CSF[C9 + neurofilament light] + serum[chitinase-3-like-1 + CCL27 + vitamin D binding protein + C1inh/C1s] (concordance 0.98; P < 0.001; Fig. 3, Additional file 1: Table S10).

Fig. 3
figure 3

Plot of the linear risk predictor from the optimum combined biomarker model and time to EDSS 6. The optimum combined model contained age, sex and CSF[C9 + neurofilament light] + serum[chitinase-3-like-1 + CCL27 + vitamin D binding protein + C1inh/C1s]

Discussion

Biomarkers hold promise to expedite diagnosis and stratify the risk of adverse multiple sclerosis prognosis, which would aid personalisation of therapy. Candidate multiple sclerosis biomarkers have been proposed but are often studied in isolation, and comparing studies that have used different clinical outcomes is challenging [12]. Given the complexity of multiple sclerosis biology, a combination of markers may improve predictions versus any marker in isolation. Here we have shown that combinations of up to six biomarkers measured in CSF or serum provide incremental improvements in the prediction of diagnosis and prognosis in multiple sclerosis (Fig. 2). Importantly, we found that predictions using combinations of blood-based biomarkers rivalled those using CSF markers, offering promise for less invasive approaches to monitoring and personalising multiple sclerosis care.

Neurofilament-light is a well-established candidate biomarker in multiple sclerosis and in our study was among the optimum single biomarkers to distinguish multiple sclerosis from control cases in both serum and CSF. Meta-analysis and systematic review have both highlighted evidence supporting neurofilament light as a marker to predict future multiple sclerosis disease activity, brain atrophy and to monitor treatment response to disease modifying therapies [17, 55]. Our model confirmed the predictive power of neurofilament light in CSF to distinguish multiple sclerosis from control cases (AUC=0.85). With the addition of further CSF markers, CSF NfL did not feature in our final predictive model (C9, chitinase-3-like-1, TNFR1, sCD27; AUC to 0.94) although this may have been affected by the relationship of NfL with age, since CSF NfL did feature in a model adjusted for the effects of age and sex within controls (Fig. 2, Table 3, Additional file 1: Table S8). The limitations of ELISA in detecting low concentrations of neurofilament light in blood have been overcome using highly sensitive SiMoA technology, and some clinics are starting to incorporate neurofilament light into practice [8]. In our study, while serum neurofilament light was one of the top biomarkers to distinguish multiple sclerosis from control cases, it did not feature in the final serum model of diagnosis. Some studies have demonstrated utility in combining neurofilament light with GFAP [18], which can also be measured using SiMoA. We used an ELISA assay for GFAP and were unable to adequately detect GFAP in CSF; serum GFAP detected using ELISA did not feature in any of our predictive models.

Chitinase-3-like-1 featured in CSF predictive models of multiple sclerosis versus control status in our cohort, and also as a serum marker in the combined CSF/serum model predicting time to EDSS 6. Chitinase-3-like-1 is a marker of macrophage and astrocyte activation and is potentially neurotoxic [19]. It has been found to be elevated in other studies of multiple sclerosis CSF [20,21,22], and has shown promise in predicting conversion from CIS to multiple sclerosis and future multiple sclerosis disease activity [22, 23]. Other markers that predicted multiple sclerosis versus control status in our study were chemokines/cytokines known to promote a pro-inflammatory milieu such as osteopontin, MCP1, CCL27, TNFR1 and soluble CD27. Osteopontin is a cytokine produced by a range of cells including macrophages, lymphocytes and dendritic cells, providing cross talk between the innate and adaptive immune system. Osteopontin regulates the differentiation of pro-inflammatory lymphocytes, inhibits apoptosis of inflammatory cells [24], and has a role in microglia-mediated synaptic engulfment [25]. CCL27 and MCP1 (CCL2) are both chemokines involved in regulating migration of macrophages/monocytes [26]. MCP1 is abundantly expressed by microglia located at the active rim of multiple sclerosis lesions [27], and is thought to play a role in multiple sclerosis pathogenesis [28], although we and others found it to be significantly lower in serum and CSF of people with multiple sclerosis than controls [29]. CCL27 induces the homing of memory T cells to sites of inflammation and has been found by others to be elevated in the serum of people with multiple sclerosis [30]. Tumour necrosis factor alpha (TNFα) is recognised as a key function in autoimmune disease, where excessive activation of TNFα mediates cytotoxic and pro-inflammatory responses via TNFR1 [31]. CD27 is a T-cell activation marker, whose soluble form has been shown by others to be significantly elevated in people with multiple sclerosis [32], and predictive of future multiple sclerosis disease activity [33], including transition from CIS to multiple sclerosis [34].

CSF biomarkers that combined to predict time to next relapse also included neurofilament light and MCP1. Vitamin D binding protein, a regulator of the distribution, stability and bioavailability of vitamin D, was one of the markers whose expression differed significantly between multiple sclerosis and control cases in serum and also featured in the CSF and serum models predicting relapse. There is evidence for an immunomodulatory role of vitamin D in multiple sclerosis [35]. Other studies have demonstrated some utility of vitamin D binding protein in distinguishing people with multiple sclerosis versus controls [36, 37], or risk of developing multiple sclerosis [38]. Vitamin D binding protein has been shown to be expressed on spinal cord neurons, pia mater and grey matter within the brains of people with multiple sclerosis, and in an animal model of multiple sclerosis, high vitamin D binding protein appeared to mitigate beneficial effects of vitamin D3 supplementation and inhibit recovery [37]. CXCL12, a chemoattractant protein for T cells as well as monocytes appeared in the CSF and serum models predicting relapse, and the CSF model predicting time to EDSS 6. CXCL12 has been found highly expressed in active multiple sclerosis lesions and appears to play a role in enhancing the inflammatory response in multiple sclerosis [39].

Our models also demonstrated evidence of complement activation and consumption. In our cohort, significantly differences in the concentration of several complement proteins between the CSF and serum of multiple sclerosis versus controls suggested dysregulation of this pathway in multiple sclerosis. The complement molecules TCC (terminal complement protein) and iC3b (activation product) featured in our final model to predict multiple sclerosis status. C1inh/C1s complex also featured prominently in the prediction of relapses and disability in our cohort. This is an indication of activation of the classical pathway of the complement cascade [40], and in line with our earlier work suggesting this is relevant to multiple sclerosis biology [16]. C3 and Factor H both featured in our serum model predicting relapses. We have previously shown both to be present at high levels in multiple sclerosis lesions [41], while others have shown these complement proteins to be elevated in the blood and CSF of people with multiple sclerosis versus controls, and the be relevant in predictions of disability outcomes [42, 43]. Overall, our data support previous findings that imply ongoing local and systemic complement dysregulation in multiple sclerosis [16, 44].

Others have explored the potential utility of combining more than one protein biomarker to predict multiple sclerosis diagnosis or prognosis. Lucchini et al. measured a panel of 8 protein candidates in the CSF of people with multiple sclerosis or controls. They found the combination of chitinase-3-like-1, CXCL10, CXCL12, CXCL13 increased the AUC for predicting conversion to clinically definite multiple sclerosis after first attack above any single biomarker in isolation [22]. Bielekova et al. [45] also found that the combination of CSF IL-12/IL-23p40, CXCL13 and IL-8 in CSF was more predictive of multiple sclerosis versus control status than any of these markers in isolation.

There are some practical considerations around combining biomarkers in clinical practice. Some biomarkers were not statistically significant in univariate analyses, but still provided a small differential addition to the final model. Conversely, some of the biomarkers that were most discriminatory in univariate analysis autocorrelated to some extent and therefore did not all appear in final combination models. There may be a balance between selecting the model that was statistically optimal versus a choosing smaller combination of markers that gives similar predictive value but is simpler to measure in multiplex (Additional file 1: Table S7 and S9). While optimum models sometimes combined CSF and serum markers, there are practical advantages in using a single sample type, ideally serum, which can be serially sampled more easily than CSF.

The present study is subject to some limitations or caveats. Using research operating procedures, we aimed for all blood and CSF samples to reach the freezer within 2 h. Variations in sample handling in routine clinical practice may introduce the risk of degradation of small molecules such as complement [46]. The demographics of our multiple sclerosis and control group differed, which could affect the results we are observing, even though we adjusted for sex and age in all analyses. Body mass index (BMI) data were unavailable for this cohort, so were not included, despite emerging evidence that correction of some biomarkers for BMI improves correlations [47]. Our cohort was too small to explore biomarker signatures of different multiple sclerosis subtypes. We included all people with multiple sclerosis in the analysis of time to relapse, even though people with progressive disease are somewhat less likely to experience relapses. While over-fitting of statistical models is a potential limitation, our test and train AUCs were very similar which suggests this was not the case in our study. However, the main caveat to our results is that this represents “discovery” work, which means that conclusions on the potential utility of these combinations could only be recognised with validation studies in independent cohorts.

In conclusion, using well-optimised assays for 24 candidate protein biomarkers in the blood and 20 CSF we have demonstrated that combination models showed better prediction of multiple sclerosis diagnosis and prognosis than single biomarkers. We also demonstrated for the first time that combination serum models rivalled those of CSF, holding promise for a non-invasive approach. This study supports the premise that combining several biomarkers in a single test will aid diagnostic and prognostic accuracy in the future. The utility of our biomarker models can only be established by robust validation in different and varied cohorts.

Availability of data and materials

The data that support the findings of this study, including code used in analyses, are available in the following repository: https://github.com/EmmaTallantyre/Biomarkers_MS_2024.

References

  1. Shah AA, Corboy JR. Discontinuing disease-modifying therapies in multiple sclerosis. Pract Neurol. 2022.

  2. Merkel B, Butzkueven H, Traboulsee AL, Havrdova E, Kalincik T. Timing of high-efficacy therapy in relapsing-remitting multiple sclerosis: a systematic review. Autoimmun Rev. 2017;16(6):658–65.

    Article  PubMed  Google Scholar 

  3. Tüzün E. Immunopathological factors associated with disability in multiple sclerosis. Arch Neuropsychiatry. 2018;55(Suppl 1):S26.

    Google Scholar 

  4. Confavreux C, Vukusic S, Moreau T, Adeleine P. Relapses and progression of disability in multiple sclerosis. N Engl J Med. 2000;343(20):1430–8.

    Article  CAS  PubMed  Google Scholar 

  5. Sellebjerg F, Royen L, Soelberg Sørensen P, Oturai AB, Jensen PEH. Prognostic value of cerebrospinal fluid neurofilament light chain and chitinase-3-like-1 in newly diagnosed patients with multiple sclerosis. Mult Scler J. 2019;25(11):1444–51.

    Article  CAS  Google Scholar 

  6. Kuhle J, Barro C, Andreasson U, Derfuss T, Lindberg R, Sandelius Å, et al. Comparison of three analytical platforms for quantification of the neurofilament light chain in blood samples: ELISA, electrochemiluminescence immunoassay and Simoa. Clin Chem Lab Med (CCLM). 2016;54(10):1655–61.

    Article  CAS  PubMed  Google Scholar 

  7. Huang J, Khademi M, Fugger L, Lindhe Ö, Novakova L, Axelsson M, et al. Inflammation-related plasma and CSF biomarkers for multiple sclerosis. Proc Natl Acad Sci. 2020;117(23):12952–60.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  8. Ferreira-Atuesta C, Reyes S, Giovanonni G, Gnanapavan S. The evolution of neurofilament light chain in multiple sclerosis. Front Neurosci. 2021;15: 642384.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Thebault S, Booth RA, Rush CA, MacLean H. Freedman multiple sclerosis. Serum neurofilament light chain measurement in multiple sclerosis: hurdles to clinical translation. Front Neurosci. 2021;15: 654942.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Comabella M, Fernandez M, Martin R, Rivera-Vallve S, Borras E, Chiva C, et al. Cerebrospinal fluid chitinase 3-like 1 levels are associated with conversion to multiple sclerosis. Brain. 2010;133(4):1082–93.

    Article  PubMed  Google Scholar 

  11. Börnsen L, Khademi M, Olsson T, Sørensen PS, Sellebjerg F. Osteopontin concentrations are increased in cerebrospinal fluid during attacks of multiple sclerosis. Mult Scler J. 2011;17(1):32–42.

    Article  Google Scholar 

  12. Castle D, Wynford-Thomas R, Loveless S, Bentley E, Howell OW, Tallantyre EC. Using biomarkers to predict clinical outcomes in multiple sclerosis. Pract Neurol. 2019;19(4):342–9.

    Article  PubMed  Google Scholar 

  13. Harding K, Williams O, Willis M, Hrastelj J, Rimmer A, Joseph F, et al. Clinical outcomes of escalation vs early intensive disease-modifying therapy in patients with multiple sclerosis. JAMA Neurol. 2019;76(5):536–41.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Kurtzke JF. Rating neurologic impairment in multiple sclerosis: an expanded disability status scale (EDSS). Neurology. 1983;33(11):1444.

    Article  CAS  PubMed  Google Scholar 

  15. Ingram G, Colley E, Ben-Shlomo Y, Cossburn M, Hirst CL, Pickersgill T, et al. Validity of patient-derived disability and clinical data in multiple sclerosis. Mult Scler J. 2010;16(4):472–9.

    Article  CAS  Google Scholar 

  16. Zelek WM, Fathalla D, Morgan A, Touchard S, Loveless S, Tallantyre E, et al. Cerebrospinal fluid complement system biomarkers in demyelinating disease. Mult Scler J. 2020;26(14):1929–37.

    Article  CAS  Google Scholar 

  17. Williams T, Zetterberg H, Chataway J. Neurofilaments in progressive multiple sclerosis: a systematic review. J Neurol. 2021;268:3212–22.

    Article  PubMed  Google Scholar 

  18. Barro C, Healy BC, Liu Y, Saxena S, Paul A, Polgar-Turcsanyi M, et al. Serum GFAP and neurofilament light levels differentiate subsequent progression and disease activity in patients with progressive multiple sclerosis. Neurol-Neuroimmunol Neuroinflamm. 2023;10(1).

  19. Matute-Blanch C, Calvo-Barreiro L, Carballo-Carbajal I, Gonzalo R, Sanchez A, Vila M, et al. Chitinase 3-like 1 is neurotoxic in primary cultured neurons. Sci Rep. 2020;10(1):7118.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  20. Burman J, Raininko R, Blennow K, Zetterberg H, Axelsson M, Malmeström C. YKL-40 is a CSF biomarker of intrathecal inflammation in secondary progressive multiple sclerosis. J Neuroimmunol. 2016;292:52–7.

    Article  CAS  PubMed  Google Scholar 

  21. Floro S, Carandini T, Pietroboni AM, De Riz MA, Scarpini E, Galimberti D. Role of Chitinase 3-like 1 as a biomarker in multiple sclerosis: a systematic review and meta-analysis. Neurol-Neuroimmunol Neuroinflamm. 2022;9(4).

  22. Lucchini M, De Arcangelis V, Piro G, Nociti V, Bianco A, De Fino C, et al. CSF CXCL13 and chitinase 3-like-1 levels predict disease course in relapsing multiple sclerosis. Mol Neurobiol. 2023;60(1):36–50.

    Article  CAS  PubMed  Google Scholar 

  23. Gil-Perotin S, Castillo-Villalba J, Cubas-Nuñez L, Gasque R, Hervas D, Gomez-Mateu J, et al. Combined cerebrospinal fluid neurofilament light chain protein and chitinase-3 like-1 levels in defining disease course and prognosis in multiple sclerosis. Front Neurol. 2019;10:1008.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Clemente N, Raineri D, Cappellano G, Boggio E, Favero F, Soluri MF, et al. Osteopontin bridging innate and adaptive immunity in autoimmune diseases. J Immunol Res. 2016;2016.

  25. De Schepper S, Ge JZ, Crowley G, Ferreira LS, Garceau D, Toomey CE, et al. Perivascular cells induce microglial phagocytic states and synaptic engulfment via SPP1 in mouse models of Alzheimer’s disease. Nat Neurosci. 2023;26(3):406–15.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Deshmane SL, Kremlev S, Amini S, Sawaya BE. Monocyte chemoattractant protein-1 (MCP-1): an overview. J Interferon Cytokine Res. 2009;29(6):313–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Tanuma N, Sakuma H, Sasaki A, Matsumoto Y. Chemokine expression by astrocytes plays a role in microglia/macrophage activation and subsequent neurodegeneration in secondary progressive multiple sclerosis. Acta Neuropathol. 2006;112:195–204.

    Article  CAS  PubMed  Google Scholar 

  28. Mahad DJ, Ransohoff RM, editors. The role of MCP-1 (CCL2) and CCR2 in multiple sclerosis and experimental autoimmune encephalomyelitis (EAE). Semin Immunol; 2003: Elsevier.

  29. Scarpini E, Galimberti D, Baron P, Clerici R, Ronzoni M, Conti G, et al. IP-10 and MCP-1 levels in CSF and serum from multiple sclerosis patients with different clinical subtypes of the disease. J Neurol Sci. 2002;195(1):41–6.

    Article  CAS  PubMed  Google Scholar 

  30. Khaiboullina SF, Gumerova AR, Khafizova IF, Martynova EV, Lombardi VC, Bellusci S, et al. CCL27: novel cytokine with potential role in pathogenesis of multiple sclerosis. BioMed Res Int. 2015;2015.

  31. Jang D-I, Lee A-H, Shin H-Y, Song H-R, Park J-H, Kang T-B, et al. The role of tumor necrosis factor alpha (TNF-α) in autoimmune disease and current TNF-α inhibitors in therapeutics. Int J Mol Sci. 2021;22(5):2719.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Hintzen RQ, van Lier RW, Kuijpers KC, Baars PA, Schaasberg W, Lucas CJ, et al. Elevated levels of a soluble form of the T cell activation antigen CD27 in cerebrospinal fluid of multiple sclerosis patients. J Neuroimmunol. 1991;35(1–3):211–7.

    Article  CAS  PubMed  Google Scholar 

  33. Komori M, Blake A, Greenwood M, Lin YC, Kosa P, Ghazali D, et al. Cerebrospinal fluid markers reveal intrathecal inflammation in progressive multiple sclerosis. Ann Neurol. 2015;78(1):3–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Mescheriakova JY, Runia TF, Jafari N, Siepman TA, Hintzen RQ. Soluble CD27 levels in cerebrospinal fluid as a prognostic biomarker in clinically isolated syndrome. JAMA Neurol. 2017;74(3):286–92.

    Article  PubMed  Google Scholar 

  35. Galoppin M, Kari S, Soldati S, Pal A, Rival M, Engelhardt B, et al. Full spectrum of vitamin D immunomodulation in multiple sclerosis: mechanisms and therapeutic implications. Brain Commun. 2022;4(4):fcac171.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Maghbooli Z, Omidifar A, Varzandi T, Salehnezhad T, Sahraian MA. Reduction in circulating vitamin D binding protein in patients with multiple sclerosis. BMC Neurol. 2021;21(1):1–8.

    Article  Google Scholar 

  37. Yang M, Qin Z, Zhu Y, Li Y, Qin Y, Jing Y, et al. Vitamin D-binding protein in cerebrospinal fluid is associated with multiple sclerosis progression. Mol Neurobiol. 2013;47:946–56.

    Article  CAS  PubMed  Google Scholar 

  38. Grut V, Biström M, Salzer J, Stridh P, Lindam A, Alonso-Magdalena L, et al. Free vitamin D3 index and vitamin D-binding protein in multiple sclerosis: a presymptomatic case–control study. Eur J Neurol. 2022;29(8):2335–42.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Calderon TM, Eugenin EA, Lopez L, Kumar SS, Hesselgesser J, Raine CS, et al. A role for CXCL12 (SDF-1α) in the pathogenesis of multiple sclerosis: regulation of CXCL12 expression in astrocytes by soluble myelin basic protein. J Neuroimmunol. 2006;177(1–2):27–39.

    Article  CAS  PubMed  Google Scholar 

  40. Hurler L, Toonen EJ, Kajdácsi E, van Bree B, Brandwijk RJ, de Bruin W, et al. Distinction of early complement classical and lectin pathway activation via quantification of C1s/C1-INH and MASP-1/C1-INH complexes using novel ELISAs. Front Immunol. 2022:6384.

  41. Ingram G, Loveless S, Howell OW, Hakobyan S, Dancey B, Harris CL, et al. Complement activation in multiple sclerosis plaques: an immunohistochemical analysis. Acta Neuropathol Commun. 2014;2(1):1–15.

    Article  Google Scholar 

  42. Aeinehband S, Lindblom RP, Al Nimer F, Vijayaraghavan S, Sandholm K, Khademi M, et al. Complement component C3 and butyrylcholinesterase activity are associated with neurodegeneration and clinical disability in multiple sclerosis. PLoS ONE. 2015;10(4): e0122048.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Ingram G, Hakobyan S, Loveless S, Robertson N, Morgan BP. Complement regulator factor H in multiple sclerosis. J Cell Biochem. 2011;112(10):2653–4.

    Article  CAS  PubMed  Google Scholar 

  44. Ingram G, Hakobyan S, Robertson NP, Morgan BP. Elevated plasma C4a levels in multiple sclerosis correlate with disease activity. J Neuroimmunol. 2010;223(1–2):124–7.

    Article  CAS  PubMed  Google Scholar 

  45. Bielekova B, Komori M, Xu Q, Reich DS, Wu T. Cerebrospinal fluid IL-12p40, CXCL13 and IL-8 as a combinatorial biomarker of active intrathecal inflammation. PLoS ONE. 2012;7(11): e48370.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  46. Huang J, Khademi M, Lindhe Ö, Jönsson G, Piehl F, Olsson T, et al. Assessing the preanalytical variability of plasma and cerebrospinal fluid processing and its effects on inflammation-related protein biomarkers. Mol Cell Proteomics. 2021;20.

  47. Yalachkov Y, Schäfer JH, Jakob J, Friedauer L, Steffen F, Bittner S, et al. Effect of estimated blood volume and body mass index on GFAP and neurofilament light levels in the serum and CSF of patients with multiple sclerosis. Neurol-Neuroimmunol Neuroinflamm. 2023;10(1).

  48. Marastoni D, Magliozzi R, Bolzan A, Pisani AI, Rossi S, Crescenzo F, et al. CSF levels of CXCL12 and osteopontin as early markers of primary progressive multiple sclerosis. Neurol Neuroimmunol Neuroinflamm. 2021;8(6): e1083.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Mungan S, Öztekin NS, Güzel I, Uysal S. The predictive role of CXCL12 (SDF-1 alpha) in multiple sclerosis. Turkiye Klinikleri J Neurol. 2018;13:62–6.

    Article  Google Scholar 

  50. Khademi M, Kockum I, Andersson ML, Iacobaeus E, Brundin L, Sellebjerg F, et al. Cerebrospinal fluid CXCL13 in multiple sclerosis: a suggestive prognostic marker for the disease course. Mult Scler J. 2011;17(3):335–43.

    Article  CAS  Google Scholar 

  51. Stampanoni Bassi M, Iezzi E, Drulovic J, Pekmezovic T, Gilio L, Furlan R, et al. IL-6 in the cerebrospinal fluid signals disease activity in multiple sclerosis. Front Cell Neurosci. 2020;14:120.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Wei Y, Chang H, Feng H, Li X, Zhang X, Yin L. Low serum interleukin-10 is an independent predictive factor for the risk of second event in clinically isolated syndromes. Front Neurol. 2019;10:604.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Petereit H, Pukrop R, Fazekas F, Bamborschke S, Röpele S, Kölmel H, et al. Low interleukin-10 production is associated with higher disability and MRI lesion load in secondary progressive multiple sclerosis. J Neurol Sci. 2003;206(2):209–14.

    Article  CAS  PubMed  Google Scholar 

  54. Kallaur AP, Oliveira SR, Delicato de Almeida ER, Kaminami Morimoto H, Lopes J, de Carvalho Jennings Pereira WL, et al. Cytokine profile in relapsing-remitting multiple sclerosis patients and the association between progression and activity of the disease. Mol Med Rep. 2013;7(3):1010–20.

    Article  CAS  PubMed  Google Scholar 

  55. Khaibullin T, Ivanova V, Martynova E, Cherepnev G, Khabirov F, Granatov E, et al. Elevated levels of proinflammatory cytokines in cerebrospinal fluid of multiple sclerosis patients. Front Immunol. 2017;8:531.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Soilu-Hänninen M, Koskinen J, Laaksonen M, Hänninen A, Lilius E-M, Waris M. High sensitivity measurement of CRP and disease progression in multiple sclerosis. Neurology. 2005;65(1):153–5.

    Article  PubMed  Google Scholar 

  57. Fettah E, Demir A. C-reactive protein/albumin ratio in patients with multiple sclerosis and its relationship with disease subtype and disability. J Surg Med. 2020;4(11):974–7.

    Article  Google Scholar 

  58. Yalachkov Y, Anschuetz V, Jakob J, Schaller-Paule MA, Schaefer JH, Reilaender A, et al. C-reactive protein levels and gadolinium-enhancing lesions are associated with the degree of depressive symptoms in newly diagnosed multiple sclerosis. Front Neurol. 2021;12.

  59. Kassubek R, Gorges M, Schocke M, Hagenston VA, Huss A, Ludolph AC, et al. GFAP in early multiple sclerosis: a biomarker for inflammation. Neurosci Lett. 2017;657:166–70.

    Article  CAS  PubMed  Google Scholar 

  60. Abdelhak A, Hottenrott T, Morenas-Rodríguez E, Suárez-Calvet M, Zettl UK, Haass C, et al. Glial activation markers in CSF and serum from patients with primary progressive multiple sclerosis: potential of serum GFAP as disease severity marker? Front Neurol. 2019;10:280.

    Article  PubMed  PubMed Central  Google Scholar 

  61. Hinsinger G, Galéotti N, Nabholz N, Urbach S, Rigau V, Demattei C, et al. Chitinase 3-like proteins as diagnostic and prognostic biomarkers of multiple sclerosis. Mult Scler J. 2015;21(10):1251–61.

    Article  CAS  Google Scholar 

  62. Orsi G, Cseh T, Hayden Z, Perlaki G, Nagy SA, Giyab O, et al. Microstructural and functional brain abnormalities in multiple sclerosis predicted by osteopontin and neurofilament light. Mult Scler Relat Disord. 2021;51: 102923.

    Article  CAS  PubMed  Google Scholar 

  63. Wen S-R, Liu G-J, Feng R-N, Gong F-C, Zhong H, Duan S-R, et al. Increased levels of IL-23 and osteopontin in serum and cerebrospinal fluid of multiple sclerosis patients. J Neuroimmunol. 2012;244(1–2):94–6.

    Article  CAS  PubMed  Google Scholar 

  64. Carrieri P, Provitera V, Perrella M, Tartaglia G, Busto A, Perrella O. Possible role of transforming growth factor-β in relapsing-remitting multiple sclerosis. Neurol Res. 1997;19(6):599–600.

    Article  CAS  PubMed  Google Scholar 

  65. Moharami S, Nourazarian A, Nikanfar M, Laghousi D, Shademan B, Joodi Khanghah O, et al. Investigation of serum levels of orexin-A, transforming growth factor β, and leptin in patients with multiple sclerosis. J Clin Lab Anal. 2022;36(1): e24170.

    Article  CAS  PubMed  Google Scholar 

  66. Mahad D, Howell S, Woodroofe M. Expression of chemokines in the CSF and correlation with clinical disease activity in patients with multiple sclerosis. J Neurol Neurosurg Psychiatry. 2002;72(4):498–502.

    CAS  PubMed  PubMed Central  Google Scholar 

  67. Naegelin Y, Saeuberli K, Schaedelin S, Dingsdale H, Magon S, Baranzini S, et al. Levels of brain-derived neurotrophic factor in patients with multiple sclerosis. Ann Clin Transl Neurol. 2020;7(11):2251–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Hohnoki K, Inoue A, Koh C-S. Elevated serum levels of IFN-γ, IL-4 and TNF-α/unelevated serum levels of IL-10 in patients with demyelinating diseases during the acute stage. J Neuroimmunol. 1998;87(1–2):27–32.

    Article  CAS  PubMed  Google Scholar 

  69. Tahani S, Dehghani L, Jahanbani-Ardakani H, Shaygannejad V, Fazli A, Hamidavi A, et al. Elevated serum level of IL-4 in neuromyelitis optica and multiple sclerosis patients. J Immunoassay Immunochem. 2019;40(5):555–63.

    Article  CAS  PubMed  Google Scholar 

  70. Martynova E, Goyal M, Johri S, Kumar V, Khaibullin T, Rizvanov AA, et al. Serum and cerebrospinal fluid cytokine biomarkers for diagnosis of multiple sclerosis. Mediat Inflamm. 2020;2020.

  71. Beck J, Rondot P, Catinot L, Falcoff E, Kirchner H, Wietzerbin J. Increased production of interferon gamma and tumor necrosis factor precedes clinical manifestation in multiple sclerosis: do cytokines trigger off exacerbations? Acta Neurol Scand. 1988;78(4):318–23.

    Article  CAS  PubMed  Google Scholar 

  72. Matejčíková Z, Mareš J, Sládková V, Svrčinová T, Vysloužilová J, Zapletalová J, et al. Cerebrospinal fluid and serum levels of interleukin-8 in patients with multiple sclerosis and its correlation with Q-albumin. Mult Scler Relat Disord. 2017;14:12–5.

    Article  PubMed  Google Scholar 

  73. Magliozzi R, Pezzini F, Pucci M, Rossi S, Facchiano F, Marastoni D, et al. Changes in cerebrospinal fluid balance of TNF and TNF receptors in naïve multiple sclerosis patients: early involvement in compartmentalised intrathecal inflammation. Cells. 2021;10(7):1712.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  74. Mikova O, Yakimova R, Bosmans E, Kenis G, Maes M. Increased serum tumor necrosis factor alpha concentrations in major depression and multiple sclerosis. Eur Neuropsychopharmacol. 2001;11(3):203–8.

    Article  CAS  PubMed  Google Scholar 

  75. Nicoletti F, Di Marco R, Mangano K, Patti F, Reggio E, Nicoletti A, et al. Increased serum levels of interleukin-18 in patients with multiple sclerosis. Neurology. 2001;57(2):342–4.

    Article  CAS  PubMed  Google Scholar 

  76. Rinaldi AO, Sanseverino I, Purificato C, Cortese A, Mechelli R, Francisci S, et al. Increased circulating levels of vitamin D binding protein in multiple sclerosis patients. Toxins. 2015;7(1):129–37.

    Article  PubMed  PubMed Central  Google Scholar 

  77. Mashayekhi F, Salehi Z. Expression of leukemia inhibitory factor in the cerebrospinal fluid of patients with multiple sclerosis. J Clin Neurosci. 2011;18(7):951–4.

    Article  CAS  PubMed  Google Scholar 

  78. Ingram G, Hakobyan S, Hirst CL, Harris CL, Loveless S, Mitchell JP, et al. Systemic complement profiling in multiple sclerosis as a biomarker of disease state. Mult Scler J. 2012;18(10):1401–11.

    Article  CAS  Google Scholar 

  79. Tatomir A, Talpos-Caia A, Anselmo F, Kruszewski AM, Boodhoo D, Rus V, et al. The complement system as a biomarker of disease activity and response to treatment in multiple sclerosis. Immunol Res. 2017;65(6):1103–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Li Y, Qin Z, Yang M, Qin Y, Lin C, Liu S. Differential expression of complement proteins in cerebrospinal fluid from active multiple sclerosis patients. J Cell Biochem. 2011;112(7):1930–7.

    Article  CAS  PubMed  Google Scholar 

  81. de Flon P, Laurell K, Sundström P, Blennow K, Söderström L, Zetterberg H, et al. Comparison of plasma and cerebrospinal fluid neurofilament light in a multiple sclerosis trial. Acta Neurol Scand. 2019;139(5):462–8.

    Article  PubMed  Google Scholar 

  82. Ning L, Wang B. Neurofilament light chain in blood as a diagnostic and predictive biomarker for multiple sclerosis: a systematic review and meta-analysis. PLoS ONE. 2022;17(9): e0274565.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would like to acknowledge Zoltán Prohászka, (Department of Internal Medicine and Hematology, Semmelweis University, Budapest, 1085, Hungary, Research Group for Immunology and Hematology, Semmelweis University-Eötvös Loránd Research Network (Office for Supported Research Groups), Budapest, 1085, Hungary) and Erik J.M. Toonen (R&D Department, Hycult Biotech, 5405 PB, Uden, The Netherlands). For the purpose of open access, the author has applied a CC BY public copyright licence to any Author Accepted Manuscript version arising.

Funding

This work was supported by the MRC Confidence in Concept scheme (Cardiff University, 2019). EK’s salary and consumables were supported by the MRC Confidence in Concept scheme (Cardiff University, 2019). WMZ is supported by Race Against Dementia Alzheimer’s Research UK Fellowship. SL receives salary support by the Health and Care Research Wales funded BRAIN Unit.

Author information

Authors and Affiliations

Authors

Contributions

SL, ECT, NPR and WMZ conceived the project. EK, SL and WMZ performed laboratory work. WJW performed the statistical analysis and KLK contributed to interpretation and contextualisation. AR, VA, SL and ECT collated clinical data on the cohort. LH provided assay(s) and advice for the laboratory analysis. EK and ECT drafted the manuscript and all authors reviewed and edited the final draft.

Corresponding author

Correspondence to Emma C. Tallantyre.

Ethics declarations

Ethics approval and consent to participate

All participants provided written informed consent. This study has Research Ethics Committee approval REC Ref: 05/ WSE03/111, 19/WA/0289.

Competing interests

The authors report no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

An overview of the assays used for this study. BDNF: brain-derived neurotrophic factor; CCL27: CC motif chemokine ligand 27; CRP: C-reactive protein; CXCL: C-X-C motif chemokine ligand; GFAP: glial fibrillary acidic protein; IL: interleukin; LIF: leukaemia inhibitory factor; MCP1: monocyte chemoattractant protein 1; TGF-β: transforming growth factor beta; TNFα: tumour necrosis factor alpha; TNFR1: tumour necrosis factor receptor 1; TCC: terminal complement complex, VDBP: vitamin D binding protein. Table S2. Rates of missing or imputed data in assay results. Table S3. Characteristics of 4 test/ train cohorts. Table S4. Comparison of CSF biomarkers between multiple sclerosis and non-multiple sclerosis groups. Mann–Whitney comparisons, corrected for multiple comparisons using Benjamini–Hochberg procedure. Table S5. Comparison of serum biomarkers between multiple sclerosis and non-multiple sclerosis groups Mann–Whitney comparisons corrected for multiple comparisons using Benjamini–Hochberg procedure. Table S6. A progression from single through combinations of multiple biomarkers to predict multiple sclerosis versus non-multiple sclerosis status. Table S7. Breakdown of the Train / Test results for the combined CSF & serum modelling of MS versus non-MS status. Mean AUC values were ordered from lowest to highest, and the optimum model was selected when addition of a further analyte resulted in an AUC increase < 0.01. Table S8. Sensitivity analysis: all biomarker concentrations were corrected for age and sex according to a linear model generated in control samples. A progression from single through combinations of multiple biomarkers to predict multiple sclerosis versus non-multiple sclerosis status. Table S9. Concordance of biomarkers in predicting time to next relapse in univariate analysis (adjusted for sex and age). Table S10. A progression from single through combinations of multiple biomarkers to predict time to relapse and time to disability, adjusted for age and sex. Table S11. Concordance of biomarkers in predicting time to EDSS 6 in univariate analysis (adjusted for sex, age and disease modifying therapy). Fig. S1. The plots show the 3 models below, run on 1000 random selection of (approx 75%) Train and (approx. 25%) Test data. The left hand plot is the range of AUC’s from the Train data when used also as Test and the right-hand plot shows the range of AUC’s when Test data are used as test.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kodosaki, E., Watkins, W.J., Loveless, S. et al. Combination protein biomarkers predict multiple sclerosis diagnosis and outcomes. J Neuroinflammation 21, 52 (2024). https://doi.org/10.1186/s12974-024-03036-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12974-024-03036-4

Keywords