MR Elastography-Based Assessment of Matrix Remodeling at Lesion Sites Associated With Clinical Severity in a Model of Multiple Sclerosis

Magnetic resonance imaging (MRI) with gadolinium based contrast agents (GBCA) is routinely used in the clinic to visualize lesions in multiple sclerosis (MS). Although GBCA reveal endothelial permeability, they fail to expose other aspects of lesion formation such as the magnitude of inflammation or tissue changes occurring at sites of blood-brain barrier (BBB) disruption. Moreover, evidence pointing to potential side effects of GBCA has been increasing. Thus, there is an urgent need to develop GBCA-independent imaging tools to monitor pathology in MS. Using MR-elastography (MRE), we previously demonstrated in both MS and the animal model experimental autoimmune encephalomyelitis (EAE) that inflammation was associated with a reduction of brain stiffness. Now, using the relapsing-remitting EAE model, we show that the cerebellum—a region with predominant inflammation in this model—is especially prone to loss of stiffness. We also demonstrate that, contrary to GBCA-MRI, reduction of brain stiffness correlates with clinical disability and is associated with enhanced expression of the extracellular matrix protein fibronectin (FN). Further, we show that FN is largely expressed by activated astrocytes at acute lesions, and reflects the magnitude of tissue remodeling at sites of BBB breakdown. Therefore, MRE could emerge as a safe tool suitable to monitor disease activity in MS.


INTRODUCTION
Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system (CNS) that represents the most common cause of non-traumatic disability in young adults. MS is considered to be an autoimmune disease in which selfreactive immune cells gain access to the CNS leading to myelin destruction and neuronal damage and the subsequent formation of multifocal lesions (1). These pathological hallmarks are also characteristic for experimental autoimmune encephalomyelitis (EAE), which is the prototypical model for MS (2).
In MS, clinical relapses correlate with the development of perivenular inflammatory lesions inside the CNS. Lesion monitoring is commonly required to make an accurate diagnosis and to monitor disease progression and response to treatment. Typically, new lesions are visualized using gadolinium-based contrast agents (GBCA) on MRI. During active inflammatory processes, GBCA cross the leaking blood-brain barrier (BBB), enter the CNS parenchyma, and alter the magnetic properties of the tissue, reducing the T1 relaxation time. Generally, the more severe the inflammatory activity, the greater the burden of GBCA-enhancement on post-contrast T1-weighted scans (3). However, GBCA-MRI has limitations due to safety reasons and lack of sensitivity. We previously showed in EAE that certain lesions may remain undetectable by GBCA (4,5). Furthermore, there are emerging concerns that GBCA may deposit in the tissue after repeated applications (6).
Magnetic resonance elastography (MRE) could represent a promising alternative to the use of GBCA. MRE provides information about the mechanical properties of tissues (7,8) by analyzing their response to oscillatory shear stress (9). Using MRE, we reported on reduced brain viscoelasticity in patients with clinically isolated syndrome (10) as well as patients with established relapsing-remitting (11) and chronic-progressive (12) MS. Brain tissue softening was also observed in patients with a neuromyelitis optical spectrum disorder (NMOSD) (13). However, the mechanisms underlying brain softening in patients remain elusive. Our aim is therefore to dissect these mechanisms using animal models of neuroinflammation. Application of MRE in the cuprizone mouse model of demyelination demonstrated that brain softening was associated with demyelination (14). The observed reduction of brain viscoelasticity in this model was not related to tissue inflammation. However, in the EAE model of MS, we previously showed that inflammation was associated with reduced brain elasticity on MRE (15,16).
In the present study, we evaluated the capacity of MRE to reveal acute inflammatory activity and disease severity, comparing sagittal MRE with conventional GBCA-MRI in the relapsing-remitting EAE model. Furthermore, we aimed to clarify the nature of the tissue alterations at lesion sites by examining changes to the extracellular matrix (ECM). For this we concentrated attention on the proteoglycan fibronectin (FN), as it has been shown to accumulate within perivascular lesions, and correlate with the degree of inflammation (17,18). We thereby endeavor to improve our understanding of pathological changes detected by MRE but not by conventional GBCA-MRI.

