Progression of Microstructural Damage in Spinocerebellar Ataxia Type 2: A Longitudinal DTI Study

BACKGROUND AND PURPOSE: The ability of DTI to track the progression of microstructural damage in patients with inherited ataxias has not been explored so far. We performed a longitudinal DTI study in patients with spinocerebellar ataxia type 2. MATERIALS AND METHODS: Ten patients with spinocerebellar ataxia type 2 and 16 healthy age-matched controls were examined twice with DTI (mean time between scans, 3.6 years [patients] and 3.3 years [controls]) on the same 1.5T MR scanner. Using tract-based spatial statistics, we analyzed changes in DTI-derived indices: mean diffusivity, axial diffusivity, radial diffusivity, fractional anisotropy, and mode of anisotropy. RESULTS: At baseline, the patients with spinocerebellar ataxia type 2, as compared with controls, showed numerous WM tracts with significantly increased mean diffusivity, axial diffusivity, and radial diffusivity and decreased fractional anisotropy and mode of anisotropy in the brain stem, cerebellar peduncles, cerebellum, cerebral hemisphere WM, corpus callosum, and thalami. Longitudinal analysis revealed changes in axial diffusivity and mode of anisotropy in patients with spinocerebellar ataxia type 2 that were significantly different than those in the controls. In patients with spinocerebellar ataxia type 2, axial diffusivity was increased in WM tracts of the right cerebral hemisphere and the corpus callosum, and the mode of anisotropy was extensively decreased in hemispheric cerebral WM, corpus callosum, internal capsules, cerebral peduncles, pons and left cerebellar peduncles, and WM of the left paramedian vermis. There was no correlation between the progression of changes in DTI-derived indices and clinical deterioration. CONCLUSIONS: DTI can reveal the progression of microstructural damage of WM fibers in the brains of patients with spinocerebellar ataxia type 2, and mode of anisotropy seems particularly sensitive to such changes. These results support the potential of DTI-derived indices as biomarkers of disease progression.

S pinocerebellar ataxia type 2 (SCA2) is the second most frequent autosomal dominant inherited ataxia worldwide, after SCA3. 1 It is caused by expansion in excess of 32 CAG repeats in the gene encoding the Ataxin-2 protein, which mainly targets several pontine neurons and Purkinje cells in the cerebellum, 2 and it is associated with a pathologic pattern of pontocerebellar atrophy. 1,3 MR T1-weighted imaging enables in vivo detection of brain stem and cerebellar atrophy in patients with SCA2 in crosssectional 4-6 and longitudinal 7 studies.
We performed a longitudinal DTI study in 10 patients with SCA2 and 16 age-matched healthy controls to explore the ability of DTI to detect and map the progression of microstructural damage reflecting advance of neurodegeneration. In particular, we analyzed several DTI-derived indices, including mean diffusivity (MD), axial diffusivity (AD), radial diffusivity (RD), fractional anisotropy (FA), and mode of anisotropy (MO), by using tractbased spatial statistics (TBSS), which enable a robust and unbiased voxelwise whole-brain analysis of the main white matter tracts. 19,21,31,32

Subjects
Ten patients (4 women, 6 men; mean age, 47.5 Ϯ 12.7 years) with a genetically determined diagnosis of SCA2 33 gave their consent to participate in this longitudinal study, which was approved by the local ethics committee. The cutoff number of triplet repeat expansions that qualified for diagnosis of SCA2 was 34 CAG repeats on 1 allele, and the mean number of abnormal triplets was 40.6 Ϯ 1.4. Each patient underwent MR imaging twice on the same scanner and with the same acquisition protocol at an average of 3.6 Ϯ 0.7 years apart (range, 2.2-4.1 years). At the time of both MR imaging examinations, each patient underwent evaluation by the same clinician (A.G.), who assessed the neurologic deficit by using the Inherited Ataxia Clinical Rating Scale. 34 This is a 0 -38 scale (38 corresponds to maximum clinical deficit) that evaluates signs and symptoms related to ataxia and also pyramidal tract dysfunction and impaired vibration or position sense, which are frequently observed in SCA2. 1 At the time of baseline MR imaging, the mean disease duration from clinical onset in patients with SCA2 was 12.8 Ϯ 7.3 years (range, 2-23 years), and the mean Inherited Ataxia Clinical Rating Scale score was 17.2 Ϯ 4.3 (range, 9 -25). Demographic, genetic, and clinical data of the patients with SCA2 are summarized in Table 1.
Sixteen age-and sex-matched healthy subjects (7 women, 9 men; mean age, 50.3 Ϯ 18.8 years) were recruited as controls and provided written consent to participate in the study. They had no history of neurologic or psychiatric dysfunction, and results of their neurologic examinations at baseline and follow-up were negative. Control subjects also underwent MR imaging twice with the same scanner and acquisition protocol at an of average of 3.3 Ϯ 1.0 years apart (range, 1.9 -4.7 years).
The TBSS results from the baseline examinations of 9 of the patients with SCA2 and 10 of the healthy controls enrolled in the present study were reported previously. 19

