Article Text


A voxel based morphometry study on mild cognitive impairment
  1. C Pennanen1,3,
  2. C Testa6,8,
  3. M P Laakso1,2,3,6,
  4. M Hallikainen1,
  5. E-L Helkala4,
  6. T Hänninen1,
  7. M Kivipelto1,
  8. M Könönen2,
  9. A Nissinen5,
  10. S Tervo1,3,
  11. M Vanhanen1,
  12. R Vanninen2,
  13. G B Frisoni6,7,
  14. H Soininen1,3
  1. 1Department of Neurology, Kuopio University Hospital, Kuopio, Finland
  2. 2Department of Clinical Radiology, Kuopio University Hospital, Kuopio, Finland
  3. 3Department of Neuroscience and Neurology, University of Kuopio, Kuopio, Finland
  4. 4Department of Public Health, University of Kuopio, Kuopio, Finland
  5. 5The National Institute of Public Health, Helsinki, Finland
  6. 6The Laboratory of Epidemiology & Neuroimaging, IRCCS San Giovanni di Dio-FBF, Brescia, Italy
  7. 7AFaR- Associazione Fatebenefratelli per la Ricerca, Rome, Italy
  8. 8The Machine Vision Laboratory, Department of Mathematics and Computer Science, University of Udine, Udine, Italy
  1. Correspondence to:
 Dr H Soininen
 Department of Neurology, Bldg. 5, Kuopio University Hospital, POB 1777, 70211 Kuopio, Finland;


Background: Mild cognitive impairment (MCI) is the most widely used concept in classifying cognitive impairment in the elderly who do not fulfil the criteria for dementia. MCI is considered to confer an increased risk of progressing to dementia and most often Alzheimer’s disease (AD). Various approaches such as imaging of the brain have been applied to predict the conversion of MCI to dementia. A number of volumetric magnetic resonance imaging (MRI) studies have detected atrophy of the medial temporal lobe in subjects with MCI, but for the other cerebral regions the results have been inconsistent.

Objective: To study the pattern of brain atrophy in MCI.

Methods: Thirty two controls and 51 individuals with MCI deriving from population based cohorts were studied by MRI using voxel based morphometry. The threshold of t maps was set at p<0.001.

Results: Individuals with MCI had significant unilateral atrophy in the medial temporal lobe on the right side. Less extensive atrophy was found elsewhere—for example, in the temporal lobe, left superior parietal lobule, left anterior cingulate gyrus, and bilaterally in the thalami.

Conclusions: The MRI findings in MCI resemble those seen in early AD.

  • AD, Alzheimer’s disease
  • CDR, clinical dementia rating
  • CSF, cerebrospinal fluid
  • GM, grey matter
  • HERA, hemispheric encoding/retrieval asymmetry
  • MCI, mild cognitive impairment
  • MRI, magnetic resonance imaging
  • VBM, voxel based morphometry
  • WM, white matter
  • MRI
  • mild cognitive impairment
  • voxel based method

Statistics from

Mild cognitive impairment (MCI) is currently the most widely used concept in classifying the cognitive impairment in the elderly who do not fulfil the criteria for dementia. This is important and challenging from at least two points of view. Firstly, in 2001 the Report of the Quality Standards Subcommittee of the American Academy of Neurology reviewed a number of studies and, even though these studies had used various criteria for MCI, they pointed to an annual conversion rate of 6%–25% from MCI to Alzheimer’s disease (AD), a rate greatly exceeding that seen in the normal age matched population.1 Thus, investigating MCI might provide a means to study AD in its earliest phases. Secondly, if MCI were indeed to be a significant predictor of future AD, this might prove beneficial in terms of preventive or interventive measures for AD. However, some recent longitudinal population based studies have cast some doubt on the concept of MCI. Ritchie et al suggested that MCI is a poor predictor of dementia; that the group is unstable, with cases changing category almost yearly, and called for modifications to the current criteria.2 In the study by Larrieu et al, an annual conversion rate of 8.3% was observed during a five year period, but again the cases had a tendency to fluctuate, and far more than 40% of MCI cases reverted to normal instead of progressing to dementia during follow up.3