Animals
All procedures were approved by the Animal Welfare Department of the State Office of Health and Social Affairs Berlin (LAGeSo), in accordance with national and international guidelines to minimize discomfort to animals (86/609/EEC). Experimental SJL mice (Charles River Laboratories, Sulzfeld, Germany) were housed in the central animal facility of the Charité-Universitätsmedizin Berlin. All animals were kept in a temperature-and humidity-controlled colony room and maintained on a light/dark cycle of 12/12 h with ad libitum access to food and water. Over multiple experiments, mice were divided into two groups, EAE (n = 25) and non-manipulated healthy control (n = 7).

In vivo Scans
In vivo MRE and MRI scans were performed as described previously (16) on a 7 T Bruker Pharmascan 70/16 rodent MR scanner (Bruker Biospin, Ettlingen, Germany), running Paravision 5.1 software, with a 20 mm RF quadrature volume head coil (RAPID Biomedical GmbH, Rimpar, Germany). Mice were anesthetized with 1.5-2.0% isoflurane in 30% O 2 and 70% N 2 O administered via face mask, with continuous respiration monitoring using a pressure-sensitive pad placed on the thorax (Small Animal Instruments Inc., Stony Brook, NY, USA). The animals were placed on a bed with circulating heated water to maintain constant body temperature. In those animals investigated by both MRE and MRI, the MRE measurements were acquired first, followed by GBCA-MRI after 24 h.

Magnetic Resonance Elastography (MRE)
The MRE images were acquired in one 2 mm midsagittal slice as described previously (16). Mechanical vibration was generated by an air-cooled electromagnetic Lorentz coil in the fringe field of the MRI scanner. Vibrations were initiated by a trigger pulse from the control unit of the MRI scanner, and transferred to the animal through a carbon fiber piston, which was connected to the bite bar transducer. The transducer was gimbaled through a rubber bearing and retaining bracket at the temperaturecontrolled mouse bed. A plastic disk held up the entire setup in the center of the magnet bore (16).
The timing of the vibration was defined and recorded by a fast low-angle shot (FLASH) sequence especially made for MRE measurement. The direction of the motion sensitizing gradient (MSG), with a strength of 285 mT/m, a frequency of 900 Hz, and 9 periods, was maintained parallel to the principal axis of the magnetic field (maximum amplitude of the mechanical driver was 10 microns). To compensate for the static phase contributions, phase difference images were calculated from two images differing in the sign of the MSG. Frequency amplitude and the number of cycles were controlled by a waveform generator connected via an audio amplifier to the driving coil. Additional scan parameters were as follows: TE = 14.3 ms; TR array = 166.0 ms; slice thickness = 2.0 mm; matrix = 128, FOV = 25 mm; two averages; eight dynamic scans over a vibration period and an acquisition time of 12 min.

MRE Data Analysis
Complex wave images corresponding to the harmonic drive frequency were extracted by temporal Fourier transformation of the unwrapped phase-difference images. To reduce noise, in addition to the Butterworth band pass filter, a spatiotemporal directional filter was applied to the wave images (20). The spatiotemporal filter filtered waves that were propagating from bottom-to-top in the sagittal slice. A 2D-Helmholtz inversion was performed to the filtered data, yielding the complex shear modulus G * and the magnitude modulus |G * | = abs (G * ). The calculated spatially averaged G * -values were represented by the real part of the complex shear modulus G' = Re (G * ), known as the storage modulus that represents tissue elasticity, and the imaginary part G ′ = Im (G * ), which is the loss modulus representing tissue viscosity. The magnitude, storage and loss moduli were expressed in pascals (Pa). The loss factor calculated as the phase angle phi = arctan (G ′′ / G ′ ) represents the fluidity of the tissue, which is the degree of viscosity relative to elasticity and is interpreted as being sensitive to the architecture of viscoelastic networks in biological tissues (21). In addition to calculating values for the storage and loss moduli for the entire sagittal slice, the brain was separated into two regions of interest divided at the junction between the cerebrum and the cerebellum.

