Abstract
Background:
Spinal cord volume loss is associated with clinical disability in multiple sclerosis (MS). Aerobic capacity may mitigate the impact of central nervous system (CNS) damage accumulation, exerting beneficial effects on MS-related disability.
Objectives:
We investigated whether aerobic capacity could moderate the association between spinal cord atrophy and clinical disability in MS.
Methods:
In this cross-sectional analysis, expanded disability status scale (EDSS), peak of oxygen consumption (VO2peak), brain volumetric measures, and the normalized mean upper cervical cord area (nMUCCA) were collected from 51 MS patients and 33 healthy controls (HCs). Low aerobic capacity was defined as having a VO2peak z-score less than –1.64 standard deviations. In MS patients, we explored whether the association between nMUCCA and EDSS is moderated by the level of aerobic capacity.
Results:
The relationship between nMUCCA and EDSS was moderated by aerobic capacity, with a significant nMUCCA × aerobic capacity interaction (β = −0.099, 95% bootstrapped confidence interval [CI] = [−0.172; −0.014], p = 0.012). Lower nMUCCA was significantly associated with higher EDSS score in MS patients with low aerobic capacity (β = −0.073, p < 0.001), but not in those with high aerobic capacity (β = 0.026, p = 0.417).
Conclusions:
In MS patients with mild disability, higher aerobic capacity can potentially mitigate the negative impact of spinal cord damage on clinical disability.
Introduction
Atrophy of the spinal cord occurs from the earliest phases of multiple sclerosis (MS) and is strongly associated with clinical disability and disease progression. 1 Improvements in magnetic resonance imaging (MRI) technologies have allowed to better evaluate the relevance of spinal cord damage in explaining MS-related disability. 1 Although spinal cord atrophy is typically evaluated using dedicated MRI sequences, 1 normalized mean upper cervical cord area (nMUCCA) has been proposed as a reliable measure of spinal cord atrophy, obtainable from high-resolution brain three-dimensional (3D) T1-weighted MRI scans. 2
Disability can be considered as the result of a complex interplay between an individual’s health condition and contextual factors, such as personal and environmental influences. 3 Accordingly, several factors could influence the relation between spinal cord damage and disability in MS patients. For instance, psychological personal aspects like depression, or the level of self-efficacy can potentially have an impact on the disability level of MS patients.4,5 In addition to psychological aspects, aerobic capacity may represent another important element able to impact the relation between spinal cord neurodegeneration and disability. Aerobic capacity can be viewed as a personal physical resource for MS patients, reflecting their overall health status and potentially mitigating or protecting against the progression of motor impairment associated with MS pathology. 6 Indeed, aerobic capacity, which represents cardiovascular and respiratory fitness, refers to the ability of the body to utilize oxygen efficiently during prolonged physical activity. The peak of maximum oxygen consumption (VO2peak) obtained during a cardiopulmonary exercise testing (CPET) represents the gold standard measure for aerobic capacity. In MS patients, aerobic activity has emerged as a potential lifestyle factor with positive effects on overall fitness, disability, and brain integrity.7,8 Some exploratory studies have suggested several hypotheses to explain the mechanisms behind these positive effects, including neuroplasticity, 9 anti-inflammatory effects,10,11 and improved immune regulation,12,13 although with conflicting evidence.
Here, we explored whether aerobic capacity may moderate the association between spinal cord atrophy and clinical disability in a well-characterized group of mildly disabled MS patients.
Methods
Ethics approval
Approval was received from the institutional ethical standards committee on human experimentation of IRCCS Ospedale San Raffaele for any experiments using human subjects (Protocol number 59/INT/2015). Written informed consent was obtained from all subjects prior to study participation according to the Declaration of Helsinki.
Study design and population
This study is a secondary retrospective cross-sectional analysis conducted using baseline data obtained from an ongoing interventional rehabilitative protocol at our Institute (ClinicalTrials.gov ID NCT04097418). From the Neuroimaging Research Unit, IRCCS San Raffaele Scientific Institute (Milan, Italy) database, we retrospectively selected 51 consecutive patients with a diagnosis of MS according to the 2017 McDonald criteria and 33 healthy controls (HCs), meeting the following inclusion criteria: 18–65 range of age; no neurological (apart MS), psychiatric, orthopedic, or rheumatological disorders; no concomitant antidepressants, myorelaxants, beta-blockers, psychoactive, or steroids therapies; no drug or alcohol abuse history; no MRI or CPET contraindications. 14 In addition, patients with MS had to be relapsed and steroid free in the 3 months preceding study enrollment, with a stable disease-modifying treatment for 1 month or more and with an expanded disability status scale (EDSS) score ⩽ 6.0. In accordance with the criteria set by the original study, we excluded all participants engaged in structured physical activities or sport sessions three or more days per week, to eliminate the inclusion of well-trained individuals. The flow diagram of the study selection process is shown in Figure 1.

