Background Corticospinal tract (CST) degeneration and cortical atrophy are consistent features of amyotrophic lateral sclerosis (ALS). We hypothesised that neurite orientation dispersion and density imaging (NODDI), a multicompartment model of diffusion MRI, would reveal microstructural changes associated with ALS within the CST and precentral gyrus (PCG) ‘in vivo’.
Methods 23 participants with sporadic ALS and 23 healthy controls underwent diffusion MRI. Neurite density index (NDI), orientation dispersion index (ODI) and free water fraction (isotropic compartment (ISO)) were derived. Whole brain voxel-wise analysis was performed to assess for group differences. Standard diffusion tensor imaging (DTI) parameters were computed for comparison. Subgroup analysis was performed to investigate for NODDI parameter differences relating to bulbar involvement. Correlation of NODDI parameters with clinical variables were also explored. The results were accepted as significant where p<0.05 after family-wise error correction at the cluster level, clusters formed with p<0.001.
Results In the ALS group NDI was reduced in the extensive regions of the CST, the corpus callosum and the right PCG. ODI was reduced in the right anterior internal capsule and the right PCG. Significant differences in NDI were detected between subgroups stratified according to the presence or absence of bulbar involvement. ODI and ISO correlated with disease duration.
Conclusions NODDI demonstrates that axonal loss within the CST is a core feature of degeneration in ALS. This is the main factor contributing to the altered diffusivity profile detected using DTI. NODDI also identified dendritic alterations within the PCG, suggesting microstructural cortical dendritic changes occur together with CST axonal damage.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Statistics from Altmetric.com
In amyotrophic lateral sclerosis (ALS) degeneration of the corticospinal system, cranial and spinal anterior horn motor neurons is a consistent clinicopathological feature. However, the mechanisms underlying degeneration are not well understood.1 It is unclear whether degeneration of corticospinal neurons represents a ‘dying back’ or ‘dying forward’ process,2 or whether these categories proposed from histological (ie, end stage) analyses are relevant to the dynamics of cell damage and death in vivo. Understanding these processes is particularly relevant in the context of disease heterogeneity, since there may be common pathways linking the many genetic factors that contribute to pathogenesis.
In this context MRI has provided important insights into the localisation and extent, but not the fine detail, of cellular pathology in ALS. Arguably the most promising application has been diffusion tensor imaging (DTI), a quantitative MRI technique able to indirectly detect microstructural abnormalities due to its sensitivity to water molecule displacement within the tissue.3 4
However, DTI has limitations. First, changes to the DTI-derived indices of fractional anisotropy (FA) and mean diffusivity (MD) are non-specific.5 Second, their interpretation in areas of complex axonal or dendritic architecture is not straightforward.6 7 Although indices such as the mode of anisotropy can provide some useful information in regions of crossing or fanning,8 a quantification of neurite density in such areas is challenging with DTI.
Recent developments in diffusion MRI have addressed some of the limitations of standard DTI and advanced capacity to characterise tissue microstructure. Neurite morphology evaluated by diffusion MRI is correlated with histological analyses of neuropil in postmortem brain tissue9 and with the orientation of neurites in Golgi-impregnated brain tissue in vivo.10 Neurite orientation dispersion and density imaging (NODDI) provides a simplified yet sophisticated model of diffusion MRI, which separates the signal arising from three different tissue compartments: intraneurite water, extraneurite water and cerebrospinal fluid (CSF).11 12 This model, recently validated with histological analyses in multiple sclerosis spinal cord,13 facilitates the analysis of neurite morphology and focuses on quantifying the altered architecture of these multifaceted structures. NODDI estimates the density and fanning of neurites, as well as the partial volume contamination from CSF, using standard MRI scanners within a clinically feasible time frame.
NODDI has been employed in a range of imaging studies, investigating the normal changes observed with ageing,14 as well as alterations associated with neurodegenerative diseases,15 16 white matter (WM) reorganisation following stroke17 and assessment of focal cortical dysplasia.18 NODDI has also been used to assess WM abnormalities in first-episode psychosis.19
We used NODDI to analyse the mechanisms underlying the neurodegenerative changes in ALS.20 21 Our principal aim was to test the hypothesis that cortical pathology, such as loss or altered complexity of dendrites in the neuropil of the precentral gyrus (PCG), is occurring in ALS. We hypothesised that NODDI would detect reduced neurite density and orientation dispersion, accompanied by an increase in free water compartment. An additional aim was to clarify the microstructural changes within the WM, where we hypothesised that axonal degeneration within the corticospinal tract (CST) would be detected as a loss of neurite density using NODDI.
We used NODDI in parallel with DTI to test the hypothesis that, by removing the confounding effect of orientation dispersion, NODDI would detect more axonal damage within the WM tracts than conventional DTI.
NODDI parameters were used to assess whether a combination of bulbar and limb involvement clinically would reveal differences when compared with limb-confined disease, to test the hypothesis that more widespread clinical involvement would relate to more extensive microstructural damage. Finally correlation of the NODDI parameters with clinical measures of ALS severity was explored with the anticipation that as clinical measures of ALS severity reflect the underlying pathological processes, there may be correlation between these and the NODDI parameters.
This cross-sectional study using NODDI provides the first microstructural analysis of altered density and complexity of neurites in ALS. Analysing these characteristics and how they relate to the disease process ‘in vivo’ provides an exciting opportunity to advance the understanding of ALS pathogenesis.
Participants with ALS (23; 16 men and 7 women; median age 67 years, range 45–73 years) diagnosed with El Escorial definite, probable or laboratory-supported ALS were recruited from Brighton and Sussex University Hospitals NHS Trust. All had vital capacity (VC) values of over 60% predicted and capacity to give fully informed consent. Exclusion criteria were any contraindication to MRI or inability to tolerate lying flat for 1 hour due to either severe illness or respiratory compromise. Healthy controls (23; 14 men and 9 women, median age 64 years, range 43–76 years) were included in the study for comparison. Genetic screening for known ALS mutations was not carried out for ethical reasons.
Clinical assessments of participants with ALS
Participants with ALS underwent a standardised neurological examination, the Revised ALS Functional Rating Scale (ALSFRS-R)22 and the Edinburgh Cognitive and Behavioural ALS Screen (ECAS).23 The clinical assessment was used to stratify participants with ALS according to the site of ALS onset, disease distribution at the time of assessment, disease duration, Medical Research Council (MRC)24 manual muscle testing power score and upper motor neuron (UMN)3 score. The ALSFRS-R was used to calculate the monthly rate of change in ALSFRS-R (∆ALSFRS-R; reduction in ALSFRS-R score from 48 at baseline to the score at the time of assessment, divided by the number of months between the two scores), which we used as a substitute for the rate of progression. The ECAS was used to assess for the presence of any cognitive deficits related to the frontotemporal spectrum disorder (ALS-FTSD). These measures were all used to explore for associations between clinical phenotype and NODDI measures. All except for one of the participants with ALS were taking riluzole at the time of the clinical assessment and MRI.
MRI data were acquired on a 1.5T Siemens AVANTO scanner. The acquisition included the following:
Dual-echo turbo spin-echo (echo time (TE)=11 and 86 ms; repetition time (TR) =3040 ms; echo-train length=6; field of view (FoV)=240×210 mm2; matrix=256×224; slice thickness=5 mm).
Fast fluid-attenuated inversion recovery (TE=89 ms; TR=9720 ms; inversion time (TI)=2578; echo-train length=16; FoV=240×210 mm2; matrix=256×224; slice thickness=5 mm).
Volumetric high-resolution magnetization prepared rapid gradient echo (MPRAGE) (TE=3.57 ms; TR=2730 ms; TI=100 ms; flip angle=70; FoV=256×240 mm2; matrix=254×240; slice thickness=1 mm).
Multishell diffusion-weighted (DW) MRI was acquired with single-shot, twice-refocused pulse-gradient spin-echo (echo planar imaging; EPI),25 using 3 b-values (9 directions for b=300 smm−2, 30 directions with b=800 smm−2 and 60 diffusion directions with b=2400 smm−2), optimised for NODDI.11 12 Ten non-diffusion-weighted (b=0) volumes were acquired. The following were additional acquisition parameters: TE=99 ms; TR=8400; FoV=240×240 mm2; matrix=96×96; slice thickness=2.5 mm.
In addition, we also obtained resting-state functional MRI and quantitative magnetisation transfer data (reported elsewhere). The total acquisition time was approximately 1 hour, although the time required for the sequences detailed above was 30 min.
The DW images were first corrected for eddy-current distortions and for involuntary movement, using the following pipeline. A ‘within-b-value’ coregistration was performed, and the average b=0, b=300, b=800 and b=2400 images were created. Each of these averages was coregistered to the mean b=0 image to obtain the transformation matching them. Each DW volume was realigned to the mean DW image with the same b-value, and the final transformation matching each DW volume with the b=0 image was obtained by combining the matrices from either step. The b-matrices were rotated accordingly.26
The resulting images were skull-stripped using the FMRIB Software Library (FSL) brain extraction tool (http://fsl.fmrib.ox.ac.uk/fslwiki). Data were analysed using the software distributed by the developers of NODDI (http://www.nitrc.org/projects/noddi_toolbox), implemented in Matlab (V.2012b), which yielded maps of the neurite density index (NDI), orientation dispersion index (ODI) and free water volume fraction for each participant (isotropic compartment (ISO)).11 Camino (camino.org.uk) was used to fit the diffusion tensor to the same data using weighted least squares regression and to compute FA and MD maps.
The FA images were warped to common space to create a group-specific template, using Advanced Normalization Tools V.2.1.0 (http://stnava.github.io/ANTs). The resulting template was coregistered with the Montreal Neurological Institute (MNI) template from the FSL. The combination of these two transformations was applied to each NODDI parametric map, FA and MD. Three-dimensional Gaussian smoothing with full width at half maximum of 6 mm was performed in SPM12 (http://www.fil.ion.ucl.ac.uk/spm/).
Statistical Package for the Social Sciences (SPSS IBM V. 22 for Windows) was used for statistical analysis of the demographic data.
Voxel-wise comparisons of whole brain NDI, ODI, ISO, FA and MD maps between the ALS and control groups were performed within the framework of the general linear model of SPM12.
The ALS group was then divided into two subgroups based on the presence or absence of bulbar symptoms (11 bulbar plus limb vs 12 limb alone). Voxel-wise whole brain analysis of NDI, ODI and ISO maps between the ALS subgroups was performed using the general linear model framework within SPM12.
Finally significant associations between the NODDI parametric maps and all of the clinical measures of ALS disease severity were assessed using SPM12.
In all analyses in SPM12, the age of the participant was used as a covariate. The results were accepted as significant where p<0.05 after family-wise error correction at the cluster level, where clusters were formed using p<0.001.
The demographics of all participants and the clinical characteristics of participants with ALS are summarised in table 1. The two groups were well matched in terms of age (t=1.11, p=0.27) and gender (χ2=0.38, p=0.53).
Clinically the ALS group was phenotypically heterogeneous with a wide breadth of weakness, functional disability, rate of progression and cognitive deficit.
The majority of participants with ALS had the typical mixture of upper and lower motor neuron involvement. We recruited four ALS participants with lower motor neuron (LMN)-predominant ALS (two flail arm, two flail leg) and one participant with UMN-predominant phenotype.
None of the participants with ALS fulfilled the revised diagnostic criteria for ALS-Frontotemporal Dementia ALS-FTD.27 However five participants with ALS met the criteria for ALS-FTSD: one with ALS cognitive and behavioural impairment and four with ALS cognitive impairment. Six of the participants with ALS had abnormal subscores but did not meet the criteria for ALS cognitive impairment, and one patient had an isolated borderline language subscore.
Neurite orientation dispersion and density imaging
In the whole brain analysis of NODDI parameters, NDI was reduced in the ALS group bilaterally within a large continuous region of the CST extending from the posterior limb of the internal capsule rostrally through the corona radiata into the subcortical WM of the PCG (p<0.001) (figure 1 and table 2). The area of decreased NDI also extended into the transcallosal connection fibres between the primary motor cortices (figure 1A and table 2).
Diffusion tensor imaging
In the whole brain analysis of DTI parameters, FA was reduced in the ALS group bilaterally within the two regions of the CST: the pontine CST (p=0.006) and the superior corona radiata and subcortical WM of the PCG (right p<0.001, left p=0.002) (figure 2A and table 2).
MD was significantly increased in the ALS group within the internal capsule of the left CST (p=0.048) (figure 2B and table 2). No changes in MD were detected in the brainstem portion of the CST or the corpus callosum.
Subgroup analysis using the NODDI parameters
Whole brain analysis for differences between subgroups revealed that NDI was significantly reduced in the corona radiata and subcortical WM of the right hemisphere in the subgroup with both bulbar and limb involvement, compared with those with limb-only involvement (figure 3).
There were no significant differences in ODI or ISO between these ALS subgroups.
Correlation of NODDI parameters with disease severity
Using whole brain analysis of NODDI parameters, both ODI and ISO showed significant correlation with disease duration. ODI correlated with disease duration predominantly in the right PCG, but also with a smaller region in the medial aspect of the left PCG (see figure 4A). Further areas where ODI correlated significantly with disease duration were detected within the precuneus and the dorsolateral prefrontal cortex bilaterally. The ISO correlated with disease duration in the interhemispheric fissure, in close proximity to the paracentral lobule and precuneus, as well as the fourth ventricle (see figure 4B).
No significant correlations were found between the NODDI parameters and measures of disease severity: ALSFRS-R, ∆ALSFRS-R, UMN score, MRC power score and ECAS.
The major findings in this study using NODDI in ALS relate, firstly, to significant reduction of NDI in the CST and corpus callosum in ALS, suggesting loss of axonal density in these pathways. Secondly, the reduction of ODI within the PCG in the ALS group and the correlation of ODI with disease duration suggest a relationship between cortical pathology and the severity of the ALS disease process. Finally, our results are consistent with previous DTI findings; however, NODDI provides a more definitive interpretation of the microstructural changes underlying ALS by separating neurite density from orientation dispersion.
In relation to the WM changes, the area of reduced NDI extends from the posterior limb of the internal capsule and continues rostrally into the subcortical WM of the PCG. This supports the notion that loss of motor neuron axon density within the cerebral WM is a core feature of the neurodegenerative process associated with ALS, consistent with many pathological studies.28 29 It is possible that alterations in glial cells and myelin may also contribute to these changes in NDI. However, as the density of axons per voxel is much higher than that of glial cells, and loss of motor neuron axons is a key feature of ALS, we expect that the largest changes are explained by axonal loss.
Regarding reductions in the NDI and ODI in the PCG, we interpret these changes as compatible with the loss of dendritic density and reduced dendritic arborisation within the motor cortex, as suggested by recent histopathological analyses in humans and animal models.30 31 Although changes in cortical dendrites are observed at early stages of disease evolution in animal models of ALS,32 33 the cortical findings in our cohort using NODDI are limited to a region within the right PCG. It is possible that pathological heterogeneity within the primary motor cortex contributes to difficulties in detecting consistent alterations in cortical areas using current imaging techniques for whole brain comparisons or that loss of cortical dendrites is a relatively late phenomenon in the evolution of cerebral pathology in ALS. Nonetheless, these findings indicate that NODDI can reveal new aspects of the cellular pathology of diseases such as ALS, in keeping with the recent demonstration that changes in NODDI parameters closely reflect complex histological changes in the grey matter (GM) and in the WM of the spinal cord in multiple sclerosis.13
A direct comparison of the NDI and FA whole brain findings in our study shows that NDI detected larger areas of neurodegeneration than FA at the same statistical threshold, where the combined cluster sizes for the areas of significantly reduced NDI (29 349) were greater than those for areas of significantly reduced FA (9283) (see table 2). This suggests that NDI may have a higher sensitivity than FA for detecting ALS-related WM pathology and also indicates that the primary pathology identified by the diffusion MRI signal is alteration in axon density.
In our study FA was reduced into two separate regions of the CST in the ALS group, one within the pons and another within the corona radiata extending into the subcortical WM of the PCG. The CST has previously been implicated in DTI studies demonstrating reduced FA in the ALS group.34–36 Typically studies report reduced FA in one or several isolated regions, including the pons, cerebral peduncles, posterior limb of the internal capsule, corona radiata and the subcortical WM of the PCG.3 4 35 36 Low FA within the corpus callosum is also typically demonstrated in ALS37 but not always.4 Rarely are such extensive, contiguous changes demonstrated throughout the CST and corpus callosum, as evident with our findings of reduced NDI using NODDI.
We found increased MD in the ALS group within the left internal capsule, which is consistent with previous studies.3 Reports of increased MD have been less consistent than FA findings in ALS, which may be explained by the known dependency of MD on the specific b-value used for the acquisition.26 In this study, we used a maximum b-value of 2400 smm−2, which is higher than typical DTI protocols and may account for our findings. We could have restricted our DTI analysis to the inner shell (b=800) data to make it more compatible with other DTI studies. However this would have made the comparison with NODDI (fitted with the full data set) unfair. In general the extent of WM changes we found using DTI is somewhat reduced compared with other studies.3 4 34–37 This might be the consequence of specific methodological differences, the relatively small sample size or the clinical heterogeneity of the cohort.
NODDI facilitates the interpretation of the altered diffusivity profile in ALS observed using DTI where axonal loss, demonstrated by reduced NDI in the CST and corpus callosum, is likely to be the main contributing factor for the reduced FA and increased MD. We detected reduced ODI in the right internal capsule, a region where there is crossing of WM pathways. This could be explained by the predominant degeneration of axonal fibres in one direction, thereby reducing crossing of fibres in this region. Furthermore, this may contribute to the lack of changes in FA detected in this region as the corresponding FA may not change due to the remaining fibres becoming more coherent, despite axonal degeneration.
Increased ISO was detected within the right lateral ventricle, likely representing atrophy in ALS, leading to increased CSF spaces.
When subgroup analysis was performed, the ALS subgroup with both bulbar and limb involvement showed loss of right CST axon density significantly greater than the group with limb-confined ALS. This is in keeping with our hypothesis that more widespread clinical involvement in ALS relates to more extensive microstructural damage.
Regression analysis revealed that both ODI and ISO correlated with disease duration, where regions of the right PCG and interhemispheric fissure were involved, respectively. This suggests that loss of cortical dendrites and related cerebral atrophy are occurring alongside axonal degeneration. The correlation of ODI in the PCG with disease duration provides support for the rationale that loss of complex neuropil architecture within the motor cortex is occurring in ALS over time.
Clinical heterogeneity of the ALS participants in this study was considerable, which could explain why we failed to detect a significant group difference in more regions of the PCG. A larger study using a more homogeneous ALS phenotype may be required in order to achieve this on whole brain analysis of ODI. Although statistically significant cortical GM changes in ODI and NDI were limited to a region within the right PCG, these abnormalities could be more relevant to the site and severity of disease. Moreover, the topographical pattern of disease within the PCG in ALS is likely to have a larger spread of regions involved, with higher variability between patients, when compared with the converging axons of the CST, the common output pathway of all cortical motor fibres. This reduces the detection of significant differences on whole brain analysis, which may therefore only become more apparent in larger cohorts. Although the GM results are not as strong as those obtained in the WM, they nevertheless support the hypothesis of dendritic abnormalities paralleling disease progression.
We recognise that there are limitations related to the technique used in this study. NODDI may not separate the effects of orientation dispersion and neurite density effectively within the GM and may therefore not fully distinguish intracortical neuropil pathology in ALS.38 More recently developed quantitative diffusion MRI techniques, incorporating microscopic anisotropy mapping, may separate different signal components better and highlight differences more effectively39 40 but would require additional imaging sequences that are not yet available clinically. The NODDI technique used in this study11 models orientation dispersion isotropically and may therefore be limited in its capacity to model multiple fibre orientations arising from complex dendritic structures.7 NODDI has subsequently been developed further to incorporate a quantification of anisotropic orientation dispersion, and analysis using this more sophisticated modelling system is warranted in future studies in ALS.12
We used whole brain voxel-wise analysis, which does not require any a priori hypotheses and provides information about changes in both WM and GM. Although whole brain analysis is prone to potential biases due to interindividual anatomical variation and partial volume effects, NODDI has the advantage over DTI of explicitly modelling the free water (isotropic) component and should therefore be less susceptible to the partial volume contamination of CSF. Alternative processing approaches to whole brain analysis, such as tract-based spatial statistics41 and GM special statistics,14 have been optimised to minimise the potential sources of error and to detect changes in WM and GM, respectively. Tractometry, which maps quantitative indices (such as those derived from NODDI) to points of streamlines reconstructed by diffusion tractography, is a powerful visualisation tool and can potentially give insight into the regions of degeneration but suffers from well-known limitations.42
It should also be noted that many of the studies using NODDI14 15 were performed at 3T, while our data were performed at 1.5T. The field strength has obvious implication for the signal to noise ratio (SNR). However, the reduced signal intensity at 1.5T is partially compensated by the longer T2 relaxation time at this field strength. We used a voxel size of 2.5 mm3, which is almost twice as large as that used in the original NODDI paper11 at 3T, and should ensure a comparable SNR at the price of reduced spatial resolution.
Another limitation is the sample size, which is particularly relevant given emerging insights into disease heterogeneity in ALS.1 Nevertheless, our findings illustrate the potential for NODDI to interrogate hitherto undetectable details of cellular pathology of ALS in life.
We conclude that NODDI reveals evidence compatible with structural damage in the PCG in ALS, provides new insights into the cellular pathology of WM damage, and demonstrates that reduced density of corticospinal axons is the main contributing factor underlying the altered diffusivity profile detected by DTI. Our observations are compatible with relatively early involvement of cortical neuropil including dendrites in ALS. More extensive longitudinal studies in larger phenotypically and genetically defined cohorts will be needed to refine further studies on the cortical microstructure in ALS using NODDI and developments thereof.
The authors would like to thank all the participants who contributed to this study.
Contributors As clinical research fellow, RJB recruited and consented participants, performed all clinical assessments, supervised imaging, and processed and interpreted the data. RJB also constructed the first draft of this article and prepared all of the figures. MCG performed the data preprocessing and advised on statistical analysis. NGD designed the NODDI pipeline for the institution and provided supervision to RJB and MCG in data preprocessing. DJS, AKS, HZ and DCA contributed to the design of the study and contributed to the drafting of the paper. HZ and DCA provided the necessary skills underpinning NODDI, advised on the development of the hypothesis and contributed to writing the funding proposal. They also advised on data interpretation and on all drafts of the paper. MC contributed to the design of the study, and to developing and submitting the funding proposal, and designed and implemented the imaging protocol. MC also provided supervision for the data processing and analysis, and revised the drafts of the paper at all stages. PNL was chief investigator for this study. PNL screened the subjects with ALS recruited through the Sussex ALS clinical service. PNL and MC conceived the hypotheses and submitted the funding proposal, supervised the clinical research fellow (RJB) and oversaw all aspects of the study including ethics and research governance. PNL edited and revised the drafts of the paper at all stages. All authors reviewed and accepted the final manuscript.
Funding This study was supported by a grant from the Motor Neurone Disease Association of England and Northern Ireland (Reference: Leigh/Apr14/824/791). EPSRC grants N018702 and M020533 support DCA's work on this topic.
Competing interests None declared.
Patient consent Obtained.
Ethics approval Favourable ethical opinion was awarded by the National Health Service Research Ethics Committee London (14/LO/0195).
Provenance and peer review Not commissioned; externally peer reviewed.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.