Magnetic Resonance Imaging (MRI)
T1 maps were generated using a saturation recovery RAREVTR method, in which the repetition time (TR) was varied to acquire a series of axial T1 weighted images, from which the T1 map was produced. Scan parameters were as follows: TE = 8.3 ms; TR array = 230, 460, 1061, 1485, 2080, 3080 and 7500 ms; flip angle = 90 • /180 • ; RARE factor = 2; slice thickness = 1.0 mm; matrix = 128, FOV = 1.92 cm; NA = 1, 10 slices, scan time = 17 min 44 sec. After acquiring the precontrast T1 map, the animals were administered 0.2 mmol/kg gadopentate dimeglumine (Gd-DTPA, Magnevist, Bayer Vital GmbH, Leverkusen, Germany) by intravenous injection. After 5 min, the post-contrast T1 maps were acquired using the same parameters as above. Data acquisition was done with ParaVision 5.1 (Bruker Biospin, Germany). The raw data files were exported as NIFTI image files, and analyzed in ImageJ v. 1.51 (NIH, open source). A region of interest (ROI) defining the brain was manually traced for all 10 slices, and the mean T1 value of each ROI calculated. The mean T1 from all 10 slices was determined for each animal, and the post-contrast mean was subtracted from the pre-contrast mean, to yield the difference-delta T1. The delta T1 was used for the statistical analysis.
In addition, T1-weigthed images were acquired as described previously (22). As a complimentary method, we also calculated the T1 signal intensity change directly from the T1-weighted images. An ROI defining the brain was manually traced for 20 slices and the mean signal intensity (SI) value of each ROI was calculated. The mean SI from all 20 slices was determined for each animal in both pre-and post-contrast, calculated as: signal intensity change (SI%) = [(SI post-contrast-SI pre-contrast)/SI pre-contrast] * 100.

Tissue Processing
Mice were sacrificed 1 day after MRE and MRI measurements. Animals were deeply anesthetized with ketamine/xylazine, then transcardially perfused with PBS. Twenty-five brains were extracted and 19 of them were cut sagittally in two symmetrical halves. Half of the brain was postfixed in 4% paraformaldehyde (PFA) overnight at 4 • C and then PBS washed, followed by 30% sucrose in PBS, soaking until the tissue sunk to the bottom. Meanwhile, the other half of the brain was reserved for RNA extraction. For histological analysis, brain tissue was embedded in O.C.T., frozen in methylbutane with dry ice, and stored at −80 • C. For analysis, frozen tissues were cut into 12 µm coronal or sagittal cryosections and stored at 4 • C.
Quantitative Reverse-Transcription PCR The brain tissue for PCR was divided into two portions, the anterior (cerebral) and posterior (cerebellar) region, according to the MRE scanning regions. Total RNA was extracted from the tissue by using the Trizol method. The RNA was reverse transcribed, and quantitative PCR (qPCR) carried out as described previously, using an ABI Prism 7000 SequenceDetection System (Applied Biosystems, Darmstadt, Germany) (16). Primers and probes were from Eurofins MWG Operon (Ebersberg, Germany), and the sequences used were as followed: 18s, served as the endogenous reference, forward 5'-TTCGAACGTCTGCCCTATCAA-3 ′ , reverse 5 ′ -TCCC CGTCACCCATGGT-3 ′ , probe 5 ′ FAM-TGATGTTTATTGACAA CACGCTTTACTTTATACCTGAAGA-3 ′ TAMRA. We used the 2 − C T method to analyze the results.

Statistical Analysis
Data were analyzed by unpaired two-tailed t-test, paired twotailed t-test or repeated-measures analysis of variance (ANOVA), as appropriate. The non-parametric Spearman correlation was used to assess correlation between MRE parameters and EAE score. Pearson correlation was used to assess correlation between imaging parameters, and between MRE parameters and the PCR or immunostaining quantification data. Analysis was done using GraphPad Prism v.5.01. (GraphPad software, La Jolla, CA, USA). * p < 0.05, * * p < 0.01, * * * p < 0.001.