Flow diagram of the study selection process.
Clinical assessment and CPET procedure
Within three days from MRI acquisition, a neurological examination, with EDSS scoring and clinical phenotype definition (relapsing-remitting or progressive MS) was performed by an experienced neurologist blinded to MRI findings. At the same moment, all subjects performed the timed 25-foot walk test (T25FWT) to assess walking speed and the nine-hole peg test (9-HPT) to assess manual dexterity.
CPET was conducted by a cardiologist specialized in cardiac rehabilitation on a recumbent stepper (NuStep T5XR, Ann Arbor, MI, USA) to obtain peak of maximum oxygen consumption (VO2peak), considered the gold standard measure of aerobic capacity. 15 The validity of all CPET exams was carefully evaluated as previously described. 9 Body mass index (BMI) was calculated for each subject.
Low aerobic capacity was defined as a VO2peak z-score at least 1.64 standard deviations below the mean value of HC. Although no specific VO2peak thresholds tailored to our study population were identified in the literature, a threshold of –1.64 SD corresponds approximately to the fifth percentile, a commonly used criterion for defining abnormal cases. In this context, this threshold served as a practical and reliable criterion with adequate specificity for classifying patients with low aerobic capacity relative to the physiological levels observed in HCs.
MRI acquisition
Using a 3.0 Tesla Philips Ingenia CX scanner (Philips Medical Systems), the following brain MRI sequences were acquired with standardized procedures for subjects positioning (receiving Coil = dS-Head-32): (a) axial dual-echo turbo spin echo (TSE), field of view (FOV) = 240 × 240 mm, pixel size = 0.94 × 0.94 mm, 50 slices, 3 mm thick, matrix = 256 × 256, repetition time (TR) = 2599 ms, echo time (TE) = 16/80 ms; (b) 3D T1-weighted magnetization-prepared rapid gradient-echo (MPRAGE), FOV = 256 × 256, pixel size = 1 × 1 mm, 204 slices, 1 mm thick, matrix = 256 × 256, TR = 7 ms, TE = 3.2 ms, TI = 1000 ms, flip angle = 8°.
MRI analysis
T2-hyperintense white matter (WM) lesion volumes were measured using a local thresholding segmentation technique (Jim 8.0).
After lesion refilling, normalized brain volume (NBV), gray matter volume (NGMV), WM volume (NWMV), and thalamic volume were quantified on 3D T1-weighted images using FSL-SIENAx software and FSL-FIRST tool, applying the scaling factor derived from FSL-SIENAx.
Using the Jim8.0 software, the active surface (AS) method was applied on 3D T1-weighted scans of the brain to derive MUCCA. Briefly, the cord center was manually identified on axial oriented 3D T1-weighted scans by placing landmarks at the extremes of the studied cord region (i.e. the cranial limit of the odontoid process and the caudal border of C2) and approximately every 10 mm between these landmarks. These markers were used by an automated procedure to define cord contours. The mean MUCCA was obtained as the total cord volume divided by the cord length. Normalized MUCCA (nMUCCA) was obtained adjusting for FSL SIENAx brain scaling factor, inversely proportional to the total intracranial volume, according to the formula previously reported. 16
Statistical analysis
Demographic, clinical, and MRI data differences between patients and HCs, between MS patients with low or high aerobic capacity, and according to disease phenotype were tested using Chi-square, two-sample t, Mann–Whitney test, and one-way analysis of variance (ANOVA) test with Bonferroni post hoc group comparisons, as appropriate. T2-hyperintense WM LV were log-transformed. The normality of clinical and MRI variables was assessed using the Shapiro–Wilk test. In MS patients, we performed Spearman’s rank correlations between clinical measures (EDSS score, T25FWT and 9-HPT) and brain volumetric MRI measures and nMUCCA. Correlation analyses were also performed separately for MS patients with low and high aerobic capacity. Finally, we conducted a 5000 bootstrap samples moderation analysis corrected for age and sex considering nMUCCA as independent variable, EDSS as dependent variable and aerobic capacity as moderator. Analyses were performed using SPSS software (version 28.0, PROCESS Macro SPSS). The p-values < 0.05 were deemed statistically significant.
Results
Demographic and MRI findings
Compared to HCs, MS patients were significantly older (47 vs 44 years of age, p < 0.001) whereas sex and BMI were not significantly different (p ⩾ 0.560; see Table 1). MS patients had significantly worse T25FWT and 9-HPT performance (p < 0.001), and they had significantly higher brain T2-hyperintense WM LV (p < 0.001), as well as significantly lower NBV, NGMV, NWMV, thalamic volume, and nMUCCA (all p < 0.001; see Table 1).
Comparison of demographic, clinical, and MRI characteristics between HCs and MS patients, as well as between MS with low and high aerobic capacity.
9-HPT: nine hole peg test; BMI: body mass index; DMT: disease-modifying therapy; EDSS: expanded disability status scale; HC: healthy control; IQR: interquartile range; kg: kilograms; min: minutes; mL: milliliters; mm2: square millimeters; MRI: magnetic resonance imaging; MS: multiples sclerosis; n: number; NBV: normalized brain volume; NGMV: normalized gray matter volume; nMUCCA: normalized mean upper cervical cord area; NWMV: normalized white matter volume; P: progressive phenotype; RR: relapsing-remitting phenotype; SD: standard deviation; T25FWT: timed 25-foot walk test; VO2peak: peak of maximal oxygen consumption; WM LV: white matter lesion volume.
t-test.
Pearson’s Chi-square test.
Mann–Whitney test.
Compared to relapsing-remitting MS patients, those with a progressive disease were more frequently males, had significantly higher EDSS score (p = 0.012), worse T25FWT and 9-HPT performance (p ⩽ 0.024). VO2peak values or MRI data were not significantly different between the two groups (see Supplementary Table 1).
CPET findings
CPET reliability criteria have been proposed within the context of MS population to ensure consistent and reproducible results during exercise testing especially in patients with mild disability, considering the physiological and functional challenges associated with this disease. 17
In all HC subjects considered for this secondary analysis, a minimum of three out of four CPET criteria were consistently fulfilled. Specifically, the VO2 plateau criterion was achieved in 75% of the sample, the respiratory exchange ratio (RER) ⩾ 1.10 criterion in 100%, the rating of perceived exertion (RPE) measured on the Borg Scale ⩾ 17 criterion in 100%, and the HR criterion in 92%.
In all MS patients considered for this secondary analysis, a minimum of two out of four CPET criteria were consistently fulfilled. Specifically, the VO2 plateau criterion was achieved in 53% of the sample, the RER ⩾ 1.10 criterion in 63%, the RPE ⩾ 17 criterion in 87%, and the HR criterion in 17%. Compared to HC, MS patients had a lower value of zVO2peak (p < 0.001). Thirty-one (61%) patients were characterized by low aerobic capacity. Compared to MS patients with high aerobic capacity, those with low aerobic capacity showed significantly higher BMI values (p ⩽ 0.003), no other demographic, clinical and MRI measures differences were found (p ⩾ 0.065; see Table 1).
Associations between clinical and MRI measures
In MS patients, a higher EDSS score was significantly associated with lower nMUCCA (r = −0.353, p = 0.011; see Figure 2, Table 2). Such an association remained statistically significant only for MS patients with low aerobic capacity (r = −0.536, p = 0.002; see Supplementary Table 2). Worse 9-HPT performance was significantly associated with lower NBV, NGMV, and thalamic volume (r values from −0.294 to −0.466, p ⩽ 0.027), whereas the association with nMUCCA was not statistically significant (p = 0.154; see Table 2). T25FWT performance showed no significant association with brain and cord MRI findings (see Table 2).