Image Processing and TBSS
Preliminarily visual evaluations of all B0 and diffusion-weighted images revealed no motion artifacts in any patient with SCA2 or control subject. Diffusion-weighted images were corrected for head motion and eddy current distortions by using the FMRIB Diffusion Toolbox (2.0; http://www.fmrib.ox.ac.uk/fsl/fdt/index. html), part of the FMRIB Software Library (FSL 5.0.2; http:// www.fmrib.ox.ac.uk/fsl), 35 after which brain tissue was segmented by using the FSL Brain Extraction Tool (http://fsl.fmrib. ox.ac.uk/fsl/fslwiki/BET). The b-matrix was subsequently reoriented by applying the rotational part of the affine transformation used in the eddy/motion-correction step. 36 A tensor model was then fitted to the raw data by using a constrained nonlinear leastsquares procedure implemented in the software package Camino (www.camino.org.uk), 37 and residual nonpositive definite tensors (in regions in which the nonlinear algorithm failed to converge) were removed by tensor interpolation in the log-Euclidean domain. 38 Diffusion tensor invariants, MD, AD, RD, FA, and MO, were then computed from the estimated tensor field by using the DTI ToolKit software package (http://software.incf.org/software/dtitoolkit-dti-tk). 39 The MO is an adimensional tensor invariant (range, Ϫ1 to 1) that characterizes anisotropy type. 40 In particular, the MO specifies the type of anisotropy as a continuous measure reflecting differences in the shape of the diffusion tensor ranging from planar (eg, in regions of crossing fibers from 2 fiber populations of similar attenuation or regions of "kissing" fibers [low MO value]) to linear (eg, in regions in which 1 fiber population orientation predominates [high MO value]). 41 Voxelwise statistical analyses of all tensor invariant data were carried out by using TBSS, 32 also part of the FMRIB Software Library, which involves 1) nonlinear registration of every FA image to every other one, 2) identification of the "most representative image" as the one that requires the least warping to align every other image to it, 3) affine alignment of this target image into Montreal Neurological Institute 152 (MNI152) space, 4) transformation of every other image into MNI space by combining the single nonlinear transforms into the target image with the affine transform into MNI space, 5) creation of a mean FA image, and 6) thinning to create a mean FA skeleton that represents the centers of all tracts common to the group. Each subject's aligned FA data were then projected onto this skeleton, and previously computed warps were also successively applied to MD, AD, RD, and MO images, thereby transforming all images into MNI space for subsequent skeletonization as described above. The resulting data were then fed into voxelwise intersubject and intrasubject statistics. Each analysis included full correction for multiple comparisons over space by using permutation-based nonparametric inference within the framework of the general linear model. Specifically, a first general linear model included only group ("between" factor) to investigate the effect of group (patients versus controls) at baseline. A second general linear model included both group and time ("within" factor), which enabled us to investigate the grouptime interaction (patient versus control differences in follow-up versus baseline changes). In another general linear model, we investigated within-patient correlation between changes in DTI-derived indices and changes of clinical variables. All the general linear model designs included age and sex as nuisance covariates. Also, with the aim of taking into account potential confounders resulting from longitudinal atrophy, 7 every subject's follow-up T1 scan was registered to the corresponding baseline T1 scan by using high-dimensional nonlinear registration, after which a Jacobian determinant image (which quantifies regional volume change) was computed for the inverse transformations. The baseline scans were all registered to a population-specific T1 template, which in turn was registered to the MNI152 T1 template. The Jacobian determinant image was then transformed into MNI152 space through a single interpolation step in which all the aforementioned transformations were combined. Full details about this procedure can be found in Mascalchi et al 7 and Tessa et al. 42 After masking out CSF, the whole-brain average of the Jacobian determinant image was inserted into the longitudinal general linear models as a subject-specific nuisance covariate of global atrophy.
To improve the confidence intervals on resulting P-value maps, in each comparison we used 50,000 permutations (as opposed to the commonly used value of 5000). P values were calculated and corrected for multiple comparisons by using the 2D parameter settings with threshold-free cluster enhancement, thereby avoiding the use of an arbitrary threshold for the initial cluster formation. 43

Baseline Evaluations
TBSS showed changes in all DTI-derived indices in several WM tracts of the brain in patients with SCA2 that were significant compared with those of the controls. These tracts included the cerebral and cerebellar WM, the corpus callosum, brain stem, cerebellar peduncle, and thalamus (On-line Fig 1) and consisted of increased MD and RD and decreased FA and MO. AD was increased in the brain stem, cerebellum, and cerebellar peduncles, and the WM was underlying the right primary sensorimotor cortex, corpus callosum, and thalami. However, AD was decreased in the right internal capsule and frontal WM. There were no significant decreases in MD or RD or significant increases in FA or MO.

Longitudinal Evaluations
TBSS revealed longitudinal AD and MO changes in patients with SCA2 that were significant compared with those of the controls (Fig 1).
Longitudinal AD changes were significantly greater (ie, increased) in patients with SCA2 than in controls in the right cerebral hemisphere and corpus callosum (Fig 1). The mean longitudinal changes in these regions are listed in Table 2.
Longitudinal MO changes were significantly lower (ie, decreased) in patients with SCA2 than in controls in many WM tracts, including bilateral hemispheric cerebral WM, corpus callosum, internal capsule, cerebral peduncle, pons, and left inferior, middle, and superior cerebellar peduncle tracts and in the WM of the left paramedian vermis (Fig 1). The mean longitudinal changes in these regions are listed in Table 2.
No WM tract showed significant differences in longitudinal MD, RD, or FA changes in patients with SCA2 versus those in controls.

MR Imaging Clinical Correlations
No significant correlations between clinical progression (modification of the Inherited Ataxia Clinical Rating Scale) and changes in the DTI-derived indices were observed.

DISCUSSION
In patients with SCA2, histopathologic studies have shown widespread neuronal loss, which is prominent in the cerebellar cortex and pontine nuclei, midbrain, medulla, and motor cortex. 1,3,44 WM damage in patients with SCA2 consists of a loss of myelinated fibers and gliosis affecting the transverse pontine fibers, central cerebellum and cerebellar folia, middle and inferior cerebellar peduncles, medial lemniscus and trigeminal tracts, fasciculus gracilis and cuneatus, and spinocerebellar tracts. 1 Gliosis has been observed also in the thalamus, globus pallidus, and subthalamic region. 1 MR imaging shows loss of bulk of the brain stem and cerebellum in presymptomatic and early symptomatic SCA2 gene carriers 4-6 and circumscribed atrophy of the frontal, parietal, and temporal cortex, the thalamus, and frontal and temporal WM in more advanced cases. 19,45,46 MR imaging can also show diffuse high signal intensity of the brain stem and cerebellar WM with sparing of the corticospinal tracts in proton-attenuation and T2-weighted images. 47 Two studies evaluated patients with SCA2 with DTI using regions of interest 16 or MD histogram analysis and TBSS. 19 Both studies revealed a widespread decrease in FA and an increase in the apparent diffusion coefficient, MD, and RD, which are consistent with microstructural damage of WM tracts in the brain stem, cerebellar peduncles, and cerebellum and in the corpus callosum and several WM tracts in the cerebral hemispheres, including the corticospinal tracts, inferior longitudinal fasciculus, and inferior fronto-occipital fasciculus.
To the best of our knowledge, this is the first DTI study to have explored longitudinal changes of the brain regional microstructure in patients with SCA2 compared with those of healthy age-matched controls. In particular, patients with SCA2 showed progression of WM damage in the hemispheric cerebral WM, corpus callosum, pons, and left cerebellar peduncles and in the WM of the left paramedian vermis. We noted a hemispheric asymmetry of some longitudinal DTI changes in our patients with SCA2, which may reflect either the relatively low sample size or a truly uneven distribution of the neurodegenerative process. Asymmetry of damage was previously reported from MR imaging studies of patients with Huntington disease, familial Alzheimer disease, and Friedreich ataxia 21,[48][49][50][51] and is not uncommon in pathologic studies of neurodegenerative diseases. 52 In our opinion, the issue of the symmetric versus asymmetric distribution of microstructural changes revealed by DTI needs further investigation with larger sample sizes.
We believe that a thorough analysis of the cross-sectional and longitudinal changes of the various DTI-derived indices that were used in this study, and their interplay in a typical neurodegenerative disease such as SCA2, can be of interest. At baseline, TBSS showed widespread changes in DTI-derived indices in patients with SCA2 compared with those in the controls. These changes include increases in MD and RD and decreases in FA and MO. AD was increased in the brain stem, cerebellum and cerebellar peduncles, thalami, and WM underlying the (right) primary sensory cortex and decreased in the right internal capsule and frontal WM. The increases in MD and RD and decreases in FA are in line with those reported so far in cross-sectional studies of various neurodegenerative diseases. 5,[8][9][10][11][13][14][15][16][17][18][19][20][21][22][23][24][25][27][28][29]53 In our study, we observed areas of increased AD that are in line with those shown in cross-sectional studies of patients with Alzheimer disease, [53][54][55][56] Friedreich ataxia, 20 and Huntington disease. 13,55 It is notable that an increase in AD per se is not inconsistent with a decrease in FA. In fact, it is sufficient for the increased AD to be accompanied by a more pronounced increase in RD to result in a decrease of FA in degenerating WM fibers. 20 Also, we observed an area of decreased AD in the right internal capsule and frontal WM. The coexistence of WM areas exhibiting increased and decreased AD was reported from a study of Alzheimer disease 53 and supports the view that AD may show time-dependent behavior related to a variable combination of factors determining either a decrease caused by Wallerian degeneration and gliosis or an increase caused by axonal damage and atrophy of WM fiber tracts with expanded extracellular space. In our opinion, the complex interplay of these processes should be addressed in animal-model studies. 57 To our knowledge, so far, only 1 cross-sectional study in patients with mild cognitive impairment and mild Alzheimer disease addressed MO changes in patients with a neurodegenerative disease. 41 An initially counterintuitive increase in the MO was observed in a region of crossing fibers (ie, planar anisotropy) in the centrum semiovale of patients with mild cognitive impairment versus controls and was explained by selective preservation of motor-related projection fibers (at an early stage of the disease) crossing the degenerating association fibers of the superior longitudinal fasciculus, which results in more linear anisotropy. 41 In our study, we observed at baseline a widespread decrease of the MO in patients with SCA2 versus that of controls, which included areas with crossing fibers such as cerebellar  hemispheric WM, centrum semiovale, lobar WM, and areas without crossing fibers such as the cerebellar peduncles, brain stem longitudinal tracts, external and internal capsules, and corpus callosum. Given that a decrease in the MO indicates a reduction of the degree of linear anisotropy, our results suggest a general disorganization of WM fibers in patients with SCA2. The longitudinal data in our study reveal significantly more pronounced modifications of AD and MO in patients with SCA2 than in the controls. In the patients with SCA2, the increase in AD was circumscribed to the right frontal WM and corpus callosum, whereas the decrease in the MO involved symmetrically and extensively the hemispheric cerebral WM and internal capsules, the corpus callous, the cerebral peduncles, the pons, all 3 of the cerebellar peduncles, and the left paramedian vermis WM. No areas of significant difference in changes in MD, RD, or FA in patients with SCA2 were detected compared with those in controls.
The increase in AD at follow-up in the patients with SCA2 is in line with the data from a longitudinal study of symptomatic patients with Huntington disease and was interpreted as the reflection of increased extra-axonal space, resulting from reduced axonal caliber, which allows faster water-molecule movement parallel to axons. 55 The decrease in the MO at follow-up in patients with SCA2 is in line with the progressive WM fiber disorganization expected in a neurodegenerative disease. It is notable that the far more extensive longitudinal changes in the MO compared with the modifications in AD, including WM tracts in the brain stem and cerebellum, may support the potential use of the MO as an additional valuable and sensitive descriptor of microstructural changes. 41 We submit the following non-mutually exclusive explanations for the lack of significant differences in longitudinal MD, FA, and RD modification in our group of patients with SCA2 compared with those in the control group: 1) the occurrence of a "ceiling" effect; 2) low statistical power (␤ error) (in this regard, a gain in the signal-to-noise ratio, which can be obtained at higher field strength [eg, 3T] or by using a head coil with a greater number of channels [eg, 32] for advanced neuroimaging) and the use of high angular resolution diffusion-weighted imaging acquisitions would have improved the sensitivity of our study in detecting any white matter changes in the brain stem and/or cerebellum.
Finally, we preliminarily explored whether longitudinal DTI changes were correlated with clinical progression in patients with SCA2. The lack of a significant correlation between longitudinal DTI changes and clinical deterioration in our study (which has been observed also for morphometry findings in patients with SCA2 7 and those with SCA1, SCA3, or SCA6 58 ) might reflect the small sample size, but in our opinion, supports the role of DTI metrics as potential biomarkers of the progression of disease-related microstructural alterations that may precede the evolution of clinical disease manifestations.
However, we did not obtain a structured assessment of cognitive or behavioral changes in our patients. We therefore could not explore the possibility that DTI modifications are correlated with the evolution of nonmotor features. Given that a recent study indicated that a simple clinical examination can track disease progression in paucisymptomatic SCA2 gene carriers, 59 we feel that a tight comparative evaluation of clinical and MR techniques including DTI and morphometry in monitoring neurodegenera-tion is advisable for ascertaining the optimal strategy in view of future trials, especially in patients in the early stages of disease. 60 CONCLUSIONS DTI is capable of tracking the progression of microstructural brain damage and therefore represents a candidate biomarker of the progression of neurodegeneration in patients with SCA2.