RESULTS
During Relapsing-Remitting EAE, Brain Viscoelasticity Is Altered Particularly in the Cerebellum To assess the sensitivity of the sagittal MRE protocols in the SJL EAE model (16), we monitored mechanical brain changes in SJL mice immunized with PLP. Mice developed a typical relapsingremitting disease course, showing first clinical signs by 9-10 days after immunization, reaching peak disease 3-4 days later (day 12-14 p.i.; Figure 1A). MRE measurements were performed at day 14 post-immunization, when clinical signs of disease were well-established and coinciding with the expected disease peak. MRE data were acquired in one 2 mm midsagittal slice and confirmed our previous MRE data in relapsing-remitting EAE acquired in coronal slices (15). We observed that at day 14-15 p.i. the overall viscoelasticity of the tissue (|G * | magnitude modulus) as well as the storage modulus G' (elasticity) and the loss modulus G" (viscosity) were significantly diminished in EAE mice compared to controls (Figures 1B-D). No alteration of the phase angle (loss factor) was observed (Figure 1E), indicating that the overall architecture of the brain tissue was not affected by acute inflammation. We previously showed in C57/BL6 mice that different brain regions also show distinct viscoelastic properties (16). However, it remained unclear whether this was also applicable in the  SJL EAE model, and whether different brain regions may also differ in the MRE values during relapsing-remitting EAE. A representative example of a midsagittal slice from a SJL EAE mouse is shown in Figure 2A, illustrating the magnitude image A (1), the wave deflection image A (2) and the magnitude of the complex shear modulus |G * | A (3). The acquisition of MRE data in midsagittal slices permits the separate analysis of the anterior region (cerebrum) and the posterior brain region (cerebellum) (ROIs in Figure 2A-1). The results show that during EAE, the cerebellum shows a striking reduction of the magnitude, storage and loss moduli, when compared with the cerebrum or the whole brain ( Figure 2B). The phase angle remained stable in the investigated regions. Furthermore, to estimate the effect of disease in these two different brain regions, we calculated the viscoelastic changes of the EAE tissue in relation to sex and age matched healthy controls. Figure 2C demonstrates that viscoelastic values at the time of the expected peak of EAE were decreased in both regions when normalized to healthy mice, as indicated by mean difference |G * | < 0, indicating that during EAE both regions undergo a "softening" of the tissue. Nevertheless, the magnitude modulus decreased more strikingly in the cerebellum, consistent with the fact that this region is more affected by inflammatory pathology than the cerebrum in the SJL EAE model.

Sagittal MRE Measurements Show a Correlation Between Mechanical Brain Properties and Clinical Disability in the Relapsing-Remitting EAE Model
It is well-established in both MS patients (23) and EAE (24), that MRI measurements do not always correlate with clinical disability. Using MRE acquired in the sagittal plane, we observed statistically significant correlations between the magnitude and storage moduli and the clinical score on the day of MRE acquisition, p = 0.0117, r = −0.6306 and p = 0.0373, r = −0.5411, respectively (Figures 3A,B). No significant correlations between loss modulus, phase angle and  Figures 3C,D). These results indicate that there is an association between brain softening and more severe clinical signs in EAE mice.

MRE Does Not Correlate With Gadolinium Enhancement in Acute EAE
Contrast agent-based MRI represents the standard MR tool in MS to detect acute inflammatory lesions (25). Therefore, we asked whether the intensity of gadolinium (Gd) enhancement in brain MRI also correlated with clinical disability and MRE values in the SJL EAE model. We acquired axial and coronal T1-weighted images before and after intravenous injection of 0.2 mmol/kg Gd contrast agent, 24 h subsequent to the MRE scans. Brain lesions in active EAE are unevenly distributed in space, and are highly diffuse, lacking clear boarders. In order to better evaluate the extent of these diffuse lesions, we generated T1 maps to yield a quantitative measure of brain tissue contrast changes resulting from GBCA leakage across the blood-brain barrier. The T1 values were averaged from ROIs defining the entire brain, in order to obtain a global metric of Gd enhancement in the whole brain. Representative images illustrating the T1 maps preand post-contrast are shown Figure 4A. As expected, the postcontrast T1 values for all animals were significantly reduced, compared to the pre-contrast values (data not shown). Contrary to expectations, there was no statistically significant correlation between the delta T1 and the magnitude modulus (p = 0.8797, r = 0.05181, Figure 4B) or between T1 signal intensity changes directly obtained from the whole brain T1-weighted images and MRE values (p = 0.3402, r = 0.2315, Figure 4C). Accordingly, no significant correlation was observed between T1 signal intensity changes and EAE scores of the mice (p = 0.2325, r = −0.4195, Figure 4D).

Reduction of Viscoelasticity Correlates With an Increased Expression of FN
In view of the correlation between acute inflammatory events and loss of brain stiffness as measured by MRE, we considered whether viscoelastic changes may reflect processes of molecular remodeling that occur during the formation of brain lesions. Here we found that reduced viscoelasticity during EAE was significantly correlated with increased gene expression of FN in brain tissue, p = 0.0041, r = −0.9473 (Figure 5A). This correlation was robust, and was also confirmed in a separate analysis using frozen tissue from our previous study in SJL EAE acquiring MRE data in the coronal orientation (15), p = 0.0163, r = −0.6737 (Figure 5B).
To assess whether this association may involve changes of other components of the ECM, we investigated the gene expression of other key ECM proteoglycans, associated with the BBB (eg. glypican 5), or with the perineuronal and interstitial matrix (e.g., neurocan and brevican). No significant correlation was observed between brain viscoelasticity and the expression of these key proteoglycans, p = 0.9755 and r = −0.01632, p = 0.3062 and r = −0.5056, p = 0.8671 and r = −0.08886, respectively ( Figure 5C).

