Evolution of Brain Volume Loss Rates in Early Stages of Multiple Sclerosis
Citation Manager Formats
Make Comment
See Comments

Abstract
Objective To describe the dynamics of brain volume loss (BVL) at different stages of relapsing-remitting multiple sclerosis (RRMS), to describe the association between BVL and clinical measures, and to investigate an effect of treatment escalation on the rate of BVL.
Methods Together, 1903 patients predominantly with RRMS from the Avonex-Steroids-Azathioprine cohort (N = 166), the study of early IFN-β1a treatment cohort (N = 180), and the quantitative MRI cohort (N = 1,557) with ≥2 MRI scans and ≥1-year of follow-up were included. Brain MRI scans (N = 7,203) were performed using a single 1.5-T machine. Relationships between age or disease duration and global and tissue-specific BVL rates were analyzed using mixed models.
Results Age was not associated with the rate of BVL (β = −0.003; Cohen f2 = 0.0005; adjusted p = 0.39). Although disease duration was associated with the rate of BVL, its effect on the BVL rate was minimal (β = −0.012; Cohen f2 = 0.004; adjusted p = 4 × 10−5). Analysis of association between tissue-specific brain volume changes and age (β = −0.019 to −0.011; adjusted p = 0.028–1.00) or disease duration (β = −0.028 to −0.008; adjusted p = 0.16–0.96) confirmed these results. Although increase in the relapse rate (β = 0.10; adjusted p = 9 × 10−9), Expanded Disability Status Scale (EDSS; β = 0.17; adjusted p = 8 × 10−5), and EDSS change (β = 0.15; adjusted p = 2 × 10−5) were associated with accelerated rate of BVL, their effect on the rate of BVL was minimal (all Cohen f2 ≤ 0.007). In 94 patients who escalated therapy, the rate of BVL decreased following treatment escalation by 0.29% (β = −0.29; Cohen f2 = 0.133; p = 5.5 × 10−8).
Conclusions The rate of BVL is relatively stable throughout the course of RRMS. The accelerated BVL is weakly associated with concurrent higher disease activity, and timely escalation to high-efficacy immunotherapy helps decrease the rate of BVL.
Glossary
- ASA=
- Avonex-Steroids-Azathioprine;
- BVL=
- brain volume loss;
- CIS=
- clinically isolated syndrome;
- DMT=
- disease-modifying treatment;
- EDSS=
- Expanded Disability Status Scale;
- MS=
- multiple sclerosis;
- QMRI=
- Quantitative MRI;
- RRMS=
- relapsing-remitting MS;
- SET=
- Study of Early IFN-β1a Treatment;
- SPMS=
- secondary progressive MS
Accelerated whole-brain volume loss (BVL) in patients with multiple sclerosis (MS) is associated with an increased risk of neurologic and cognitive impairment. Assessment of BVL may help to evaluate early treatment response and improve the stratification of MS patients who are at highest risk of disease progression.1,–,5 Although BVL is clinically relevant, our understanding of BVL dynamics throughout the disease course is limited. Previous studies have reported contradictory findings. Although some studies have suggested constant BVL rate throughout the MS disease course,6,7 other studies have reported accelerated BVL rate in early8 or advanced2,9 disease stages.
It is widely accepted that patients reaching secondary progressive phase have in general greater disease severity than age-matched patients with clinically isolated syndrome (CIS) or relapsing-remitting MS (RRMS). However, the previous studies investigating BVL dynamics in MS compared BVL rates between CIS or RRMS and secondary progressive MS (SPMS). Comparison of patients with different MS phenotypes but also with different disease severity may introduce bias into longitudinal evaluation of BVL trajectory. Only 1 study examined evolution of BVL rate on a continuous scale, but was limited by having a relatively small number of time points over a limited follow-up duration.8
A thorough understanding of the average BVL trajectory throughout the course of MS is important for balanced recruitment of patients into clinical studies where BVL is used as an outcome measure. Furthermore, an expectation of higher rate of BVL in certain stages of MS would increase clinical awareness during such periods and might lead to timely intensified therapeutic intervention in patients exhibiting higher disease activity, on the proviso that such therapies can mitigate the rate of BVL. Given that the frequency of relapses is highest at young age and during early disease,10,–,12 and that highly inflammatory MS phenotype is associated with faster BVL,13 we hypothesized that BVL is accelerated in early stages of MS. In this study, we investigated the dynamics of BVL and tissue-specific brain volume changes in patient with different ages and disease durations. Given that previous studies showed a spatiotemporal pattern of gray and white matter volume changes during aging14,15 and MS progression,2,8,14 we also investigated an evolution of tissue-specific brain volume changes through disease course, such as total, cortical and deep gray matter, white matter, thalamus, and lateral ventricle. The secondary aims were to describe association between BVL and clinical measures and to investigate an effect of treatment change on the rate of BVL.
Methods
Patients
We analyzed a combined cohort of adult (aged >18 years), predominantly patients with RRMS, who participated in the Avonex-Steroids-Azathioprine (ASA) trial (N = 166), Study of Early IFN-β1a Treatment (SET) (N = 180), and in a cohort of patients from the quantitative MRI (QMRI) project (N = 1557).
ASA Cohort
The original ASA study was a 2-year randomized, double-blind, placebo-controlled trial with a 5-year open-label extension, investigating clinical and MRI outcomes among patients with RRMS treated with IFNß-1a (enrollment from April 1999 to December 2003).16
SET Cohort
The original SET study was a 2-year investigator initiated, prospective, observational clinical study of IFNß-1a treatment in patients after the first demyelinating event suggestive of MS involving 8 centers (enrollment from October 2005 to June 2009; clin.gov #NCT01592474).17
QMRI Cohort
All patients with MS routinely followed at the MS Centre at the General University Hospital in Prague, and with brain MRI from March 2000 to November 2015, were included in the QMRI program.18
Clinical Evaluation
After the SET and ASA were completed, patients' standard clinical follow-up, consisting of 3- to 6-monthly visits, was aligned with the QMRI cohort. In this study, we included clinical data collected before November 2015. Confirmed disability worsening was defined as worsening by an increase in the Expanded Disability Status Scale (EDSS) of 1.5 points if the baseline EDSS score was 0; 1.0 point if the baseline EDSS score was between 1.0 and 5.0; and 0.5 point for baseline EDSS scores of 5.5 or higher. We required confirmation of disability worsening after 12 months.19 Additional details about study cohorts and clinical evaluations are provided in supplemental material 1, links.lww.com/NXI/A448.
Selection Criteria
We excluded MRI scans performed: within 1 month after clinical relapse or treatment with high-dose steroids, within 10 months (MRI monitoring every 12 ± 2 months) after previous MRI scan, in pregnancy, or within 12 months after delivery. In addition, MRI scans performed within first 12 months after disease-modifying treatment (DMT) initiation or DMT escalation were excluded to prevent confounding effect of pseudo-atrophy.20,21 DMT escalation was defined as switching from low-efficacy DMT (dimethyl fumarate, glatiramer acetate, interferon β, and teriflunomide) to high/moderate-efficacy DMT (fingolimod, mitoxantrone, and natalizumab). Finally, we included only patients with ≥1-year (±2 months) follow-up and ≥2 MRI scans.
Ethics Statement
The study protocols were approved by the medical ethics committees, and all patients gave their written informed consent.
MRI Acquisition and Analysis
All MRI scans analyzed in the study were performed on the same scanner (1.5-T Gyroscan; Philips Medical Systems, Best, the Netherlands) that did not undergo a major hardware upgrade. The standardized protocol consisted of 2 sequences: fluid-attenuated inversion recovery and T1-weighted 3-dimensional turbo field echo (T1-WI/TFE 3D). Volumetric assessment of T2 lesions and brain volumes was performed using semi-automated ScanView software. Lesion filling on T1-WI 3D images was used to reduce the impact of T1 hypointensities on tissue segmentation.22 For longitudinal changes of the BVL, we applied the SIENA technique.23 In MRI scans from the ASA and SET cohorts, we performed validation analysis of tissue-specific brain volumes by FMRIB's Integrated Registration and Segmentation Tool24 and SIENAX multi-time-point.25 Details of MRI acquisition and analysis are described in supplemental material 2, links.lww.com/NXI/A449.
Statistical Analysis
All analyses were performed using the R statistical software (R-project.org). We analyzed annualized % volume changes (BVL rate) relative to the preceding MRI examination.
The Relationship Between Age or Disease Duration as Continuous Variables and BVL Rate
The relationship between age or disease duration (independent variables) and annualized BVL rate (dependent variable) was analyzed using linear mixed models with a random intercept specified for each patient. We estimated effect size (Cohen f2) for the effect of age and disease duration on the BVL rate. Cohen f2 ≥ 0.02, f2 ≥ 0.15, and f2 ≥ 0.35 represent small, medium, and large effect sizes, respectively.26 Regression coefficients (β) for the age or disease duration of the models were interpreted as average change of the BVL rate for each 1-year increase in patient age or disease duration.
Mixed models with an interaction term were used to investigate additive effect of age at time point and disease duration at baseline (or vice versa) on evolution of MRI measures. Variance inflation factors were used to examine multicollinearity among independent variables.
The Relationship Between Age or Disease Duration as Categorical Variables and BVL Rate
For each patient, we calculated the mean annual BVL rate at different available ages and disease durations (figure 1) and compared them by mixed models. For each patient, we also fitted individual slopes (e.g., calculated regression coefficients) of the BVL rates at different available age and disease duration intervals using linear regression. These fitted individual BVL slopes (e.g., regression coefficients) were compared by mixed models to investigate differences in BVL trajectories throughout disease course.
Relationship between annualized whole-brain volume loss rate and age (A) and disease duration as continuous variables (B) with a smooth curve fitted by loess regression technique; relationship between mean annualized whole-brain volume loss rates at different ages (C) and disease duration intervals (D); and relationship between mean fitted individual slopes (regression coefficients) of annualized whole-brain volume loss rate trajectories at different ages (E) and disease duration intervals (F). The dotted black horizontal line represents no change in annualized whole-brain volume loss rate, and the red dotted horizontal line represents the reference cutoff value (0.4% per year) for pathologic whole-brain volume change.
The Relationship Between Clinical Disease Activity and the Rate of BVL
We then analyzed the relationships of annualized relapse rate, EDSS score, and annualized EDSS change (all independent log+1 transformed variables) with BVL rates using mixed models. Mixed models were also used to investigate interactions between age (or disease duration) and measures of clinical disease activity and their associations with MRI measures.
In addition, we analyzed associations between the BVL during the initial 2 years of follow-up and time to confirmed disability worsening using adjusted Cox proportional hazard models. The cutoff point for BVL rate was based on the third quartile. Association between annualized EDSS change over the entire follow-up and age or disease duration at baseline was analyzed by multivariable linear regression.
The Effect of Treatment Change on BVL Rate
Finally, in a subgroup of patients who had completed at least 2 MRI scans before and 2 MRI scans after DMT escalation, we compared BVL rates before and after DMT escalation using mixed models.
Sensitivity Analyses
Sensitivity analyses were completed (a) in a subgroup of patients treated with only low-efficacy DMTs to eliminate the effect of treatment change and (b) in trial patients from the ASA and the SET cohorts to analyze association between age and disease duration with tissue-specific brain % volume changes.
Adjustment of Multivariable Models
Unless otherwise indicated, multivariable mixed and Cox proportional hazard models were adjusted for sex, study cohort, T2 lesion volume at time point, normalized brain volume (brain parenchymal fraction) at baseline, proportion of time spent on low-efficacy DMTs, and proportion of time spent on high/moderate-efficacy DMTs.
False Discovery Rate Control
The Holm procedure with q = 0.05 was used to control false discovery rate of dependent tests. The number in [brackets] following adjusted p refers to the number of comparisons.
Data Availability
Anonymized data not published within this article will be made available by reasonable request from any qualified investigator for purposes of replicating study results.
Results
Baseline and Follow-up Clinical Characteristics
Of the 3,415 patients (20,101 MRI scans), 1,903 patients (7,203 MRI scans) fulfilled the inclusion criteria and were included in this study. Age, disease duration, and EDSS at baseline differed among the cohorts (small to median effect sizes). The mean rate of BVL was 0.48% ± 0.93% among all patients; 0.45% ± 0.94% among patients without DMTs; 0.50% ± 0.94% among patients treated with low-efficacy DMTs; and 0.43% ± 0.82% among patients treated with high/moderate-efficacy DMTs. Table 1 describes the demographic, clinical, and MRI characteristics of the sample.
Sample Characteristics
Brain Volume Loss and Age
When we explored age as continuous variable, we did not find a significant association between younger age and higher BVL rate (β = −0.003; 95% CI = −0.006 to 0.0002; Cohen f2 = 0.0005; p = 0.065; adjusted p[7] = 0.39) (table 2 and figure 1A). Similar results were seen when different DMTs groups were analyzed separately (table 2). To eliminate an effect of disease duration on this association, we analyzed patients with disease duration at baseline, of 1, 2, or 3 years, but again we did not find an effect of age on BVL rate (all Cohen f2 < 0.001) (table 2).
Association Between Age (or Disease Duration) and Rates of Whole-Brain Volume Loss in Subgroups of Patients Defined by Their Age/Disease Duration at Baseline or by Their DMT
When age was analyzed as categorical variable, we found similar BVL rates in different age groups (table 3).
Rates of Whole-Brain Volume Loss in Subgroups of Patients Defined by Their Age and Disease Duration
We did not find differences in means of linear regression coefficients of BVL rates among age intervals (β = 0.047; 95% CI = −0.06 to 0.15; Cohen f2 < 0.001; p = 0.40) (figure 1E).
Brain Volume Loss and Disease Duration
Although BVL rate tended to decrease with disease duration, the effect of disease duration on BVL rate was minimal (β = −0.012; 95% CI = −0.017 to −0.007; Cohen f2 = 0.004; p = 6 × 10−6; adjusted p[7] = 4 × 10−5) (table 2 and figure 1B). Among patients on low-efficacy DMTs, each year of disease duration was associated with the mean reduction of 0.009 BVL rate (β = −0.009; 95% CI = −0.016 to −0.003; Cohen f2 = 0.002; p = 0.005; adjusted p[7] = 0.024). To eliminate a potential confounding effect of age on this association, we separately analyzed patients with baseline age of 25–30 years, 30–35 years, or 35–40 years, but we again did not find any effect of disease duration on BVL rate (all Cohen f2 < 0.007) (table 2).
When age analyzed as categorical variable, we found similar BVL rates across different disease duration intervals (table 3).
We did not find differences in means of linear regression coefficients of the BVL rates among disease duration interval subgroups (β = 0.013; 95% CI = −0.03 to 0.06; Cohen f2 < 0.001; p = 0.59) (figure 1F).
We did not find differential effect of disease duration across different ages on BVL rate (interaction between age at time point and disease duration at baseline: β = −0.00002; p = 0.41; adjusted p[2] = 0.41; interaction between age at baseline and disease duration at time point: β = −0.000001; p = 0.16; adjusted p[2] = 0.32).
Association Between Brain Volume Loss and Relapse Activity
The relapse rate (e.g., number of relapses per year) declined with increasing age and disease duration of the studied patients (figure e-1, links.lww.com/NXI/A443). Increase in relapse rate by 1.0 was associated with the mean BVL rate increase by 0.10% (β = 0.10; 95% CI = 0.07 to 0.13; Cohen f2 = 0.007; p = 3 × 10−9; adjusted p[3] = 9 × 10−9) (figure e-2, links.lww.com/NXI/A444). Younger age was associated with a stronger association between relapse rate and BVL rate (interaction between age and relapse rate: β = 0.012; p = 0.008; adjusted p[3] = 0.016). Disease duration did not further modify the association between relapse rate and BVL rate (interaction between disease duration and relapse rate: β = 0.003; p = 0.61; adjusted p[3] = 0.61).
Association Between Brain Volume Loss and Disability
The BVL rate was associated with the EDSS score (β = 0.17; 95% CI = 0.09 to 0.24; Cohen f2 = 0.003; p = 2 × 10−5; adjusted p[5] = 8 × 10−5) and annualized EDSS absolute change (β = 0.15; 95% CI = 0.09–0.21; Cohen f2 = 0.006; p = 3 × 10−6; adjusted p[5] = 2 × 10−5). The association between BVL rate and annualized EDSS absolute change was stronger among older patients (interaction between age and BVL rate: β = 0.003; p = 0.00015; adjusted p[5] = 0.0005). Disease duration did not further modify the association between BVL rate and annualized EDSS absolute change (interaction between disease duration and BVL rate: β = 0.001; p = 0.28; adjusted p[5] = 0.28).
The mean BVL rate during the first 2 years of the follow-up (cutoff < −0.87%; third quartile) was associated with a higher risk of confirmed disability worsening over the long term (OR = 1.4; 95% CI = 1.1–1.7; p = 0.005; adjusted p[5] = 0.010) (figure e-3, links.lww.com/NXI/A445). Annualized absolute change in EDSS did not differ among patients with different ages and disease durations (both p = 0.38–0.41) (e-figure 4, links.lww.com/NXI/A446).
Association Between Brain Volume Loss and Lesion Burden
An average increase in T2 lesion volume was 6.5% per year (95% CI = 6.0–7.1; p < 1 × 10−10; adjusted p[8] < 1 × 10−10). We found a weak association between the higher rate of T2 lesion volume accumulation (mL) and younger age or shorter disease duration. On average, the rate of T2 lesion volume increase decelerated by 0.006 mL for every year of age (β = −0.006; 95% CI = −0.010 to −0.002; p = 0.002; adjusted p[8] = 0.008) and by 0.014 mL for every year of disease duration (β = −0.014; 95% CI = −0.02 to −0.008; p = 4 × 10−6; adjusted p[8] = 2 × 10−5; figure 2).
Mean annualized T2 lesion volume absolute changes in different age (A) and disease duration (B) intervals. Statistical analysis was performed on continuous data.
The higher rate of T2 lesion volume accumulation (mL) was associated with a higher BVL rate (β = −0.15; 95% CI = −0.17 to −0.12; Cohen f2 = 0.028; p < 1 × 10−10; adjusted p[8] < 1 × 10−10), higher annualized change in EDSS (β = 0.02; 95% CI = 0.01–0.03; Cohen f2 = 0.003; p = 0.0002; adjusted p[8] = 0.0008), and higher relapse rate (β = 0.20; 95% CI = 0.05–0.34; Cohen f2 = 0.001; p = 0.007; adjusted p[8] = 0.017). The strength of the association between the rate of accumulation of T2 lesion volume and BVL rate was similar among patients of different ages and disease durations (interaction between age and T2 lesion volume accumulation rate: β = 0.0003; p = 0.83; adjusted p[8] = 0.83; interaction between disease duration and T2 lesion volume accumulation rate: β = 0.004; p = 0.027; adjusted p[8] = 0.054).
Tissue-Specific Brain Volume Analysis
The sensitivity analysis that included the 348 trial patients (1,569 scans) from the ASA and the SET cohorts confirmed the absence of relationship between age and BVL rate (β = −0.001; p = 0.75; adjusted p[8] = 1.00). Similarly, the rate of segmented volume loss was constant across patient ages: white matter (β = −0.011; p = 0.014; adjusted p[8] = 0.70), total gray matter (β = −0.016; p = 0.011; adjusted p[8] = 0.066), deep gray matter (β = −0.019; p = 0.046; adjusted p[8] = 0.18), thalamic (β = −0.017; p = 0.10; adjusted p[8] = 0.30), and lateral ventricle volume changes (β = 0.015; p = 0.56; adjusted p[8] = 1.00; in all Cohen f2 <0.005). Only on 1 occasion, a weak association between age and cortical gray matter % volume loss was seen (β = −0.019; p = 0.004; Cohen f2 = 0.006; adjusted p[8] = 0.028). This association was confirmed in analysis of patients treated only with low-efficacy DMTs (β = −0.021; Cohen f2 = 0.008; p = 0.002; adjusted p[8] = 0.016).
The magnitude of the association between disease duration and the BVL rate was similar to the primary analysis, but did not reach the level of formal statistical significance due to smaller sample size (β = −0.009; p = 0.33; adjusted p[7] = 0.96). Similarly, we did not find any evidence for associations of disease duration with segmented brain volumes: white matter (β = −0.016; p = 0.32; adjusted p[7] = 0.96), total gray matter (β = −0.024; p = 0.094; adjusted p[7] = 0.38), cortical gray matter (β = −0.028; p = 0.045; adjusted p[7] = 0.27), deep gray matter (β = −0.014; p = 0.51; adjusted p[7] = 0.27), thalamic (β = −0.013; p = 0.023; adjusted p[7] = 0.16), and lateral ventricle volume changes (β = 0.008; p = 0.73; adjusted p[7] = 0.96; in all Cohen f2 <0.005) (figure e-5, links.lww.com/NXI/A447).
An Effect of Disease-Modifying Treatment on Brain Volume Loss
In a subgroup of 94 patients who escalated DMTs (609 MRI scans), the mean BVL rates substantially decreased following treatment escalation (before: −0.45; 95% CI = −0.54 to −0.37 vs after: −0.10; 95% CI = −0.13 to −0.07). Differences in BVL rates following treatment escalation were confirmed in adjusted mixed models, where treatment escalation resulted in a mean reduction of the BVL rate by 0.29% (β = −0.29; 95% CI = −0.40 to −0.19; Cohen f2 = 0.133; p = 5.5 × 10−8) (figure 3). An average time between the cessation of low-efficacy DMT and the initiation of high/moderate-efficacy DMT was 1.9 months (range 0–11.4 months).
Evolution of annualized whole-brain volume loss after escalation from the low-efficacy on high/moderate-efficacy disease-modifying treatments. The treatment change occurred at time 0.
Discussion
In this large longitudinal study of patients with RRMS, the rate of BVL declined with longer disease duration but not with older age. Patients with shorter disease duration showed marginally higher rates of BVL than patients with longer disease history; however, the effect of disease duration on the BVL rate was minimal. Similarly, higher rate of BVL was weakly associated with the higher frequency of recent relapses and greater disability, although the magnitude of these associations was relatively small. On the other hand, a stronger association was found between the rate of BVL and concurrent enlargement of cerebral T2 lesions or time to future confirmed disability worsening. Importantly, BVL was ameliorated by escalation of MS therapy. We did not find any clinically relevant associations between time and the rate of tissue-specific BVL.
In a study of 140 patients with MS, Andorra et al.8 showed that the rate of BVL is accelerated during the first 5 years after disease onset, followed by its subsequent stabilisation. Similarly, we have demonstrated that the rate of BVL is slightly higher in patients earlier after disease onset. Each year of disease duration was associated with a mean reduction in BVL rate of 0.013% (or 0.009% in patients treated exclusively with low-efficacy DMTs). However, the decline in the BVL rate was constant throughout the disease course.
Another large study of 963 untreated patients (1,926 time points) reported that the BVL rate was independent of MS subtype, age, or disease duration, when adjusting for baseline brain volume.6 A recent study of 1,417 patients (3,604 time points) showed accelerated cortical gray matter atrophy in SPMS compared with RRMS.2 These observations are complementary. First, in SPMS, the disease is more progressed and patients reaching secondary progressive phase have in general greater disease severity than age-matched patients with CIS or RRMS. Second, SPMS is usually associated with older age and longer disease history, and other determinants such as aging15 and comorbidities27 start to contribute to BVL.
The rates of BVL found in our study were comparable with those reported in previous studies that included treated MS patients,6,8,14,28 In the present study, we focused on relatively young cohorts of patients with short disease duration and with relatively more active disease from its onset. In keeping with a published study,13 our work showed that the frequency of relapses is associated with the rate of BVL; however, the magnitude of the association between recent relapse activity and the rate of BVL was marginal. The rate of recent increase in the cerebral T2 lesion burden was moderately associated with the rate of BVL.13,29,30 Of interest, the weak association between relapse frequency and the BVL rate was relatively stronger in younger patients but not specifically during early disease. These observations may suggest that the rate of BVL is closely related to the episodes of neuroinflammation that underlie formation of new or reactivation of existing cerebral lesions that presents clinically as MS relapses. This is supported by the demonstrated strong relationship between acute lesional inflammation and early axonal disruption as well as subsequent degeneration of chronically demyelinated axons resulting in brain tissue loss.31,–,34 In addition, we found a relationship between early BVL and future disability progression over long-term follow-up, which is in agreement with previous research and supports clinical relevance of BVL.1,–,5,9,13,17 As expected, we found a somewhat stronger relationship between BVL and disability worsening in older patients.35 Taken together, these observations may suggest that although formation of new lesions is a strong driver of BVL among younger patients, the loss of brain tissue in turn determines impaired neurologic function at older age. Given that the effect of age on the association between the rate of BVL and relapse frequency or disability worsening was only minimal, clinical relevance of these findings requires further validation.
Establishing the definition of BVL that is considered normal in the context of normal aging is a key step toward clinical interpretation of BVL in MS. The constant rate of BVL across different age groups of patients with MS simplifies the clinical use of this metric. However, the drivers of BVL may differ over time. In elderly patients, the aging process and comorbid diseases are likely to contribute more to the overall BVL than in younger patients, and therefore, a faster BVL in older patients may not always represent accelerated brain atrophy due specifically to MS.14,36 On the contrary, in early MS, higher BVL rates are more likely to indicate high level of disease activity requiring more aggressive therapy. Considering the relationship between early BVL and future development of disability,2,4,9,13 a neurologist should be particularly vigilant in screening for accelerated BVL in patients recently diagnosed with MS. This has become of particular importance as the new and emerging DMTs hold the promise of decelerating the rate of BVL and the loss of neurologic capacity.37,38 Our present study has also shown that escalation of DMT can help normalize the increased rate of BVL.28,39
Limitations
Our study included mostly patients treated with DMTs; therefore, the effect of pseudo-atrophy may be present.20,21 To minimize this, we analyzed only MRI scans performed at least 12 months after DMT initiation. Importantly, to minimize a potential effect of DMT escalation, we replicated our findings in a subgroup of patients who were treated exclusively with low-efficacy DMTs.
Moreover, the majority of patients started MRI monitoring at the time of disease onset or at presentation with relapses; therefore, the regression-to-the-mean phenomenon might play a role in detecting marginally increased rates of BVL in early disease stages. In addition, patients with RRMS who remained untreated over the long term are probably those with a favorable disease phenotype. On the other hand, untreated patients during early disease stages represent a more heterogeneous population.
It would be expected that the rate of BVL becomes accelerated in older age.2,9,14,36 In practice, older patients with more severe disability and during the transition to SPMS often receive less frequent MRI monitoring due to limited therapeutic options. Such patients might be slightly underrepresented in this study and explain the absence of accelerated rate of BVL in older patients.
Although BVL rates showed in the study are in agreement with previous research,6,8,14,28 the differences in volumetric outcomes are often driven by differences in MRI acquisition and volumetric technique. Therefore, absolute changes of BVL rates should not be directly compared with the results from other studies, in particular where different methodology was used.
In the available cohorts, almost 70% of time points were recorded during treatment with low-efficacy DMTs. Our results need to be replicated among patients treated with more contemporary, high/moderate-efficacy DMTs.
Finally, the study did not include patients with progressive MS, which may limit generalizability of our finding to relapsing MS only. We did not study healthy controls and therefore cannot comment on the contribution of age vs MS to the total observed BVL.14,36
In this longitudinal study of a large sample of patients with RRMS, we did not find a clinically relevant relationship between BVL rate and age or disease duration. Accelerated BVL was weakly associated with concurrent higher disease activity and eventually led to more profound neurologic disability. Evidence of higher rate of BVL should prompt the consideration of escalation of MS therapy to prevent neurologic disability.
Study Funding
Study funded by Sanofi (GZ-2017-11718), the Czech Ministry of Education project Progress Q27/LF1, and Czech Ministry of Health project RVO-VFN64165 and NV18-04-00168.
Disclosure
The authors of the study have been consultants or received support for conference travel or research from pharmaceutical companies producing immune-modulatory drugs for escalation treatment in MS. T. Uher received financial support for conference travel and speaker honoraria from Biogen Idec, Novartis, Roche, Genzyme, and Merck Serono and support for research activities from Biogen Idec and Sanofi (GZ-2017-11718). J. Krasensky received financial support for conference travel and research activities from Biogen Idec. C.B. Malpas and N. Bergsland have nothing to disclose. M. Dwyer received personal compensation from Claret Medical and Keystone Heart Ltd. and research grant support from Novartis, Bristol Myers Squibb, and Mapi Pharma. E.K. Havrdova received speaker honoraria and consultant fees from Biogen Idec, Merck Serono, Novartis, Genzyme, Teva, Actelion, and Receptos and support for research activities from Biogen Idec and Merck Serono. M. Vaneckova received compensation for travel, speaker honoraria and consultant fees from Biogen Idec, Novartis, Merck Serono, Genzyme, Roche and Teva, and support for research activities from Biogen Idec. D. Horakova received compensation for travel, speaker honoraria and consultant fees from Biogen Idec, Novartis, Merck Serono, Bayer Shering, and Teva, and support for research activities from Biogen Idec. R. Zivadinov received personal compensation from EMD Serono, Sanofi, Novartis, Bristol-Myers Squibb, and Keystone Heart for speaking and consultant fees. He received financial support for research activities from Sanofi, Novartis, Bristol-Myers Squibb, Mapi Pharma, Protembis, V-VAWE Medical, and Keystone Heart. Dr. Zivadinov serves on the editorial board of J Alzh Dis, BMC Med, BMC Neurol, Vein and Lymphatics, and Clinical CNS Drugs. Tomas Kalincik served on scientific advisory boards for Roche, Genzyme-Sanofi, Novartis, Merck, and Biogen, steering committee for Brain Atrophy Initiative by Genzyme, received conference travel support and/or speaker honoraria from WebMD Global, Novartis, Biogen, Genzyme-Sanofi, Teva, BioCSL, and Merck, and received research support from Biogen. Go to Neurology.org/NN for full disclosures.
Acknowledgment
The authors thank statistician Vaclav Capek, MS, PhD, from the Department of Neurology and Center of Clinical Neuroscience, Charles University in Prague for help with developing a data analysis plan, review of statistics, and help with statistical analysis and interpretation. The authors also thank patients who participated in this study.
Appendix Authors

Footnotes
Go to Neurology.org/NN for full disclosures. Funding information is provided at the end of the article.
The Article Processing Charge was funded by Sanofi grant GZ-2017-11718.
- Received January 15, 2020.
- Accepted in final form January 5, 2021.
- Copyright © 2021 The Author(s). Published by Wolters Kluwer Health, Inc. on behalf of the American Academy of Neurology.
This is an open access article distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND), which permits downloading and sharing the work provided it is properly cited. The work cannot be changed in any way or used commercially without permission from the journal.
References
- 1.↵
- 2.↵
- 3.↵
- Azevedo CJ,
- Pelletier D
- 4.↵
- 5.↵
- 6.↵
- 7.↵
- 8.↵
- Andorra M,
- Nakamura K,
- Lampert EJ, et al.
- 9.↵
- 10.↵
- 11.↵
- 12.↵
- 13.↵
- 14.↵
- Azevedo CJ,
- Cen SY,
- Jaberzadeh A,
- Zheng L,
- Hauser SL,
- Pelletier D
- 15.↵
- 16.↵
- 17.↵
- 18.↵
- Uher T,
- Vaneckova M,
- Krasensky J, et al.
- 19.↵
- 20.↵
- 21.↵
- 22.↵
- 23.↵
- 24.↵
- 25.↵
- 26.↵
- 27.↵
- 28.↵
- Arnold DL,
- Fisher E,
- Brinar VV, et al.
- 29.↵
- Varosanec M,
- Uher T,
- Horakova D, et al.
- 30.↵
- Chard DT,
- Brex PA,
- Ciccarelli O, et al.
- 31.↵
- 32.↵
- 33.↵
- 34.↵
- 35.↵
- 36.↵
- Battaglini M,
- Gentile G,
- Luchetti L, et al.
- 37.↵
- He A,
- Merkel B,
- Brown JWL, et al.
- 38.↵
- Brown JWL,
- Coles A,
- Horakova D, et al.
- 39.↵
- Gaetano L,
- Haring DA,
- Radue EW, et al.
Letters: Rapid online correspondence
REQUIREMENTS
You must ensure that your Disclosures have been updated within the previous six months. Please go to our Submission Site to add or update your Disclosure information.
Your co-authors must send a completed Publishing Agreement Form to Neurology Staff (not necessary for the lead/corresponding author as the form below will suffice) before you upload your comment.
If you are responding to a comment that was written about an article you originally authored:
You (and co-authors) do not need to fill out forms or check disclosures as author forms are still valid
and apply to letter.
Submission specifications:
- Submissions must be < 200 words with < 5 references. Reference 1 must be the article on which you are commenting.
- Submissions should not have more than 5 authors. (Exception: original author replies can include all original authors of the article)
- Submit only on articles published within 6 months of issue date.
- Do not be redundant. Read any comments already posted on the article prior to submission.
- Submitted comments are subject to editing and editor review prior to posting.
You May Also be Interested in
Hastening the Diagnosis of Amyotrophic Lateral Sclerosis
Dr. Brian Callaghan and Dr. Kellen Quigg
► Watch
Related Articles
- No related articles found.
Topics Discussed
Alert Me
Recommended articles
-
Article
Correlation between brain volume loss and clinical and MRI outcomes in multiple sclerosisErnst-Wilhelm Radue, Frederik Barkhof, Ludwig Kappos et al.Neurology, January 28, 2015 -
Articles
A reassessment of the plateauing relationship between T2 lesion load and disability in MSMaria Pia Sormani, Marco Rovaris, Giancarlo Comi et al.Neurology, September 30, 2009 -
Article
Teriflunomide slows BVL in relapsing MSA reanalysis of the TEMSO MRI data set using SIENAErnst-Wilhelm Radue, Till Sprenger, Laura Gaetano et al.Neurology: Neuroimmunology & Neuroinflammation, August 09, 2017 -
Article
Spinal cord volume lossA marker of disease progression in multiple sclerosisCharidimos Tsagkas, Stefano Magon, Laura Gaetano et al.Neurology, June 27, 2018