Original Research ARTICLE
Assessing age-related gray matter decline with voxel-based morphometry depends significantly on segmentation and normalization procedures
- 1Movement Control and Neuroplasticity Research Group, Department of Kinesiology, KU Leuven, Belgium
- 2CNRS, INCIA, UMR 5287, University of Bordeaux, Talence, France
- 3Department of Electrical Engineering – ESAT – PSI & iMinds - Future Health Department, KU Leuven, Belgium
- 4Neural Control of Movement Laboratory, Health Sciences and Technology, ETH Zurich, Switzerland
Healthy ageing coincides with a progressive decline of brain gray matter (GM) ultimately affecting the entire brain. For a long time, manual delineation-based volumetry within predefined regions of interest (ROI) has been the gold standard for assessing such degeneration. Voxel-Based Morphometry (VBM) offers an automated alternative approach that, however, relies critically on the segmentation and spatial normalization of a large collection of images from different subjects. This can be achieved via different algorithms, with SPM5/SPM8, DARTEL of SPM8 and FSL tools (FAST, FNIRT) being three of the most frequently used. We complemented these voxel based measurements with a ROI based approach, whereby the ROIs are defined by transforms of an atlas (containing different tissue probability maps as well as predefined anatomic labels) to the individual subject images in order to obtain volumetric information at the level of the whole brain or within separate ROIs. Comparing GM decline between 21 young subjects (mean age 23) and 18 elderly (mean age 66) revealed that volumetric measurements differed significantly between methods. The unified segmentation/normalization of SPM5/SPM8 revealed the largest age-related differences and DARTEL the smallest, with FSL being more similar to the DARTEL approach. Method specific differences were substantial after segmentation and most pronounced for the cortical structures in close vicinity to major sulci and fissures. Our findings suggest that algorithms that provide only limited degrees of freedom for local deformations (such as the unified segmentation and normalization of SPM5/SPM8) tend to overestimate between-group differences in VBM results when compared to methods providing more flexible warping. This difference seems to be most pronounced if the anatomy of one of the groups deviates from custom templates, a finding that is of particular importance when results are compared across studies using different VBM methods.
The human brain undergoes continuous structural changes due to development and aging. Post mortem studies have shown that even healthy aging is accompanied by notable cortical atrophy and loss of brain weight from the sixth life decade onwards (Skullerud, 1985). This brain tissue loss is further accelerated by most neurodegenerative or neuropsychiatric disorders such that magnetic resonance imaging (MRI) in combination with morphometry offers an interesting diagnostic tool for differentiating between healthy and pathological brain changes (for an overview, see Raz and Rodrigue, 2006; Mueller et al., 2012a,b). One common approach is to evaluate the volume of gray matter (GM), white matter (WM) and cerebrospinal fluid (CSF) within structures of interest determined through manual delineation. Even though this approach has long been considered the “gold standard” it is very labor-intensive and requires high anatomical expertise such that investigations were frequently restricted to a limited set of predefined anatomical regions of interest (ROI). Alternatively, semi-automated procedures were incorporated for detecting region boundaries and reducing observer-dependence (Mu et al., 1999; Jernigan et al., 2001; Tisserand et al., 2002; Raz et al., 2003; Allen et al., 2005a; Walhovd et al., 2005, 2010; Greenberg et al., 2008; Hasan and Pedraza, 2009). During the last decade, volumetry of predefined brain structures has been complemented by whole-brain, voxel-based methods like voxel-based morphometry (VBM; Good et al., 2001; Ashburner and Friston, 2000, 2001; Tisserand et al., 2004; Gonoi et al., 2010), or estimating voxel-based cortical thickness (Hutton et al., 2008). Particularly VBM has been frequently used for investigating GM decline in the context of neurodegenerative pathologies (e.g., Good et al., 2002; Testa et al., 2004; Chetelat et al., 2005; Giuliani et al., 2005; Douaud et al., 2006; Keller and Roberts, 2008; Davies et al., 2009; Ibarretxe-Bilbao et al., 2009; Ferreira et al., 2011) and healthy ageing (Good et al., 2001; Van Laere and Dierckx, 2001; Resnick et al., 2003; Tisserand et al., 2004; Grieve et al., 2005; Jernigan and Gamst, 2005; Lemaıtre et al., 2005; Lehmbeck et al., 2006; Fjell et al., 2009; Kalpouzos et al., 2009; Kennedy et al., 2009; Ziegler et al., 2012). In a nutshell, VBM assumes that each voxel of a T1-weighted high resolution MR image contains a mixture of GM, WM, and CSF. The image is segmented to yield the content or probability of each tissue class at the voxel level. After transforming each individual brain image into the same stereotactic space (thus correcting for global differences in position, size, and shape across individuals) these voxel-based tissue probabilities can be statistically compared between different populations to determine local alterations of brain structure.
VBM has been frequently applied to show that healthy aging is accompanied by GM tissue loss. Even though changes are distributed across the brain there is particular vulnerability with respect to the frontal lobes (Good et al., 2001; Van Laere and Dierckx, 2001; Tisserand et al., 2004; Grieve et al., 2005; Lehmbeck et al., 2006; Kennedy et al., 2009; Kalpouzos et al., 2009). This general finding was further confirmed by manual delineation and semi-automated ROI-based studies (Jernigan et al., 2001; Raz et al., 2004; Allen et al., 2005a; Walhovd et al., 2005, 2010). However, it is worth noting that when VBM and manual delineations were applied to the same data set, results were only partly congruent. Investigating (pre)frontal subregions, Tisserand et al. (2002) consistently identified pronounced age-related GM reduction in fronto-lateral areas with both VBM and manual volumetry. Method-dependent differences, however, were revealed for orbital portions: whereas these exhibited age-effects when using manual delineation, the voxel-based approach identified only cingulate subregions. Allen et al. (2005b) assessed VBM against manual volumetry for insular cortex and detected a different rate of age-induced degeneration: while VBM indicated the pace of insular decline to exceed that of frontal cortex, the manual delineation results instead attributed a slower pattern of decline. In a seminal study, Kennedy et al. (2009) compared VBM to manual tracings of predefined ROIs in a large sample of 200 subjects. They demonstrated that VBM might reveal reliable information when aggregated within meaningful anatomical regions. However, when interpreted at the voxel level, VBM tended to overestimate age related GM differences, particularly for voxels bordering WM or CSF.
One important factor that may contribute to inconsistencies across studies is that VBM appears to be very sensitive to small methodological variations. Thus, even minor differences in spatial transformations or smoothing procedures might yield widely divergent results as has been demonstrated for longitudinal VBM studies in GM changes (Thomas et al., 2009; Thomas and Baker, 2013) and the analysis of WM structures (Jones et al., 2007). Here we revisit this issue for the comparison of brain structure and, particularly, GM volume between different age groups. The motivation for our study is that image analysis methods are constantly modified and improved, while it is rarely the case that published results are reanalyzed using the new algorithms. VBM relies critically on the segmentation and normalization of the structural images and different algorithms have been implemented to solve these problems. Two of the most frequently used software packages are SPM (Wellcome Trust Centre for Neuroimaging, London, UK) and FSL (Smith et al., 2004; Douaud et al., 2006), and both offer VBM data processing pipelines. Within SPM two versions are available: one uses the unified segmentation and normalization routines as implemented in standard SPM5/SPM8. This method determines GM, WM, and CSF tissue classes in an iterative procedure that warps the T1- weighted images of each individual subject to the ICBM152 space (i.e., MNI space) and segments them using SPM5/SPM8 standard prior tissue probability maps (Ashburner and Friston, 2005). One disadvantage of the SPM5/SPM8 procedure is that local, non-linear transformations are rather limited because warping is achieved through a restricted set of spatial basis functions (Ashburner and Friston, 2004) that are regularized by minimizing the bending energy of the deformation fields. In the widely used SPM8 release a new segmentation procedure has been implemented which determines additional tissue classes to reduce misclassification and uses prior tissue probability maps based on a large sample with a broad age range. Moreover, an alternative normalization procedure, known as DARTEL, uses a diffeomorphic algorithm to warp individual subject images into a common space that is defined in an iterative procedure to determine one common template optimized for the population under study (Ashburner, 2007; Ashburner and Friston, 2009). DARTEL is believed to allow more local, non-linear deformations than the standard SPM5/SPM8 normalization and to reveal better results in populations with deviant anatomy. FSL uses a segmentation procedure independent of prior information to determine the different tissue classes. Instead it applies a hidden Markov random field model together with an associated Expectation-Maximization algorithm to associate each intensity value from the anatomical image with a specific mixture of GM, WM and CSF probabilities (Zhang et al., 2001). Normalization is achieved by combining affine with nonlinear registration (Andersson et al., 2007a,b). The use of a b-spline representation of the registration warp field (Rueckert et al., 1999) allows application of highly local deformations in order to match individual anatomy to a common template.
Here we address the question of how far methodological differences regarding the segmentation (i.e., partitioning anatomical images in GM, WM, and CSF probability maps) and normalization of MRI images (i.e., transforming the individual anatomical images to a template) influence VBM results. We use the aging brain as a model for a significant reduction in GM volume that has been consistently demonstrated by a large number of previous studies using various methods. Moreover, we compared two groups (young vs. elderly subjects) at the extreme spectrum of the age range to compare how the different methods influence the effect size of age-related structural changes.
In addition to those software packages freely available, we complemented the VBM results with an automated ROI approach using an independent algorithm to determine volumetric measurements within predefined anatomical structures (D'Agostino et al., 2003, 2006a,b; Stiers et al., 2005; Machilsen et al., 2007).
Materials and Methods
Thirty-nine subjects, 21 young (23 ± 1.7 years; range 20–27; 8 male/13 female) and 18 elderly (66.2 ± 3.4 years; range 62–73; 9 male/9 female), participated in the study. All subjects were strongly right-handed, as determined by the Edinburgh handedness questionnaire (Oldfield, 1971) and none of the participants reported a history of neurological disease, sensorimotor dysfunction or used psycho- or vasoactive medication. General functions were assessed using the Mini Mental State examination and all scores ranged within normal limits (score ≥26). Participants were informed about the experimental procedures and provided written informed consent. The informed consent and the study design were approved by the local Ethics Committee of Biomedical Research at KU Leuven, in accordance with the ethical standards laid down in the 1964 Declaration of Helsinki.
Whole-brain T1-weighted images were collected on a 3-T Intera MRI scanner (Philips, Best, Netherlands), using a six element SENSE head coil (MRI Devices Corp., Waukesha, WI). Each volume consisted of 182 slices, with echo time (TE) 4.6 ms, repetition time (TR) 9.68/9.58 ms, inversion time (TI) 1100 ms, field of view (FOV) 250 mm, matrix 256 × 256, in plane resolution of 0.98 × 0.98 mm2, slice thickness 1.2 mm, and SENSE factor 2. None of the participants exhibited abnormalities in brain structure.
We analyzed the data with four different approaches: (1) the unified segmentation/normalization algorithm “SPM8 segment” as implemented in SPM5/SPM8 (subsequently called SPM VBM), (2) DARTEL using the “SPM8 new segment” and normalization procedure implemented in SPM8 (subsequently called DARTEL VBM), (3) tools of the FSL package (subsequently called FSL VBM), and (4) a procedure combining Intensity-Based Segmentation and atlas-to-image Non-Rigid Registration (IBSNRR) with volumetry for predefined ROIs. For each method, we followed the procedures advised by the software developers to generate GM, WM, and CSF segments in individual subject space. For the first three methods we also determined a set of warping parameters defining the transformation from individual subject space to MNI ICBM152 template space.
All SPM analyses were performed with SPM8 (v4667, Wellcome Trust Centre for Neuroimaging, London, UK) run in Matlab Version 18.104.22.1684 (R2008a). Segmentation of the images was performed by the unified segmentation routine “SPM8 segment” as implemented in SPM5/SPM8 (Ashburner and Friston, 2005). Briefly, this procedure segments the T1-weighted image in GM, WM, and CSF tissue classes, accomplished by iterative registration to ICBM152 space (i.e., MNI template space) and segmentation using SPM5/SPM8 standard prior tissue probability maps. It is worth noting, however, that the standard SPM5/SPM8 probability maps were derived based on a sample of young adults only and, thus, might introduce a bias when examining age effects (Figure 1A). For each participant, the segmentation procedure yields GM, WM, and CSF probability maps in individual subject space as well as deformation fields describing the transformation of the tissue priors to the native subject images. The GM probability maps were warped to standard ICBM152 space by applying the inverse of the deformation fields determined during segmentation.
Figure 1. GM probability maps used by the unified segmentation and normalization procedure of SPM5/SPM8 (A) and the IBSNRR approach (B). Note that the IBSNRR atlas provides much more anatomical detail than the GM probability map in (A). Labels of the IBSNRR atlas are shown in (C). Anatomical structures are numbered in accordance to the label names (numbers are provided in brackets).
Segmentation was performed using the unified segmentation (Ashburner and Friston, 2005) as implemented in the “New segment” function of SPM8. This segmentation determines more tissue classes than standard SPM5/SPM8 including GM, WM, CSF, soft tissue, skull, and out-of-brain regions of the image. Prior probability maps were derived from a large sample of healthy adults across the lifespan (Good et al., 2001) that are available in ICBM152 space. First, a customized template was created based on the GM segmentations using a diffeomorphic method known as DARTEL (Ashburner, 2007; Ashburner and Friston, 2009). DARTEL uses an iterative process to progressively refine warping parameters describing the transformation from individual subject space to a common space. During this procedure flow fields were calculated and subsequently applied to the native GM image of each subject. Next, the affine transformation from the customized template space to ICBM152 space was calculated and applied to the individual GM segmentations.
Using the FSL 4.1.5 tools (Smith et al., 2004), all individual T1 weighted images were brain-extracted using BET (Smith, 2002) and visually inspected. Then tissue-type segmentation was carried out using FAST v4.1 (Zhang et al., 2001). Importantly, FAST does not use prior information to determine the different tissue classes but applies a hidden Markov random field model together with an associated Expectation-Maximization algorithm to associate each intensity value from the anatomical image with a specific mixture of GM, WM, and CSF probabilities. Then a study specific template was created. Therefore, the resulting GM probability images were normalized to ICBM152 standard space using the standard FSL template and the affine registration tool FLIRT (Jenkinson and Smith, 2001; Jenkinson et al., 2002), followed by nonlinear registration executed by FNIRT (Andersson et al., 2007a,b), which uses a b-spline representation of the registration warp field (Rueckert et al., 1999). The resulting images were flipped and averaged to create a symmetric, study-specific template. In a second iteration, the native GM probability maps were then non-linearly re-registered to the customized template.
IBSNRR: Intensity-based Segmentation and Non-rigid Registration
The general idea of IBSNRR is that tissue class priors of anatomical ROIs are predefined by a statistical brain atlas that gets transformed to individual subject space. As a result the shape of each structural ROI is adjusted to individual anatomy and GM, WM, CSF and non-brain tissue probability maps are obtained within predefined structures in native subject space (Stiers et al., 2005; Machilsen et al., 2007). Subsequently, tissue volumes are quantified within each ROI. As such, this method implements automated ROI volumetry rather than a VBM approach. IBSNRR uses tissue class priors for segmentation derived from a statistical brain atlas constructed from 20 manually labeled brain images of 10 males and 10 females (median age 31 years) acquired at the Imperial College, Hammersmith Hospital (London), at an average voxel size of 0.937 mm3 (Hammers et al., 2002, 2007; Wang et al., 2005). This atlas consists of a T1-weighted image, with corresponding probability maps for global tissue classes of GM, WM, and CSF (Figure 1B) and for 49 predefined anatomical structures (Figure 1C). Initial affine atlas-to-image registration was performed based upon the atlas mean gray-scale image and the T1-weighted subject images, using mutual information (Maes et al., 1997). Next, segmentation occurred within a unified framework that alternates between an intensity-based segmentation step using a Gaussian mixture model (analogous to Van Leemput et al., 1999) and a non-rigid registration step that matches the tissue class maps of the atlas to the segmented subject images using a viscous-fluid model (D'Agostino et al., 2003, 2006a,b). This algorithm yields atlas-to-image deformation fields allowing for local non-linear transformations that outperform direct atlas-to-image registrations. By applying these deformation fields to the tissue probability maps of the atlas anatomical structures, a segmentation solution was obtained. Segmentation results were visually inspected by two independent observers and no major misclassifications were detected.
Data Analysis and Statistics
Whole brain volumes
All four methods yielded GM, WM, and whole brain CSF probability maps in native space. From these images, total GM, WM, and CSF volumes were determined by thresholding each tissue probability map at p ≥ 0.2, multiplying the tissue probability with the volume of the voxel and integrating these values across the brain. Subsequently, these volumes were subjected to an analysis of variance for repeated measurements (repeated measures ANOVA) with the between factor AGE (young, old) and METH (SPM, DARTEL, FSL, IBSNRR). One important confound for between-group analyses of brain volume is the average head size of the participants. For example, if our elderly subject sample had larger brains on average, it would not be surprising if they also had more GM volume, since brain size and GM volume are expected to be proportionally related. To prevent this potential bias we estimated the intracranial volume (ICV) for each participant. We aimed to have an ICV estimate independent of the 4 evaluated segmentation methods and therefore used the dilated brain mask of the ICBM152 template, available from the FSL software package. This mask was transformed to native subject space based on the inverted, non-linear warping parameters derived from the FSL normalization procedure and visually inspected (see Supplementary Table 1 for additional justification and analyses motivating this specific approach). Subsequently, ICVwas estimated by integrating the volume across all non-zero voxels of the mask and ICV was used when whole brain GM, WM, and CSF volumes were statistically compared either by adding it as a covariate of no interest, or by proportionally scaling GM, WM, and CSF volumes to the ICV volume [i.e., GM (%) = GM (l)/ICV (l) *100] prior to the statistical comparisons.
Next we asked whether the distribution of GM, WM, and CSF probabilities would differ between methods. Distributions were determined for each method, age group and tissue class by counting voxels with p ≥ 0.05 and p ≤ 1 in bins of 0.05, and expressing them as a percentage of ICV. Repeated measures ANOVAs with the factor BIN (0.05–0.1, …., 0.95–1) and METH (SPM, DARTEL, FSL, IBSNRR) were calculated separately for each tissue class and Tukey post-hoc tests were used to test for method specific differences within each bin. We then quantified age specific differences in the GM probability distribution. Therefore, the mean volume of the young subjects was calculated for each bin and the GM volume of each elderly participant was expressed as a fraction of this value. These values were subjected to a repeated measures ANOVA with the within factors BIN and METH.
For all analyses alpha was set to 0.05 and Tukey post-hoc tests were applied when appropriate. All statistical analyses above were performed with Statistica 10 (StatSoft Inc., Tulsa USA).
GM volumetry: interactions between age and segmentation method
GM segmentations derived by each of the four methods were thresholded at p ≥ 0.2, a minimum probability value often recommended for SPM VBM analyses. We wanted to quantify whether the segmentation method affects age specific GM differences either with respect to the localization or the size of the effect. Two analyses were carried out to answer this question. First, the thresholded GM probability maps derived by each of the four methods were normalized to ICBM152 space by applying the same non-rigid transformation determined by SPM8 to all images of one subject. Normalized images were resampled to 2 × 2 × 2 mm voxels using a nearest neighbor interpolation method and no modulation for resizing voxels was applied. Next the normalized GM probability maps were smoothed with a 10 mm FWHM Gausian kernel and subjected to a repeated measures ANOVA in SPM with the between factor AGE (young, old) and the within factor METH (SPM, DARTEL, FSL, ISBNRR). Finally we calculated all AGE × METH interaction effects and the resulting statistical parametric maps (Friston et al., 1994) were thresholded at puncor < 0.001 (i.e., uncorrected for multiple comparisons) and at pFWE < 0.05 (i.e., after applying family wise error correction). This analysis revealed where in the brain age-related GM differences depend significantly on the segmentation method. Note however, that normalization into the common ICBM 152 space might have been suboptimal for comparing young vs. elderly subjects.
Therefore, we performed a second complementary analysis determining mean GM probability within 21 bilateral ROIs derived by the IBSNRR method (for a list of these areas see Supplementary Table 1). More specifically, GM probability values were multiplied with the voxel volume and these values were integrated across all voxels within a given ROI and with p > 0.2 and subsequently divided by ICV. This procedure was performed for each individual and each method. Next, for each ROI and each method, data were z-transformed (i.e., the mean across all subjects was subtracted and this value was divided by the standard deviation). This allows the identification of the effect size of age specific differences while removing variability introduced by the segmentation method per se. For each ROI and each segmentation method, z-transformed GM volume of the young and the elderly were subjected to independent t-tests and a Bonferroni correction was applied to the resulting p-values. Additionally, Cohen's d effect sizes were calculated.
GM VBM: the influence of segmentation and normalization methods on age specific GM differences
Finally, GM probability maps derived from SPM5/SPM8, DARTEL, and FSL segmentation were thresholded at p > 0.2 and transformed to ICBM152 space by applying the warpings as obtained with SPM5/SPM8, DARTEL and FSL normalization, respectively (resulting in nine normalized GM data sets in total). Additionally, a modulation step was included to keep the total GM amount constant despite local expansion or contraction introduced by image normalization. For each of the 3 segmentation × 3 normalization data sets, a separate second level model was specified in SPM, calculating independent t-tests between young and elderly subjects while GM data were proportionally scaled for ICV. The contrast young > elderly was calculated and thresholded at puncor < 0.001 as well as PFWE < 0.05. These results were qualitatively compared to describe differences introduced by the various normalization procedures.
GM, WM, and CSF: Cross-method Comparison of Detecting Ageing Effects on Overall Volumes
GM, WM, CSF, and ICV are shown separately for the young and the elderly group in Table 1. For all tissue classes, there was a significant AGE × METH interaction irrespective of whether ICV was used as a covariate (upper part of the Table 1) or as proportional scaling factor (lower part of Table 1). More specifically, GM volumes were generally smaller in the elderly than in the young subjects. However, the extent of this age specific effect differed across methods. Testing the AGE × METH interaction pairwise between methods (and applying a Bonferroni correction to the p-values) revealed that SPM segmentation resulted in significantly larger age related differences in GM volume than all other methods (p < 0.005). By contrast, DARTEL segmentation resulted in significantly smaller age effects than all other methods (P < 0.005). These results were driven by the fact that SPM yielded exceptionally small GM volumes for the old group.
Table 1. Whole brain volumes (Vol) estimated based on tissue segmentations derived from SPM, DARTEL, FSL, or ISBNRR.
WM volumes were estimated much more similarly across methods. Small, but significant differences, were revealed because SPM segmentation indicated a larger increase in WM volume for the elderly as compared to the young than in the other methods. Accordingly, the age effect was significantly larger for SPM than for DARTEL or FSL segmentations (p < 0.05).
CSF volumes were larger in the elderly than in the young group, however this age effect was much more pronounced when estimated with SPM than with all other methods (p < 0.0003) because CSF estimates yielded by SPM were exceptionally large in the elderly.
In summary, we found large differences across segmentation methods when quantifying the effect of ageing on changes in overall tissue class volumes. With respect to GM, SPM tended to over-estimate age specific differences, while DARTEL tended to under-estimate these structural changes when compared to the other methods.
One might be concerned that age related effects on brain volume might be driven by partial volume effects, i.e., by voxels with a relatively low probability of belonging to a given tissue class. Therefore we determined the distributions of GM, WM and CSF probabilities as derived by either method and displayed them separately for each age group (Figure 2A). This analysis revealed several insights: First, for all tissue classes and all segmentation methods, by far the largest amount of voxels has a probability of p ≥ 0.9 (except for CSF derived by DARTEL). Second, distributions differ between segmentation methods. For example, with the FSL segmentation, more voxels fall into the probability range between 0.2 and 0.8 as compared to the other methods that use prior probability maps. Third, statistics revealed significant BIN × METH interactions. In particular, all methods differ significantly from each other regarding the number of voxels assigned p > 0.95, and this result was found for all tissue classes.
Figure 2. (A) Distributions of GM, WM and CSF probabilities in the young (open bars) and the elderly group (filled bars). Volumes per bin were normalized to ICV (Vnorm) and expressed as a percentage. Distributions were derived from the segmentation solutions provided by the unified segmentation/normalization of SPM5/SPM8 (upper panel), DARTEL (second panel), FSL (third panel) and IBSNRR (lower panel). Results of the statistics are symbolized by the dot diagram on top. For each bin, Tukey post-hoc tests, i.e., pairwise comparisons between all methods, were performed and significant results (p < 0.05) were indicated by a dot (S, SPM; D, DARTEL; F, FSL; I, IBSNRR). (B) Age specific differences in the GM probability distribution are shown for each method (SPM: yellow, DARTEL: red, FSL: green, ISBNRR: blue). For each bin, the GM volume of the elderly participants was expressed as a fraction of the mean volume of the young group. Significant between-method differences as tested by Tukey post-hoc tests (p < 0.05) are expressed by the dot diagram above the graph.
Next we focused on GM estimates and determined how the age effect is influenced by different segmentation methods (Figure 2B). Figure 2B shows GM volume of the elderly group as a fraction of the GM volume of the young. Thus, values below 1 indicate that elderly have less voxels than the average of the young group within a given probability bin, while values larger than 1 indicate that elderly have more voxels. Again, FSL differed from all other methods, as it revealed consistently lower GM volumes in the elderly than in the young for nearly all probability ranges. By contrast, SPM, DARTEL, and IBSNRR segmentation suggested a shift of GM distribution with age, such that the elderly had a relative larger number of voxels with low GM probabilities and a relatively lower number of voxels with high GM probabilities than the young. Statistics revealed a significant BIN × METH interaction further supporting our observation (see Figure 2 for statistics summary).
In summary, it is unlikely that age specific differences in GM volume are driven by voxels with relatively low GM probabilities. Instead, SPM, FSL, and IBSNRR segmentation find pronounced volume differences between elderly and young for voxels with high GM probability (p ≥ 0.9). Only DARTEL differs in this respect since the GM reduction in the elderly was mainly observed for voxel probabilities ranging between 0.6 and 0.9
GM Volumetry: Interactions between Age and Segmentation Method
After previous analyses have established that the different segmentation methods have a significant influence on quantifying age specific changes in overall GM volumes, we asked where these method-specific differences are located in the brain. To this end we transformed all GM probability maps from native space to ICBM152 space and, importantly, applied for each individual the same transformation parameters to all GM probability maps. Subsequently we determined the AGE × METH interaction. Note that this will not reveal general differences in GM probability across methods or areas exhibiting GM decline with age. Instead this analysis was performed to identify where in the brain the age related GM decline is over- or under-estimated when compared across segmentation methods. Comparing SPM to DARTEL segmentation, SPM revealed significantly larger GM differences between young and elderly throughout the surface of the brain (Figure 3A). When SPM was compared to FSL, differences in the aging effect were most pronounced for dorsal cortical region, particularly around the central sulcus (Figure 3B). By contrast when comparing SPM to IBSNRR (Figure 3C), differences tended to be located more inferiorly, for example around the sylvian fissure. DARTEL segmentation revealed also a significantly smaller age effect on GM when compared to FSL (Figure 3D) and IBSNRR (Figure 3E): differences with both methods were particularly located around the interhemispheric and the sylvian fissure even though this effect was much more pronounced when compared to FSL than to IBSNRR. Finally, also FSL and IBSNRR segmentation exhibited minor differences when quantifying age related GM decline, such that FSL tended to reveal larger age effects for inferior located cortical surface areas (particularly around the sylvian fissure, Figure 3F) and lower age related effects for the dorsal cortical areas (Figure 3G). In summary, method specific differences seem to arise mainly from differential segmentation results for the cortical surface and particularly close to large sulci and fissures. By contrast, differences were only minor for subcortical gray matter structures.
Figure 3. Method × Age interaction of GM probabilities at the voxel level after common normalization to MNI space (see Methods section for details). Voxels exhibiting significantly larger GM differences between young versus elderly group when analyzed with SPM than with DARTEL (A), than with FSL (B) and than with ISBNRR (C). Also, FSL (D) and ISBNRR (E) revealed for some voxels larger GM differences between young versus elderly subjects than DARTEL. There were only minor differences in age-specific GM decline between FSL and ISBNRR (F,G). Statistical parametric maps were either thresholded at p < 0.001 uncorrected (red) or p < 0.05 FWE corrected (yellow).
One has to keep in mind though, that GM segmentations were transformed from native to ICBN152 space. Even though this transformation was identical across segmentation methods, it differed across subjects and mis-registrations might have exaggerated or diminished age related effects at the voxel level. Therefore, we performed a complementary ROI-based analysis for the anatomical structures predefined in the IBSNRR atlas. This analysis used GM estimates for large anatomical structures and is performed in native space. GM differences were estimated by independent t-tests calculated for each ROI and segmentation method. T-values and effect sizes are reported in Table 2. It can be seen that SPM revealed significantly lower GM volume for elderly than for young for nearly all cortical ROIs containing large surface area and that the size of these effects was usually large (> 0.08). FSL and IBSNRR revealed age related effects for frontal lobe, cingulum, insula, and parietal lobe, but less so for temporal and occipital ROIs. Moreover, the size of the age related effects were often smaller than when estimated based on SPM segmentations. For the DARTEL segmentations none of the cortical ROIs reach statistical significance. However, this might be due to our very strict correction for multiple comparisons. When inspecting effect sizes, also the DARTEL segmentation revealed large age related effects regarding GM loss in frontal and parietal lobe (Cohen's d > 0.95), as well as for nucleus accumbens and caudate nucleus.
In summary, when age related GM decline is used as a model to compare different segmentation methods, SPM is the most liberal approach revealing the largest GM differences, whereas DARTEL revealed the most conservative results (i.e., smallest effect size). This difference is particularly pronounced for cortical surface areas neighboring large fissures or large sulci.
As a final control, we warped the GM segmentations to ICBM152 space using the warpings derived with SPM5/SPM8, DARTEL, and FSL. Figure 4 shows the VBM results in form of statistical parametric maps when GM differences are contrasted between the young and the elderly sample. Even though this analysis reveals purely qualitative information, it is apparent that the general results pattern is similar, but also that some differences persist after normalization, particularly when more stringent statistical thresholding is used. Figures 4A,E,I show the results revealed by the SPM5/SPM8 VBM, DARTEL VBM, and FSL VBM pipeline, respectively. Representative slices show that all methods reveal age-related differences around large sulci, like the lateral fissure, but that extent and peak location differ. In accordance to the segmentation results, DARTEL VBM revealed less age-related GM differences than the other methods. SPM5/SPM8 VBM and FSL VBM revealed age-related changes of similar extent but results of FSL VBM were located more medially and superior to those of SPM5/SPM8 VBM. Furthermore, we applied each normalization method to each segmentation result. Overall, SPM5/SPM8 normalization (Figures 4A,D,G) seems to further enhance age-related differences in comparison to DARTEL (Figures 4B,E,H) while the FSL provides an intermediate solution, particularly when thresholded at pFWE < 0.05 (note though that the FSL developers advise to use non-parametric statistics for VBM analyses). In summary, when age-related decline is estimated with VBM, method specific differences persist also after normalization.
Figure 4. GM segmentations derived by SPM (A–C), DARTEL (D–F), and FSL (G–I) were normalized using the warping and modulation parameters of SPM (A,D,G), DARTEL (B,E,H) and FSL (C,F,I). Statistical parametric maps were either thresholded at p < 0.001 uncorrected (red) or p < 0.05 FWE corrected (yellow).
VBM is a largely automated approach to quantify structural changes of the brain and an interesting tool for studying large populations. In the present study we asked whether cross-sectional VBM studies are significantly influenced by segmentation and normalization algorithms provided by different software packages frequently used by the neuroimaging community. We used structural changes due to aging as a well-established model of pronounced GM reduction particularly in the frontal lobe. Our findings show that mainly segmentation but also normalization procedures have a profound influence on VBM results. Overall age-specific differences in GM estimates were largest using standard SPM VBM as implemented in SPM5/SPM8 and smallest using DARTEL VBM. Even with a rather small sample as tested here, differences between methods reached significance which is important to consider when comparing studies that used different software packages for VBM analyses.
Ageing Effects of Total GM, WM, and CSF Volume
Irrespective of which segmentation method was used, we found lower GM and higher CSF volume in the elderly as compared to the young (after correction for ICV). This result was highly expected based on previous literature (Good et al., 2001; Van Laere and Dierckx, 2001; Resnick et al., 2003; Tisserand et al., 2004; Grieve et al., 2005; Lemaıtre et al., 2005; Lehmbeck et al., 2006; Kalpouzos et al., 2009; Kennedy et al., 2009). However, regarding GM volume, we found large differences across segmentation methods when quantifying the effect size of the age-related decline: SPM tended to over-estimate age specific differences, while DARTEL tended to under-estimate these structural changes when compared to the other methods. This pronounced difference between SPM5/SPM8 and DARTEL segmentation has been reported before based on a much larger sample including data across the lifespan (see supplementary material in Peelle et al., 2012). One has to note however, that we tested two groups at the extreme end of the age range and this selection contributed to the large difference between SPM and DARTEL segmentation since the analysis of Peelle et al. (2012) indicated that the methodological differences regarding GM volume estimation are smaller when less extreme age groups are compared. Our analysis suggests that methods using a-priory probability maps during the segmentation algorithm (i.e., SPM, DARTEL, and IBSNRR) revealed a slightly shifted distribution of GM probabilities in the elderly: in the elderly we found more voxels with low GM probability and less with high GM probability than in the young. However, the largest part of the volume loss was driven by a profound reduction of high GM probability voxels.
All approaches revealed an increased CSF volume in the elderly vs. young subjects. SPM5/SPM8 segmentation however revealed distinctly more variable and notably higher CSF estimates than the other methods. This large CSF overestimation generated by SPM5/SPM8 is possibly derived from misclassification of non-CSF voxels in the subarachnoid space (within the dura, Mueller et al., 1998). Although CSF volume is generally not at the core interest of VBM studies, this overestimation might indirectly influence the results when incorporated for the calculation of a covariate such as ICV (Buckner et al., 2004). Typically this covariate is required when assessing tissue volumes relative to global brain volume and the advised method is to sum up thresholded GM, WM and CSF probability maps. One common workaround is to approximate total intracranial volume by summing GM and WM volume whilst omitting CSF or by using other covariates depending on which aspect of age related GM changes is investigated (Peelle et al., 2012).
The Influence of Different Segmentation Methods on VBM
We showed that method specific differences resulted from the segmentation and were particularly pronounced for the cortical surface close to large sulci and fissures. The largest differences were found between SPM5/SPM8 and DARTEL whereas FSL and IBSNRR segmentation resulted in intermediate sensitivity regarding age-related GM changes. Method specific deviations were not only observed for voxel wise comparisons but also when anatomical regions of interest were used.
SPM5/SPM8, DARTEL, and IBSNRR made use of prior probability maps and our results suggest that segmentation depends critically on how well the anatomical image of the individual subject aligns with this prior information. Warping between the individual anatomy and a template requires non-linear, local deformations particularly around ventricles, sulci and fissures that are typically enlarged in elderly (Blatter et al., 1995; Mu et al., 1999; Raz et al., 2004; Resnick et al., 2007). This presents a challenge for normalization procedures. Comparing the unified segmentation/normalization algorithm provided by SPM5/SPM8 to DARTEL suggests that DARTEL provides a better normalization for deviant anatomy than the unified segmentation/normalization algorithm provided by SPM5/SPM8. SPM5/SPM8 results may be additionally biased because prior tissue probability maps were based on a sample consisting of relatively young, healthy subjects whereas DARTEL probability maps provide a better representation of brain anatomy across the lifespan.
One important difference between the SPM5/SPM8 and the DARTEL procedure is that the latter creates a custom template based on the population under study, thus minimizing spatial transformation across the sample. In the past it has been argued that the unified segmentation/normalization algorithm of SPM5/SPM8 eliminates the necessity for a custom template, even for anatomically deviant populations (Thomas et al., 2009). However, it might be advisable to use a customized template when extreme age-groups are compared, as in our study.
IBSNRR revealed age-related differences that were smaller than the SPM5/SM8 solution but larger than the DARTEL solution. These inconsistencies result most likely from different warping algorithms since IBSNRR uses an integrated segmentation and registration on the basis of a different deformation model, allowing larger local non-linear deformations (more degrees of freedom) than SPM5/SPM8. In agreement with our results, Klein et al. (2009) have argued that normalization solutions tend to be better when methods allow more degrees of freedom of the deformation. Alternatively, advanced segmentation methods can be used to avoid an age–related bias (Ziegler et al., 2012).
Moreover, IBSNRR uses an atlas that was created by a non-linear registration procedure so that the resulting probability maps contained much anatomical detail regarding cortical sulci (see Figure 1 for a comparison with the SPM5/SPM8 GM probability map).
FSL VBM was the only method that revealed a segmentation solution that did not rely on prior tissue probability maps and age-related GM differences were very similar to those revealed by IBSNRR. However, when FSL normalization was qualitatively compared to SPM5/SPM8 and DARTEL normalization, localization of age-related GM differences was slightly different. This is one possible explanation why some VBM results obtained with SPM5/SPM8 (Draganski et al., 2004) could not be reproduced when FSL was used (Scholz et al., 2009), even though these studies tested longitudinal rather than cross sectional GM changes.
In summary, differences in segmentation algorithms (that often rely on an initial spatial normalization step) seem to be a major source for between-method differences when VBM is used to compare the brain structure of young vs. elderly subjects (Courchesne et al., 2000; Bookstein, 2001; Davatzikos, 2004). We suggest that this effect is particularly pronounced when subjects have deviant anatomy (as it is typically the case for elderly individuals) and normalization employs only limited degrees of freedom for local deformations. Therefore, differences in the normalization algorithm might significantly influence the segmentation step when prior tissue probability maps are used, but also the localization of significant GM differences when individuals are compared within a common template space.
Limitations of the Study
In this study we investigated whether different VBM procedures have a significant influence on GM estimates of young vs. elderly subjects. Even though we found substantial differences between methods we cannot infer which VBM method is the most correct one since it is difficult to determine the ground truth for voxel-based methods. Typically this important step for method validation is done by the developers and published when a new method is introduced (see for example Ashburner and Friston, 2011). However, it is rarely the case that previous datasets are reanalyzed what makes comparison across different studies difficult. Here we report that age-specific GM differences are substantial across methods and should be considered when different findings are compared. This is particularly important when early VBM studies are considered which were mostly performed using SPM5/SPM8. Consequently, a-priory hypotheses should be formulated with care considering that SPM5/SPM8 VBM tends to overestimate age-related GM decline in relation to DARTEL, particularly when groups deviate from the anatomical SPM5/SPM8 templates and probability maps.
Another caveat is the fact that our sample was relatively small and drawn from the extremes of the age spectrum. Method specific differences are probably less pronounced when a smaller age range is considered (Peelle et al., 2012). Finally, we applied parametric statistics whilst non-parametric approaches might be generally more appropriate for VBM studies (Douaud et al., 2006). One has to keep in mind, though, that we were not interested in reporting age-specific GM differences per se, but rather whether these results depended on the algorithms provided by software packages frequently used by the neuroimaging community.
Morphometric measurements are increasingly applied for the detection of GM changes in healthy ageing as well as in neurodegenerative disease. Due to its automated and near hypothesis-free character, VBM has gained popularity as a substitute for manual demarcations of GM within volumes of interest. Here, we used the aging brain as a well-known model for structural atrophy and asked whether comparing GM between young and elderly by VBM depends on methodological differences between commonly used software packages, i.e., the unified segmentation/normalization procedure provided by SPM5/SPM8, DARTEL, and FSL. We showed that VBM revealed the largest age-related GM differences when using the SPM5/SPM8 and the smallest when using DARTEL for segmentation. These method-specific differences reached significance when tested at various levels of description (total brain volume, regions of interest, voxel based). We argue that the segmentation procedure can have a major influence on cross-sectional VBM results, particularly when anatomical deviations are more outspoken in one group than the other. Methods that provide only limited degrees of freedom for local deformations (such as SPM5/SPM8) might overestimate between-group differences in VBM results. This is important when results are compared across studies using different VBM methods and particularly when a-priori hypotheses are derived from early VBM studies that were often performed using SPM5/SPM8 segmentation.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We would like to thank Prof. Patrick Dupont for his valuable comments. This work was supported by the Research Council of the KU Leuven (CREA/07/037) and the Flanders Fund for Scientific Research (FWO Project G.0758.10).
The Supplementary Material for this article can be found online at: http://www.frontiersin.org/journal/10.3389/fnagi.2014.00124/abstract
Allen, J. S., Bruss, J., Brown, C. K., and Damasio, H. (2005a). Normal neuroanatomical variation due to age: the major lobes and a parcellation of the temporal region. Neurobiol. Aging 26, 1245–1260. doi: 10.1016/j.neurobiolaging.2005.05.023
Allen, J. S., Bruss, J., Brown, C. K., and Damasio, H. (2005b). Methods for studying the aging brain: volumetric analyses versus VBM (Commentary). Neurobiol. Aging 26, 1275–1278. doi: 10.1016/j.neurobiolaging.2005.05.017
Andersson, J. L. R., Jenkinson, M., and Smith, S. (2007a). “Non-linear optimisation,” in FMRIB Technical Report TR07JA1. Available online at: www.fmrib.ox.ac.uk/analysis/techrep
Andersson, J. L. R., Jenkinson, M., and Smith, S. (2007b). “Non-linear registration, aka Spatial normalisation,” in FMRIB Technical Report TR07JA2. Available online at: www.fmrib.ox.ac.uk/analysis/techrep
Blatter, D. D., Bigler, E. D., Gale, S. D., Johnson, S. C., Anderson, C. V., Burnett, B. M., et al. (1995). Quantitative volumetric analysis of brain MR: normative database spanning 5 decades of life. AJNR Am. J. Neuroradiol. 16, 241–251.
Buckner, R. L., Head, D., Parker, J., Fotenos, A. F., Marcus, D., Morris, J. C., et al. (2004). A unified approach for morphometric and functional data analysis in young, old, and demented adults using automated atlas-based head size normalization: reliability and validation against manual measurement of total intracranial volume. Neuroimage 23, 724–738. doi: 10.1016/j.neuroimage.2004.06.018
Chetelat, G., Landeau, B., Eustache, F., Mezenge, F., Viader, F., de la Sayette, V., et al. (2005). Using voxel-based morphometry to map the structural changes associated with rapid conversion in MCI: a longitudinal MRI study. Neuroimage 27, 934–946. doi: 10.1016/j.neuroimage.2005.05.015
Courchesne, E., Chisum, H. J., Townsend, J., Cowles, A., Covington, J., Egaas, B., et al. (2000). Normal brain development and aging: quantitative analysis at in vivo MR imaging in healthy volunteers. Radiology 216, 672–682. doi: 10.1148/radiology.216.3.r00au37672
D'Agostino, E., Maes, F., Vandermeulen, D., and Suetens, P. (2003). A viscous fluid model for multimodal non-rigid image registration using mutual information. Med. Image. Anal. 7, 565–575. doi: 10.1016/S1361-8415(03)00039-2
D'Agostino, E., Maes, F., Vandermeulen, D., and Suetens, P. (2006a). An information theoretic approach for non-rigid image registration using voxel class probabilities. Med. Image. Anal. 10, 413–431. doi: 10.1016/j.media.2005.03.004
Davies, R. R., Scahill, V. L., Graham, A., Williams, G. B., Graham, K. S., and Hodges, J. R. (2009). Development of an MRI rating scale for multiple brain regions: comparison with volumetrics and with voxel-based morphometry. Neuroradiology 51, 491–503. doi: 10.1007/s00234-009-0521-z
Douaud, G., Gaura, V., Ribeiro, M. J., Lethimonnier, F., Maroy, R., Verny, C., et al. (2006). Distribution of gray matter atrophy in Huntington's disease patients: a combined ROI-based and voxel-based morphometric study. Neuroimage 32, 1562–1575. doi: 10.1016/j.neuroimage.2006.05.057
Ferreira, L. K., Diniz, B. S., Forlenza, O. V., Busatto, G. F., and Zanetti, M. V. (2011). Neurostructural predictors of Alzheimer's disease: a meta-analysis of VBM studies. Neurobiol. Aging 32, 1733–1741. doi: 10.1016/j.neurobiolaging.2009.11.008
Fjell, A. M., Westlye, L. T., Amlien, I., Espeseth, T., Reinvang, I., Raz, N., et al. (2009). High consistency of regional cortical thinning in aging across multiple samples. Cereb. Cortex 19, 2001–2012. doi: 10.1093/cercor/bhn232
Friston, K. J., Holmes, A., Worsley, K. J., Poline, J. B., Frith, C. D., and Frackowiak, R. S. (1994). Statistical parametric maps in functional imaging: a general linear approach. Hum. Brain Mapp. 2, 189–210. doi: 10.1002/hbm.460020402
Giuliani, N. R., Calhoun, V. D., Pearlson, G. D., Francis, A., and Buchanan, R. W. (2005). Voxel-based morphometry versus region of interest: a comparison of two methods for analyzing gray matter differences in schizophrenia. Schizophr. Res. 74, 135–147. doi: 10.1016/j.schres.2004.08.019
Gonoi, W., Abe, O., Yamasue, H., Yamada, H., Masutani, Y., Takao, H., et al. (2010). Age-related changes in regional brain volume evaluated by atlas-based method. Neuroradiology 52, 865–873. doi: 10.1007/s00234-009-0641-5
Good, C. D., Johnsrude, I. S., Ashburner, J., Henson, R. N., Friston, K. J., and Frackowiak, R. S. (2001). A voxel-based morphometric study of ageing in 465 normal adult human brains. Neuroimage 14, 21–36. doi: 10.1006/nimg.2001.0786
Good, C. D., Scahill, R. I., Fox, N. C., Ashburner, J., Friston, K. J., Chan, D., et al. (2002). Automatic differentiation of anatomical patterns in the human brain: validation with studies of degenerative dementias. Neuroimage 17, 29–46. doi: 10.1006/nimg.2002.1202
Greenberg, D. L., Messer, D. F., Payne, M. E., Macfall, J. R., Provenzale, J. M., Steffens, D. C., et al. (2008). Aging, gender, and the elderly adult brain: an examination of analytical strategies. Neurobiol. Aging 29, 290–302. doi: 10.1016/j.neurobiolaging.2006.09.016
Hammers, A., Chen, C. H., Lemieux, L., Allom, R., Vossos, S., Free, S. L., et al. (2007). Statistical neuroanatomy of the human inferior frontal gyrus and probabilistic atlas in a standard stereotaxic space. Hum. Brain Mapp. 28, 34–48. doi: 10.1002/hbm.20254
Hammers, A., Koepp, M. J., Free, S. L., Brett, M., Richardson, M. P., Labbe, C., et al. (2002). Implementation and application of a brain template for multiple volumes of interest. Hum. Brain Mapp. 15, 165–174. doi: 10.1002/hbm.10016
Hasan, K. M., and Pedraza, O. (2009). Improving the reliability of manual and automated methods for hippocampal and amygdala volume measurements. Neuroimage 48, 497–498. doi: 10.1016/j.neuroimage.2009.05.004
Jenkinson, M., Bannister, P.R., Brady, J. M., and Smith, S. M. (2002). Improved optimisation for the robust and accurate linear registration and motion correction of brain images. Neuroimage, 17, 825–841. doi: 10.1006/nimg.2002.1132
Jernigan, T. L., Archibald, S. L., Fennema-Notestine, C., Gamst, A. C., Stout, J. C., Bonner, J., et al. (2001). Effects of age on tissues and regions of the cerebrum and cerebellum. Neurobiol. Aging 22, 581–594. doi: 10.1016/S0197-4580(01)00217-2
Jones, D. K., Chtines, X. A., Job, D., Khong, P. L., Leung, L. T., Marenco, S., et al. (2007). “What happens when nine different groups analyze the same DT-MRI dataset using voxel-based methods?,” in Proceeds of ISMRM 15th Annual Meeting (Berlin).
Kalpouzos, G., Chetelat, G., Baron, J. C., Landeau, B., Mevel, K., Godeau, C., et al. (2009). Voxel-based mapping of brain gray matter volume and glucose metabolism profiles in normal aging. Neurobiol. Aging 30, 112–124. doi: 10.1016/j.neurobiolaging.2007.05.019
Kennedy, K. M., Erickson, K. I., Rodrigue, K. M., Voss, M. W., Colcombe, S. J., Kramer, A. F., et al. (2009). Age-related differences in regional brain volumes: a comparison of optimized voxel-based morphometry to manual volumetry. Neurobiol. Aging 30, 1657–1676. doi: 10.1016/j.neurobiolaging.2007.12.020
Klein, A., Andersson, J., Ardekani, B. A., Ashburner, J., Avants, B., Chiang, M. C., et al. (2009). Evaluation of 14 nonlinear deformation algorithms applied to human brain MRI registration. Neuroimage 46, 786–802. doi: 10.1016/j.neuroimage.2008.12.037
Lehmbeck, J. T., Brassen, S., Weber-Fahr, W., and Braus, D. F. (2006). Combining voxel-based morphometry and diffusion tensor imaging to detect age-related brain changes. Neuroreport 17, 467–470. doi: 10.1097/01.wnr.0000209012.24341.7f
Lemaıtre, H., Crivello, F., Grassiot, B., Alpérovitch, A., Tzourio, C., and Mazoyer, B. (2005). Age- and sex-related effects on the neuroanatomy of healthy elderly. Neuroimage 26, 900–911. doi: 10.1016/j.neuroimage.2005.02.042
Machilsen, B., D'Agostino, E., Maes, F., Vandermeulen, D., Hahn, H. K., Lagae, L., et al. (2007). Linear normalization of MR brain images in pediatric patients with periventricular leukomalacia. Neuroimage 35, 686–697. doi: 10.1016/j.neuroimage.2006.12.037
Maes, F., Collignon, A., Vandermeulen, D., Marchal, G., and Suetens, P. (1997). Multimodality image registration by maximization of mutual information. IEEE Trans. Med. Imaging 16, 187–198. doi: 10.1109/42.563664
Mu, Q., Xie, J., Wen, Z., Weng, Y., and Shuyun, Z. (1999). A quantitative MR study of the hippocampal formation, the amygdala, and the temporal horn of the lateral ventricle in healthy subjects 40 to 90 years of age. AJNR Am. J. Neuroradiol. 20, 207–211.
Mueller, E. A., Moore, M. M., Kerr, D. C., Sexton, G., Camicioli, R. M., Howieson, D. B., et al. (1998). Brain volume preserved in healthy elderly through the eleventh decade. Neurology 51, 1555–1562. doi: 10.1212/WNL.51.6.1555
Mueller, S., Keeser, D., Reiser, M. F., Teipel, S., and Meindl, T. (2012a). Functional and structural mr imaging in neuropsychiatric disorders, part 1: imaging techniques and their application in mild cognitive impairment and Alzheimer Disease. AJNR Am. J. Neuroradiol. 33, 1845–1850. doi: 10.3174/ajnr.A2799
Mueller, S., Keeser, D., Reiser, M. F., Teipel, S., and Meindl, T. (2012b). Functional and structural mr imaging in neuropsychiatric disorders, part 2: application in schizophrenia and autism. AJNR Am. J. Neuroradiol. 33, 2033–2037. doi: 10.3174/ajnr.A2800
Peelle, J. E., Cusack, R., and Henson, R. N. (2012). Adjusting for global effects in voxel-based morphometry: gray matter decline in normal aging. Neuroimage 60, 1503–1516. doi: 10.1016/j.neuroimage.2011.12.086
Raz, N., Gunning-Dixon, F., Head, D., Rodrigue, K. M., Williamson, A., and Acker, J. D. (2004). Aging, sexual dimorphism, and hemispheric asymmetry of the cerebral cortex: replicability of regional differences in volume. Neurobiol. Aging 25, 377–396. doi: 10.1016/S0197-4580(03)00118-0
Raz, N., Rodrigue, K. M., Kennedy, K. M., Head, D., Gunning-Dixon, F., and Acker, J. D. (2003). Differential aging of the human striatum: longitudinal evidence. Am J. Neuroradiol. 24, 1849–1856. doi: 10.1136/jnnp.2007.142786
Resnick, S. M., Lamar, M., and Driscoll, I. (2007). Vulnerability of the orbitofrontal cortex to age-associated structural and functional brain changes. Ann. NY Acad. Sci. 1121, 562–575. doi: 10.1196/annals.1401.027
Rueckert, D., Sonoda, L. I., Hayes, C., Hill, D. L. G., Leach, M. O., and Hawkes, D. J. (1999). Non-rigid registration using free-form deformations: application to breast MR images. IEEE Transac. Med. Imaging, 18, 712–721. doi: 10.1109/42.796284
Skullerud, K. (1985). Variations in the size of the human brain. Influence of age, sex, body length, body mass index, alcoholism, Alzheimer changes, and cerebral atherosclerosis. Acta Neurol. Scand. Suppl. 102, 1–94.
Smith, S. M., Jenkinson, M., Woolrich, M. W., Beckmann, C. F., Behrens, T. E. J., Johansen-Berg, H., et al. (2004). Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 23, 208–219. doi: 10.1016/j.neuroimage.2004.07.051
Stiers, P., Swillen, A., De Smedt, B., Lagae, L., Devriendt, K., D'Agostino, E., et al. (2005). Atypical neuropsychological profile in a boy with 22q11.2 Deletion Syndrome. Child Neuropsychol. 11, 87–108. doi: 10.1080/09297040590911220
Testa, C., Laakso, M. P., Sabattoli, F., Rossi, R., Beltramello, A., Soininen, H., et al. (2004). A comparison between the accuracy of voxel-based morphometry and hippocampal volumetry in Alzheimer's disease. J. Magn. Reson. Imaging 19, 274–282. doi: 10.1002/jmri.20001
Thomas, A. G., Marrett, S., Saad, Z. S., Ruff, D. A., Martin, A., and Bandettini, P.A. (2009). Functional but not structural changes associated with learning: an exploration of longitudinal voxel-based morphometry (VBM). Neuroimage 48, 117–125. doi: 10.1016/j.neuroimage.2009.05.097
Thomas, C., and Baker, C. I. (2013). Teaching an adult brain new tricks: A critical review of evidence for training-dependent structural plasticity in humans. Neuroimage 73, 225–236. doi: 10.1016/j.neuroimage.2012.03.069
Tisserand, D. J., Pruessner, J. C., Sanz Arigita, E. J., van Boxtel, M. P., Evans, A. C., Jolles, J., et al. (2002). Regional frontal cortical volumes decrease differentially in aging: an MRI study to compare volumetric approaches and voxel-based morphometry. Neuroimage 17, 657–669. doi: 10.1006/nimg.2002.1173
Tisserand, D. J., van Boxtel, M. P., Pruessner, J. C., Hofman, P., Evans, A. C., and Jolles, J. (2004). A voxel-based morphometric study to determine individual differences in gray matter density associated with age and cognitive change over time. Cereb. Cortex 14, 966–973. doi: 10.1093/cercor/bhh057
Van Laere, K. J., and Dierckx, R. A. (2001). Brain perfusion SPECT: age- and sex-related effects correlated with voxel-based morphometric findings in healthy adults. Radiology 221, 810–817. doi: 10.1148/radiol.2213010295
Van Leemput, K., Maes, F., Vandermeulen, D., and Suetens, P. (1999). Automated model-based bias field correction of MR images of the brain. IEEE Trans. Med. Imaging 18, 885–896. doi: 10.1109/42.811268
Walhovd, K. B., Fjell, A. M., Reinvang, I., Lundervold, A., Dale, A. M., Eilertsen, D. E., et al. (2005). Effects of age on volumes of cortex, white matter and subcortical structures. Neurobiol. Aging 26, 1261–1270. doi: 10.1016/j.neurobiolaging.2005.05.020
Walhovd, K. B., Westlye, L. T., Amlien, I., Espeseth, T., Reinvang, I., Raz, N., et al. (2010). Consistent neuroanatomical age-related volume differences across multiple samples. Neurobiol. Aging 32, 916–932. doi: 10.1016/j.neurobiolaging.2009.05.013
Wang, Q., Seghers, D., D'Agostino, E., Maes, F., Vandermeulen, D., Suetens, P., et al. (2005). Construction and validation of mean shape atlas templates for atlas-based brain image segmentation. Inf. Process Med. Imaging 19, 689–700. doi: 10.1007/11505730_57
Zhang, Y., Brady, M., and Smith, S. (2001). Segmentation of brain MR images through a hidden Markov random field model and the expectation maximization algorithm. IEEE Trans. Med. Imaging, 20, 45–57. doi: 10.1109/42.906424
Keywords: ageing, gray matter, atrophy, brain mapping, volumetry, MRI
Citation: Callaert DV, Ribbens A, Maes F, Swinnen SP and Wenderoth N (2014) Assessing age-related gray matter decline with voxel-based morphometry depends significantly on segmentation and normalization procedures. Front. Aging Neurosci. 6:124. doi: 10.3389/fnagi.2014.00124
Received: 10 February 2014; Accepted: 27 May 2014;
Published online: 23 June 2014.
Edited by:Rodrigo Orlando Kuljiš, Zdrav Mozak Limitada, Chile
Reviewed by:Rongqiao He, Chinese Academy of Sciences, China
Olof Erik Lindberg, Karolinska institute, Sweden
Copyright © 2014 Callaert, Ribbens, Maes, Swinnen and Wenderoth. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Nicole Wenderoth, Neural Control of Movement Laboratory, Health Sciences and Technology, ETH Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland e-mail: email@example.com