FN Deposits Are Predominantly Found in Perivascular Areas, and Reduced Viscoelasticity Correlates With Increased FN Protein Expression
To confirm that FN gene expression is associated with protein deposits, we performed immunohistochemical analysis in the corresponding tissue, and demonstrated that at the time of the expected peak of EAE severity FN was highly expressed in perivascular areas ( Figure 6A). In agreement with the gene expression data, reduced viscoelasticity correlated with increased FN immunofluorescence intensity ( Figure 6B). In particular, a significant correlation was observed with the magnitude modulus and storage modulus (p = 0.0083, r = −0.8447 and p = 0.0127, r = −0.8201, respectively); whereas, no correlation of FN immunostaining with the loss modulus or phase angle could be detected (p = 0.1766, r = −0.5301 and p = 0.2889, r = −0.4290, respectively).
FN present in the lesions may arise either from the circulating plasma Fn (pFn), deposited following BBB disruption, or synthesized locally by glial cells (cellular cFn). These two FN variants contain alternatively spliced domains, EIIIA and EIIIB. To clarify the source of the FN observed in EAE, we used an anti EIIIA-FN antibody (IST9), which only recognizes cellular FN, in combination with the broad-spectrum anti-FN antibody. The staining of cellular FN showed a consistently overlapping pattern with the total FN, indicating that the main source of the perivascular FN deposits was indeed cellular, rather than circulating pFn (Figure 6C). The co-staining of FN with the astrocytic marker GFAP further indicated that not only astroglia but also endothelial cells may be major sources of FN production in the lesion sites ( Figure 6D).

