Objectives: To test the applicability of an automated method of magnetic resonance image analysis (voxel based morphometry) to detect presence and severity of regional grey matter density reduction—a proxy of atrophy—in Alzheimer's disease.
Methods: Twenty nine probable Alzheimer's patients and 26 non-demented controls (mini-mental state examinations mean (SD) 21 (4) and 29 (1)) underwent high resolution 3D brain magnetic resonance imaging. Spatial normalisation to a stereotactic template, segmentation into grey matter, white matter, and cerebrospinal fluid, and smoothing of the grey matter were carried out based on statistical parametric mapping (SPM99) algorithms. Analyses were carried out: (a) contrasting all Alzheimer's patients with all controls (p<0.05 corrected for multiple comparisons); (b) contrasting the three Alzheimer's patients with mini-mental state of 26 and higher with all controls (p<0.0001 uncorrected); and (c) correlating grey matter density with mini-mental state score within the Alzheimer's group (p<0.0001 uncorrected).
Results: When all Alzheimer's patients were compared with controls, the largest atrophic regions corresponded to the right and left hippocampal/amygdalar complex. All parts of the hippocampus (head, body, and tail) were affected. More localised atrophic regions were in the temporal and cingulate gyri, precuneus, insular cortex, caudate nucleus, and frontal cortex. When the mildest Alzheimer's patients were contrasted with controls, the hippocampal/amygdalar complex were again found significantly atrophic bilaterally. The mini-mental state score correlated with grey matter density reduction in the temporal and posterior cingulate gyri, and precuneus, mainly to the right.
Conclusions: Voxel based morphometry with statistical parametric mapping is sensitive to regional grey matter density reduction in mild Alzheimer's disease.
- Alzheimer's disease
- voxel based morphometry
- hippocampus, atrophy
- AD, Alzheimer's disease
- MMSE, mini-mental state examination
- SPM, statistical parametric mapping
- MR, magnetic resonance
Statistics from Altmetric.com
- AD, Alzheimer's disease
- MMSE, mini-mental state examination
- SPM, statistical parametric mapping
- MR, magnetic resonance
Of the biological indicators of Alzheimer's disease (AD), those detecting medial temporal lobe atrophy based on three dimensional T1 weighted magnetic resonance (MR) imaging are among the most sensitive. Hippocampal and entorhinal cortex volumetry and surface measures of the entorhinal cortex are believed to have high sensitivity and specificity in the detection of AD from non-demented elderly controls1,2 and can predict conversion of mild cognitive impairment to AD.3 However, the variability of the measurements because of the human tracer has so far limited direct comparisons of the results of different research groups. The development of observer independent tools might be a significant advancement.
To date, some computer based tools are available that can detect volume or shape changes in the regions of interest to AD.4–7 Among these, one of the most user independent is statistical parametric mapping (SPM).8 Initially developed for brain activation studies with positron emission tomography,9 it has been applied to different image analysis issues. In the case of morphometry, SPM permits the comparison of the density of grey matter on a voxel by voxel basis in a group of patients with that of a group of controls. The result of the comparison is a three dimensional map of significant regional differences of grey matter density reduction, a close proxy of grey matter atrophy. The process is largely automatic and human intervention is comparatively low. Three studies have been published so far on AD,10–12 but inconsistent results, different study populations and analysis protocols, and the present availability of more sophisticated image processing protocols13 argue for an independent confirmation.
The aim of this study is to test the sensitivity of voxel based morphometry with SPM in the detection of grey matter density reduction and its association with global disease severity in patients with mild to moderate AD. Moreover, the sensitivity to the very mild stages of AD will be assessed in those AD patients with highest mini-mental state score (26 and above).14
This study comprised 29 patients who met NINCDS-ADRDA criteria for probable AD15 and 26 non-demented controls recruited between September 1993 and December 1994 for a study on linear measures of brain atrophy in AD and have been described in previous reports.16 AD were outpatients seen at the Alzheimer's Unit, Brescia, Italy. Routine dementia assessment including standardised history, laboratory examinations, physical and neurological examination, neuropsychological assessment, and routine computed tomography was carried out in all.16 Control subjects were patients' relatives (mostly spouses) without detectable cognitive deficit. All were given the mini-mental state examination (MMSE) and clinical dementia rating scale17 and were judged not to be demented by a neurologist and a psychologist involved in the evaluation of the patients. Uncontrolled physical diseases (heart or renal failure, cirrhosis, etc) were considered exclusion criteria, while risk factors such as hypertension and diabetes were accepted. Of the original 46 patients and 31 controls,16 17 and 4 could not be included in this study as some slices were missing. One additional control was excluded for the development of multiple system atrophy a few years after scanning.
Table 1 shows that patients were about five years older (p=0.04) than controls, while sex and education were not different (p>0.09). Dementia severity was of mild to moderate degree. All patients and controls were right handed.
Three AD patients were in the highest decile of the MMSE distribution and were selected to investigate atrophy in very mild AD. Age, sex, and MMSE were: 66, F, 26; 81, F, 27; and 79, M, 27.
Written informed consent was obtained by both patients and control subjects or their primary caregivers after discussion of risks and benefits of participation. No compensation was provided. The study was approved by the local ethics committee.
MR images acquisition
For each subject, a high resolution sagittal T1 weighted volumetric MR scan was acquired at the Radiology Department, University of Verona, using a 1.5 tesla Magnetom scanner (Siemens, Erlangen, Germany), with a gradient echo 3D technique: TR=10 ms, TE=4 ms, TI=300 ms, flip angle=10°, field of view=250 mm, acquisition matrix 160×256, slice thickness 1.3 mm, 128 slices. The same parameters were used for every subject.
After removing the voxels below the cerebellum with MRIcro (www.psychology.nottingham.ac.uk/staff/cr1/mricro.html),18 MR scans were analysed with SPM99 (www.fil.ion.ucl.ac.uk/spm) running under Matlab 6.0 (Mathworks, Sherborn, MA, USA) on a Sun Sparc Ultra 30 workstation (Sun Microsystem, Mountain View, CA). Before carrying out group comparisons, three steps need to be carried out with SPM99: creation of a grey matter template, normalisation of the grey matter to the template, and smoothing. A complete description of the preprocessing procedures can be found elsewhere.6
Grey matter template
This was made through creation of a customised template of the whole brain, normalisation of the original images to the template, extraction of the grey matter, smoothing with a 12 mm filter, and averaging. The customised template of the whole brain was obtained through normalisation of the images of all patients and controls to the T1 stereotactic template of SPM9919 using a 12 parameter affine transformation followed by smoothing with an 8 mm isotropic Gaussian kernel and averaging of the smoothed normalised images. Starting estimates for the registration were assigned by specifying the position of the anterior commissure.20
The normalisation of the original images to the template was obtained through a 12 parameter affine and non-linear transformation. Grey matter was then extracted through segmentation of images into grey matter, white matter, and cerebrospinal fluid with a modified mixture model cluster algorithm, smoothed with a 12 mm filter, and averaged thus obtaining a stereotactic customised grey matter template. Customised priors were computed from this template that will be used to improve the accuracy of the extraction of the grey matter of each participant in the following normalisation step.13
The grey matter of patients was extracted with the same procedure as described above for controls. The grey matter of patients and controls was normalised onto the grey matter template with affine and non-linear transformations, medium regularisation, reslicing 2×2×2 mm, and no masking.10 The resulting images were visually inspected one by one to exclude macroscopical segmentation errors.
The normalised grey matter images were smoothed with an isotropic Gaussian filter with a full width at half maximum of 8 mm.
Two statistical models were built to compare grey matter density between all patients and controls and between the three very mild patients and controls with the “Compare populations–AnCova” procedure. A separate model was built to detect the regions where grey matter density was associated with dementia severity with the “Single subject–Covariates only” procedure. Proportional scaling was not used and grey matter thresholding was set at 0.8.
In the first two models the contrast of interest was that of decreased grey matter density of patients compared with controls and age was entered as a nuisance covariate. In the first model (all patients versus controls) the significance threshold of the t statistics was set at 0.05 corrected for multiple comparisons. In the second model (the three very mild patients versus controls) the significance threshold was set at a higher level (0.0001 uncorrected) to compensate for the small sample size. The anatomical localisation of the t statistics peaks was made based on visual examination with the aid of a brain atlas.21 In the third model, MMSE was entered as the covariate of interest. The significance threshold of the t statistics was set at 0.0001 uncorrected.
Table 2 and figure 1 show the results of the comparison between all patients and controls. The table shows that the largest regions of significant grey matter density reduction were located bilaterally in the medial temporal lobe structures (hippocampus, amygdala, and uncus). Grey matter density reduction was significant in all parts of the hippocampus (head, body, and tail). Figure 1 shows that the largest cluster located in the head of the right hippocampus and amygdala tended to involve the two structures separately. Other smaller areas of grey matter density reduction were located in the temporal lobe (superior and middle temporal gyri and fusiform gyrus), parietal lobe (precuneus), mesial aspect of the brain (anterior and posterior cingulate), and insular cortex. Scattered but significant voxels were detected in the head of the caudate nucleus and frontal gyri. Overall, the structures on the right side were more heavily involved than those on the left.
The results were checked in a subgroup of 17 patients and 17 controls individually matched (±2 years) for age (mean (SD) ages 72.5 (8.7) and 73.0 (8.4) respectively). The largest and most significant clusters were confirmed in the medial temporal lobe (left hippocampal head/posterior amygdala at −26, −8, −14 with z=5.07 and right anterior amygdala at 22, 4, −20 with z=4.96).
Table 3 and figure 2 show the results of the comparison between the three patients with the highest MMSE and controls. Medial temporal regions on both sides were again the most significant or largest regions. Smaller regions of grey matter density reduction were detected in the temporal and frontal gyri. One voxel resulted significant in the precentral gyrus. To detect false positive voxels that might arise from the comparatively low significance threshold, we checked which voxels were negative in the first experiment, where a much stricter p threshold was used (p=0.05 corrected). We found that the only surviving voxels were located in the hippocampal/amygdalar region.
The results were checked with an age matched group of 14 controls with age range (66 to 81 years) equal to that of the three very mild AD patients (mean (SD) ages 74.2 (5.5) and 75.3 (6.7) respectively). The largest and most significant cluster was detected in the left amygdala (at −30, 4, −24 with z=4.89).
The areas where grey matter density reduction was correlated with dementia severity (table 4 and fig 3) corresponded to the temporoparietal cortex mainly to the right involving the superior and inferior temporal gyri, parietal lobule, and precuneus. A comparatively large region was also detected in the left posterior cingulate. The grey matter density values of the most significant voxels in the right and left regions showed a significant correlation with MMSE scores (fig 4).
All the analyses were re-run by covarying by sex, but the results did not change appreciably.
We have shown that voxel based morphometry with SPM: (1) is sensitive to medial temporal grey matter density reduction in AD and that this is the case even in the very mild patients; (2) demonstrates grey matter density reduction outside the medial temporal lobe confined to areas in the temporal gyri, precuneus, insular and cingulate cortex, and caudate nucleus; (3) is sensitive to the temporoparietal atrophy associated with the progression of the disease. These results support the use of this tool to study regional atrophy in AD.
Three studies have performed voxel based morphometry in AD. In seven comparatively young AD patients (mean age of 65 years) with mild to moderate AD compared with seven controls, Rombouts et al11 found atrophy of the hippocampus, insula, and, unexpectedly, caudate nucleus. In the hippocampus and insula, atrophy tended to be greater on the right. The authors used a moderately restrictive criterion for atrophy detection (p<0.0005 uncorrected for multiple comparisons). Baron et al10 have studied 19 mild to moderate AD patients and compared them with 16 controls finding atrophy in the medial temporal structures and cingulate, insular, and temporoparietal cortex. Medial temporal atrophy was centred in the amygdalae on both sides although more marked on the right, where also the posterior hippocampus was involved. The criterion for atrophy detection was p<0.10 corrected for multiple comparisons. Ohnishi et al12 compared 26 mild with moderate AD patients to 92 controls and found atrophy remarkably confined to the whole length of the hippocampus bilaterally but more marked on the right side. The SPM parameters used to obtain these results are not described in sufficient detail to permit comparisons with other works. All the studies have used as template brain the one computed at the Montreal Neurological Institute resulting from averaging non-atrophic brains of 152 young adults (mean age 25 years).19
In this study, we have assessed a patient group representative of the population of sporadic AD (mean age 74 years) and used a preprocessing protocol (customised template and customised priors computed from our study groups) that should optimise the normalisation and segmentation of atrophic brains. We have confirmed all of the previous findings (atrophy of head body and tail of the hippocampus, amygdala, insula, lateral temporal and cingulate cortex, precuneus, and caudate nucleus) and have found evidence of involvement of other regions demonstrated to be affected in pathological and structural studies of AD patients of mild to moderate severity (uncus and fusiform gyrus).22,23
The finding of atrophy of the anterior cingulate gyrus was unexpected as most studies have found involvement of the posterior cingulate gyrus.24,25 However, the anterior cingulate gyrus is involved in complex attentional tasks,26 and attention is known to be the first non-memory domain to be affected in AD, before deficits in language and visuospatial functions appear.27 With the use of structural and functional imaging, other investigators have found that anterior cingulate atrophy is predictive of conversion to AD in memory impaired individuals,28,29 suggesting that this region might be affected relatively early in the course of the disease.
Caudate atrophy in AD has been reported by other investigators with voxel based morphometry10,11 and has been explained in view of the demonstrated involvement of striatial nuclei by neurofibrillary changes.30 Recently, in 21 AD patients with a mean MMSE of 16, Barber and colleagues have found a 9%–14% reduction of the volume of the caudate nucleus.31 Further evidence that the caudate nucleus is abnormal in AD comes from MR studies showing an increase of iron lad comparable to that of Huntington's disease.32
The findings in the subgroup of three mildest AD patients are largely but not completely overlapping with those in the whole group. The comparatively looser significance threshold that we used (p<0.0001 uncorrected versus p<0.05 corrected), might lead to the hypothesis that some significant voxels might be false positives. However, with the exception of one voxel in the right precentral gyrus and two in the left medial frontal gyrus, all other positive voxels were located in regions found atrophic also in the experiment with the larger group. Moreover, the voxels in the right precentral and left medial frontal gyrus were those with the lowest significance values (table 3). Lastly, when a voxel by voxel comparison was performed between the results of the two analyses, voxels significant in both analyses were located in the medial temporal region. This supports the view that voxel based morphometry with SPM might be sensitive to atrophy even in the mildest phases of AD, although present algorithms are still unsatisfactorily specific. Moreover, these data will need to be confirmed in a study where significant voxels survive correction for multiple comparisons, which will probably require large number of mild AD cases.
Although the regions that we found correlated with global cognitive severity are those expected based on current knowledge, we found a striking and unexpected asymmetry (right greater than left). This is not in agreement with the notion that the MMSE is heavily weighted on verbal rather than non-verbal tasks and with most previous studies with functional33 and traditional structural imaging.34 However, others have indicated that the MMSE is more affected by right than left focal hemispheric lesions35 and that structural markers of severity in AD are more strictly correlated with right than left temporoparietal hypometabolism.36 To our knowledge, this is the first study of the correlation between MMSE and atrophy with voxel based morphometry, and replication is mandatory.
The main limitation of this study relates to the arbitrary setting of the significance threshold for the detection of grey matter density reduction. However, this is an intrinsic limitation of SPM and studies published to date have all used different thresholds.10,11,12,37 For each voxel SPM computes a t statistics denoting the difference of mean grey matter density between two groups. As is the case of the traditional t statistics, probability is sensitive to group size in that experiments with larger groups have more power to detect effect differences than experiments with smaller groups. In experiments with smaller groups or effect size, the lower power can be accounted for by setting a lower probability threshold. This increases sensitivity, but at the expense of lower specificity. This problem is particularly cogent to the correlation between MMSE and atrophy, where the poor specificity is an unsolved problem. The implementation of Bayesian analysis into SPM will help investigate this issue.38
The quest for observer independent tools to assess brain changes in degenerative brain diseases has just begun. The increasing use of voxel based morphometry with SPM in brain disorders39–43 is attributable to the high independence from the human operator in that once the preprocessing analysis parameters and the significance threshold of the analysis have been set, the process is largely automatic. However, its advantages and limitations are only beginning to be elucidated and more work will be necessary. The development of other techniques such as deformation based and tensor based4–6 as alternatives to or integrated into the voxel based morphometry approach6 might provide additional information that will hopefully prove useful in clinical and research settings.
We are indebted to Dr Maria Luisa Gorno Tempini for valuable comments on SPM analysis.
Funding: this study was funded by the Research Council for Health of the Academy of Finland, and EVO grants 5510, 5772722 from the Kuopio University Hospital.
Competing interests: none declared.