Segmentation of MUCCA measure on (a) 3D T1-weighted image and (b) scatter plot within the MS group highlighting the association between EDSS and nMUCCA.
Spearman rank correlations between clinical and MRI measures in MS patients.
9-HPT: nine-hole peg test; EDSS: expanded disability status scale; NBV: normalized brain volume; NGMV: normalized gray matter volume; nMUCCA: normalized mean upper cervical cord area; NWMV: normalized white matter volume; r: r-value; T25FWT: timed 25-foot walk test; WM LV: white matter lesion volume.
Statistically significant associations are shown in bold.
Moderation analysis
In MS patients, VO2peak and nMUCCA were normally distributed, as confirmed by the Shapiro–Wilk test (p = 0.206 and p = 0.486, respectively), while the EDSS was not normally distributed (p = 0.006). Therefore, we decided to conduct the moderation analysis using bias corrected and accelerated 5000-sample bootstrapping.
In MS patients, the interaction between nMUCCA and aerobic capacity was significant (β = −0.099, 95% bootstrapped confidence interval [CI] = [−0.172; −0.014], p = 0.012), indicating that the relationship between nMUCCA and EDSS was moderated by aerobic capacity. Lower nMUCCA was significantly associated with higher EDSS score only in MS patients with low aerobic capacity (β = −0.073, 95% CI = [−0.115; −0.032], p < 0.001), whereas no significant association was found in MS patients with high aerobic capacity (β = 0.026, 95% CI = [−0.037; 0.089], p = 0.417; see Figure 3).