DISCUSSION
In this study, we investigated the capacity of MRE to detect acute inflammatory events in the mouse brain during relapsingremitting EAE. MRE investigations were performed together with T1-weighted MRI using GBCA and were correlated with clinical disability. Furthermore, we assessed the association between MRE values and molecular aspects of extracellular matrix remodeling at the sites of inflammatory lesions.
Studies using MRE to show changes in brain tissue stiffness in the context of pathology should consider the baseline values of normal brain, though it is crucial to take into consideration technical differences among various research groups. For example, Murphy et al. used a higher stimulation frequency (1,500 Hz), yielding higher mouse brain stiffness values (25.0 kPa) (26). Other studies using stimulation frequencies similar to the setup in the current study show that the values we obtained are within the range of other published values [reviewed in Table 1 of Bertalan et al. (27)]. Consistent with our previous studies using coronal MRE (15), we show here that MRE using sagittal slices in SJL EAE revealed a decrease in the overall brain viscoelasticity at the time of the expected peak of disease, compared to the corresponding healthy controls. Furthermore, confirming our data in C57/BL6 mice (16), we demonstrated also in SJL mice that the cerebellum is softer than the cerebrum, and that these regional differences were maintained during EAE. In the SJL model, it is well-documented that the cerebellum is especially susceptible to blood-brainbarrier disruption and lesion formation (28). This together with the significant softening of the cerebellum suggests that areas of high inflammatory activity are especially well-identified by MRE. However, confirming the sensitivity of the sagittal MRE scans for visualizing whole brain inflammation, we observed a significant correlation between whole-brain brain viscoelasticity and EAE clinical score. No correlation was established between cerebellar viscoelasticity and EAE score. Thus, although changes in the cerebellum are relevant and pronounced, they do not seem to reflect all the processes that ultimately determine the clinical EAE signs.
To date, GBCA-MRI represents the standard approach to identify BBB breakdown in MS patients and also in animal models (29-31). However, there is increasing concern from recent studies and case reports pointing out potential toxic side effects of GBCA (32), thus motivating the desire for alternative MRI methods that can both better detect and improve understanding of the nature of the pathology during disease. In the active EAE model, GBCA-enhancing lesions are numerous, but rather small and diffusely distributed, and are thus inherently difficult to quantify using T1-weighted images. Therefore, we used T1 mapping to yield an unbiased, quantitative readout of the burden of lesion activity in the entire brain-a strategy that was recently applied in a mouse brain tumor model (33). Using the delta T1 averaged over the entire brain as a metric of BBB disruption also avoids potential sampling bias that might occur when attempting to quantify lesions in selected brain slices, given the uneven distribution of brain lesions in SJL EAE. Contrary to our expectations, there was no significant correlation between the overall magnitude of GBCA-enhancing lesions and the wholebrain MRE parameters. This underscores the complex nature of multiple pathological processes that occur simultaneously in EAE. It may be that the relationship between viscoelastic changes and GBCA-enhancement cannot be detected with a simple linear correlation. This may reflect that the disease processes that lead to viscoelastic changes have different kinetics from those of acute BBB disruption and active lesions. This remains a topic for further investigations.
Inflammatory lesions are not only defined by enhanced endothelial permeability and disruption of the BBB, but also by astroglial activation and remodeling of the extracellular matrix components (34). In particular FN-mediated signaling seems to promote vascular remodeling during demyelinating disease (35). Therefore, we investigated the relationship between alteration of mechanical properties of the brain and expression of FN, as an indicator for neurovascular junction remodeling during lesion formation. Our data indicate that reduction of brain viscoelasticity is indeed associated with an overall increase of FN expression in the brain tissue (36, 37).
FN is a multidomain glycoprotein binding to cell-surfacereceptors, mostly integrins, and to the ECM and appears in two forms: plasma and cellular FN. Cellular FN is secreted as a soluble covalent dimer, and in a complicated process in which the FN molecule undergoes different conformational changes, is assembled into a stable matrix (38). Our data point to a cellular source of the FN deposits detected in the perivascular areas. FN fibers were in close proximity to the reactive astrocytes, confirming previous reports that glia cells are major producers of the protein (39). However, further studies are required to determine the sources of FN. So far, we could not exclude other cellular sources such as brain endothelial cells. Additionally, we observed that FN deposits characterized disturbed BBB with enhanced perivascular space and astrocyte endfoot detachment. Thus, although the FN fibers are extremely elastic (40), and might be expected to contribute to enhanced tissue elasticity, in this context their presence might reflect an assembly of processes involved in disruption of the neurovascular unit at the lesion sites. This may lead to tissue softening, due to enhanced endothelial permeability, alteration of blood flow, enlargement of the perivascular spaces, inflammation or astrocytic endfeed detachment.
Importantly, the correlation of MRE data with FN deposition could also be confirmed in brain tissue from our previous MRE study in SJL animals, in which MRE data was acquired using a coronal slice (15). In contrast, no correlation was found between MRE values and other ECM components such neurocan, brevican, or glypican, which are reported to be relevant for neuroinflammation, but which do not directly reflect the process of acute lesion formation (36, 37). Furthermore, we found that deposits of FN are prominent around inflamed vessels and that the presence of FN also correlated with the overall softening of the brain tissue at acute EAE.
Altogether our study sheds light on the mechanism of brain softening due to inflammation. Previous studies in mouse models and patients have demonstrated in vivo that brain tissue becomes softer during the progression of neuroinflammatory processes. Our data demonstrate that some of these observations might be explained by mechanisms involving the enhancement of perivascular spaces and astrocyte endfoot detachment leading to weaker couplings between the neuronal-vascular networks, and reduced tissue stiffness. This might also explain the high sensitivity of MRE to processes of acute inflammation and lesion formation in the brain when compared to GBCA-MRI. The application of MRE in mouse models of brain disease is an emerging field. Future experiments to evaluate how well MRE can predict disease development in the pre-onset phase, or to identify tissue changes and accurately reflect disease severity during clinical relapse and remission, hold promise. Further technical developments that will allow faster acquisition times to obtain MRE data covering the entire mouse brain with reasonable scan times, as well as improvements in spatial resolution will be essential to advance our understanding of the complex processes occurring during neuroinflammation. The results of the present study lay a foundation for such upcoming studies to investigate in more detail the temporal and spatial changes of tissue mechanics in relation with lesion development.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the corresponding author, without undue reservation, to any qualified researcher.

ETHICS STATEMENT
This study was carried out in accordance with the recommendations of national and international guidelines to minimize discomfort to animals (86/609/EEC. The protocol was approved by the Animal Welfare Department of the State Office of Health and Social Affairs Berlin (LAGeSo).