One way of addressing the issue could be to use imaging studies, which do seem to have some potential to detect early AD. Indeed, a number of volumetric magnetic resonance imaging (MRI) studies have detected temporal atrophy (the most common MRI finding in AD) to some degree in cases with MCI as well; however, for the other cerebral regions the results have been rather inconsistent.4,5 In addition to the traditional volumetric methods, recent advances in image analysis have produced new tools to map the entire brain instead of being restricted to single regions. These techniques can be done with relatively little individual effort and have high reproducibility.6–8 Perhaps the most widely used of these is voxel based morphometry (VBM), which has not been applied to MCI, with one exception9 where 22 patients with amnestic MCI were compared with 22 healthy controls and found to have significant grey matter loss predominantly affecting the hippocampal region and anterior cingulate gyri, and extending into the temporal neocortex. Therefore, the aim of our study was to try to replicate the findings of the preliminary study of MCI in a larger sample, by using more recent techniques related to VBM, and by examining a group of cases with MCI who fulfil probably the most dominant criteria (the first version of the Mayo criteria) for MCI.



The study included 51 individuals with MCI and 32 controls, derived from two population based cohorts in which cognitive functions of the elderly had been evaluated in the Department of Neurology, Kuopio University Hospital, Kuopio, Finland.10,11 In both cohorts, the evaluation consisted of a structured interview including clinical dementia rating (CDR) scale and an extensive neuropsychological assessment described in detail elsewhere.12 The scoring of the CDR was independent of the scores obtained from neuropsychological tests.

The controls showed no impairment in the cognitive tests, and had no history of neurological or psychiatric diseases. MCI was diagnosed using the criteria proposed by Mayo Clinic Alzheimer’s Disease Research Center. Recently, these criteria have been modified, but when our study was conducted the criteria required: (1) memory complaint by patient, family, or physician; (2) normal activities of daily living; (3) normal global cognitive function; (4) objective impairment in memory or in one other area of cognitive function as evident by scores >1.5 SD below the age appropriate mean; (5) CDR score of 0.5, and (6) absence of dementia.13,14 The ApoE genotype was determined from blood leukocytes as described previously.15

The study was approved by the local ethics committee, and all the participants gave informed consent for their participation in the study.

Magnetic resonance imaging

Patients and controls underwent a high resolution MRI at the Kuopio University Hospital, Kuopio, Finland. The MR images were acquired using a 1.5 Tesla Vision Magnetom scanner (Siemens, Erlangen, Germany), with a 3D gradient echo technique (TR = 9.7 ms, TE = 4 ms, flip angle  = 10°, field of view  = 250 mm, and the acquisition matrix was 256×256 resulting in 128 2.0 mm thick T1 weighted coronal slices). T2 weighted images were also acquired, and used to study the possibility of vascular pathology in the cases. None of the subjects was excluded from the study on the basis of T2 weighted images. To control for differences in head size, intracranial area measured at the level of the anterior commissure was used.12

Voxel based morphometry preprocessing