Moderation analysis corrected for age and sex performed in MS patients considering EDSS as outcome variable, nMUCCA as predictor and aerobic capacity as moderator.
Discussion
We found that, compared to HC, MS patients showed a significant lower nMUCCA,18–20 which was significantly associated with more severe disability,18–20 thus supporting the relevance of spinal cord volume loss in determining clinical impairment. 1 Moreover, in line with previous studies, 21 we found that MS patients had significantly lower aerobic capacity compared to HC, suggesting a state of physical deconditioning.
Our results suggest that aerobic capacity may influence the association between spinal cord atrophy and clinical disability. We found a significant interaction between nMUCCA and aerobic capacity, along with the evidence that nMUCCA was negatively associated with EDSS scores only in MS patients with low aerobic capacity. Considering that no significant differences were found in demographic, clinical, and MRI measures between MS patients with low or high aerobic capacity, except for a significantly lower BMI in those with higher aerobic capacity, these findings suggest that higher aerobic capacity may mitigate the detrimental impact of irreversible spinal cord volume loss on MS-related clinical disability. In contrast, in MS patients with lower aerobic capacity, the severity of spinal cord damage has a more pronounced impact on clinical disability. Therefore speculating, a higher aerobic capacity may represent a personal physical resource useful for reducing the impact of spinal cord damage. Indeed, better cardiovascular fitness seems to be associated with better brain health through multiple mechanisms. Although no univocal evidence has been published within the MS field, aerobic exercise (which usually results in better aerobic capacity) may lead to improvements in the cerebral blood flow,22,23 to the promotion of neuroplasticity9,24 and stimulation of brain-derived neurotrophic factors (BDNFs). 25 In addition, better aerobic capacity may positively influence mood and fatigue, as well as reduce stress, anxiety, and depressive symptoms.12,26–28
Moreover, aerobic exercise may have potential direct effect not only on the brain but also on the spinal cord itself. 29 Some experimental studies seem to suggest that exercise may reduce immune cell infiltration, including T cells, macrophages, and B cells, and decreases markers of inflammation, oxidative stress, and axonal damage. It also preserves motor neurons in the lumbar spinal cord and enhances myelin integrity. These changes are linked to exercise-induced modulation of inflammatory cytokines, such as interferon gamma (IFN-γ) and interleukin-17 (IL-17), and reduced blood–brain barrier permeability. 29 Moreover, exercise influences neurotrophic factors, such as BDNF and nerve growth factor (NGF), which support neural plasticity and repair in the spinal cord and brain. 29 Further investigations are needed to better understand the underling mechanisms and potential therapeutic targets of aerobic exercise and active lifestyle.
Unexpectedly, 19 we found no significant association between brain volumetric MRI measures and EDSS score, or T25FWT performance. Similar results were obtained when the analysis was stratified by MS patients with low and high aerobic capacity. These findings may be explained by the evaluation of MS patients with relatively mild disability, limited variability in EDSS scores, and a narrow range of walking speed, combined with a global rather than regional assessment of brain tissue loss. On the contrary, 9-HPT performance was associated with brain volumetric measures, potentially representing a more appropriate clinical measure of neurodegeneration in these early phases of the disease. However, we confirmed the more relevant role of spinal cord rather than brain damage in explaining disability.1,30 Nonetheless, since our study included MS patients with mild disability, our results should be interpreted with caution. Accordingly, they cannot be generalized to all MS patients and require confirmation in future studies involving larger cohorts, particularly those with higher levels of disability.
Our study has several limitations. Although age has been shown to influence VO2peak values, with a gradual decline across the lifespan similarly observed in both HC and MS patients, 31 we did not apply any age-related correction when calculating zVO2peak for determining aerobic capacity levels. Unfortunately, previously published data utilized a different protocol compared to our study for determining VO2peak values, and our small cohort of HC did not allow for a reliable estimation of the age effect on VO2peak. As a result, some older MS patients may have been wrongly classified as having low aerobic capacity, despite being comparable to age-matched HC. This could explain the relatively high prevalence (61%) of MS patients classified as having low aerobic capacity, despite their relatively low disability level (median EDSS score = 2.5). Due to the design of the study, HC and MS patients who regularly engage in sports or physical activity were not included. Moreover, detailed information regarding the amount of sport participation, exercise habits, and lifestyle was not collected. Therefore, this study mainly explored the impact of an active lifestyle, rather than the effects of physical activity. Future studies that incorporate a comprehensive assessment of physical activity and lifestyle, as well as MS patients who regularly engage in physical activity, could provide further insights into the beneficial effects of cardiovascular fitness in mitigating the detrimental impact of spinal cord damage on clinical disability.
Finally, the cross-sectional nature of our study precluded direct assessment of whether longitudinal changes in aerobic capacity (perhaps obtained directly through an aerobic training period) could impact the relationship between nMUCCA and EDSS score, as would be possible with an interventional longitudinal design. Moreover, due to the cross-sectional design of this study, causality between a high level of aerobic capacity and low disability, as suggested by our data, cannot be established. MS patients with lower levels of disability are more likely to engage in physical activity, which could, in turn, contribute to higher aerobic capacity. However, in the direct comparison, we found that while progressive MS patients had significantly higher EDSS scores compared to relapsing-remitting MS patients, no significant between-group difference in VO2peak values was found. These additional findings suggest that aerobic capacity may be, at least partially independent from the severity of clinical disability in MS patients. Nevertheless, further studies, particularly longitudinal ones, are required to better explore this association.
In conclusion, our study showed that a higher aerobic capacity may act as a personal physical resource able to reduce the impact of spinal cord damage on clinical disability in MS patients. For this reason, our study supports the importance of aerobic exercise and regular physical activity aimed at improving aerobic capacity in MS patients.32,33 This approach can be a rewarding strategy complementing specific pharmacologic therapies to improve MS patients’ management and their global health.32,33
Supplemental Material
sj-docx-1-msj-10.1177_13524585251318647 – Supplemental material for Aerobic capacity moderates the association between cervical cord atrophy and clinical disability in mildly disabled multiple sclerosis patients
Supplemental material, sj-docx-1-msj-10.1177_13524585251318647 for Aerobic capacity moderates the association between cervical cord atrophy and clinical disability in mildly disabled multiple sclerosis patients by Matteo Albergoni, Paolo Preziosa, Alessandro Meani, Chiara Dallari, Paola Valsasina, Maria A Rocca and Massimo Filippi in Multiple Sclerosis Journal
Supplemental Material
sj-docx-2-msj-10.1177_13524585251318647 – Supplemental material for Aerobic capacity moderates the association between cervical cord atrophy and clinical disability in mildly disabled multiple sclerosis patients
Supplemental material, sj-docx-2-msj-10.1177_13524585251318647 for Aerobic capacity moderates the association between cervical cord atrophy and clinical disability in mildly disabled multiple sclerosis patients by Matteo Albergoni, Paolo Preziosa, Alessandro Meani, Chiara Dallari, Paola Valsasina, Maria A Rocca and Massimo Filippi in Multiple Sclerosis Journal
Supplemental Material
sj-docx-3-msj-10.1177_13524585251318647 – Supplemental material for Aerobic capacity moderates the association between cervical cord atrophy and clinical disability in mildly disabled multiple sclerosis patients
Supplemental material, sj-docx-3-msj-10.1177_13524585251318647 for Aerobic capacity moderates the association between cervical cord atrophy and clinical disability in mildly disabled multiple sclerosis patients by Matteo Albergoni, Paolo Preziosa, Alessandro Meani, Chiara Dallari, Paola Valsasina, Maria A Rocca and Massimo Filippi in Multiple Sclerosis Journal
Footnotes
Acknowledgements
The authors acknowledge the role of Dr. Francesco Maranta in performing the CPET.
Author Contributions
Study concept and design was done by M.A., M.A.R., and M.F. Data acquisition and analysis were contributed by M.A., P.P., A.M., C.D., P.V., and M.A.R. Manuscript was drafted by M.A. and P.P. Intellectual content and manuscript editing was done by M.A., P.P., A.M., C.D., P.V., M.A.R., and M.F. Guarantor was M.F.
Data Availability
The dataset used and analyzed during this study is available from the corresponding author on reasonable request.
Declaration of Conflicting Interests
The author(s) declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: The authors declare that they have no competing interests in relation to this work.
Potential conflicts of interest outside the submitted work are as follows:
M. Albergoni, C. Dallari, A. Meani, and P. Valsasina have nothing to disclose. P. Preziosa received speaker honoraria from Roche, Biogen, Novartis, Merck, Bristol Myers Squibb, Genzyme, Horizon, and Sanofi. He has received research support from Italian Ministry of Health and Fondazione Italiana Sclerosi Multipla. M. Filippi is Editor-in-Chief of the Journal of Neurology, Associate Editor of Human Brain Mapping, Neurological Sciences, and Radiology; received compensation for consulting services from Alexion, Almirall, Biogen, Merck, Novartis, Roche, Sanofi; speaking activities from Bayer, Biogen, Celgene, Chiesi Italia SpA, Eli Lilly, Genzyme, Janssen, Merck-Serono, Neopharmed Gentili, Novartis, Novo Nordisk, Roche, Sanofi, Takeda, and TEVA; participation in Advisory Boards for Alexion, Biogen, Bristol-Myers Squibb, Merck, Novartis, Roche, Sanofi, Sanofi-Aventis, Sanofi-Genzyme, Takeda; scientific direction of educational events for Biogen, Merck, Roche, Celgene, Bristol-Myers Squibb, Lilly, Novartis, Sanofi-Genzyme; he receives research support from Biogen Idec, Merck-Serono, Novartis, Roche, the Italian Ministry of Health, the Italian Ministry of University and Research, and Fondazione Italiana Sclerosi Multipla. M.A. Rocca received consulting fees from Biogen, Bristol Myers Squibb, Eli Lilly, Janssen, Roche; and speaker honoraria from AstraZaneca, Biogen, Bristol Myers Squibb, Bromatech, Celgene, Genzyme, Horizon Therapeutics Italy, Merck Serono SpA, Novartis, Roche, Sanofi and Teva. She receives research support from the MS Society of Canada, the Italian Ministry of Health, the Italian Ministry of University and Research, and Fondazione Italiana Sclerosi Multipla. She is Associate Editor for Multiple Sclerosis and Related Disorders and Co-Associate Editor for Europe and Africa for Multiple Sclerosis Journal.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study has been partially supported by the Italian Ministry of Health (GR-2019-12369599).
Ethical Considerations
Approval was received from the institutional ethical standards committee on human experimentation of IRCCS Ospedale San Raffaele for any experiments using human subjects (Protocol number 59/INT/2015). Written informed consent was obtained from all subjects prior to study participation according to the Declaration of Helsinki.
Supplemental Material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
