- Open Access
Ultra-High Field Diffusion MRI Reveals Early Axonal Pathology in Spinal Cord of ALS mice
Translational Neurodegenerationvolume 7, Article number: 20 (2018)
Amyotrophic lateral sclerosis (ALS) is a disease characterized by a progressive degeneration of motor neurons leading to paralysis. Our previous MRI diffusion tensor imaging studies detected early white matter changes in the spinal cords of mice carrying the G93A-SOD1 mutation. Here, we extend those studies using ultra-high field MRI (17.6 T) and fluorescent microscopy to investigate the appearance of early structural and connectivity changes in the spinal cords of ALS mice.
The spinal cords from presymptomatic and symptomatic mice (80 to 120 days of age) were scanned (ex-vivo) using diffusion-weighted MRI. The fractional anisotropy (FA), axial (AD) and radial (RD) diffusivities were calculated for axial slices from the thoracic, cervical and lumbar regions of the spinal cords. The diffusion parameters were compared with fluorescence microscopy and membrane cellular markers from the same tissue regions.
At early stages of the disease (day 80) in the lumbar region, we found, a 19% decrease in FA, a 9% decrease in AD and a 35% increase in RD. Similar changes were observed in cervical and thoracic spinal cord regions. Differences between control and ALS mice groups at the symptomatic stages (day 120) were larger. Quantitative fluorescence microscopy at 80 days, demonstrated a 22% reduction in axonal area and a 22% increase in axonal density. Tractography and quantitative connectome analyses measured by edge weights showed a 52% decrease in the lumbar regions of the spinal cords of this ALS mice group. A significant increase in ADC (23.3%) in the ALS mice group was related to an increase in aquaporin markers.
These findings suggest that the combination of ultra-high field diffusion MRI with fluorescent ALS mice reporters is a useful approach to detect and characterize presymptomatic white matter micro-ultrastructural changes and axonal connectivity anomalies in ALS.
Amyotrophic lateral sclerosis (ALS) involves progressive deterioration of upper and lower motor neurons within the brainstem, corticospinal tracts, and anterior horn areas of the spinal cord (SC) . It has been shown that genetic mutations in ALS patients lead to changes in molecular pathways and neuronal degeneration in selective groups of cells subsequently promoting abnormalities in genomic expressions  and axonal function . Among the ALS transgenic animal models available, the G93A-SOD1 model is the most widely used to study the neurological deterioration associated with this disease .
Unfortunately, other than an increase in the number of available animal models [5,6,7], little progress has been made toward the early detection of ALS to improve the outcome of therapeutic interventions for these patients. However, a better way to understand the early changes occurring in cellular structures is to include transgenic mouse lines with constitutively expressed fluorescent protein. Thus, a specific neuronal population can be visualized with quantitative fluorescence techniques to assess the morphological aspects of neuronal structures and to understand the role of axonal connectivity in neurological diseases [8,9,10]. One of the most commonly used fluorescent probes in research is the Thy1-YFP mouse model, expressing yellow fluorescent protein (YFP) in a subset of neurons across different brain structures and axon trajectories in specific layers of the cerebral cortex and spinal cord (SC) . Although the study axonal degeneration with transgenic fluorescent ALS mice has been described in the past , such studies were mainly performed at the end stage of the disease with limited focus on changes in the SC white matter (WM) microstructure.
On the other hand, since the average life expectancy of ALS patients from the time of diagnosis ranges from two to five years, any meaningful improvement in the ALS survival rate would depend on the establishment of methods for the detection of early and critical pathogenic events. To this end, magnetic resonance imaging (MRI) provides the best non-invasive way to assess the appearance of neurodegenerative diseases at early presymptomatic stages [13, 14]. In addition, diffusion tensor imaging (DTI) provides substantial information on the organization of neural tissues . As an example, the radial and axial diffusion parameters provide an index of axonal fiber tract integrity and fractional anisotropy gives a measure of brain tissue organization, which is important in the context of ALS. In the search for early ALS biomarkers, recent studies using tractography reconstructions related fiber track changes with a selective degeneration in WM axonal populations [16, 17]. Ultimately, improvements in the sensitivity and specificity of MRI diffusion methods have introduced new imaging biomarkers that could able to identify the changes occurring in the susceptible WM tissues caused by ALS .
Although several reports have shown that DTI methods can detect early alterations in WM, and extensive neuropathological evidence points towards early white matter structural alterations in ALS [19,20,21,22,23], the underlying link between the diffusion and the underlying axonal neuropathological process at early stages of the disease has not been fully elucidated. Thus, our study is designed to validate DTI markers obtained by ultra-high field MRI (UHF-MRI) and to characterize early changes in axonal ultrastructure and WM connectivity using an ex vivo tissue and a fluorescent ALS mouse.
All procedures used to obtain tissues followed an approved protocol from the animal care committee (ACC) at the University of Illinois in Chicago (UIC). C57BJ6 mice, overexpressing the SOD1 transgene with the G93A mutation, were obtained from the Jackson Laboratory (JAX # 004435). The G93A-SOD1 mice have been extensively characterized as an animal model for ALS, developing motor symptoms at approximately 110 days of age and dying around 160 days [24, 25]. Based on these finding, we considered two groups of animals for this work: a presymptomatic group at postnatal day 80 (P80) and a symptomatic group at postnatal day 120 (P120). To evaluate morphologic axonal anomalies in the context of ALS, mice encoding a yellow fluorescent protein (YFP) transgene specifically associated with a neuronal Thy1 promoter was chosen (JAX#003709). Thus, we generated double transgenic mice (YFP, G93A-SOD1) and littermates carrying only the YFP transgene used as a control group. For this study, a total 18 animals (n=5 SCs per P80 and P120 YFP, G93A-SOD1 groups (10 SCs) and n=4 SCs per P80 and 120 YFP control groups (8 SCs), were used for scanning and further histological analysis. Mice had easy access to food and water and were checked daily to assess their level of well-being and health. In any situation of animal distress or pain, animals were euthanized in carbon dioxide using standard protocols.
Animal preparation for MRI imaging
Animals were rendered unconscious with CO2 inhalation, then transcardiac perfused with a PBS and 4% paraformaldehyde (PFA) solution. A laminectomy was performed and the spinal cords were extracted intact then immersed in PFA (> 48 hours). Prior to imaging, the cords were soaked overnight in PBS to removed free fixative, then placed in individual 5 mm diameter NMR tubes (New Era #NEML5-7, 300-400 MHz) filled with fluorocarbon oil (Fluorinert®, 3M, Maplewood, MN). Each set of 9, 5-mm tubes (each containing a spinal cord) was positioned in a 20-mm diameter NMR tube (New Era # NE-L25-7) using a custom-made plastic tube holder. Images were acquired with a 17.6T vertical-bore Avance II scanner using a 20-mm RF coil, Micro-2.5 gradients, and Paravision 6.0 software (Bruker, Karlsruhe, Germany).
Diffusion weighted imaging and data processing
For each set of 9 spinal cords (2 imaging sessions) a total of 60 MRI slices were acquired in blocks of 20 slices, centered at the cervical, thoracic, and lumbar levels each, and oriented along the rostral-caudal axis of the spinal cord. Diffusion weighted images were acquired using a spin echo sequence with TR= 4000 ms and TE =28 ms, interleaved 0.15 mm thick slices, field of view = 20 x 20 x 3 mm3 in each block of slices, in-plane acquisition matrix = 133 x 133, for an isotropic image resolution of 150 μm. For connectomics calculations (), two images were acquired with b = 0 s/mm2 and diffusion weighting was applied with b = 700 s/mm2 in 12 directions, and b = 2500 s/mm2 in 64 directions , with 3.5 ms gradient pulses and 17.5 ms separation. This acquisition was averaged twice for a total acquisition time of 19 hrs. per set.
Diffusion data processing was performed using FSL  to calculate fractional anisotropy (FA), axial diffusivity (AD), and radial diffusivity (RD) . The average apparent diffusion coefficient (ADC) over the 64 directions was calculated using b = 0 and b = 2500 s/mm2. In addition, the calculation of the water displacement probability density function  was used to estimate white matter tracts. White matter regions-of-interest (ROIs), defined within slices at the top and bottom of each of spinal cord segments (cervical, thoracic, and lumbar), were manually outlined following the anterolateral distribution of the white matter on each slice (Fig. 1) using ImageJ software (NIH, Bethesda).
White matter connectivity
To determine white matter structural connectivity in the three regions of the spinal cords, deterministic tractography was performed using in-house software using 125 seeds per voxel, an angular threshold of 50o, and a step size of 75 μm. Fiber tracks were visualized with Trackvis (Version 0.6.1, Massachusetts General Hospital, Boston, MA). Using a network approach to define white matter integrity along the spinal cord, the top and bottom slices in each segment were used to define network nodes. White matter tracts were used to define network edges and the connectivity quantified using a dimensionless, scale-invariant edge weight (EW) .
Histology and immuno-fluorescence analysis
Although it is feasible to histologically examine all segments of the spinal cord (cervical and thoracic and lumbar), results from this and previous MRI diffusion ex-vivo studies  have shown predominant alterations in MRI diffusion at the lumbar segments of ALS mice. Thus, this study is focused in the histological analysis of white matter anomalies in the lumbar regions, specifically between the third and fifth lumbar segments (Fig. 2).
After MRI scanning, oil media was removed and spines were placed in progressive solutions of sucrose [5-30 %] for an additional 24 hrs for cryo-protection. After embedding in optical cutting temperature (OCT) polymer compound (Tissue Tek, Sakura, Finetek, cat #4583), 50 μm thick spinal cord sections were obtained using a microtome (Leica cryostat CM 1850 Cryostat, Buffalo Grove, IL). Based in our previous finding to detect early SC white matter (WM) changes , only coronal section of lumbar spinal cord were used for histological analysis. Spinal cord sections were mounted on slides (Fisher-brand Superforst, cat# 12-550-15) and dried for 15 minutes. Then, the OCT was removed by washing three times with Tris base buffer (TBS). Sections were permeabilized with Triton-X100 0.25 % for 10 minutes and blocked with 5 % goat serum for an hour in TBS. Spinal cord sections (50μm thickness) were mounted on slides. Slides were dried and mounted in Vecta-Shield mounting media (Vector Laboratories, Burlingame, CA). Images were acquired by confocal microscopy (Leica LMS-710 confocal microscope, Germany). Each coronal section in the lumbar spinal cord was selected using similar anatomical reference following a spinal cords stereotaxic coordinates (The Spinal Cord: A Christopher and Dana Reeve Foundation Text and Atlas, 1st Ed, Watson & Paxinos, 2008). Confocal microscopy images and z stack images for three-dimensional reconstruction were obtained by background subtraction using negative controls samples without primary antibody and collected by two independent channels: 534 nm channel for the YFP yellow signal and 647 nm channel to detect fluorescent emission from antibodies from other markers. To evaluate the role of water permeability in ALS, we used anti aquaporin-4 staining (AQP4) (StressMarq Bioscence Inc., Victoria, BC) (Cat #SPC-505D, 1:400). Quantitative measures were obtained by counting the mean pixel value per equal picture area using auto-threshold methods and the pixel aggregates of each figure compiled and tabulated for analysis. Briefly, the procedure divides the image into objects and background by taking an initial threshold. Averages of the number of pixels at, below, or above the threshold were computed and subsequent averages of these two values were used.
Quantitative data were tabulated and analyzed using GraphPad Prism 6 software (La Jolla, CA). Based on the results from pilot experiments, the group size of animals per experimental group were established using power analysis and sample size calculations. For quantitative analysis of YFP, and AQP4 fluorescence levels and MRI diffusion values (ADC, FA, AD and RD), one-way ANOVA and Tukey’s post hoc tests were used to determine statistical differences among experimental animal groups (P80 and P120). A value of p < 0.05 was used to demonstrate statistical significance. Results were replicated by application of non-parametric statistical tools (Mann–Whitney test). Error bars in all the figures represent standard error of the mean (SEM).
MRI diffusion demonstrate presymptomatic spinal cord changes of YFP, G93A-SOD1 mice
Previous clinical and animal diffusion imaging studies ALS SCs have shown that regional changes in parameters were more affected at the distal than proximal regions [14, 33]. Hence, we centered our MRI diffusion studies on the anterolateral white matter funiculi of three SC regions (cervical, thoracic and lumbar segments) (Fig. 1a). Manual segmentation of WM ROIs across diffusion maps from five SC slices were considered in this analysis, (Fig. 1b). Two time-points were considered; a presymptomatic stage (P80) and a symptomatic stage (P120). At presymptomatic stages (P80), we observed a significant reduction of FA in the cervical region (p < 0.05) in the YFP, G93A-SOD1 mice (YFP mice =0.60 +/- 0.01 versus YFP, G93A-SOD1 mice = 0.57 +/- 0.01) (- 6.6%) as well as in the thoracic regions (p < 0.01) (YFP mice = 0.66 +/- 0.01 versus YFP, G93A-SOD1 mice = 0.58 +/- 0.01) (- 12.6%) and lumbar region (p < 0.001) (YFP mice = 0.66 +/- 0.01 versus YFP, G93A-SOD1 mice = 0.53 +/- 0.03) (- 19.7%). A significant decrease in axial diffusion (AD) in the lumbar segment was seen at P80 (YFP mice = 6.7 +/- 0.1x10-4 mm2/s versus YFP, G93A-SOD1 mice = 6.3 +/- 0.05x10-4 mm2/s) (p < 0.01) (-8.7 %). These changes were decreased further on each spinal cord segment during the symptomatic stage. Presymptomatic results in radial diffusion (RD) demonstrate a significant (p < 0.01) increase not only in the lumbar levels (YFP mice = 2.0 +/- 0.1x10-4 mm2/s versus YFP,G93A-SOD1 mice = 2.7 +/- 0.5x10-4 mm2/s) (p < 0001) (+35 %) in the dorsal region (YFP mice = 2.5 +/- 0.1x10-4 mm2/s versus YFP,G93A-SOD1 mice = 3.1 +/- 0.1x10-4 mm2/s) (p < 0001) (+24.8%) and cervical levels (YFP mice = 2.2 +/- 0.2x10-4 mm2 /s versus YFP, G93A-SOD1 mice = 2.52 +/- 0.9x10-4 mm2/s) (p<001) (+16.1 %). Overall, these DTI findings demonstrated that lumbar changes in diffusion across axonal structures in this ALS mice can be detected before symptoms manifest (Fig. 4a).
Presymptomatic changes in white matter diffusion are linked to axonal structural anomalies
Alterations in axonal features from spinal cords (SC) in G93A-SOD1 mice associated to a fluorescent reporter has been previously reported . However, the results from these studies only carried descriptive findings and no specific quantitative analysis has been done. To determine the structural changes underlying the WM alterations in diffusion observed in previous MRI scans, we focused our histological analysis on the anterolateral region at the lumbar SC section from YFP and YFP, G93A-SOD1 mice (Fig. 2a). To make results equivalent to the voxel size performed during the MRI sessions (voxel size = 100 x100 μm2) we performed structural analysis using confocal fluorescence images with similar size. Specifically, for morphological evaluations we manually registered WM ROIs for each animal and group following similar topography and guided by stereotaxic coordinates . Measurements from each groups showed a significant (p < 0.01) reduction (-22.1 %) in axonal areas in the YFP, G93A-SOD1 mice at P80 (YFP mice = 2.9 +/- 0.06 μm2 versus. YFP, G93A-SOD1 mice= 2.3 +/-0.05 μm2) and -41.9% P120 (YFP mice = 3.69 +/-0.09 μm2 versus YFP, G93A-SOD1 mice = 2.4 +/-0.05 μm2) (n > 2000 axons per group). Nonetheless, we also observed a significant (p < 0.05) increase in the number YFP-positive axons delimited within each ROI (Axonal Density) in the YFP, G93ASOD1 mice group (+21.6 %) at P80 (YFP mice 1024 +/- 61 axons/100μm2 vs YFP, G93A-SOD1 mice 1272 +/- 48 axons/100μm2) and more significant (p < 0.01) increase (+42.4 %) at P120 (Fig. 2b). All together, these results point towards a significant remodeling of the SC axonal structures at early stage of the disease.
Presymptomatic structural changes in axonal fibers are associated with alterations in axonal connectivity
Using tractography reconstructions and histological confocal reconstructions we have demonstrated a critical impact of ALS towards axonal connectivity. Specifically, tractography methods showed early anomalies in fiber organization in the YFP, G93A-SOD1 mice and histological reconstructions evaluated the specific structural anomalies across individual axons (Fig. 3a and b). Among new techniques has been developed for quantitative connectome analysis , edges weights (EW) analysis has been introducing the concept of the adjustment and quantitation of connection strengths  as calculated from our data set (Fig. 4b). Results from this analysis showed an early significant decrease of this parameters (p < 0.01) at P80, at the cervical levels (YFP mice = 0.024 +/- 0.002 versus YFP, G93A-SOD1 mice = 0.0136 +/- 0.004) (p < 0.05) (-43.5 %), as well as in the SC lumbar regions (YFP mice = 0.039 +/- 0.0001 versus YFP, G93A-SOD1 mice = 0.023 +/- 0.0004) (p<0.01) (-40.9 %), (Fig. 4c).
Early ADC anomalies in the spinal cord of the ALS mice are coupled with changes in water transport protein membranes
The role of ADC in SC ALS tissue has been described before but it changes in the early stages of the disease remain unknown [16, 37]. Measurements from our experimental model have shown early changes of this parameter particularly at the thoracic (YFP mice = 3.2 +/- 0.1x10-4 mm2 /s versus YFP, G93A-SOD1 mice = 3.8 +/- 0.1x10-4 mm2/s) (p<001), (+16.4 %) and lumbar level (YFP mice = 3.2 +/- 0.1x10-4 mm2 /s versus YFP, G93A-SOD1 mice = 3.53 +/- 0.9x10-4 mm2/s) (p<001) (+23.3 %) (Fig. 5a). The complex nature of diffusional water exchange in diverse neuropathological conditions is associated to changes in membrane permeability [38, 39]. One of the factors regulating this water diffusion process is determined by the transmembrane proteins (aquaporin channels) . Aquaporin 4 (AQP4) has shown a high prevalence in the central nervous system (CNS). Quantitative IHC analysis revealed a significant (p < 0.001) increase of AQP4 in the ALS mice SC at the presymptomatic stage (P80) (YFP mice = 541.6 +/- 131.5 a.u. vs YFP, G93A-SOD1 mice = 5403 +/-176 a.u.) (1-2 fold increase) predominantly in the extra-axonal compartment (Fig. 5b and c).
ALS is characterized by a selective loss of motor neurons in the brain and spinal cord . Although the underlying mechanism of the disease is unknown, pathological observations and experimental data establish that alterations in synaptic and neuronal function occur well before neuronal death , supporting the principal role of axonal degeneration in the neuropathological process . The earlier compromise of motor axons, and alterations in water diffusion in the spinal cords from an animal model of ALS have been increasingly gaining attention [42,43,44]. An unproven hypothesis of early axonal structural changes has been proposed to explain the preponderance of early phenotype and symptoms in ALS and other neurodegenerative diseases [10, 45]. Furthermore, our previous MRI studies demonstrated that early axonal injury in the G93A-SOD1 mice model occurred at the lower level of the spinal cord and that such microstructural changes could be monitored using DTI . Moreover, in vivo UHF-MRI diffusion studies have a critical advantage for longitudinal studies acquired with high spatial resolution . From the medical perspective, using DTI to gain an understanding of the ultrastructural changes occurring during the early stage of ALS may improve the detection and treatment of this disease at earlier stages [10, 47].
Previous work G93ASOD1 mice were not able to demonstrate presymptomatic changes in DTI diffusion parameters using 7T and 9.4T MRI instruments in combination with conventional histology [18, 48,49,50,51,52] (Table 1). To the best of our knowledge, this study is the first to use UHF-MRI at 17.6T to interrogate the WM microstructure of the spinal cord from an ALS mouse model. The combination of higher signal to noise ratios (SNR) available at higher magnetic fields [53, 54], and the histological detail provided by an ALS transgenic mice expressing a neuronal specific endogenous fluorescent protein (YFP), enhanced the detection of early changes in DTI parameters and WM microstructure (Figs. 1b and 4a). Thus, we have identified presymptomatic changes in MRI diffusion are based on alterations
in of morphological axonal features, such as a reduction in axonal areas and increase in axonal density (Fig. 2).
Besides the increased number of mechanisms responsible for the alteration of axonal function in ALS [3, 55, 56], the microstructural alterations causing the changes the diffusion signal in ALS are still not well-known. Changes in axonal are widely used in diverse scientific and clinical scenarios [57,58,59]. Nonetheless, such tractography features are an under representation of the real neuropathological changes occurring inside a single voxel (Fig. 3). Thus, such non-quantitative parameters have to be perfected and validated before can be used to make further decisions in clinical scenarios [60,61,62,63].
Although the micro-anatomical anomalies displayed in this work could explain the early changes in diffusion in ALS, the analysis of axonal connectivity could unveil a further insight in the real axonal disconnection problem observed in many neurodegenerative diseases [57, 59]. Obtained by diffusion path probabilities, such techniques can estimate the connection strengths across different WM regions (nodes) enabling the reconstruction of connectomes, consistent with CNS networks properties, mapped by previous imaging modalities and post-mortem brain studies . Using quantitative techniques measuring the diffusion properties along edge voxels, such as edge weight (EW) analysis [26, 36], we evaluated the early axonal disconnection in this animal model (Fig. 4b). Specifically, we addressed how structural information across spinal cord axons was linked and early impacted by ALS. These results have pointed to an early disconnection predominantly in cervical and lumbar regions, as observed in two thirds of ALS patients with upper and lower limbs symptoms (spinal form) .
In addition to the microstructural changes observed during early stages, another point of attention are the events related to alteration in cellular and membrane water exchange and their relationship with changes in MRI diffusion signals [66,67,68]. One of the mechanisms to maintain the osmotic balance across cellular membranes in the central nervous system is regulated by an intricate mechanism of water permeably controlled by membrane channels called Aquaporin (AQP) [69,70,71,72,73]. Although many structural alterations can be the cause of changes in water diffusion across cellular membranes, recent work has shown a significant link between AQPs and the MRI apparent diffusion coefficient (ADC) . Although such parameter has been found highly impaired in ALS models , changes in ADC has not be fully characterized in early stages of ALS. In that regard, our studies have demonstrated that changes in AQP4 and ADC are present at earlier stages of the disease (Fig. 5a and b). Moreover, based in our histological finding such increase in AQP4 channel expression could reveal additional anomalies across different cellular types located in the non- axonal compartment (Fig. 5c). Nevertheless, the specific role of AQP channels and their contribution to the water exchange and MRI diffusion in the early stages of ALS remains to be determined [75, 76].
Overall, this work proposes new insights to understand the biological nature of changes in MRI diffusion parameters during early stages of ALS (Fig. 6). Although previous MRI work has demonstrated a decrease in SC volume in ALS patients at symptomatic stages [77, 78], one of the limitations by our current MRI and histological techniques of axonal quantification is the exclusive assessment of ROIs without accounting for the volume reduction of the entire CNS structure (SC) interrogated. Moreover, recent findings using alternative diffusion techniques have demonstrated a similar increase in axonal densities during early stages of neurodegenerative diseases  pointing to an increase in the total number of axons per volume unit of WM tissues at early stage among different neurodegenerative diseases. Among the different mechanisms to explain the early structural changes, one is related to the specific susceptibility of subsets of neurons by the mutated G93A-SOD1 protein, producing an axonal dying-back mechanism [3, 41, 80]. Another theory is the relative impairment in the glial cell population , and both are considered the subject of future DTI in vivo studies.
Another limitation in the co-registration of UHF-MRI and histological techniques is the mismatch between the resolution of the MRI signal (100 μm) and the maximum optical resolution given by the laser confocal systems and fluorescent probe used in our biological preparations (1/2 the wavelength of the excitation light, where YFP = 0.534 microns). Hence, results from a single reconstructed tract from DTI data (Fig. 3a) represents an area of approx. 100-200 axons/voxel (Fig. 3b). Moreover, considering the overall WM axonal sizes (1-20 microns) and the mosaic expression of the YFP tag in this particular fluorescent animal (only 10-30% of the neuronal population) the current model clearly underrepresent the real number of axonal elements per voxel area (approx. 400-1000 axons/voxel).
Our studies have shown that UHF-MRI and DTI techniques are useful tools to detect early changes in the microstructural features of ALS (Table 2). Yet, our current techniques have some limitations due to the small range of b-values studied. In addition, the involvement of different cell populations as described in this work points towards the existence of an anomalous water diffusion process in a porous biological material marking the need for higher b-values and the use of more complex diffusion models [82,83,84] to evaluate the complex WM neuroanatomical changes occurring during the development of ALS . Finally, more extensive neurobiological work should be done in order to validate each bioimaging marker during the early stages of this disease.
This work demonstrates that the combination of UHF-MRI and fluorescent mouse reporters is a useful approach to detect and characterize microstructural changes and axonal connectivity anomalies in ALS. Using a transgenic mouse fluorescent reporter, our studies were able to identify changes in axonal size, density and connectivity in the spinal cord of ALS mice before the disease was fully manifest. We found that anomalies in diffusion markers captured by DTI are possibly connected with alterations in cellular and molecular markers observed in our fluorescent ALS mice. We described early WM changes using connectomics as an additional method to evaluate early axonal connectivity anomalies in ALS. We believe this investigation is a step in the characterization of early bioimaging markers to detect and eventually monitor future treatments in patients with ALS.
Amyotrophic lateral Sclerosis
Diffusion Tensor Imaging
transgenic mice with the overexpression of human mutant gene copper zinc superoxide dismutase identified in familiar forms of ALS patients
ratio of the inner axonal diameter to the total outer diameter
Ultra-High Field MRI
Weight edge connectomics
Mehta P, Antao V, Kaye W, Sanchez M, Williamson D, Bryan L, Muravov O, Horton K, Division of T, Human Health Sciences AfTS, et al. Prevalence of amyotrophic lateral sclerosis - United States, 2010-2011. MMWR Suppl. 2014;63:1–14.
de Oliveira GP, Alves CJ, Chadi G. Early gene expression changes in spinal cord from SOD1(G93A) Amyotrophic Lateral Sclerosis animal model. Front Cell Neurosci. 2013;7:216.
Morfini GA, Bosco DA, Brown H, Gatto R, Kaminska A, Song Y, Molla L, Baker L, Marangoni MN, Berth S, et al. Inhibition of fast axonal transport by pathogenic SOD1 involves activation of p38 MAP kinase. PLoS One. 2013;8:e65235.
Gurney ME, Pu H, Chiu AY, Dal Canto MC, Polchow CY, Alexander DD, Caliendo J, Hentati A, Kwon YW, Deng HX, et al. Motor neuron degeneration in mice that express a human Cu,Zn superoxide dismutase mutation. Science. 1994;264:1772–5.
Fogarty MJ, Klenowski PM, Lee JD, Drieberg-Thompson JR, Bartlett SE, Ngo ST, Hilliard MA, Bellingham MC, Noakes PG. Cortical synaptic and dendritic spine abnormalities in a presymptomatic TDP-43 model of amyotrophic lateral sclerosis. Sci Rep. 2016;6:37968.
Liu Y, Pattamatta A, Zu T, Reid T, Bardhi O, Borchelt DR, Yachnis AT, Ranum LP. C9orf72 BAC Mouse Model with Motor Deficits and Neurodegenerative Features of ALS/FTD. Neuron. 2016;90:521–34.
Nolan M, Talbot K, Ansorge O. Pathogenesis of FUS-associated ALS and FTD: insights from rodent models. Acta Neuropathol Commun. 2016;4:99.
Oglesby E, Quigley HA, Zack DJ, Cone FE, Steinhart MR, Tian J, Pease ME, Kalesnykas G. Semi-automated, quantitative analysis of retinal ganglion cell morphology in mice selectively expressing yellow fluorescent protein. Exp Eye Res. 2012;96:107–15.
Bannerman PG, Hahn A. Enhanced visualization of axonopathy in EAE using thy1-YFP transgenic mice. J Neurol Sci. 2007;260:23–32.
Gatto RG, Chu Y, Ye AQ, Price SD, Tavassoli E, Buenaventura A, Brady ST, Magin RL, Kordower JH, Morfini GA. Analysis of YFP(J16)-R6/2 reporter mice and postmortem brains reveals early pathology and increased vulnerability of callosal axons in Huntington's disease. Hum Mol Genet. 2015;24:5285–98.
Porrero C, Rubio-Garrido P, Avendano C, Clasca F. Mapping of fluorescent protein-expressing neurons and axon pathways in adult and developing Thy1-eYFP-H transgenic mice. Brain Res. 2010;1345:59–72.
King AE, Blizzard CA, Southam KA, Vickers JC, Dickson TC. Degeneration of axons in spinal white matter in G93A mSOD1 mouse characterized by NFL and alpha-internexin immunoreactivity. Brain Res. 2012;1465:90–100.
Wang Y, Liu L, Ma L, Huang X, Lou X, Wang Y, Wu N, Liu T, Guo X. Preliminary study on cervical spinal cord in patients with amyotrophic lateral sclerosis using MR diffusion tensor imaging. Acad Radiol. 2014;21:590–6.
Nair G, Carew JD, Usher S, Lu D, Hu XP, Benatar M. Diffusion tensor imaging reveals regional differences in the cervical spinal cord in amyotrophic lateral sclerosis. Neuroimage. 2010;53:576–83.
Underwood CK, Kurniawan ND, Butler TJ, Cowin GJ, Wallace RH. Non-invasive diffusion tensor imaging detects white matter degeneration in the spinal cord of a mouse model of amyotrophic lateral sclerosis. Neuroimage. 2011;55:455–61.
Evans MC, Serres S, Khrapitchev AA, Stolp HB, Anthony DC, Talbot K, Turner MR, Sibson NR. T(2)-weighted MRI detects presymptomatic pathology in the SOD1 mouse model of ALS. J Cereb Blood Flow Metab. 2014;34:785–93.
Ong HH, Wehrli FW. Quantifying axon diameter and intra-cellular volume fraction in excised mouse spinal cord with q-space imaging. Neuroimage. 2010;51:1360–6.
Marcuzzo S, Bonanno S, Figini M, Scotti A, Zucca I, Minati L, Riva N, Domi T, Fossaghi A, Quattrini A, et al. A longitudinal DTI and histological study of the spinal cord reveals early pathological alterations in G93A-SOD1 mouse model of amyotrophic lateral sclerosis. Exp Neurol. 2017;293:43–52.
Boillee S, Vande Velde C, Cleveland DW. ALS: a disease of motor neurons and their nonneuronal neighbors. Neuron. 2006;52:39–59.
Figini M, Scotti A, Marcuzzo S, Bonanno S, Padelli F, Moreno-Manzano V, Garcia-Verdugo JM, Bernasconi P, Mantegazza R, Bruzzone MG, Zucca I. Comparison of Diffusion MRI Acquisition Protocols for the In Vivo Characterization of the Mouse Spinal Cord: Variability Analysis and Application to an Amyotrophic Lateral Sclerosis Model. PLoS One. 2016;11:e0161646.
Kim JH, Wu TH, Budde MD, Lee JM, Song SK. Noninvasive detection of brainstem and spinal cord axonal degeneration in an amyotrophic lateral sclerosis mouse model. NMR Biomed. 2011;24:163–9.
Ozdinler PH, Benn S, Yamamoto TH, Guzel M, Brown RH Jr, Macklis JD. Corticospinal motor neurons and related subcerebral projection neurons undergo early and specific neurodegeneration in hSOD1G(9)(3)A transgenic ALS mice. J Neurosci. 2011;31:4166–77.
Saberi S, Stauffer JE, Schulte DJ, Ravits J. Neuropathology of Amyotrophic Lateral Sclerosis and Its Variants. Neurol Clin. 2015;33:855–76.
Heiman-Patterson TD, Deitch JS, Blankenhorn EP, Erwin KL, Perreault MJ, Alexander BK, Byers N, Toman I, Alexander GM. Background and gender effects on survival in the TgN(SOD1-G93A)1Gur mouse model of ALS. J Neurol Sci. 2005;236:1–7.
Durand J, Amendola J, Bories C, Lamotte d'Incamps B. Early abnormalities in transgenic mouse models of amyotrophic lateral sclerosis. J Physiol Paris. 2006;99:211–20.
Colon-Perez LM, Spindler C, Goicochea S, Triplett W, Parekh M, Montie E, Carney PR, Price C, Mareci TH. Dimensionless, Scale Invariant, Edge Weight Metric for the Study of Complex Structural Networks. PLoS One. 2015;10:e0131493.
Jones DK. The effect of gradient sampling schemes on measures derived from diffusion tensor MRI: a Monte Carlo study. Magn Reson Med. 2004;51:807–15.
Jenkinson M, Beckmann CF, Behrens TE, Woolrich MW, Smith SM. Fsl. Neuroimage. 2012;62:782–90.
Basser PJ, Mattiello J, LeBihan D. MR diffusion tensor spectroscopy and imaging. Biophys J. 1994;66:259–67.
Callaghan P, Eccles CD, Xia Y. NMR microscopy of dynamic displacements: k-space and q-space imaging. Journal of Physics E: Scientific Instruments. 1988;21:820–2.
Gatto RG, Li W, Magin RL. Diffusion tensor imaging identifies presymptomatic axonal degeneration in the spinal cord of ALS mice. Brain Res. 2018;1679:7.
Tallon C, Russell KA, Sakhalkar S, Andrapallayal N, Farah MH. Length-dependent axo-terminal degeneration at the neuromuscular synapses of type II muscle in SOD1 mice. Neuroscience. 2015;312:179–89.
Agosta F, Rocca MA, Valsasina P, Sala S, Caputo D, Perini M, Salvi F, Prelle A, Filippi M. A longitudinal diffusion tensor MRI study of the cervical cord and brain in amyotrophic lateral sclerosis patients. J Neurol Neurosurg Psychiatry. 2009;80:53–5.
Watson CP, G.; Kayalioglu, G.: The Spinal Cord : A Christopher and Dana Reeve Foundation Text and Atlas. 1st Edition edn. London: Academic Press 2008; 2008.
Preti MG, Bolton TA, Van De Ville D. The dynamic functional connectome: State-of-the-art and perspectives. Neuroimage. 2016;160:41–54.
Colon-Perez LM, Couret M, Triplett W, Price CC, Mareci TH. Small Worldness in Dense and Weighted Connectomes. Front Phys. 2016;4
Niessen HG, Angenstein F, Sander K, Kunz WS, Teuchert M, Ludolph AC, Heinze HJ, Scheich H, Vielhaber S. In vivo quantification of spinal and bulbar motor neuron degeneration in the G93A-SOD1 transgenic mouse model of ALS by T2 relaxation time and apparent diffusion coefficient. Exp Neurol. 2006;201:293–300.
Ye AQ, Hubbard Cristinacce PL, Zhou FL, Yin Z, Parker GJ, Magin RL. Diffusion tensor MRI phantom exhibits anomalous diffusion. Conf Proc IEEE Eng Med Biol Soc. 2014;2014:746–9.
Nilsson M, Latt J, van Westen D, Brockstedt S, Lasic S, Stahlberg F, Topgaard D. Noninvasive mapping of water diffusional exchange in the human brain using filter-exchange imaging. Magn Reson Med. 2013;69:1573–81.
Agre P. The aquaporin water channels. Proc Am Thorac Soc. 2006;3:5–13.
Fischer LR, Culver DG, Tennant P, Davis AA, Wang M, Castellano-Sanchez A, Khan J, Polak MA, Glass JD. Amyotrophic lateral sclerosis is a distal axonopathy: evidence in mice and man. Exp Neurol. 2004;185:232–40.
Kawamura Y, Dyck PJ, Shimono M, Okazaki H, Tateishi J, Doi H. Morphometric comparison of the vulnerability of peripheral motor and sensory neurons in amyotrophic lateral sclerosis. J Neuropathol Exp Neurol. 1981;40:667–75.
Bradley WG, Good P, Rasool CG, Adelman LS. Morphometric and biochemical studies of peripheral nerves in amyotrophic lateral sclerosis. Ann Neurol. 1983;14:267–77.
Hammad M, Silva A, Glass J, Sladky JT, Benatar M. Clinical, electrophysiologic, and pathologic evidence for sensory abnormalities in ALS. Neurology. 2007;69:2236–42.
Murmu RP, Li W, Holtmaat A, Li JY. Dendritic spine instability leads to progressive neocortical spine loss in a mouse model of Huntington's disease. J Neurosci. 2013;33:12997–3009.
Brennan FH, Cowin GJ, Kurniawan ND, Ruitenberg MJ. Longitudinal assessment of white matter pathology in the injured mouse spinal cord through ultra-high field (16.4 T) in vivo diffusion tensor imaging. Neuroimage. 2013;82:574–85.
Kamagata K, Kerever A, Yokosawa S, Otake Y, Ochi H, Hori M, Kamiya K, Tsuruta K, Tagawa K, Okazawa H, et al. Quantitative Histological Validation of Diffusion Tensor MRI with Two-Photon Microscopy of Cleared Mouse Brain. Magn Reson Med Sci. 2016;15:416–21.
Caron I, Micotti E, Paladini A, Merlino G, Plebani L, Forloni G, Modo M, Bendotti C. Comparative Magnetic Resonance Imaging and Histopathological Correlates in Two SOD1 Transgenic Mouse Models of Amyotrophic Lateral Sclerosis. PLoS One. 2015;10:e0132159.
El Mendili MM, Cohen-Adad J, Pelegrini-Issac M, Rossignol S, Morizot-Koutlidis R, Marchand-Pauvert V, Iglesias C, Sangari S, Katz R, Lehericy S, et al. Multi-parametric spinal cord MRI as potential progression marker in amyotrophic lateral sclerosis. PLoS One. 2014;9:e95516.
Cohen-Adad J, El Mendili MM, Morizot-Koutlidis R, Lehericy S, Meininger V, Blancho S, Rossignol S, Benali H, Pradat PF. Involvement of spinal sensory pathway in ALS and specificity of cord atrophy to lower motor neuron degeneration. Amyotroph Lateral Scler Frontotemporal Degener. 2013;14:30–8.
Valsasina P, Agosta F, Benedetti B, Caputo D, Perini M, Salvi F, Prelle A, Filippi M. Diffusion anisotropy of the cervical cord is strictly associated with disability in amyotrophic lateral sclerosis. J Neurol Neurosurg Psychiatry. 2007;78:480–4.
Rossi C, Boss A, Lindig TM, Martirosian P, Steidle G, Maetzler W, Claussen CD, Klose U, Schick F. Diffusion tensor imaging of the spinal cord at 1.5 and 3.0 Tesla. Rofo. 2007;179:219–24.
Behr VC, Weber T, Neuberger T, Vroemen M, Weidner N, Bogdahn U, Haase A, Jakob PM, Faber C. High-resolution MR imaging of the rat spinal cord in vivo in a wide-bore magnet at 17.6 Tesla. MAGMA. 2004;17:353–8.
Pallebage-Gamarallage M, Foxley S, Menke RAL, Huszar IN, Jenkinson M, Tendler BC, Wang C, Jbabdi S, Turner MR, Miller KL, Ansorge O. Dissecting the pathobiology of altered MRI signal in amyotrophic lateral sclerosis: A post mortem whole brain sampling strategy for the integration of ultra-high-field MRI and quantitative neuropathology. BMC Neurosci. 2018;19:11.
Shang Y, Huang EJ. Mechanisms of FUS mutations in familial amyotrophic lateral sclerosis. Brain Res. 2016;1647:65–78.
Sharma A, Lyashchenko AK, Lu L, Nasrabady SE, Elmaleh M, Mendelsohn M, Nemes A, Tapia JC, Mentis GZ, Shneider NA. ALS-associated mutant FUS induces selective motor neuron degeneration through toxic gain of function. Nat Commun. 2016;7:10465.
Agosta F, Pagani E, Petrolini M, Caputo D, Perini M, Prelle A, Salvi F, Filippi M. Assessment of white matter tract damage in patients with amyotrophic lateral sclerosis: a diffusion tensor MR imaging tractography study. AJNR Am J Neuroradiol. 2010;31:1457–61.
Iwata NK, Aoki S, Okabe S, Arai N, Terao Y, Kwak S, Abe O, Kanazawa I, Tsuji S, Ugawa Y. Evaluation of corticospinal tracts in ALS with diffusion tensor MRI and brainstem stimulation. Neurology. 2008;70:528–32.
Steinbach R, Loewe K, Kaufmann J, Machts J, Kollewe K, Petri S, Dengler R, Heinze HJ, Vielhaber S, Schoenfeld MA, Stoppel CM. Structural hallmarks of amyotrophic lateral sclerosis progression revealed by probabilistic fiber tractography. J Neurol. 2015;262:2257–70.
Campbell JS, Pike GB. Potential and limitations of diffusion MRI tractography for the study of language. Brain Lang. 2014;131:65–73.
Thomas C, Ye FQ, Irfanoglu MO, Modi P, Saleem KS, Leopold DA, Pierpaoli C. Anatomical accuracy of brain connections derived from diffusion MRI tractography is inherently limited. Proc Natl Acad Sci U S A. 2014;111:16574–9.
Daducci A, Dal Palu A, Descoteaux M, Thiran JP. Microstructure Informed Tractography: Pitfalls and Open Challenges. Front Neurosci. 2016;10:247.
Mukherjee P, Chung SW, Berman JI, Hess CP, Henry RG. Diffusion tensor MR imaging and fiber tractography: technical considerations. AJNR Am J Neuroradiol. 2008;29:843–52.
Sotiropoulos SN, Zalesky A. Building connectomes using diffusion MRI: why, how and but. NMR Biomed. 2017;
Wijesekera LC, Leigh PN. Amyotrophic lateral sclerosis. Orphanet J Rare Dis. 2009;4:3.
Obata T, Kershaw J, Kuroiwa, D.; Shibata, S.; , Yoichiro, A.; , Yasui, M.; Aoki I.: Effect of Cell Membrane Water Permeability on Diffusion-Weighted MR Signal: A Study Using Expression-Controlled Aquaporin4 Cells. In International Society for Magnetic Resonance in Medicine 20th Annual Meeting Melbourne, Australia. Proc. Intl. Soc. Mag. Reson. Med; 2012: 1830.
Mukherjee A, Wu D, Davis HC, Shapiro MG. Non-invasive imaging using reporter genes altering cellular water permeability. Nat Commun. 2016;7:13891.
Beaulieu C. The basis of anisotropic water diffusion in the nervous system - a technical review. NMR Biomed. 2002;15:435–55.
Amiry-Moghaddam M, Ottersen OP. The molecular basis of water transport in the brain. Nat Rev Neurosci. 2003;4:991–1001.
Foglio E, Rodella LF. Aquaporins and neurodegenerative diseases. Curr Neuropharmacol. 2010;8:112–21.
Oklinski MK, Skowronski MT, Skowronska A, Rutzler M, Norgaard K, Nieland JD, Kwon TH, Nielsen S. Aquaporins in the Spinal Cord. Int J Mol Sci. 2016;17:2050.
Vitellaro-Zuccarello L, Mazzetti S, Bosisio P, Monti C, De Biasi S. Distribution of Aquaporin 4 in rodent spinal cord: relationship with astrocyte markers and chondroitin sulfate proteoglycans. Glia. 2005;51:148–59.
Oshio K, Binder DK, Yang B, Schecter S, Verkman AS, Manley GT. Expression of aquaporin water channels in mouse spinal cord. Neuroscience. 2004;127:685–93.
Nicaise C, Soyfoo MS, Authelet M, De Decker R, Bataveljic D, Delporte C, Pochet R. Aquaporin-4 overexpression in rat ALS model. Anat Rec (Hoboken). 2009;292:207–13.
Dai J, Lin W, Zheng M, Liu Q, He B, Luo C, Lu X, Pei Z, Su H, Yao X. Alterations in AQP4 expression and polarization in the course of motor neuron degeneration in SOD1G93A mice. Mol Med Rep. 2017;16:1739–46.
Bataveljic D, Nikolic L, Milosevic M, Todorovic N, Andjus PR. Changes in the astrocytic aquaporin-4 and inwardly rectifying potassium channel expression in the brain of the amyotrophic lateral sclerosis SOD1(G93A) rat model. Glia. 2012;60:1991–2003.
de Albuquerque M, Branco LM, Rezende TJ, de Andrade HM, Nucci A, Franca MC Jr. Longitudinal evaluation of cerebral and spinal cord damage in Amyotrophic Lateral Sclerosis. Neuroimage Clin. 2017;14:269–76.
Grolez G, Kyheng M, Lopes R, Moreau C, Timmerman K, Auger F, Kuchcinski G, Duhamel A, Jissendi-Tchofo P, Besson P, et al. MRI of the cervical spinal cord predicts respiratory dysfunction in ALS. Sci Rep. 2018;8:1828.
Benitez A, Jensen JH, Helpern JA. Axonal Density and Myelin Integrity in Cognitice Decline: A diffusional Kurtosis Imaging Study. In: Alzheimer's Association International Conference; July 17, vol. 2017. London: Alzheimer & Dementia; 2018. p. P774–5.
Dadon-Nachum M, Melamed E, Offen D. The "dying-back" phenomenon of motor neurons in ALS. J Mol Neurosci. 2011;43:470–7.
Lasiene J, Yamanaka K. Glial cells in amyotrophic lateral sclerosis. Neurol Res Int. 2011;2011:718987.
Liang Y, Ye AQ, Chen W, Gatto RG, Colon-Perez L, Mareci TH, Magin RL. A fractal derivative model for the characterization of anomalous diffusion in magnetic resonance imaging. Commun Nonlinear Sci Numer Simul. 2016;39:529–37.
Magin RL, Akpa BS, Neuberger T, Webb AG. Fractional Order Analysis of Sephadex Gel Structures: NMR Measurements Reflecting Anomalous Diffusion. Commun Nonlinear Sci Numer Simul. 2011;16:4581–7.
Magin RL, Ingo C, Colon-Perez L, Triplett W, Mareci TH. Characterization of Anomalous Diffusion in Porous Biological Tissues Using Fractional Order Derivatives and Entropy. Microporous Mesoporous Mater. 2013;178:39–43.
Fukui Y, Hishikawa N, Sato K, Nakano Y, Morihara R, Shang J, Takemoto M, Ohta Y, Yamashita T, Abe K. Detecting spinal pyramidal tract of amyotrophic lateral sclerosis patients with diffusion tensor tractography. Neurosci Res. 2017;
Rasoanandrianina H, Grapperon AM, Taso M, Girard OM, Duhamel G, Guye M, Ranjeva JP, Attarian S, Verschueren A, Callot V. Region-specific impairment of the cervical spinal cord (SC) in amyotrophic lateral sclerosis: A preliminary study using SC templates and quantitative MRI (diffusion tensor imaging/inhomogeneous magnetization transfer). NMR Biomed. 2017;30
Budrewicz S, Szewczyk P, Bladowska J, Podemski R, Koziorowska-Gawron E, Ejma M, Slotwinski K, Koszewicz M. The possible meaning of fractional anisotropy measurement of the cervical spinal cord in correct diagnosis of amyotrophic lateral sclerosis. Neurol Sci. 2016;37:417–21.
Iglesias C, Sangari S, El Mendili MM, Benali H, Marchand-Pauvert V, Pradat PF. Electrophysiological and spinal imaging evidences for sensory dysfunction in amyotrophic lateral sclerosis. BMJ Open. 2015;5:e007659.
Romano A, Guo J, Prokscha T, Meyer T, Hirsch S, Braun J, Sack I, Scheel M. In vivo waveguide elastography: effects of neurodegeneration in patients with amyotrophic lateral sclerosis. Magn Reson Med. 2014;72:1755–61.
We would like to especially acknowledge Dr. Gerardo Morfini for providing chemicals and materials used in our experiments and Ehsan Tavassoli for proofreading the manuscript.
This study was supported in part by a Chicago Biomedical Consortium (CBC) postdoctoral fellowship grant (Award #085740) to RG at the University of Illinois in Chicago. Data collection was supported by the Magnetic Laboratory Visiting Scientist Program (Award VSP #278) of the National High Magnetic Field Laboratory (NHMFL) and Advanced Magnetic Resonance Imaging and Spectroscopy (AMRIS) to RG. A portion of this work was performed in the McKnight Brain Institute at the National High.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
This study was carried out in accordance with the recommendations of the Animal Care Committee (ACC) guidelines at the University of Illinois at Chicago. The protocol was approved by the Institutional ACC at the University of Illinois at Chicago.
Consent for publication
The authors declare that they have no competing interests.
About this article
- Amyotrophic Lateral Sclerosis
- Spinal Cord
- Ultra-high Field MRI
- Diffusion Tensor Imaging
- Yellow Fluorescent Protein
- G93A-SOD1 mice
- Axonal Degeneration