After removing the voxels below the cerebellum with MRIcro software ( the MR images were analysed with SPM99 ( running under Matlab 5.3 (Mathworks, Sherborn, MA, USA).

The MR images were preprocessed following a protocol which includes (1) generation of a customised template, (2) generation of customised prior probability maps, and (3) the main VBM steps: normalisation of the original MR images, segmentation of normalised images, cleaning of grey matter images, modulation of grey matter images, and smoothing of modulated images.

The advantages of this procedure over the traditional simple protocol are that (1) the use of template and prior probability maps computed from the population under study reduces error in the normalisation and segmentation steps, (2) the cleaning step permits removal non-brain voxels erroneously classified as grey matter,17 and (3) modulation permits preservation of the original volume of grey matter within each voxel.7,8,18

Customised template

The customised template was obtained by normalising the images of the 51 MCI patients and 32 controls to the Montreal Neurological Institute template19 of SPM99 using a 12 parameter affine transformation, smoothing the normalised images with an 8 mm isotropic Gaussian kernel and averaging the smoothed images. The anterior commissure was manually set as the origin of the spatial coordinates and images were reoriented coronally perpendicular to the intercommissural line. The normalisation procedure uses a bilinear interpolation algorithm to reslice images to voxel size of 2×2×2 mm. This voxel size is used in the following processing and analysis.

Customised prior probability maps

Customised prior probability maps were computed by segmenting the normalised images into grey matter (GM), white matter (WM), and cerebrospinal fluid (CSF), then smoothing with an 8 mm Gaussian filter, and finally averaging the segmented images, thus obtaining the customised prior probability maps specific for GM, WM, and CSF.17

Main VBM steps

Original images were normalised to the customised template through affine and non-linear transformations, medium regularisation, reslicing 2×2×2 mm, and no masking.20 The normalised images were segmented into GM, WM, and CSF using the customised prior probability maps. The Xbrain routine, based on erosions and dilatations, was used to remove voxels of non-brain tissue from the segmented images, thus obtaining a brain mask to clean the GM images by intersection with the mask.

In the modulation step, voxel values of the cleaned GM images were multiplied by the measure of relative volumes of warped and unwarped structures derived from the non-linear step of spatial normalisation (Jacobian determinant).17 The modulated GM images were smoothed with an 12 mm isotropic Gaussian kernel.

The final output is a 3D matrix where the three indices are the spatial x, y, and z coordinates of voxels in the reference space, and each value of the matrix is proportional to the volume of GM within each voxel. It should be emphasised that the output of each stage of the analysis was visually checked to make sure that the algorithms had actually carried out the expected changes.

Statistical analysis

A “single subject: conditions and covariates” model was used to compare the volume of GM between MCI and controls on a voxel by voxel basis. The analysis was controlled for apolipoprotein E (apoE) genotype including genotype (ε3/3, ε3/4ε, and ε4/4) as a dummy variable. Intracranial area, age, and sex were included as nuisance covariates. The resulting t map was thresholded at p<0.001, uncorrected. The demographic and clinical data were analysed using one way ANOVA or χ2 test, the level of significant was set at p<0.05.


There was no difference in age (F = 1.96, p = 0.17), sex (χ2 = 0.45, p = 0.50), education (F = 0.20, p = 0.66), or frequency of the apoE ε4 allele (χ2 = 2.5, p = 0.11) between the controls and cases with MCI. As expected, the Mini-Mental State Examination scores were significantly lower in cases with MCI compared with the controls (F = 24.8, p<0.0001) (table 1).

Table 1

 Clinical features of the MCI* cases and controls

Compared with controls, the cases with MCI had the greatest atrophy in the right hippocampus-amygdala region, and in the right hippocampal tail and thalamus (table 2, fig 1). The medial temporal structures on the left did not quite achieve the level of statistical significance. Other regions with significant GM atrophy were found in the right superior temporal gyrus, the left thalamus, the left inferior temporal gyrus, and the left anterior cingulate gyrus (table 2).

Table 2

 Atrophic regions of MCI cases compared with controls (p<0.001, uncorrected)

Figure 1

 Medial temporal involvement in MCI patients compared with controls with p<0.001, uncorrected. Significant voxels are superimposed on MR images of a single control subject. Atrophy affects the amygdala and the head and tail of the hippocampus.


In this study we observed that individuals with MCI have unilateral atrophy in the medial temporal lobe on the right side. Less extensive areas of atrophy were found elsewhere in the temporal lobe, left superior parietal lobule, left cingulate gyrus and, notably, in the thalami bilaterally.

To some degree these findings fit comfortably into the widely accepted theory of a hierarchical distribution of AD pathology—that is, the pathology first starts on the medial temporal lobe, and then spreads to cortical (parietal) association areas.21 Our observation, indicating a spread from a predominantly atrophied right medial temporal lobe to the left parietal lobes is in accordance, to some extent, with the so called functional HERA (hemispheric encoding/retrieval asymmetry) model.22 Although the HERA model is usually considered to involve mostly the medial temporal lobe and the prefrontal cortices, functional imaging studies have also revealed parietal involvement.23,24 Thus, as far as the medial temporal lobe and the parietal cortex are considered, these data suggest that MCI does indeed represent preclinical AD. However, we must remember that MCI is a heterogeneous condition and our cohort may include individuals who do not progress to dementia.

To the best of our knowledge, there is only one published study using VBM in MCI,9 but with a slightly smaller sample and different operational criteria for MCI. Their data are in accordance with ours in that they found bilateral thalamic atrophy and some atrophy in the temporal neocortex. Moreover, we found a small area of atrophy in the anterior cingulate cortex. This is of interest because in animal models of AD it is claimed entorhinal lesions that may be the initial site of AD pathology, lead to hypometabolism of the cingulate.25 In addition, in humans, it should be noted that hypoperfusion of the cingulate has been proposed to predict conversion from MCI to AD.26,27 The minor differences between the two studies can be explained by the participants studied, the definition of MCI, and differences in image analysis. However, it would seem that the results between the two studies, in general, are more convergent than divergent.

With respect to the temporal lobe in MCI (and AD), the involvement of regions other than the medial temporal lobe remains a matter of dispute, at least in volumetric studies.4,5 It also should be emphasised that the data obtained by volumetry versus VBM can be highly divergent.17,28,29 Thus, a direct comparison between volumetric and VBM studies may be difficult. In fact, this proviso also applies to comparison between different VBM studies. A good example of this is the thalamic involvement found in this study. In our previous VBM study on AD, no thalamic involvement was detected,30 nor has it been noted in some other VBM studies on AD,17,31,32 though others have detected thalamic involvement.20,29,33 A recent PET study on MCI did hint at thalamic hypoperfusion even when controlled for atrophy seen in MR images.34 The pathology of AD has been suggested to also affect the thalamus,35 yet this issue, and the factors underlying the discrepancy, remain to be resolved in future studies. One possibility that is worth mentioning is the bias inherent in the method. As Karas et al stated, the thalamus is one of the structures which “lie geographically in locations which are hard to evaluate with computational models”.29

Indeed, it should be kept in mind that VBM is a very sensitive and complex method that can be greatly affected by a number of incidental factors. These might explain some discrepancies between our own and other studies such as low quality scans, the reduction of grey/white matter contrast of the elderly or diseased people, and the greater frequency of motion effects in these groups. In addition to these factors, spatial normalisation still produces registration errors with regional specificity: small structures (such as the hippocampus) are highly variable particularly in diseased groups, which might reduce the sensitivity of VBM to detect changes. Methods such as tensor based mapping6 might overcome this problem.


In this study we found a unilateral medial temporal atrophy in individuals with MCI. This is a consistent finding in volumetric studies on MCI and invariably detected in volumetric and VBM studies on AD, where the atrophy is bilateral. Some other regions were also found to be atrophic in MCI, but these were relatively small, and it might be premature to draw definitive conclusions based on those findings. However, this study further emphasises the importance of detecting medial temporal atrophy in AD. It also suggests that, as a group, individuals with MCI, defined according to the first version of the Mayo criteria, are more likely to progress to dementia and AD than their age matched peers. Nonetheless, caution must be exercised when interpreting the results due to the heterogeneous nature of MCI. A follow up of the participants of this study is hoped to further clarify this issue.


The study was supported by the Research Council for Health of the Academy of Finland, and an EVO grant 5772720 from the Kuopio University Hospital.


View Abstract


  • Competing interests: none declared

Request permissions

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.