Original Research ARTICLE
Selective Brain Network and Cellular Responses Upon Dimethyl Fumarate Immunomodulation in Multiple Sclerosis
- 1Department of Neurology, Focus Program Translational Neuroscience (FTN), Rhine-Main Neuroscience Network (rmn2), University Medical Center of the Johannes Gutenberg University Mainz, Mainz, Germany
- 2Department of Neurology, Institute of Emergency Medicine, Chisinau, Moldova
- 3Laboratory of Neurobiology and Medical Genetics, Nicolae Testemiţanu State University of Medicine and Pharmacy, Chisinau, Moldova
- 4Department of Neuroradiology, University of Kiel, Kiel, Germany
- 5Department of Neurology With Institute of Translational Neurology, University of Münster, Münster, Germany
Background: Efficient personalized therapy paradigms are needed to modify the disease course and halt gray (GM) and white matter (WM) damage in patients with multiple sclerosis (MS). Presently, promising disease-modifying drugs show impressive efficiency, however, tailored markers of therapy responses are required. Here, we aimed to detect in a real-world setting patients with a more favorable brain network response and immune cell dynamics upon dimethyl fumarate (DMF) treatment.
Methods: In a cohort of 78 MS patients we identified two thoroughly matched groups, based on age, disease duration, disability status and lesion volume, receiving DMF (n = 42) and NAT (n = 36) and followed them over 16 months. The rate of cortical atrophy and deep GM volumes were quantified. GM and WM network responses were characterized by brain modularization as a marker of regional and global structural alterations. In the DMF group, lymphocyte subsets were analyzed by flow cytometry and related to clinical and MRI parameters.
Results: Sixty percent (25 patients) of the DMF and 36% (13 patients) of the NAT group had disease activity during the study period. The rate of cortical atrophy was higher in the DMF group (−2.4%) compared to NAT (−2.1%, p < 0.05) group. GM and WM network dynamics presented increased modularization in both groups. When dividing the DMF-treated cohort into patients free of disease activity (n = 17, DMFR) and patients with disease activity (n = 25, DMFNR) these groups differed significantly in CD8+ cell depletion counts (DMFR: 197.7 ± 97.1/μl; DMFNR: 298.4 ± 190.6/μl, p = 0.03) and also in cortical atrophy (DMFR: −1.7%; DMFNR: −3.2%, p = 0.01). DMFR presented reduced longitudinal GM and WM modularization and less atrophy as markers of preserved structural global network integrity in comparison to DMFNR and even NAT patients.
Conclusions: NAT treatment contributes to a reduced rate of cortical atrophy compared to DMF therapy. However, patients under DMF treatment with a stronger CD8+ T cell depletion present a more favorable response in terms of cortical integrity and GM and WM network responses. Our findings may serve as basis for the development of personalized treatment paradigms.
Loss of structural integrity driven by inflammation, demyelination, and degeneration in multiple sclerosis (MS) involves white matter (WM) and gray matter (GM) compartments, the latter playing a key role in disability and disease progression (1–4). Existing data suggest that CD4+ and CD8+ T cells contribute to the damage of the cortical GM (5, 6). Evaluation of magnetic resonance imaging (MRI)-derived parameters of GM structural alterations was incorporated into studies to track the responses to disease-modifying drugs (DMDs) (7) which have been proven to slow, to various extents, the rate of GM tissue loss. We have recently shown that advanced measures of brain network architecture closely mirror the disease course and clinical impairment (8, 9). Only an exact longitudinal quantification of local and global GM and WM tissue properties enables the development of precise disease course models, thereby creating the basis for personalized therapeutic decisions.
Presently, for relapsing-remitting MS (RRMS) several promising DMDs are available but the long-term benefits of the therapeutic algorithms are still unclear. No unambiguous personalized solutions to halt or ideally reverse the disease course exist; however, first avenues for very efficient immunomodulatory remedies arise. At the same time, markers predicting the favorable response to a specific DMD are under meticulous development but are not yet validated for clinical pathways. Mainly, MRI parameters are regarded as surrogate measures of treatment response to DMDs, although cerebrospinal fluid (CSF) or peripheral blood immune response may also be a valuable biomarker (10). In this respect, treatment response to dimethyl fumarate (DMF) was reflected by reduced counts of CD4+ and CD8+ T cells in patients without disease activity (11–13). Clinically and radiologically stable patients under DMF therapy showed a more pronounced CD8+ than CD4+ T cell reduction compared to active patients (11, 13).
Here, we have postulated that the clinical and brain structural response is tightly linked to the immune cell dynamics under DMF treatment that could be reliably monitored by analyzing peripheral blood lymphocyte subsets. Thus, we aimed to recognize patients with no disease activity and structural deterioration as mirrored by both cortical and subcortical integrity and brain network changes and relate these to immune cell dynamics. A second cohort of patients treated with natalizumab (NAT) served as a reference group to compare clinical, structural and network responses. To this end, we computed regional rates of cortical atrophy, constructed structural GM (from cortical thickness) and WM [from probabilistic tractography (PT)] networks and correlated the atrophy rates with longitudinal changes in lymphocyte subsets.
Materials and Methods
In this longitudinal study, 78 patients (mean age ± standard deviation (SD) 32.7 ± 8.7 years; 28 males; mean disease duration of 51.1 ± 37.8 months) were selected out of 1,156 patients recruited at the Department of Neurology at the University Medical Center of the Johannes Gutenberg University Mainz in Germany according to the following inclusion criteria: (1) patients aged 18–60 years, (2) patients diagnosed with RRMS, (3) starting DMF or NAT treatment, (4) scanned with a standardized MRI protocol (14), (5) serially acquired MRI scans at several time points, (6) no corticosteroid use within 30 days prior to MRI, (7) peripheral blood samples available at baseline and follow-up time points for DMF-treated patients. Exclusion criteria: (1) necessity in treatment escalation; (2) participation in any interventional trial during the study period; (3) serious adverse events requiring premature study termination; (4) patients with primary or secondary MS progression. Forty-two patients (34.5 ± 9.0 years; 14 males) were identified on DMF treatment (DMF group) and 36 patients (30.6 ± 8.1 years; 14 males) on NAT treatment (NAT group). All patients fulfilled the revised 2010 McDonald diagnostic criteria for RRMS (15). Prior to DMF or NAT therapy patients were exposed to interferon β-1 β, interferon β-1 α or glatiramer acetate. Patients' characteristics are included in Table 1.
Expanded Disability Status Scale (EDSS) score assessment was obtained at treatment initiation and later at 3-month intervals. MRI data were acquired at the time (1.9 ± 1.7 months) of treatment (DMF/NAT) onset, 6 months after the treatment and then on an annual basis. To track more robust structural and clinical changes, the longest follow-up MRI and EDSS performed after 16 months (15.8 ± 7.2 in the DMF group, 16.1 ± 4.3 in the NAT group) were considered for the analysis. Disease activity was evaluated during the entire course of the study and was based on MRI activity—appearance of new/enlarging T2 lesions or gadolinium-enhancing lesions, and/or on clinical activity—presence of relapse (new neurological symptom not associated with fever/infection, lasting at least 24 h) and sustained disability progression [increase in EDSS by ≥ 1.5 points if the baseline EDSS score was 0, by ≥ 1.0 point if the baseline EDSS score was 1.5, and by ≥ 0.5 points if the baseline EDSS score was > 5.0; (16)].
The study protocol was approved by institutional ethics committee and patients gave written informed consent in accordance with the Declaration of Helsinki.
In the DMF group, absolute lymphocyte counts and lymphocyte subsets (CD3+, CD4+, CD8+, CD56+, CD19+) were quantified with flow cytometry. Blood samples were collected at baseline (at treatment onset) and later repeatedly with a 6-month interval. Samples collected after almost one and a half years follow-up (15.8 ± 7.2 months) were used for the analysis. Fresh blood samples were drawn into EDTA-containing tubes and exposed to corresponding monoclonal antibodies (BD Biosciences) at room temperature. After erythrocyte lysis and double washing, absolute values of lymphocyte subsets were counted with TruCount beads (BD Biosciences).
Baseline and follow-up MRI scans were acquired in the study setting with a 32-channel head coil 3T MRI scanner (Magnetom Tim Trio, Siemens Healthcare) according to a standardized protocol (14) comprising sagittal three-dimensional (3D) T1-weighted magnetization prepared rapid gradient echo (MP-RAGE), 3D T2-weighted fluid attenuated inversion recovery (FLAIR) and diffusion tensor imaging (DTI) sequences. Acquisition parameters of applied sequences were: T1 MP-RAGE—repetition time (TR) = 1,900 ms, echo time (TE) = 2.52 ms, inversion time (TI) = 900 ms, echo train length (ETL) = 1, flip angle (FA) = 9°, matrix size = 256 × 256, field of view (FOV) = 256 × 256 mm2, slice thickness (ST) = 1 mm; T2-FLAIR – TR = 5000 ms, TE = 388 ms, TI = 1800 ms, ETL = 848, matrix size = 256 × 256, FOV = 256 × 256 mm2, ST = 1 mm; DTI – single-shot echo-planar readout, TR = 9000 ms, TE = 102 ms, 30 gradients directions with b = 900 s/mm2 and one no diffusion image with b = 0 s/mm2, matrix size = 128 × 128, FOV = 256 × 256 mm, 62 slices, in-plane resolution = 2 × 2 mm2, ST = 2.5 mm.
Sequential study pipeline is shown in Figure 1.
Figure 1. Data analysis pipeline. (A) The processed diffusion tensor images (DTI) were used for the derivation of probabilistic tractography. The number of streamlines from each region of interest (ROI, according to the Automated Anatomical Labeling (AAL) atlas) to other ROIs was calculated and a connectivity matrix for each subject was constructed. Subsequently, these white matter derived connectivity matrices were subjected to graph theoretical analysis and the network measures were used in statistical analysis. (B) Using flow cytometry, absolute counts of lymphocyte subsets in patient blood samples were estimated and their longitudinal changes correlated with the rates of cortical atrophy; differences in lymphocyte counts between the patient subgroups were analyzed. (C) From T1-weighted magnetic resonance images (MRI), cortical thickness for each ROI according to the Desikan–Killiany atlas was calculated and used to estimate the rates of cortical atrophy and construct the connectivity matrices. Subsequently, these gray matter derived connectivity matrices were subjected to graph theoretical analysis and the network measures were used in statistical analysis.
Cortical Thickness: Longitudinal Analysis
Sagittal T1-weighted images were processed using FreeSurfer (version 5.3.0, http://surfer.nmr.mgh.harvard.edu/) (17) for cortical surface reconstruction and volumetric segmentation. The longitudinal pipeline is based on creation of an unbiased within-subject template space and image, using robust inverse consistent registration. Initialization of processing steps—skull stripping, Talairach transformations, atlas registration, and parcellation—runs on the common information from the within-subject template (18–20). Cortical thickness (in mm) was quantified at each vertex of the tessellated surface as the average of the shortest distance between the GM-WM and the GM–CSF interface. All cortical surfaces and subcortical segmentations were manually checked for errors prior to the group analysis. To avoid lesion-induced tissue misclassification errors, gray matter segmentation was performed after filling of T1 hypointense lesions.
Longitudinal changes in cortical thickness between baseline and follow-up were assessed by computing the percent change: thickness at baseline MRI scan was subtracted from the thickness at follow-up MRI scan and divided by the scan interval (in years) and by the average thickness. The resulting surface maps were smoothed with a full width at half maximum (FWHM) Gaussian kernel of 10 mm.
Diffusion tensor computation and tractography analysis was performed using inbuilt functionality in FSL (ver. 5.0.8,http://www.fmrib.ox.ac.uk/fsl); details of this analysis can be found elsewhere (21, 22). In brief, the acquired diffusion data were corrected for subjects' head motion artifacts and eddy current distortions, and subjected to skull and other non-brain tissues removal. Thus, the pre-processed data were then used for computation of the tensor. The distribution of crossing fibers at each voxel of the brain for the computation of PT was estimated using BEDPOSTX (implemented in FSL) and the probability of major and secondary fiber directions was calculated (23). All images were aligned and affine-transformed into the Montreal Neurological Institute (MNI)-152 space. At each voxel a multi-fiber model was fit to the diffusion data, enabling to trace the fibers through regions of crossing or complexity. To obtain an estimate of the probability distribution of connections from each seed voxel 5,000 streamline samples were drawn. The generated tracts are volumes wherein the values at each voxel represent the number of samples (or streamlines) that passed through that particular voxel. Each tract from every seed mask in the atlas was repeatedly sampled (5,000 times) and only those tracts, which passed through at least one other seed mask were retained. For the elimination of spurious connections, tractography in individual subjects was thresholded to include only voxels through which at least 10 percent of all streamline samples had passed.
Brain Volumes and Lesion Segmentation
Quantification of GM and WM volumes at two time points was done by using voxel-based morphometry (VBM) analysis in Statistical Parametric Mapping (SPM8) software (http://www.fil.ion.ucl.ac.uk/spm). Anatomical 3D T1 and T2-FLAIR images were subjected for spatial normalization, tissue segmentation and spatial smoothing to obtain GM, WM and total brain (TB) volumes (24). Lesion segmentation tool (LST) (version 1.2.3; http://www.applied-statistics.de/lst.html) (25) was used to compute the lesion volume (LV), details of which are mentioned elsewhere (9). Briefly, 3D T2-FLAIR images were co-registered to 3D T1 images and lesion segmentation was run with 20 different thresholds for the lesion growth algorithm (25).
Graph Theoretical Analysis
GM Network Construction
The entire cerebral cortex was parcellated into 68 bilateral anatomical regions of interest (ROIs) (34 ROIs for each hemisphere) based on the Desikan-Killiany atlas (26). Cortical thickness from each cortical ROI was extracted and served for the construction of GM connectivity matrices. These connectivity matrices (size 68 × 68) for each group were obtained by computing the Pearson correlation coefficient between the anatomical regions across the group (8). Graph Analysis Toolbox (GAT) was used to threshold the matrices into multiple densities (ranged from 0.10 to 0.50) and compute the graph theoretical network measures (27).
WM Network Construction
The obtained streamlines information from PT (as described above) connecting each pair of ROIs (116—as defined in the Automated Anatomical Labeling (AAL) atlas) was used to construct the connectivity matrix for each subject (28). A more detailed description of the network construction is presented in our previous study (8). The obtained connectivity matrices were included in the graph network analysis for the computation of network properties using Brain Connectivity Toolbox (https://sites.google.com/site/bctnet/) (29).
Topological organization of GM and WM networks was assessed by computing the modularity. Modules are groups of nodes forming a distinct subnetwork, where the within module connection (correlation) is higher than the between module correlation (30). Modularity (Q) represents the strength of division of the network into modules and was calculated using the Newman's spectral algorithm (31). Since modularity is the measure of networks' segregation, higher modularity indicates more isolated subnetworks within a given network. Hence, with increasing modularity, the long-distance paths between the modules decrease and the local interconnections within the module increase.
Statistical analyses were performed using SPSS software (version 23.0; IBM, Armonk, NY, USA.). Normal distribution of the examined data was checked via Shapiro-Wilk test. For a balanced matching of subjects, a multivariate model was tested on our cohort of 1,156 patients with RRMS to select two groups of patients matched upon demographical, clinical and neuroimaging parameters at study entrance.
The H0 hypothesis: we assumed no association between clinical and brain structural responses and immune cell dynamics under DMF treatment. The between-group differences in disease activity (MRI activity, clinical relapse) over the study period were assessed by Pearson's χ2-test.
The comparison of lymphocyte counts at baseline and follow-up, and between baseline and follow-up time points in DMF-treated patients, Mann–Whitney U and Wilcoxon signed-rank tests were used respectively.
To test if there are any differences in rates of cortical atrophy between patient groups we performed general linear model (GLM) analysis on vertex-by-vertex basis, accounting for the effects of age and gender. Generated statistical parametric maps of significant group differences were corrected for multiple comparisons with Monte Carlo Z permutation cluster analysis (10.000 iterations) at a threshold of Z = 1.3 equivalent to a p-value of 0.05. In the DMF group, to assess if there is any relationship between the rate of cortical atrophy and lymphocyte subsets we applied separate GLM models, followed by Monte Carlo Z correction for multiple comparisons as described above.
To examine whether the subcortical volumes vary in time longitudinally as a function of group, a mixed-design repeated-measures ANOVA was performed with hemisphere (left and right) and time (baseline and follow-up) as within-subject factors and group as a between-subject factor with two levels (DMF and NAT), followed by Bonferroni correction.
The GM and WM network measures between baseline and follow-up were compared by paired t-test.
The multivariate analysis revealed no significant differences between the DMF and NAT groups at baseline for age [F(1, 76) = 2.75, p = 0.10], disease duration [F(1, 76) = 0.26, p = 0.60], EDSS [F(1, 76) = 1.31, p = 0.25], lesion volume [F(1, 76) = 0.24, p = 0.62], GM volume [F(1, 76) = 0.56, p = 0.45), WM volume [F(1, 76) = 0.17, p = 0.67], MRI activity [F(1, 76) = 2.55, p = 0.11] and relapses [F(1, 76) = 1.37, p = 0.24]. EDSS, GM, WM, TB and T2 lesion volumes at follow-up did not differ significantly from baseline values within the groups (p > 0.05), indicating a stable disease period of 16 months on average in these patients despite a 4-year mean disease duration. However, the subgroup analysis within the patients with disease activity showed that the lesion volume was significantly higher at follow-up (12.3 ± 8.2 mL) than at baseline (11.2 ± 8.0 mL, p = 0.04) only in the DMF group. During the 16-month follow-up both groups were homogenous in terms of MRI activity (χ2 = 2.534, p > 0.05) and clinical relapse (χ2 = 1.381, p > 0.05). In the DMF group, 60% (25 patients) had disease activity [of them 60% (15)—MRI activity, 68% (17)—clinical relapse, 28% (7)—both] and were defined as DMF non-responders (DMFNR), while 40% (17 patients) showed no signs of disease activity (DMF responders, DMFR). In the NAT group, 36% (13 patients) had disease activity [among them 19% (7) presented MRI activity, 28% (10) had clinical relapse, 11% (4) presented both] and 64% (23 patients) were without disease activity.
During the entire study period, patients with disease activity either from the DMF group or patients from the NAT group didn't have their DMD treatment escalated.
Peripheral Blood Cell Response Under DMF Treatment
The follow-up period between the baseline and last blood samples for DMFR was 14.1 ± 5.8 months and 16.1 ± 4.9 months for DMFNR. Examining the immunological profile of the DMF subgroups, lymphocyte counts did not differ between the DMFR and DMFNR at baseline. At follow-up, DMF-treated patients showed lower counts of CD3+, CD4+, CD8+, CD56+ and ALC, and a higher CD4/CD8 ratio (all p < 0.05) compared to baseline, except CD19+, which did not reach significance in DMFNR (p > 0.05) (Supplementary Table 1). However, the DMFR subgroup displayed a stronger reduction in CD8+ cells in comparison to the DMFNR (197.7 ± 97.1 /μl vs. 298.4 ± 190.6 /μl; p = 0.03) (Figure 2) and a greater change (Δ) in CD8+ cells (−206 /μl vs.-−158 /μl; p = 0.01). Similarly, a stronger drop in CD8+ cells was detected separately in patients without MRI activity compared to those with (DMFR: 202.3 ± 107.4 /μl vs. DMFNR: 312.7 ± 125.0 /μl, p = 0.02) and in patients without clinical activity compared to patients with (DMFR: 211.1 ± 98.7 /μl vs. DMFNR: 327.0 ± 144.5 /μl, p = 0.03; Figure 3).
Figure 2. Differences in CD8+, CD4+ T cell counts and CD4/CD8 ratio at follow-up between the DMF responders and DMF non-responders subgroups. DMF responders (DMFR) had lower counts of CD8+ T cells at follow-up (after 14.1 ± 5.8 months) in comparison to DMF non-responders (DMFNR); *p < 0.05, ns, not significant.
Figure 3. Differences in follow-up CD8+ T cell counts between patients with and without MRI activity, and between patients with and without clinical relapses under DMF treatment. Patients without MRI and clinical activity (DMF responders, DMFR) had lower counts of CD8+ T cells at follow-up (after 14.1 ± 5.8 months) in comparison to patients with MRI and clinical activity (DMF non-responders, DMFNR); *p < 0.05.
Cortical Atrophy Rates
The between-group comparison revealed greater rates of mean cortical atrophy in the DMF group (−2.4%) than in the NAT group (−2.1%, p < 0.05), mostly within the frontal and temporal lobes (Figure 4).
Figure 4. Comparison of cortical atrophy rates between the DMF and NAT groups. Cortical areas displaying the difference in cortical atrophy rates between the DMF and NAT groups, mapped on lateral and medial pial surfaces of the left (LH) and right (RH) hemispheres. Negative values (blue spectrum) denote cortical areas showing greater rates of cortical atrophy in the DMF group in comparison to the NAT group. Color bar indicates the significance levels in the clusters obtained from Monte Carlo simulation at p < 0.05 (Z = 1.3).
Within the DMF group, clusters of cortical atrophy were identified mainly in the frontal, temporal and parietal lobes of both hemispheres (Figure 5A). The NAT group showed clusters of regional cortical atrophy only in the right inferior parietal and rostral middle frontal areas (Figure 5B).
Figure 5. Longitudinal changes in cortical thickness and associated network measures within the DMF and NAT groups. Left: cortical areas displaying the rates of cortical atrophy (negative values, blue spectrum) and cortical thickening (positive values, red spectrum) in the (A) DMF and (B) NAT groups, mapped on lateral and medial pial surfaces of the left (LH) and right (RH) hemispheres. Color bar indicates the significance levels in the clusters obtained from Monte Carlo simulation at p < 0.05 (Z = 1.3). Right: modularity (Q) of gray matter (GM) and white matter (WM) networks at baseline (t0) and follow-up (t1) in the (A) DMF and (B) NAT groups; **p < 0.0001, ns, not significant.
DMFR patients had lower mean rates of cortical atrophy in comparison to DMFNR (−1.7% and −3.2%, respectively, p < 0.05). The areas showing a subgroup difference in atrophy rates were identified mainly in frontal, temporal and parietal lobes (Figure 6).
Figure 6. Comparison of cortical atrophy rates between the DMFR and DMFNR subgroups and associated network measures. Left: cortical areas displaying the difference in cortical atrophy rates between the DMFR and DMFNR subgroups, mapped on lateral and medial pial surfaces of the left (LH) and right (RH) hemispheres. Positive values (red spectrum) denote cortical areas showing lower rates of cortical atrophy in the DMFR subgroup in comparison to the DMFNR subgroup. Color bar indicates the significance levels in the clusters obtained from Monte Carlo simulation at p < 0.05 (Z = 1.3). Right: modularity (Q) of gray matter (GM) and white matter (WM) networks at baseline (t0) and follow-up (t1) in the DMFR and DMFNR subgroups; *p < 0.05, **p < 0.0001, ns, not significant.
Finally, we compared the DMFR subgroup with the NAT group in order to evaluate the differences in extent of structural GM loss. Lower regional rates of cortical atrophy were identified in the DMFR subgroup in the following clusters: left lingual, precuneus and right superior, and inferior parietal areas (Figure 7).
Figure 7. Comparison of cortical atrophy rates between the DMFR and NAT groups. Cortical areas displaying the difference in cortical atrophy rates between the DMFR and NAT groups, mapped on lateral and medial pial surfaces of the left (LH) and right (RH) hemispheres. Positive values (red spectrum) denote cortical areas showing lower rates of cortical atrophy in the DMFR subgroup in comparison to NAT group. Color bar indicates the significance levels in the clusters obtained from Monte Carlo simulation at p < 0.05 (Z = 1.3).
Clusters of regional cortical atrophy rates are presented in Supplementary Table 2.
Relation Between Lymphocyte Subsets and Cortical Atrophy
In the DMF group, regression analysis disclosed the associations between the cortical atrophy rates, and ΔCD4+ and ΔCD8+ cells between baseline and follow-up. The change in CD4+ cells correlated with the atrophy rate in the right superior frontal area (peak-vertex R2 = 0.383, p = 0.012). The change in CD8+ cells correlated with the atrophy rate in the left superior parietal (peak-vertex R2 = 0.490, p < 0.0001), cuneus (peak-vertex R2 = 0.583, p < 0.0001), and rostral middle frontal (peak-vertex R2 = 0.489, p < 0.0001) and right anterior cingulate (peak-vertex R2 = 0.518, p < 0.0001), lateral occipital (peak-vertex R2 = 0.531, p < 0.0001) and operculum (peak-vertex R2 = 0.557, p < 0.0001) (Figure 8 and Supplementary Table 3). The decrease in CD4+ and CD8+ cells was associated with the lower cortical atrophy rate.
Figure 8. Association between cortical atrophy rates and T cell subsets within the DMF group. Left: cortical areas displaying the correlation between cortical atrophy rates and change (Δ) in (A) CD4+ and (B) CD8+ cells between baseline and follow-up, mapped on lateral and medial pial surfaces of the left (LH) and right (RH) hemispheres. Color bar indicates the significance levels of the correlation in clusters obtained from Monte Carlo simulation at p < 0.05 (Z = 1.3). Right: scatter plots showing the association between the cortical atrophy rate in the (A) right superior frontal cluster (at peak-vertex) and ΔCD4+ cells and (B) left superior parietal cluster (at peak-vertex) and ΔCD8+ cells (dotted line represents the 95% confidence interval for the mean). The decrease in CD4+ and CD8+ T cell counts is associated with the lower rates of cortical atrophy.
By correlating the rates of cortical atrophy and absolute counts of CD4+ and CD8+ cells similar significant clusters could be obtained.
Subcortical structures showed no volumetric differences between the DMF and NAT groups (effect of group, p > 0.05). There were no differences as well in subcortical volumes accounting for time and hemisphere in both groups (time × hemisphere × group interaction, p > 0.05). The same was also true for the DMF subgroups; subcortical volumes did not differ between DMFR and DMFNR.
Gray Matter Network Measures
Longitudinal analysis of GM network measures in the DMF group uncovered increased modularity (t = 11.10, p < 0.0001) at follow-up (Figure 5A). The number of modules increased from 2 (at baseline) to 5 modules (at follow-up). In the NAT group, modularity (t = 5.73, p < 0.0001) at follow-up as well was higher than at baseline (Figure 5B) but with less modules at follow-up (2) than at baseline (3).
DMFR subgroup at follow-up displayed lower modularity (t = 5.20, p < 0.0001) in comparison to baseline (Figure 6). Within the DMFNR subgroup modularity (t = 32.11, p < 0.0001) at follow-up was higher than at baseline.
White Matter Network Measures
In both the DMF and NAT groups, WM modularity at follow-up did not differ (t = 0.77, p = 0.44 and t = 1.44, p = 0.15, respectively) from the modularity at baseline (Figures 5A,B). Within the DMF group, DMFR at follow-up exhibited lower modularity (t = 1.76, p = 0.046) than at baseline (Figure 6).
Emerging immunomodulatory therapies considerably modify the individual course of MS, possessing the potential to beneficially influence neuroinflammation and diffuse damage to the GM and WM. However, this therapeutic effectiveness comes at a price of rare, but life-threatening side effects such as the development of secondary immunologic disorders, hematopoietic diseases or progressive multifocal leukoencephalopathy (PML). Thus, individual patients' stratification to therapy response or failure is warranted. Importantly, easy to obtain and measure immunological markers of therapy response are extremely necessary to minimize tissue damage and long-term functional impairment (10, 32). Therefore, here the goal was to stratify patients as responders to DMF therapy by linking cellular responses to clinical, structural MRI and brain network dynamics.
While overall the group of patients treated with NAT displayed less cortical atrophy than patients receiving DMF, the DMF-treated patients free of disease activity with a stronger depletion of CD8+ T cells showed even less GM loss in comparison to the NAT-treated patients. Therefore, we highlight the feasible potential of CD8+ T cell subset monitoring as a marker of individual treatment response. Previous own work (11) and recent emerging data (33–35) showed that lymphocyte subsets present varying susceptibility to DMF. It is worth mentioning that the DMF-induced shifts in lymphocyte subsets cannot be definitely considered as measures of treatment response and currently the reduced counts of CD8+ T and other cells merely explain the DMF mechanism of action. Studies with larger sample size and longer follow-up are required to ascertain if depletion of CD8 (and other lymphocyte subsets) can serve as guiding markers of DMF therapy response and support treatment decisions.
Lymphocyte CD8+ counts correlating with cortical atrophy in patients upon DMF treatment points to the role of CD8+ T cells in ongoing inflammatory processes in GM. As it takes a longer time for DMF in order to achieve a complete effect, a greater loss of cortical GM and disease activity could emerge at first months after DMF treatment onset. Higher cortical atrophy rates in the DMF group could also be caused by inclusion of a higher proportion of patients with disease activity. This is likely to be due to selection bias, since our statistical method to match both patient groups was based not on propensity score methods but on a multivariate model that could possibly over-fit the model and select patients with higher disease activity. Despite a relatively high number of patients with disease activity in the DMF-treated (60%) and in the NAT-treated group (36%), the lesion volume and EDSS didn't significantly change over time. Apparently, this mismatch partly rises from the whole-group analysis because DMF non-responders still showed a significant increase in lesion volume at follow-up. In contrast, available studies evaluating the efficacy of DMF on clinical/MRI activity and brain atrophy measures report a 27% proportion of patients with new relapses, a relative reduction by 21% of disability progression and T2 lesion volume and a reduction by 21% of brain atrophy (36, 37). On the other hand, we are aware that the short follow-up period of our study precludes us to draw definite conclusions on delaying the brain atrophy under both therapies, this being one of the study limitations. As efficacy of DMDs is greater during the second and following years after DMD treatment onset (38–40), our study period of almost one and half years could be relatively enough to obtain approximate impressions on differences of cortical GM and network responses to the DMF and NAT therapy.
In order to quantify discrete structural alterations and depict local and global GM dynamics, we performed longitudinal brain network analysis. This is an emerging tool to explore disease-related reorganization processes that mirror the disease course (8, 9, 41–43). Modularity, a parameter reflecting long-range disconnection and integration of functionally interacting brain regions, is a very sensitive marker of structural integrity in patients with MS (8, 44–46). As the disease progresses in patients with RRMS, the brain circuits reorganize toward a topology of higher modularity with long-range disconnections and local structural homogeneity (9, 44).
With the aid of modularity analysis we found that the GM network dynamics is characterized by increased modularity and long-range disconnections in both DMF- and NAT-treated groups. Longitudinal brain network development toward increased modularity presumably driven by cortical reorganization processes could be an important structural fingerprint mirroring functional impairments or even transition into progressive forms of MS (9). In contrast, we observed an inverse pattern of network topology with decreasing modularity over time in DMF responders, and no differences between the DMF non-responders and NAT patients. Here, we postulate that DMF responders comprise a potentially DMF-induced slowing of neuronal damage and reversal of local and global reorganization processes in the given period.
In both NAT- and DMF-treated groups, WM network topological characteristics did not change over time. On one hand, the stability of WM network topology might be explained by immunomodulatory treatment success. The beneficial impact of NAT on the WM compartment can be assigned to positive effects on myelination, stabilization of the blood-brain barrier and less WM damage (47). On the other hand, the short study period and the methodology applied could have been insufficient to precisely track the WM pathology. The absence of structural alterations in the WM and deep GM compartments possibly suggests different effects of the studied therapies on GM and WM compartments (48). The fact that only DMF responders showed the characteristic decrease in modularity over time highlights a distinct network response of WM networks to DMF exposure as well. We must acknowledge that the conclusions drawn from the WM network analysis need to be interpreted with caution, since mapping of white matter connections by diffusion MRI tractography has inherent sources of errors, artifacts and biases that limit its anatomical accuracy and the validity of obtained estimates (49). WM lesions may lead to inaccurate tracking of termination sites of fibers or even cause deviations of fiber bundles nearby the lesions (50). The DTI protocol and computational algorithms used in this study can partly overcome the constraints posed by the microstructural complexity of WM tracts (fibers with crossing configurations, geometric distortion, folding patterns) (51).
Our results indicate that NAT therapy opposed to DMF treatment favors preservation of cortical and subcortical structural integrity but with equivalent network responses. But within the patient cohort treated with DMF, a more pronounced decline in circulating CD8+ T lymphocytes was associated with a favorable clinical outcome and advantageous structural network responses. Whether DMF could serve as a treatment strategy in NAT-incompatible patients under conditions of rigorous T cell monitoring should be further investigated.
The study protocol was approved by institutional ethics committee and patients gave written informed consent in accordance with the Declaration of Helsinki.
Conceived and designed the study: SM, MM, and SG. Acquired and analyzed patients' data: DC, FL, GG-E, NK, CR, VF, JK, and MM. Interpreted the data and drafted the manuscript: DC, FL, SB, SM, MM, and SG. Interpreted the data and revised the manuscript: VF, SM, MM, and SG.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank Cheryl Ernest for proofreading the manuscript. This study was supported by a grant from German Research Council (DFG; CRC-TR-128).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2019.01779/full#supplementary-material
2. Zipp F, Gold R, Wiendl H. Identification of inflammatory neuronal injury and prevention of neuronal damage in multiple sclerosis: hope for novel therapies? JAMA Neurol. (2013) 70:1569–74. doi: 10.1001/jamaneurol.2013.4391
3. Jacobsen C, Hagemeier J, Myhr K-M, Nyland H, Lode K, Bergsland N, et al. Brain atrophy and disability progression in multiple sclerosis patients: a 10-year follow-up study. J Neurol Neurosurg Psychiatry. (2014) 85:1109–15. doi: 10.1136/jnnp-2013-306906
4. Uher T, Horakova D, Bergsland N, Tyblova M, Ramasamy DP, Seidl Z, et al. MRI correlates of disability progression in patients with CIS over 48 months. NeuroImage Clin. (2014) 6:312–9. doi: 10.1016/j.nicl.2014.09.015
8. Muthuraman M, Fleischer V, Kolber P, Luessi F, Zipp F, Groppa S. Structural brain network characteristics can differentiate CIS from early RRMS. Front Neurosci. (2016) 10:14. doi: 10.3389/fnins.2016.00014
9. Fleischer V, Gröger A, Koirala N, Droby A, Muthuraman M, Kolber P, et al. Increased structural white and grey matter network connectivity compensates for functional decline in early multiple sclerosis. Multi Scler J. (2017) 23:432–41. doi: 10.1177/1352458516651503
10. Kroth J, Ciolac D, Fleischer V, Koirala N, Kramer J, Muthuraman M, et al. Increased cerebrospinal fluid albumin and immunoglobulin A fractions forecast cortical atrophy and longitudinal functional deterioration in relapsing-remitting multiple sclerosis. Mult Scler. (2019) 25:338–43. doi: 10.1177/1352458517748474
11. Fleischer V, Friedrich M, Rezk A, Bühler U, Witsch E, Uphaus T, et al. Treatment response to dimethyl fumarate is characterized by disproportionate CD8+ T cell reduction in MS. Multi Scler J. (2018) 24:632–41. doi: 10.1177/1352458517703799
12. Wright K, Winkler MD, Newton BD, Sormani MP, Okuda DT. Patient outcomes influenced by reduced lymphocyte counts after dimethyl fumarate initiation. Neurology. (2017) 4:e397. doi: 10.1212/NXI.0000000000000397
13. Medina S, Villarrubia N, Sainz De La Maza S, Lifante J, Costa-Frossard L, Roldan E, et al. Optimal response to dimethyl fumarate associates in MS with a shift from an inflammatory to a tolerogenic blood cell profile. Mult Scler. (2018) 24:1317–27. doi: 10.1177/1352458517717088
14. Droby A, Lukas C, Schänzer A, Spiwoks-Becker I, Giorgio A, Gold R, et al. A human post-mortem brain model for the standardization of multi-centre MRI studies. Neuroimage. (2015) 110:11–21. doi: 10.1016/j.neuroimage.2015.01.028
15. Polman CH, Reingold SC, Banwell B, Clanet M, Cohen JA, Filippi M, et al. Diagnostic criteria for multiple sclerosis: 2010 revisions to the McDonald criteria. Ann Neurol. (2011) 69:292–302. doi: 10.1002/ana.22366
16. Giovannoni G, Turner B, Gnanapavan S, Offiah C, Schmierer K, Marta M. Is it time to target no evident disease activity (NEDA) in multiple sclerosis? Multi Scler Related Disord. (2015) 4:329–33. doi: 10.1016/j.msard.2015.04.006
17. Fischl B, Salat DH, Busa E, Albert M, Dieterich M, Haselgrove C, et al. Whole brain segmentation: automated labeling of neuroanatomical structures in the human brain. Neuron. (2002) 33:341–55. doi: 10.1016/S0896-6273(02)00569-X
21. Groppa S, Herzog J, Falk D, Riedel C, Deuschl G, Volkmann J. Physiological and anatomical decomposition of subthalamic neurostimulation effects in essential tremor. Brain. (2013) 137:109–21. doi: 10.1093/brain/awt304
22. Koirala N, Fleischer V, Glaser M, Zeuner KE, Deuschl G, Volkmann J, et al. Frontal lobe connectivity and network community characteristics are associated with the outcome of subthalamic nucleus deep brain stimulation in patients with Parkinson's disease. Brain Topogr. (2018) 31:311–21. doi: 10.1007/s10548-017-0597-4
23. Behrens TE, Berg HJ, Jbabdi S, Rushworth MF, Woolrich MW. Probabilistic diffusion tractography with multiple fibre orientations: what can we gain? Neuroimage. (2007) 34:144–55. doi: 10.1016/j.neuroimage.2006.09.018
25. Schmidt P, Gaser C, Arsic M, Buck D, Förschler A, Berthele A, et al. An automated tool for detection of FLAIR-hyperintense white-matter lesions in multiple sclerosis. Neuroimage. (2012) 59:3774–83. doi: 10.1016/j.neuroimage.2011.11.032
26. Desikan RS, Ségonne F, Fischl B, Quinn BT, Dickerson BC, Blacker D, et al. An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage. (2006) 31:968–80. doi: 10.1016/j.neuroimage.2006.01.021
27. Hosseini SM, Hoeft F, Kesler SR. GAT: a graph-theoretical analysis toolbox for analyzing between-group differences in large-scale structural and functional brain networks. PLoS ONE. (2012) 7:e40709. doi: 10.1371/journal.pone.0040709
28. Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, et al. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage. (2002) 15:273–89. doi: 10.1006/nimg.2001.0978
33. Spencer CM, Crabtree-Hartman EC, Lehmann-Horn K, Cree BA, Zamvil SS. Reduction of CD8+ T lymphocytes in multiple sclerosis patients treated with dimethyl fumarate. Neurology. (2015) 2:e76. doi: 10.1212/NXI.0000000000000076
34. Li R, Rezk A, Ghadiri M, Luessi F, Zipp F, Li H, et al. Dimethyl fumarate treatment mediates an anti-inflammatory shift in B cell subsets of patients with multiple sclerosis. J Immunol. (2017) 198:691–8. doi: 10.4049/jimmunol.1601649
35. Schlöder J, Berges C, Luessi F, Jonuleit H. Dimethyl fumarate therapy significantly improves the responsiveness of T cells in multiple sclerosis patients for immunoregulation by regulatory T cells. Int J Mol Sci. (2017) 18:271. doi: 10.3390/ijms18020271
37. Gold R, Giovannoni G, Phillips JT, Fox RJ, Zhang A, Meltzer L, et al. Efficacy and safety of delayed-release dimethyl fumarate in patients newly diagnosed with relapsing–remitting multiple sclerosis (RRMS). Multi Scler J. (2015) 21:57–66. doi: 10.1177/1352458514537013
38. Tsivgoulis G, Katsanos AH, Grigoriadis N, Hadjigeorgiou GM, Heliopoulos I, Kilidireas C, et al. The effect of disease modifying therapies on brain atrophy in patients with relapsing-remitting multiple sclerosis: a systematic review and meta-analysis. PLoS ONE. (2015) 10:e0116511. doi: 10.1371/journal.pone.0116511
39. Vidal-Jordana A, Sastre-Garriga J, Rovira A, Montalban X. Treating relapsing–remitting multiple sclerosis: therapy effects on brain atrophy. J Neurol. (2015) 262:2617–26. doi: 10.1007/s00415-015-7798-0
40. Branger P, Parienti J-J, Sormani MP, Defer G. The effect of disease-modifying drugs on brain atrophy in relapsing-remitting multiple sclerosis: a meta-analysis. PLoS ONE. (2016) 11:e0149685. doi: 10.1371/journal.pone.0149685
42. Rocca MA, Absinta M, Moiola L, Ghezzi A, Colombo B, Martinelli V, et al. Functional and structural connectivity of the motor network in pediatric and adult-onset relapsing-remitting multiple sclerosis. Radiology. (2010) 254:541–50. doi: 10.1148/radiol.09090463
43. Shu N, Duan Y, Xia M, Schoonheim MM, Huang J, Ren Z, et al. Disrupted topological organization of structural and functional brain connectomes in clinically isolated syndrome and multiple sclerosis. Sci Rep. (2016) 6:29383. doi: 10.1038/srep29383
44. Gamboa OL, Tagliazucchi E, Von Wegner F, Jurcoane A, Wahl M, Laufs H, et al. Working memory performance of early MS patients correlates inversely with modularity increases in resting state functional connectivity networks. Neuroimage. (2014) 94:385–95. doi: 10.1016/j.neuroimage.2013.12.008
45. Kocevar G, Stamile C, Hannoun S, Cotton F, Vukusic S, Durand-Dubief F, et al. Graph theory-based brain connectivity for automatic classification of multiple sclerosis clinical courses. Front Neurosci. (2016) 10:478. doi: 10.3389/fnins.2016.00478
46. Fleischer V, Radetz A, Ciolac D, Muthuraman M, Gonzalez-Escamilla G, Zipp F, et al. Graph theoretical framework of brain networks in multiple sclerosis: a review of concepts. Neuroscience. (2019) 403:35–53. doi: 10.1016/j.neuroscience.2017.10.033
47. Miller D, Soon D, Fernando K, Macmanus D, Barker G, Yousry T, et al. MRI outcomes in a placebo-controlled trial of natalizumab in relapsing MS. Neurology. (2007) 68:1390–401. doi: 10.1212/01.wnl.0000260064.77700.fd
48. Deppe M, Krämer J, Tenberge JG, Marinell J, Schwindt W, Deppe K, et al. Early silent microstructural degeneration and atrophy of the thalamocortical network in multiple sclerosis. Hum Brain Mapp. (2016) 37:1866–79. doi: 10.1002/hbm.23144
50. Ciccarelli O, Catani M, Johansen-Berg H, Clark C, Thompson A. Diffusion-based tractography in neurological disorders: concepts, applications, and future developments. Lancet Neurol. (2008) 7:715–27. doi: 10.1016/S1474-4422(08)70163-7
Keywords: multiple sclerosis, structural integrity, gray matter networks, white matter networks, immunocellular response, personalized therapy
Citation: Ciolac D, Luessi F, Gonzalez-Escamilla G, Koirala N, Riedel C, Fleischer V, Bittner S, Krämer J, Meuth SG, Muthuraman M and Groppa S (2019) Selective Brain Network and Cellular Responses Upon Dimethyl Fumarate Immunomodulation in Multiple Sclerosis. Front. Immunol. 10:1779. doi: 10.3389/fimmu.2019.01779
Received: 08 January 2019; Accepted: 15 July 2019;
Published: 30 July 2019.
Edited by:Robert Weissert, University of Regensburg, Germany
Reviewed by:Pamela Ann McCombe, University of Queensland, Australia
Zsolt Illes, University of Southern Denmark, Denmark
Copyright © 2019 Ciolac, Luessi, Gonzalez-Escamilla, Koirala, Riedel, Fleischer, Bittner, Krämer, Meuth, Muthuraman and Groppa. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sergiu Groppa, firstname.lastname@example.org
†These authors have contributed equally to this work