Neuroimaging Biomarkers in Alzheimer's Disease
Neuroimaging Biomarkers in Alzheimer's Disease
Neuroimaging Biomarkers in Alzheimer's Disease
in Biomedical Research
Neuroimaging
biomarkers
inAlzheimer's
disease.
Dr. Samuel Barrack
Published by:
Internet Medical Publishing
Find more titles on our site: www.imedpub.com
Thematic Collections
in Biomedical Research
Neuroimaging
biomarkers
inAlzheimer's
disease.
Dr. Samuel Barrack
Disclaimer: This collection is based in articles published by the Public Library of Science (PLOS). PLOS applies the
Creative Commons Attribution License (CC-BY) to works we publish. Under this license, authors retain ownership of
the copyright for their content, but allow anyone to download, reuse, reprint, modify, distribute, and/or copy the
content as long as the original authors and source are cited.
II
III
IV
Contents
Neuroimaging
biomarkers: MRI
Diffusion Tensor Metrics as Biomarkers
in Alzheimers Disease
Neuroimaging
biomarkers: PET
2
Cost-Effectiveness of Magnetic
Resonance Imaging with a New
Contrast Agent for the Early Diagnosis
of Alzheimers Disease
16
Multi-Method Analysis of MRI Images
in Early Diagnostics of Alzheimers
Disease
30
CSF and Brain Structural Imaging
Markers of the Alzheimers
Pathological
Cascade
Sensitivity and Specificity of Medial
Temporal Lobe Visual Ratings
and Multivariate Regional MRI
Classification in Alzheimers Disease
MRI Markers for Mild Cognitive
Impairment: Comparisons between
White Matter Integrity and Gray
Matter Volume Measurements
65
72
46
89
97
118
Abstract
Background: Although diffusion tensor imaging has been a major research focus for Alzheimers disease in recent years, it
remains unclear whether it has sufficient stability to have biomarker potential. To date, frequently inconsistent results have
been reported, though lack of standardisation in acquisition and analysis make such discrepancies difficult to interpret.
There is also, at present, little knowledge of how the biometric properties of diffusion tensor imaging might evolve in the
course of Alzheimers disease.
Methods: The biomarker question was addressed in this study by adopting a standardised protocol both for the whole
brain (tract-based spatial statistics), and for a region of interest: the midline corpus callosum. In order to study the evolution
of tensor changes, cross-sectional data from very mild (N = 21) and mild (N = 22) Alzheimers disease patients were examined
as well as a longitudinal cohort (N = 16) that had been rescanned at 12 months.
Findings and Significance: The results revealed that increased axial and mean diffusivity are the first abnormalities to occur
and that the first region to develop such significant differences was mesial parietal/splenial white matter; these metrics,
however, remained relatively static with advancing disease indicating they are suitable as state-specific markers. In
contrast, increased radial diffusivity, and therefore decreased fractional anisotropythough less detectable earlybecame
increasingly abnormal with disease progression, and, in the splenium of the corpus callosum, correlated significantly with
dementia severity; these metrics therefore appear stage-specific and would be ideal for monitoring disease progression. In
addition, the cross-sectional and longitudinal analyses showed that the progressive abnormalities in radial diffusivity and
fractional anisotropy always occurred in areas that had first shown an increase in axial and mean diffusivity. Given that the
former two metrics correlate with dementia severity, but the latter two did not, it would appear that increased axial
diffusivity represents an upstream event that precedes neuronal loss.
Citation: Acosta-Cabronero J, Alley S, Williams GB, Pengas G, Nestor PJ (2012) Diffusion Tensor Metrics as Biomarkers in Alzheimers Disease. PLoS ONE 7(11):
e49072. doi:10.1371/journal.pone.0049072
Editor: Wang Zhan, University of Maryland, College Park, United States of America
Received June 21, 2012; Accepted October 4, 2012; Published November 7, 2012
Copyright: 2012 Acosta-Cabronero et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which
permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: Medical Research Council, UK (to PJN); National Institute for Health Research, Cambridge Biomedical Research Centre, UK (to PJN); and Alzheimers
Research UK (to JAC and GP). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing Interests: The authors have declared that no competing interests exist.
* E-mail: [email protected]
Introduction
Table 1. Demographic summary including cognitive features for Alzheimers disease patients and for a group of elderly controls.
General Demographics
Global Cognition
ACE-R Subscores
Control
(N = 26)
Alzheimers
Disease (N = 43)
Mild Alzheimers
disease (N = 22)
Gender, M:F
11:15
26:17
13:8
13:9
68 (6)
70 (6)
72 (5)
69 (6)
MMSE/30
29.1 (0.8)
23.7 (3.6)**
25.9 (1.6)**
21.7 (3.9)**
ACE-R/100
94.6 (3.0)
71.5 (11.9)**
81.4 (4.0)**
62.1 (8.8)**
17.9 (0.3)
15.4 (2.7)**
17.0 (1.4)*
13.9 (2.8)**
Memory/26
24.5 (1.9)
10.9 (4.3)**
13.7 (4.1)**
8.2 (2.4)**
Fluency/14
12.1 (1.7)
8.3 (3.2)**
10.6 (1.9)*
6.2 (2.6)**
Language/26
24.9 (0.9)
23.0 (2.6)**
24.6 (1.2)
21.5 (2.7)**
Visuospatial/16
15.4 (0.9)
13.7 (2.8)**
15.1 (1.2)
12.3 (3.1)**
Disease severity (as measured by ACE-R) enabled a median split of the patient cohort into very mild and mild Alzheimers disease subgroups.
Where appropriate, group values are given as mean (SD).
MMSE/30 = Mini-mental state examination score out of 30-point total; ACE-R/100 = Addenbrookes cognitive examination-revised score out of 100-point total.
Wilcoxon rank-sum significance levels: *P,0.01 (Alzheimers disease worse than controls, two-tailed); **P,0.05 (Alzheimers disease worse than controls, two-tailed,
Bonferroni-corrected on n = 28 tests); P,0.05 (Mild worse than very mild Alzheimers disease, two-tailed, Bonferroni-corrected on n = 28).
doi:10.1371/journal.pone.0049072.t001
Methods
Alzheimers Disease
(Longitudinal, N = 16)
Ethics Statement
General
Demographics
12 Months
8:8
25.1 (2.0)
22.6 (4.9)1
77.5 (7.3)
70.3 (13.6)**
14.8 (3.5)
Memory/26
13.5 (3.9)
Fluency/14
9.6 (2.5)
9.8 (3.4)**
8.5 (3.8)
Language/26
23.6 (1.8)
23.1 (2.5)
14.0 (2.9)
Subjects
Cross-sectional cohorts. Forty-three patients with earlystage probable Alzheimers disease according to criteria from the
National Institute of Neurological and Communicative Disorders
and Stroke and the Alzheimers Disease and Related Disorders
Association (NINCDS-ADRDA) [17] were recruited from the
memory clinic at Addenbrookes Hospital (Cambridge, UK). For
cross-sectional comparisons, 26 matched controls were also
Gender, M:F
Baseline
Imaging
Figure 2. Corpus callosum subdivision. Depiction of the semiautomated callosal subdivision into splenium, truncus and genu (top),
and their intersection with the mean FA skeleton inferred from N = 69
subjectsN = 43 Alzheimers disease patients and N = 26 matched
controls (bottom).
doi:10.1371/journal.pone.0049072.g002
Figure 3. Cross-sectional study of very mild Alzheimers disease. TBSS results for the very mild Alzheimers disease group compared to
controls. Statistical maps (thresholded at TFCE-P,0.05) for increased axial/radial diffusivity and reduced FA overlaid onto the mean FA skeleton and
the MNI152 template. Coronal depths are given in millimetres.
doi:10.1371/journal.pone.0049072.g003
compute the diagonal elements (l1, l2 and l3) at each brain voxel,
from which the derived metrics RD, MD and FA were also
inferred. Note that negative primary eigenvalues were deemed
unphysical and were set to 0a visual inspection of the spatial
distribution of negative eigenvalues revealed that they were located
in the periphery of white matter bundles i.e. adjacent to other
tissue types and far from tract centres.
Tract-based spatial statistics (TBSS) analysis. The TBSS
approach [25] was used to perform whole-brain statistical analyses
at white matter tract centres. Spatial normalisation was achieved
by warping all FA images to the 16161 mm3 FMRIB58_FA
standard template (FMRIB, University of Oxford, UK) in
MNI152 space (Montreal Neurological Institute, McGill University, Canada) using the FMRIBs non-linear image registration tool
Figure 4. Cross-sectional study of mild Alzheimers disease. TBSS results for the mild-stage Alzheimers disease group compared to controls.
Thresholded (TFCE-P,0.05) statistical maps for increased axial/radial diffusivity and reduced FA were overlaid onto the mean FA skeleton and the
MNI152 template. Coronal depths are given in millimetres.
doi:10.1371/journal.pone.0049072.g004
Table 3. Alzheimers disease skeletonised DTI parametric comparisons in different mid-sagittal callosal areas.
Splenium
Truncus
Genu
l1
22.24
21.39
20.22
RD
23.08**
20.76
0.53
MD
23.39**
21.08
0.30
FA
2.451
0.40
21.18
Results
Cross-sectional TBSS Study of Very Mild Alzheimers
Disease
TBSS results for the voxel-wise group contrast of very mild
Alzheimers disease versus controls are shown in Figure 3. The l1
statistical map showed significant bilateral and confluent change
that predominantly involved parietal white matter, with strongest
abnormalities along the posterior cingulum and the inter-hemispheric tracts of the caudal corpus callosum. In contrast, there was
relative sparing of caudal occipital, rostral temporal lobe and
prefrontal white matter, and, the more rostral tracts of the corpus
callosum. RD abnormalities spatially overlapped with those found
for l1 in parietal and superior temporal white matterthough only
Table 4. DTI group comparisons across disease stages and linear dependence on global cognition for all patients in the splenium.
Splenium
ACE-R Linear
Dependence (N = 43)
Cross-sectional
Very mild Alzheimers (N = 21) vs.
Control (N = 26)
l1
2.62*
1.19
1.54
0.15
RD
2.131
3.07**
20.89
20.40**
MD
3.01**
2.74*
0.21
20.21
FA
21.27
22.85*
1.4
0.44**
Group results are given as Wilcoxon rank-sum Z-statistic; statistical dependencies are given as Pearson correlation coefficient with 41 degrees of freedom.
Significance levels: 10.01,P,0.05 (Alzheimers disease worse than controls, two-tailed); *P,0.01 (Alzheimers disease worse than controls, two-tailed); **P,0.05
(Alzheimers disease worse than controls, two-tailed, Bonferroni-corrected on n = 24; or DTI measure correlated with Alzheimers disease cognitive status, two-tailed,
Bonferroni-corrected on n = 12).
doi:10.1371/journal.pone.0049072.t004
Figure 6. Cross-sectional diffusion tensor behaviour in the splenial region. Mean subject values for skeletonised DTI parameters in the
splenium as a function of cognitive status (ACE-R scores) for controls (green), very mild Alzheimers disease (blue) and mild Alzheimers disease
patients (red). The error bars represent 6 one group standard deviation. The vertical axes were scaled to 10 control standard deviations. The vertical
lines delimit the control exclusion criteria (ACE-R,88/100) and the median split (ACE-R = 74). A least-square linear fit was displayed if Pearsons
correlation coefficient was deemed statistically significant (Table 4).
doi:10.1371/journal.pone.0049072.g006
Regions of Interest
The regional analyses were consistent with the above TBSS
observations. Table 3 highlights the preservation of the midline
truncus and genu relative to the splenium when contrasting all
Alzheimers patients (N = 43) with the control group. Examining
8
Figure 7. Longitudinal study of Alzheimers disease. Whole-brain TBSS contrast on 12-month follow-up versus baseline in the longitudinal
Alzheimers disease cohort (TFCE-P,0.05) for increased radial diffusivity and reduced FA n.b. no significant results for l1 or MD were found.
doi:10.1371/journal.pone.0049072.g007
10
Discussion
This study aimed to investigate the biomarker potential of DTI
metrics in Alzheimers disease through whole-brain analyses as
well as by taking a reductionist approach in the corpus callosum to
eliminate ambiguities generated by uncertain tract visualisation
and crossing fibres. The evolution of DTI changes was examined
in both cross-sectional and longitudinal datasets; the former by
way of a median spilt in which 76% of the very mild patients were
in the mild-cognitive impairment stage of Alzheimers disease, and
the latter by contrasting 12-month follow-up scans to baseline.
Across the various analyses, a consistent picture emerged in which
an increase in l1 is the first significant change in Alzheimers
disease, but then remains relatively steady. In contrast, RDand
therefore FAbecome progressively more abnormal with disease
evolution, making them the candidate metrics for stage-specific
biomarkers. In other words, these markers could have value to
track change over time. Given that FA reduction is a function of
both RD and l1, and the early increase in l1 appeared to
attenuate with disease progression, it may be that FA is the
superior measure of decline compared to RD in longitudinal
studies. It should be emphasized, however, that the attenuation of
l1 was subtle and, unless replicated in similar studies, could simply
represent random variability in the current dataset. A further
caveat to this proposal is that the preference for FA would only
apply to areas in which the l1 abnormality has already peaked
such as the splenium in the present analyses. For instance, in
a repeated-measure longitudinal design, in areas that are initially
spared in terms of l1, the first movement of this metric would be to
increase, leading to a loss of sensitivity for FA to detect change.
This was exemplified by the cross-sectional results, in which frontal
and left temporal white matter was normal in the very mild group,
but became abnormal (mostly in terms of l1) in the mild group.
This problem would, of course, be avoided by studying RD,
whichunlike FAis independent of l1.
In contrast to RD and FA, the results indicated that l1 has no
role at all as staging biomarker for a given region of interest.
Although one could, in theory, propose it as a staging marker by
looking at the differential spatial extent of l1 abnormalities over
time, this would be difficult to implement as a marker of
progression given that some degree of heterogeneity would be
expected across subjects; furthermore, such an approach would
also be confounded if it proved correct in future studies that the l1
Longitudinal
Splenium
Baseline (N = 16) vs.
Control (N = 26)
l1
2.371
1.54
20.52
RD
1.72
2.97**
22.59*
2.71*
21.76
22.58*
22.64*
MD
2.47
FA
20.87
Baseline and 12 months versus controls comparisons are reported as Wilcoxon rank-sum Z-statistic values; paired longitudinal results are given as Wilcoxon signed-rank
Z-statistic.
Significance levels: 10.01,P,0.05 (Alzheimers disease worse than controls, two-tailed); *P,0.01 (Alzheimers disease worse than controls; or 12 months worse than
baseline, two-tailed); **P,0.05 (Alzheimers disease worse than controls, two-tailed, Bonferroni-corrected on n = 12).
doi:10.1371/journal.pone.0049072.t005
10
11
Figure 9. Longitudinal tensor behaviour in the splenium. Longitudinal pairs of mean subject skeletonised DTI parameters as a function of
cognitive status (ACE-R) for Alzheimers disease subjects at baseline (blue) and 12 months (red).
doi:10.1371/journal.pone.0049072.g009
12
11
Supporting Information
Figure S1 Increased mean diffusivity results for very
mild and mild Alzheimers disease. TBSS results for very
mild- and mild-stage Alzheimers disease groups compared to
controls. Thresholded (TFCE-P,0.05) statistical maps for increased mean diffusivity were overlaid onto the mean FA skeleton
and the MNI152 template with coronal depths given in
millimetres. Extensive, mostly bilateral distributions of DTI
abnormalities for increased mean diffusion were found in the
very mild and mild Alzheimers disease group comparisons.
Significant abnormalities were located in parietal white matter
regions including the caudal corpus callosum and the posterior
cingulum bundle, and in caudal temporal areas. All clusters of
significance found in the very mild Alzheimers disease group,
were also found in the mild group. As expected, overall MD
abnormalities were highly concordant and largely overlapped with
the spatial distribution of l1 clusters of significance shown in
Figure 3 and Figure 4 (main manuscript).
(TIF)
Figure S2 Regional analysis of the splenium in native
space. Mean subject values for DTI parameters in the native
splenial region as a function of cognitive status (ACE-R) for
12
13
(TIF)
Acknowledgments
We are grateful for the support from our patients and healthy volunteers.
Author Contributions
Conceived and designed the experiments: JAC GBW PJN. Performed the
experiments: JAC GP. Analyzed the data: JAC SA. Wrote the paper: JAC
GBW GP PJN.
References
14
1. Jack CR, Jr., Petersen RC, Xu YC, OBrien PC, Smith GE, et al. (1999)
Prediction of AD with MRI-based hippocampal volume in mild cognitive
impairment. Neurology 52: 13971403.
2. Fellgiebel A, Schermuly I, Gerhard A, Keller I, Albrecht J, et al. (2008)
Functional relevant loss of long association fibre tracts integrity in early
Alzheimers disease. Neuropsychologia 46: 16981706.
3. Mielke MM, Kozauer NA, Chan KC, George M, Toroney J, et al. (2009)
Regionally-specific diffusion tensor imaging in mild cognitive impairment and
Alzheimers disease. Neuroimage 46: 4755.
4. Takahashi S, Yonezawa H, Takahashi J, Kudo M, Inoue T, et al. (2002)
Selective reduction of diffusion anisotropy in white matter of Alzheimer disease
brains measured by 3.0 Tesla magnetic resonance imaging. Neurosci Lett 332:
4548.
5. Acosta-Cabronero J, Williams GB, Pengas G, Nestor PJ (2010) Absolute
diffusivities define the landscape of white matter degeneration in Alzheimers
disease. Brain 133: 529539.
6. Fellgiebel A, Wille P, Muller MJ, Winterer G, Scheurich A, et al. (2004)
Ultrastructural hippocampal and white matter alterations in mild cognitive
impairment: a diffusion tensor imaging study. Dement Geriatr Cogn Disord 18:
101108.
7. Huang H, Fan X, Weiner M, Martin-Cook K, Xiao G, et al. (2012) Distinctive
disruption patterns of white matter tracts in Alzheimers disease with full
diffusion tensor characterization. Neurobiol Aging 33: 20292045.
8. Shu N, Wang Z, Qi Z, Li K, He Y (2011) Multiple diffusion indices reveals white
matter degeneration in Alzheimers disease and mild cognitive impairment:
a tract-based spatial statistics study. J Alzheimers Dis 26 Suppl 3: 275285.
9. Zhang Y, Schuff N, Du AT, Rosen HJ, Kramer JH, et al. (2009) White matter
damage in frontotemporal dementia and Alzheimers disease measured by
diffusion MRI. Brain 132: 25792592.
10. Huang J, Auchus AP (2007) Diffusion tensor imaging of normal appearing white
matter and its correlation with cognitive functioning in mild cognitive
impairment and Alzheimers disease. Ann N Y Acad Sci 1097: 259264.
11. Bosch B, Arenaza-Urquijo EM, Rami L, Sala-Llonch R, Junque C, et al. (2012)
Multiple DTI index analysis in normal aging, amnestic MCI and AD.
Relationship with neuropsychological performance. Neurobiol Aging 33: 6174.
12. Zhang Y, Schuff N, Jahng GH, Bayne W, Mori S, et al. (2007) Diffusion tensor
imaging of cingulum fibers in mild cognitive impairment and Alzheimer disease.
Neurology 68: 1319.
13. Douaud G, Jbabdi S, Behrens TE, Menke RA, Gass A, et al. (2011) DTI
measures in crossing-fibre areas: increased diffusion anisotropy reveals early
white matter alteration in MCI and mild Alzheimers disease. Neuroimage 55:
880890.
14. Fellgiebel A, Muller MJ, Wille P, Dellani PR, Scheurich A, et al. (2005) Colorcoded diffusion-tensor-imaging of posterior cingulate fiber tracts in mild
cognitive impairment. Neurobiol Aging 26: 11931198.
15. Zhang J, Jones M, DeBoy CA, Reich DS, Farrell JA, et al. (2009) Diffusion
tensor magnetic resonance imaging of Wallerian degeneration in rat spinal cord
after dorsal root axotomy. J Neurosci 29: 31603171.
16. Hanyu H, Asano T, Sakurai H, Imon Y, Iwamoto T, et al. (1999) Diffusionweighted and magnetization transfer imaging of the corpus callosum in
Alzheimers disease. J Neurol Sci 167: 3744.
17. McKhann G, Drachman D, Folstein M, Katzman R, Price D, et al. (1984)
Clinical diagnosis of Alzheimers disease: report of the NINCDS-ADRDA Work
Group under the auspices of Department of Health and Human Services Task
Force on Alzheimers Disease. Neurology 34: 939944.
18. Folstein MF, Folstein SE, McHugh PR (1975) "Mini-mental state". A practical
method for grading the cognitive state of patients for the clinician. J Psychiatr
Res 12: 189198.
19. Mioshi E, Dawson K, Mitchell J, Arnold R, Hodges JR (2006) The
Addenbrookes Cognitive Examination Revised (ACE-R): a brief cognitive test
battery for dementia screening. Int J Geriatr Psychiatry 21: 10781085.
20. Reese TG, Heid O, Weisskoff RM, Wedeen VJ (2003) Reduction of eddycurrent-induced distortion in diffusion MRI using a twice-refocused spin echo.
Magn Reson Med 49: 177182.
21. Griswold MA, Jakob PM, Heidemann RM, Nittka M, Jellus V, et al. (2002)
Generalized autocalibrating partially parallel acquisitions (GRAPPA). Magn
Reson Med 47: 12021210.
22. Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TE, et al.
(2004) Advances in functional and structural MR image analysis and
implementation as FSL. Neuroimage 23 Suppl 1: S208219.
23. Jenkinson M, Smith S (2001) A global optimisation method for robust affine
registration of brain images. Med Image Anal 5: 143156.
24. Smith SM (2002) Fast robust automated brain extraction. Hum Brain Mapp 17:
143155.
25. Smith SM, Jenkinson M, Johansen-Berg H, Rueckert D, Nichols TE, et al.
(2006) Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion
data. Neuroimage 31: 14871505.
26. Smith SM, Nichols TE (2009) Threshold-free cluster enhancement: addressing
problems of smoothing, threshold dependence and localisation in cluster
inference. Neuroimage 44: 8398.
27. Hofer S, Frahm J (2006) Topography of the human corpus callosum revisited
comprehensive fiber tractography using diffusion tensor magnetic resonance
imaging. Neuroimage 32: 989994.
28. Witelson SF (1989) Hand and sex differences in the isthmus and genu of the
human corpus callosum. A postmortem morphological study. Brain 112: 799
835.
29. Snook L, Plewes C, Beaulieu C (2007) Voxel based versus region of interest
analysis in diffusion tensor imaging of neurodevelopment. Neuroimage 34: 243
252.
30. Lilliefors HW (1967) On the Kolmogorov-Smirnov test for normality with mean
and variance unknown. Journal of the American Statistical Association 62: 399
402.
31. Mann HB, Whitney DR (1947) On a Test of Whether one of Two Random
Variables is Stochastically Larger than the Other. Annals of Mathematical
Statistics 18: 5060.
32. Wilcoxon F (1946) Individual comparisons of grouped data by ranking methods.
J Econ Entomol 39: 269.
33. Pearson K (1896) Mathematical Contributions to the Theory of Evolution. III.
Regression, Heredity and Panmixia. Philosophical Transactions of the Royal
Society A 187: 253318.
34. Jack CR, Jr., Twomey CK, Zinsmeister AR, Sharbrough FW, Petersen RC, et
al. (1989) Anterior temporal lobes and hippocampal formations: normative
volumetric measurements from MR images in young adults. Radiology 172:
549554.
35. Pengas G, Pereira JM, Williams GB, Nestor PJ (2009) Comparative reliability of
total intracranial volume estimation methods and the influence of atrophy in
a longitudinal semantic dementia cohort. J Neuroimaging 19: 3746.
36. ODwyer L, Lamberton F, Bokde AL, Ewers M, Faluyi YO, et al. (2011)
Multiple indices of diffusion identifies white matter damage in mild cognitive
impairment and Alzheimers disease. PLoS One 6: e21745.
37. Oishi K, Akhter K, Mielke M, Ceritoglu C, Zhang J, et al. (2011) Multi-modal
MRI analysis with disease-specific spatial filtering: initial testing to predict mild
cognitive impairment patients who convert to Alzheimers disease. Front Neurol
2: 54.
38. Pengas G, Williams GB, Acosta-Cabronero J, Ash TWJ, Hong YT, et al. (2012)
The relationship of topographical memory performance to regional neurodegeneration in Alzheimers disease. Front Ag Neurosci 4: 17.
39. Chase TN, Foster NL, Fedio P, Brooks R, Mansi L, et al. (1984) Regional
cortical dysfunction in Alzheimers disease as determined by positron emission
tomography. Ann Neurol 15 Suppl: S170174.
40. Nestor PJ, Fryer TD, Smielewski P, Hodges JR (2003) Limbic hypometabolism
in Alzheimers disease and mild cognitive impairment. Ann Neurol 54: 343351.
41. Wyss-Coray T, Rogers J (2012) Inflammation in Alzheimer disease-a brief
review of the basic science and clinical literature. Cold Spring Harb Perspect
Med 2: a006346.
42. Di Paola M, Spalletta G, Caltagirone C (2010) In vivo structural neuroanatomy
of corpus callosum in Alzheimers disease and mild cognitive impairment using
different MRI techniques: a review. J Alzheimers Dis 20: 6795.
13
43. Correia MM, Carpenter TA, Williams GB (2009) Looking for the optimal DTI
acquisition scheme given a maximum scan time: are more b-values a waste of
time? Magn Reson Imaging 27: 163175.
44. Danielian LE, Iwata NK, Thomasson DM, Floeter MK (2010) Reliability of
fiber tracking measurements in diffusion tensor imaging for longitudinal study.
Neuroimage 49: 15721580.
15
PLOS ONE | www.plosone.org
14
Abstract
Background: Used as contrast agents for brain magnetic resonance imaging (MRI), markers for beta-amyloid deposits might
allow early diagnosis of Alzheimers disease (AD). We evaluated the cost-effectiveness of such a diagnostic test, MRI+CLP
(contrastophore-linker-pharmacophore), should it become clinically available.
Methodology/Principal Findings: We compared the cost-effectiveness of MRI+CLP to that of standard diagnosis using
currently available cognition tests and of standard MRI, and investigated the impact of a hypothetical treatment efficient in
early AD. The primary analysis was based on the current French context for 70-year-old patients with Mild Cognitive
Impairment (MCI). In alternative screen and treat scenarios, we analyzed the consequences of systematic screenings of
over-60 individuals (either population-wide or restricted to the ApoE4 genotype population). We used a Markov model of
AD progression; model parameters, as well as incurred costs and quality-of-life weights in France were taken from the
literature. We performed univariate and probabilistic multivariate sensitivity analyses. The base-case preferred strategy was
the standard MRI diagnosis strategy. In the primary analysis however, MRI+CLP could become the preferred strategy under a
wide array of scenarios involving lower cost and/or higher sensitivity or specificity. By contrast, in the screen and treat
analyses, the probability of MRI+CLP becoming the preferred strategy remained lower than 5%.
Conclusions/Significance: It is thought that anti-beta-amyloid compounds might halt the development of dementia in
early stage patients. This study suggests that, even should such treatments become available, systematically screening the
over-60 population for AD would only become cost-effective with highly specific tests able to diagnose early stages of the
disease. However, offering a new diagnostic test based on beta-amyloid markers to elderly patients with MCI might prove
cost-effective.
Citation: Biasutti M, Dufour N, Ferroud C, Dab W, Temime L (2012) Cost-Effectiveness of Magnetic Resonance Imaging with a New Contrast Agent for the Early
Diagnosis of Alzheimers Disease. PLoS ONE 7(4): e35559. doi:10.1371/journal.pone.0035559
Editor: David W. Dowdy, Johns Hopkins Bloomberg School of Public Health, United States of America
Received July 22, 2011; Accepted March 20, 2012; Published April 20, 2012
Copyright: 2012 Biasutti et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This research was financed by a grant from the Cnam, a public higher education establishment by which all authors were employed at the time. The
funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing Interests: The authors have declared that no competing interests exist.
* E-mail: [email protected]
Introduction
Alzheimers disease (AD) is the main cause of dementia in older
people, with approximately 26 million cases worldwide [1,2,3].
What is more, a major increase in this prevalence is expected in
years to come [3]. In France, while 850.000 people were
diagnosed with AD in 2004, 2.1 million may be affected by
2040 [4].
There are currently no treatments that may cure AD or halt the
course of the disease. However, in recent years, drugs such as
acetylcholinesterase inhibitors have demonstrated efficacy at
reducing the intensity of certain symptoms. Moreover, new
avenues for research are being investigated. Many scientists
believe that one of the main causes of the AD has to do with betaamyloid, microscopic protein which accumulate throughout the
cortex of Alzheimer patients. This is called the amyloid
hypothesis [2]. They believe that the destruction of brain cells
seen in AD is caused by defects in the way beta-amyloid is
16
Methods
1. Framework of the cost-effectiveness analyses
We performed cost-effectiveness analyses in order to compare
the costs and benefits of several alternative diagnostic strategies in
the French context. Costs were measured in Euros (J) and benefits
were measured in quality-adjusted life-years (QALYs), which
assign to each year of life a weight between 0 (dead) and 1 (perfect
health).
A strategy was considered dominated if another strategy had a
better or similar efficacy at a lower cost; conversely, a strategy was
considered strongly dominant when it was both more effective and
cheaper than all other strategies. We computed incremental costeffectiveness ratios (ICERs), in which changes in resource use,
compared with the next best strategy, were included in the
numerator, while additional health effects, compared with the next
best strategy, were included in the denominator. Finally, we
compared ICERs to the willingness-to-pay (WTP) for an
additional QALY, which was assumed equal to three times the
gross national product per capita in France, as recommended by
the World Health Organization Choice working group [13,14],
that is, 76.171J per QALY in 2009. If no strategy was strongly
dominant, the preferred strategy was that with the highest ICER
under the willingness-to-pay threshold.
The study was conducted from a societal perspective, meaning it
included all costs and benefits, no matter who incurred them.
Future costs and QALYs were discounted at 5% annually.
The cost-effectiveness analyses were conducted using the
TreeAge software [15]. We performed expected value analyses,
based on the computation over simulated cycles of percentages of
a hypothetical cohort in each modeled AD stage. For the
multivariate sensitivity analyses, we also performed Monte Carlo
simulations with 10.000 trials, in order to derive the distribution of
incremental cost-effectiveness ratios for the MRI+CLP strategies,
as well as acceptability curves for all strategies.
Table 1 details the compared strategies.
17
Strategy
Tests
Imaging
Standard diagnosis
None
Standard MRI
Non-enhanced MRI
MRI+CLP
doi:10.1371/journal.pone.0035559.t001
18
Figure 1. Decision tree for an individual tested for Alzheimers disease (AD). Possible outcomes of the testing procedure are depicted as a
function of the individuals health status for: (a) The primary scenario (testing of over-70 patients consulting for dementia symptoms). (b) The screen
and treat scenario (systematic screening of the over-60 population). Depending on the investigated strategy, the generic diagnostic test
mentioned in the trees may be standard diagnosis, standard MRI or MRI+CLP. When AD is diagnosed, the imaging procedures are followed by a
cognition test (MMSE) in order to determine the disease stage. No test is performed in severe AD patients, who are assumed to be diagnosed directly.
doi:10.1371/journal.pone.0035559.g001
3. Available treatments
The model assumes that all patients who receive a diagnosis of
probable Alzheimer disease receive treatment with donepezil,
memantine or a hypothetical higher-efficacy drug.
If the patient is diagnosed at a mild or moderate stage,
donepezil is prescribed, to be replaced by memantine when the
patient has evolved to a severe stage of AD. If the patient is
diagnosed when the disease is already severe, memantine is
prescribed.
In addition, a hypothetical treatment T which is efficient at
early stages of AD may be prescribed to patients who are
diagnosed with AD but have high MMSE scores in the primary
analysis, or to patients who are diagnosed with early AD in the
screen and treat analyses.
Donepezil treatment has been shown to improve significantly
memory and other cognitive functions in patients with mild to
moderate AD [21], and to reduce the annual decline in cognition
in these patients when compared with patients in a placebo group
[22]; similarly, memantine treatment has been shown to cause a
clinically noticeable reduction in deterioration over 28 weeks,
compared with placebo, in patients with moderate-to-severe AD
[23,24]. It is to be noted that actually, donepezil or memantine
help treat the symptoms of AD although there is no evidence that
they modify the underlying pathology of the disease. On the basis
of these data, we assumed that treatment modified transition
probabilities between disease stages, both reducing the speed of
AD progression and increasing the chance of symptom lessening
(Table S1). For treatment T, we only assumed a decrease in the
progression from early to mild and moderate stage AD, as we
supposed that reversal to asymptomatic AD of mild stage patients
was not possible. In the base-case analysis, these effects were
assumed to be constant throughout the duration of treatment;
long-term clinical studies suggest that treatment efficacy may last
for up to 3 years [21,25].
5. Diagnostic tests
The standard diagnostic strategy of AD was assumed to
comprise a detailed history of the patient, an assessment of
cognition and functional status using a questionnaire test such as
the MMSE, and laboratory testing. Other strategies combined
MRI (with or without a new contrast agent) for AD diagnosis with
a questionnaire test such as the MMSE aimed at determining the
stage of the disease.
In the case of severe AD, it was assumed that dementia
symptoms allowed direct diagnosis over a consultation (without
need for a diagnostic test). Hence, all severe stage AD patients
were assumed to be diagnosed and to receive treatment.
As in [10], we estimated the sensitivity of the standard
diagnostic tests at 75% and their specificity at 90%; in early stage
AD, we assumed that patients were asymptomatic and that the
standard diagnostic therefore performed as it did in non-AD
patients, declaring only 10% ( = 120.9) of them as having AD.
Based on data from a clinical study, we also hypothesized a
sensitivity of 50% in early stage, 88% in mild stage and 95% in
moderate stage AD, as well as a 96% specificity, for standard MRI
diagnosis [31].
Regarding the hypothetical diagnostic test using MRI with the
new contrast agent (MRI+CLP), we based our assumptions on
available data on PET-scan amyloid imaging. As studies on
amyloid plaques suggest that amyloid deposition reaches a plateau
by the early clinical stages of Alzheimers disease (amyloid cascade
hypothesis), we assumed the sensitivity of the MRI+CLP
diagnostic test to be independent of the disease stage. In one
study on PET tracers, the sensitivity and specificity were estimated
at 90% [32]; in another, a 95% sensitivity and 83% specificity
were found [33]. Here, we assumed a 96% sensitivity and 87%
specificity for MRI+CLP.
19
20
6.4. Costs of care. We took into account both living and care
costs. For institutionalized patients, living costs included hostel
costs, based on a French study of AD patients [18], and costs of
care included caregiver wages. In an earlier study, caregivers were
estimated to spend 517 hours over a 3 month period caring for
moderate-to-severe AD patients [38]. Here, we therefore
computed costs associated with caregivers over 6 months by
multiplying hourly wages (estimated at 13 J/h for professional
caregivers) with 1034 hours.
For patients living at home, basic costs included living and
medical expenditures [39], based on a French study of home-cared
AD patients and controls; the cost of care by unprofessional
caregivers was also estimated as an opportunity cost. In order to do
this, we valued informal care at a conservative value of 8.4 J/h,
and we assumed that these unpaid caregivers spent the same
amount of time caring for AD patients than professional
caregivers, that is, 1034 hours over a 6 month cycle.
6.5. Indirect costs. In this study, we took into account
several indirect costs of AD. First, we assigned an opportunity cost
to unprofessional caregivers who took care of AD patients leaving
at home, as described before.
Second, we took into account the burden of AD care on the
health state of unprofessional caregivers. Based on an earlier study
[38], 35% of AD caregivers take medication related to their
activity; they have high rates of depression and anxiety, as well as
high overall morbidity and mortality rates, compared to non-AD
caregiver controls [40]. We estimated the cost of this health impact
as that as of a weekly psychiatrist consultation, plus that of an
antidepressant treatment, for 35% of unprofessional caregivers.
Finally, we also evaluated the cost associated with the loss of
productivity of AD patients. This is generally not done in costeffectiveness studies of AD, since AD patients are for the most part
retired. However, recent data shows that pensioners are becoming
more and more involved in volunteer activities within nonprofit
organizations (NPOs), as well as performing informal volunteer
activities [41].
Here, based on French data [42], we estimated that 60 to 75
year-old individuals performed on average 63.8 hours of informal
volunteering activities over a 6 month period mostly within the
family sphere, such as childcare for instance. This was valued at
7.7J/hour, which was the minimum wage in France in 2009, and
multiplied by an efficiency coefficient of 0.7 to arbitrarily take
into account the reduced productivity in older individuals.
Similarly, we estimated that 50% of 60 to 75 year-old
individuals are involved in NPOs, with a mean of 12 hours of
volunteer activities per month. On average, we hence estimated
that 60 to 75 year-old individuals performed a total of 36 hours of
volunteer activities within NPOs over a 6 month period [43],
which were valued at 7.9J/hour and multiplied by the
aforementioned 0.7 efficiency coefficient.
We added the resulting estimated productivity benefit of a
French pensioner to our analyses as a cost associated to AD, in full
for moderate-to-severe patients and multiplied by 0.6 for mild AD
patients (assuming that mild AD only reduces productivity by
40%).
6.6. Effectiveness. We estimated quality-of-life weights
(QALYs) for over-60 patients without Alzheimer disease at 0.826
on a scale of 0 to 1, on the basis of the mean of time trade-off
scores for men and women aged 6584 years old published in a
study of health outcomes in the general population [44].
Quality-of-life weights for patients with Alzheimer disease at
each disease stage and care setting (institution or community) were
based on previously published Health Utilities Index Mark 2
(HUI:2) scores [17].
PLoS ONE | www.plosone.org
7. Sensitivity analysis
The ranges investigated in the Sensitivity Analysis are
summarized in Table S1, along with the data sources. We
performed both univariate and multivariate sensitivity analyses.
For the multivariate analyses, we performed a Monte-Carlo
simulation with 10.000 trials, using a priori triangular distributions
for model parameters. We then identified the most influential
parameters in the cost-effectiveness of the MRI+CLP strategy by
calculating the partial rank correlation coefficient (PRCC) between
each input parameter and the ICER of this strategy and assessing
their statistical significance.
7.1. Drug effects and prices. We investigated the impact of
the introduction of new drugs, which would have the same costs
and indications as donepezil and memantine, respectively, but
with varying efficacy: the probabilities of transitions under
treatment from mild to moderate or moderate to severe stage
were further reduced by a factor (fmM or fMS) ranging from 0.5 to
1, and the probabilities of transitions under treatment from
moderate to mild or severe to moderate stage were further
increased by a factor (fMm or fSM) ranging from 1 to 2 [17]. For
simplicity reasons, we also summarized these four avenues for
improvement of current AD treatment through a single parameter
f, assuming a linear relationship between the multiplying factors:
f~f ML ~f SM ~{2|f LM z3
Regarding the hypothetical drug efficient for early stage AD
(treatment T), we investigated both more and less efficient drugs,
with probabilities of progression from early to mild and moderate
stage AD under treatment T ranging from 0 to their values
without treatment T; fT was the reduction factor applied to these
transition probabilities due to treatment T (in [01]).
The assumed cost for a 6-month cure with treatment T was also
varied, between 0 and 1000J.
7.2. Diagnostic tests characteristics. We investigated
sensitivities from 0.9 to 1 and specificities from 0.7 to 1 in early
to moderate stage AD patients for the hypothetical diagnostic test
using MRI with the CLP contrast product. We also investigated
sensitivities ranging from 0.75 to 0.90 for the standard diagnosis in
moderate stage AD.
In the screen and treat analyses, we investigated sensitivities
of standard MRI between 0.1 and 0.5 for early stage AD.
7.3. Disease progression. To model faster or slower
progression of AD, we investigated state transition probabilities
in our Markov model ranging from 10% lower to 10% higher than
their base-case values.
7.4. Initial distribution of patients and prevalence. In
the primary analysis, we investigated AD prevalences among
consulting individuals ranging from 50 to 70%. For each fixed
prevalence, we varied the proportion of patients with mild AD
between 50 and 75% of all AD patients; the ratio between
moderate and severe AD prevalence remained the same, that is,
10 moderate stage patients for 1 severe stage patient.
In the screen and treat analyses, we investigated initial AD
prevalences among screened individuals ranging from 1 to 10%.
For each fixed prevalence, we varied the proportion of patients
with early asymptomatic AD from 30 to 75% of all AD patients;
the distribution of mild to severe stages remained the same.
7.5. Discount rate. As the screen and treat analyses
spanned a 15 year period, we assessed the impact of variations in
the assumed discount rate for costs and QALYs, from 0 to 10%
per year.
7.6. Frequency of screenings. In the screen and treat
analyses, we investigated the impact of varying durations between
screening campaigns, ranging from one to ten years.
5
Efficacy (in
QALYs)
C/E
ICER
1.7663
20 548
Dominated
Standard MRI
36 131
1.7710
20 401
MRI+CLP
36 313
1.7731
20 480
Strategy
Cost (in J)
88 439 J/QALY
1.7668
20 523
Standard MRI
36 117
1.7712
20 391
MRI+CLP
36 268
1.7737
20 447
Dominated
60 923 J/QALY
doi:10.1371/journal.pone.0035559.t002
Results
1. Primary analysis
1.1. Base Case. The first part of Table 2 summarizes the
results of the primary cost-effectiveness analysis in the base case; in
the second part of Table 2, hypothetical treatment T is offered to
all MCI patients diagnosed with AD but with high MMSE scores.
In both cases, the standard diagnosis strategy was dominated by
the standard MRI strategy (more costly and less effective).
Cost and effectiveness increased from standard MRI to
MRI+CLP strategies. In the base-case, the ICER of the
MRI+CLP strategy was higher than the French willingness-topay threshold (estimated at 76 171 J/QALY). Hence, standard
MRI was found to be the preferred strategy. However, assuming
that treatment T had been made available led to a lower ICER for
MRI+CLP, making it the preferred strategy.
1.2. Sensitivity analysis. Detailed results of the univariate
sensitivity analysis are presented in Table S2 (without treatment
Table 3. Partial rank correlation coefficients (PRCC) between input values and the ICER of the MRI+CLP strategy (compared with
the preferred strategy).
Parameter
PRCC
Parameter
PRCC
20,41
Initial prevalence of AD
20,86
Initial prevalence of AD
20,51
20,15
20,26
20,17
20,92
20,64
20,73
Discount rate
0,29
Discount rate
0,19
0,41
0,79
Cost of treatment T
0,78
0,41
20,23
0,51
All PRCCs are statistically significant at the confidence level of 5%. A higher absolute value of PRCC indicates a strong relationship between that parameter and the ICER;
a positive (resp. negative) value of PRCC implies that the value of the ICER increases (resp. decreases) when the value of the input increases.
doi:10.1371/journal.pone.0035559.t003
21
22
PLoS ONE | www.plosone.org
Figure 2. Results of the primary analysis: multivariate sensitivity analysis. The strategy with maximum net monetary benefit is depicted as
a function of the assumed sensitivity and specificity of the MRI+CLP diagnostic test, for assumed costs of the CLP contrast agent between 0 and 500J
per injection (in the absence of treatment T): (a) Cost of the CLP contrast agent at 50 J/injection. (b) Cost of the CLP contrast agent at 250 J/
injection. (c) Cost of the CLP contrast agent over 450 J/injection.
doi:10.1371/journal.pone.0035559.g002
Discussion
Here, we investigated the cost-effectiveness of a new diagnostic
tool for AD allowing early diagnosis in two different contexts. First,
we assumed that this new diagnostic tool would be made available
to the same individuals who are currently offered other diagnostic
tests. Second, we hypothesized that a treatment with significant
efficacy in early stage AD was developed, and that, as a
consequence, a national screening campaign was put into place,
using either the currently available diagnostic tests or the new
diagnostic test.
We found that in both analyses, the preferred base-case strategy
was the standard MRI strategy. However, multivariate sensitivity
analyses showed that, while in the primary analysis combining
MRI with a new contrast product could prove the preferred
strategy under a wide array of realistic scenarios, the probability of
this happening was inferior to 5% in the hypothesis of a national
screening campaign. Even assuming the availability of a low-cost
highly efficient treatment in early AD, novel contrast agents would
need to have very high specificity to be cost-effective when used for
systematic screening of the entire population.
Strategy
Efficacy (in
Cost (in J) QALYs)
C/E
ICER
Dominated
Population-wide screening
Standard diagnosis 43 559
8.0722
5 396
Standard MRI
43 009
8.0732
5 327
MRI+CLP
44 945
8.0752
5 566
Dominated
8.0377
5 563
Standard MRI
44 180
8.0386
5 496
MRI+CLP
46 075
8.0415
5 730
doi:10.1371/journal.pone.0035559.t004
23
24
PLoS ONE | www.plosone.org
Figure 3. Results of the screen and treat (population-wide screening) analysis: multivariate sensitivity analysis. The strategy with
maximum net monetary benefit is depicted as a function of the assumed efficacy and cost of the hypothetical new drug T, for assumed specificities of
the MRI+CLP diagnosis test between 0.80 and 0.99. The efficacy of treatment T is expressed as a 0-to-1 ratio between assumed probabilities of
transition from early stage AD with and without treatment T; 0 corresponds to maximum efficacy and 1 to no efficacy (in the base case, fT = 0.5: 50%
reduction). Only costs lower than the base-case cost of 500J per 6-month treatment are investigated here. (a) specificity for MRI+CLP inferior to 0.97
(including base case). (b) 0.98 specificity for MRI+CLP. (c) 0.99 specificity for MRI+CLP.
doi:10.1371/journal.pone.0035559.g003
Because the new CLP contrast agent for MRI is still under
development, assumptions had to be made regarding its sensitivity
and specificity for diagnosing AD. As mentioned in the Methods,
we based these assumptions in part on available data on PET-scan
used with amyloid markers [32,33]. We also investigated relatively
wide ranges for these characteristics in our sensitivity analyses.
Regarding the specificity, it should be noted that amyloid
plaques have been detected in 10 to 30% of otherwise apparently
normal elderly subjects. This may limit the capacity of MRI+CLP
to correctly identify healthy subjects, which is why we investigated
specificities as low as 0.70 in our sensitivity analyses. However,
several studies also suggest that the presence of amyloid plaques is
associated with declining cognitive test scores and with structural
and functional brain changes suggestive of early AD [49]. In two
recent follow-up studies, about one third of patients with mild
cognitive impairment in whom amyloid plaques were detected
converted to AD over the following 2 years [50,51]. It could
therefore be argued that at least part of these 1030% amyloidpositive individuals may not be false positive, but very early
stage AD patients.
Figure S5 illustrates the influence of false positives and false
negatives associated with the MRI+CLP diagnosis in terms of
additional cost per QALY in both our primary analysis and our
screen and treat analysis.
3. Included costs
In this study, we chose to include indirect costs in addition to
direct costs. The importance of these indirect costs, which include
costs associated with the loss of productivity of AD patients and
costs associated with informal caregiving at the homes of AD
patients, has indeed been underlined in several recent studies,
which showed that they may dominate direct costs of care in early
stages of the disease [52]. As a consequence, it has been suggested
in recent reviews that future cost-effectiveness studies should take
into account such indirect costs [48,53].
10
25
26
PLoS ONE | www.plosone.org
11
Figure 4. Results of the screen and treat (population-wide screening) analysis: multivariate sensitivity analysis. The strategy with
maximum net monetary benefit is depicted as a function of the assumed prevalence of AD in the general over-60 population (between 0 and 10%)
and specificity of the MRI+CLP diagnostic test between 0.90 and 1, for assumed costs of the CLP contrast agent between 0 and 500J per injection: (a)
Cost of the CLP contrast agent at 50 J/injection. (b) Cost of the CLP contrast agent at 250 J/injection (base-case). (c) Cost of the CLP contrast agent
at 500 J/injection.
doi:10.1371/journal.pone.0035559.g004
Supporting Information
Markov models of Alzheimers disease progression. All states are further subdivided in two, for individuals
living at home vs. inside an institution (retirement or nursing
home). A) Model for the primary analysis. B) Model for the
screen and treat analyses.
(TIF)
Figure S1
Figure S2 Multivariate sensitivity analysis of the primary cost-effectiveness study. A) Distribution of incremental
cost-effectiveness ratios (ICER) of the MRI+CLP strategy,
compared to the standard diagnosis strategy. B) Acceptability
curve: probability that each strategy either is dominant or has an
ICER inferior to the willingness-to-pay, as a function of the
willingness-to-pay threshold.
(TIF)
Figure S3 Multivariate sensitivity analysis of the primary cost-effectiveness study with treatment T. The
strategy with maximum net monetary benefit is depicted as a
function of the assumed efficacy and cost of the hypothetical new
drug T. The efficacy of treatment T is expressed as a 0-to-1 ratio
between assumed probabilities of transition from early stage AD
with and without treatment T; 0 corresponds to maximum efficacy
and 1 to no efficacy (in the base case, fT = 0.5: 50% reduction).
(TIF)
Figure S4 Multivariate sensitivity analysis of the
screen and treat (population-wide screening) costeffectiveness study. A) Distribution of incremental costeffectiveness ratios (ICER) of the MRI+CLP strategy, compared
to the standard MRI strategy. B) Acceptability curve: probability
that the MRI+CLP strategy either is dominant or has an ICER
inferior to the willingness-to-pay, as a function of the willingnessto-pay threshold.
(TIF)
6. Conclusions
Assuming that a treatment with proven efficacy in early AD
becomes available, as well as a diagnostic test allowing early
detection of the disease, the issue of screening the population will
arise. Our study suggests that, in order for this screening to be costeffective, key parameters are the specificity of the new diagnostic
test and the cost and effectiveness of the new treatment. These
preliminary results ought to be taken into account in the currently
underway research on early detection and treatment of AD,
including work on b-amyloid plaques detection and elimination.
Figure S6
12
27
compared to the reference strategy, in the primary costeffectiveness analysis, depending on the values of model
parameters.
(DOCX)
(DOCX)
Table S1 Model parameters: base-case values and
Author Contributions
References
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.
28
25. Atri A, Shaughnessy LW, Locascio JJ, Growdon JH (2008) Long-term course
and effectiveness of combination therapy in Alzheimer disease. Alzheimer Dis
Assoc Disord 22: 209221.
26. Drago V, Babiloni C, Bartres-Faz D, Caroli A, Bosch B, et al. (2011) Disease
tracking markers for Alzheimers disease at the prodromal (MCI) stage.
J Alzheimers Dis 26 Suppl 3: 159199.
27. Blum-Boisgard C, Ulmann P (2007) Medicaments specifiques de la maladie
dAlzheimer. Regime Social des Independants. pp 1130.
28. Bowen J, Teri L, Kukull W, McCormick W, McCurry SM, et al. (1997)
Progression to dementia in patients with isolated memory loss. Lancet 349:
763765.
29. Commenges D, Joly P, Letenneur L, Dartigues JF (2004) Incidence and
mortality of Alzheimers disease or dementia using an illness-death model. Stat
Med 23: 199210.
30. Katz MJ, Lipton RB, Hall CB, Zimmerman ME, Sanders AE, et al. (2011) Agespecific and Sex-specific Prevalence and Incidence of Mild Cognitive
Impairment, Dementia, and Alzheimer Dementia in Blacks and Whites: A
Report From the Einstein Aging Study. Alzheimer Dis Assoc Disord.
31. Harris GJ, Lewis RF, Satlin A, English CD, Scott TM, et al. (1998) Dynamic
susceptibility contrast MR imaging of regional cerebral blood volume in
Alzheimer disease: a promising alternative to nuclear medicine. AJNR
Am J Neuroradiol 19: 17271732.
32. Mormino EC, Kluth JT, Madison CM, Rabinovici GD, Baker SL, et al. (2009)
Episodic memory loss is related to hippocampal-mediated beta-amyloid
deposition in elderly subjects. Brain 132: 13101323.
33. Hansson O, Zetterberg H, Buchhave P, Londos E, Blennow K, et al. (2006)
Association between CSF biomarkers and incipient Alzheimers disease in
patients with mild cognitive impairment: a follow-up study. Lancet Neurol 5:
228234.
34. Ministe`re de la Sante (2009) Tarifs conventionnels des consultations de medecins
specialistes. In: Sante Mdl, ed. Paris: Assurance Maladie, Ministe`re de la Sante.
35. Ministe`re de la Sante (2009) Classification Commune des Actes Medicaux. In:
Sante Mdl, ed. Paris: Assurance Maladie, Ministe`re de la Sante.
36. Ministe`re de la Sante (2009) Nomenclature des actes de biologie medicale. In:
Sante Mdl, ed. Paris: Assurance Maladie, Ministe`re de la Sante.
37. Dictionnaire VIDAL. Paris: Editions du Vidal.
38. Feldman H, Gauthier S, Hecker J, Vellas B, Hux M, et al. (2004) Economic
evaluation of donepezil in moderate to severe Alzheimer disease. Neurology 63:
644650.
39. Rigaud AS, Fagnani F, Bayle C, Latour F, Traykov L, et al. (2003) Patients with
Alzheimers disease living at home in France: costs and consequences of the
disease. J Geriatr Psychiatry Neurol 16: 140145.
40. Elliott AF, Burgio LD, Decoster J (2010) Enhancing caregiver health: findings
from the resources for enhancing Alzheimers caregiver health II intervention.
J Am Geriatr Soc 58: 3037.
41. Reimat A (2002) Production associative et be nevolat informel : quelle
signification economique pour les activites de production des retraites ?
Innovations 15: 7398.
42. Marbot C (2008) Travailler pour des particuliers : souvent une activite
dappoint. Insee References, Les salaires en France.
43. France Benevolat (2010) La situation du benevolat en France en 2010. IFOP France Benevolat.
44. Fryback DG, Dasbach EJ, Klein R, Klein BE, Dorn N, et al. (1993) The Beaver
Dam Health Outcomes Study: initial catalog of health-state quality factors. Med
Decis Making 13: 89102.
45. Gagnon M, Rive B, Hux M, Guilhaume C (2007) Cost-effectiveness of
memantine compared with standard care in moderate-to-severe Alzheimer
disease in Canada. Can J Psychiatry 52: 519526.
46. Lejeune C, Bismuth MJ, Conroy T, Zanni C, Bey P, et al. (2005) Use of a
decision analysis model to assess the cost-effectiveness of 18F-FDG PET in the
13
47.
48.
49.
50.
51. Wolk DA, Price JC, Saxton JA, Snitz BE, James JA, et al. (2009) Amyloid
imaging in mild cognitive impairment subtypes. Ann Neurol 65: 557568.
52. Reese JP, Hessmann P, Seeberg G, Henkel D, Hirzmann P, et al. (2011) Cost
and Care of Patients with Alzheimers Disease: Clinical Predictors in German
Health Care Settings. J Alzheimers Dis.
53. Mauskopf J, Mucha L (2011) A review of the methods used to estimate the cost
of Alzheimers disease in the United States. Am J Alzheimers Dis Other Demen
26: 298309.
54. Schachter F, Faure-Delanef L, Guenot F, Rouger H, Froguel P, et al. (1994)
Genetic associations with human longevity at the APOE and ACE loci. Nat
Genet 6: 2932.
55. Caselli RJ, Walker D, Sue L, Sabbagh M, Beach T (2010) Amyloid load in
nondemented brains correlates with APOE e4. Neurosci Lett 473: 168171.
29
PLoS ONE | www.plosone.org
14
Abstract
The role of structural brain magnetic resonance imaging (MRI) is becoming more and more emphasized in the early
diagnostics of Alzheimers disease (AD). This study aimed to assess the improvement in classification accuracy that can be
achieved by combining features from different structural MRI analysis techniques. Automatically estimated MR features
used are hippocampal volume, tensor-based morphometry, cortical thickness and a novel technique based on manifold
learning. Baseline MRIs acquired from all 834 subjects (231 healthy controls (HC), 238 stable mild cognitive impairment (SMCI), 167 MCI to AD progressors (P-MCI), 198 AD) from the Alzheimers Disease Neuroimaging Initiative (ADNI) database
were used for evaluation. We compared the classification accuracy achieved with linear discriminant analysis (LDA) and
support vector machines (SVM). The best results achieved with individual features are 90% sensitivity and 84% specificity
(HC/AD classification), 64%/66% (S-MCI/P-MCI) and 82%/76% (HC/P-MCI) with the LDA classifier. The combination of all
features improved these results to 93% sensitivity and 85% specificity (HC/AD), 67%/69% (S-MCI/P-MCI) and 86%/82% (HC/
P-MCI). Compared with previously published results in the ADNI database using individual MR-based features, the presented
results show that a comprehensive analysis of MRI images combining multiple features improves classification accuracy and
predictive power in detecting early AD. The most stable and reliable classification was achieved when combining all
available features.
Citation: Wolz R, Julkunen V, Koikkalainen J, Niskanen E, Zhang DP, et al. (2011) Multi-Method Analysis of MRI Images in Early Diagnostics of Alzheimers
Disease. PLoS ONE 6(10): e25446. doi:10.1371/journal.pone.0025446
Editor: Celia Oreja-Guevara, University Hospital La Paz, Spain
Received July 12, 2011; Accepted September 5, 2011; Published October 13, 2011
Copyright: 2011 Wolz et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This project is partially funded under the 7th Framework Programme by the European Commission (http://cordis.europa.eu/ist/). Data collection and
sharing for this project was funded by the Alzheimers Disease Neuroimaging Initiative (ADNI; Principal Investigator: Michael Weiner; NIH grant U01 AG024904).
ADNI is funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering (NIBIB), and through generous contributions
from the following: Pzer Inc., Wyeth Research, Bristol-Myers Squibb, Eli Lilly and Company, GlaxoSmithKline, Merck & Co. Inc., AstraZeneca AB, Novartis
Pharmaceuticals Corporation, Alzheimers Association, Eisai Global Clinical Development, Elan Corporation, plc, Forest Laboratories, and the Institute for the Study
of Aging, with participation from the U.S. Food and Drug Administration. Industry partnerships are coordinated through the Foundation for the National Institutes
of Health. The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimers Disease
Cooperative Study at the University of California San Diego. ADNI data are disseminated by the Laboratory of Neuro Imaging at the University of California Los
Angeles. Data used in the preparation of this article were obtained from the ADNI database (www.loni.ucla.edu/ADNI). As such, the investigators within the ADNI
contributed to the design and implementation of ADNI and/or provided data but did not participate in analysis or writing of this report. ADNI investigators
include (complete listing available at http://adni.loni.ucla.edu/wp-content/uploads/how_to_apply/ADNI_Authorship_List.pdf).
Competing Interests: The authors have declared that no competing interests exist.
* E-mail: [email protected]
. These authors contributed equally to this work.
" For information on the Alzheimers Disease Neuroimaging Initiative please see the Acknowledgments section.
Introduction
30
Subjects
In the ADNI study, brain MR images were acquired at regular
intervals after an initial baseline scan from approximately 200
cognitively normal older subjects (HC), 400 subjects with mild
cognitive impairment (MCI), and 200 subjects with early AD.
Detailled inclusion/exclusion criteria used for the different subject
groups in ADNI are defined in [30]. The AD group has scores
between 2026 (inclusive) on the Mini-Mental State Examination
(MMSE) [31], and a Clinical Dementia Rating (CDR) [32] of 0.5
or 1.0. Furthermore, these subjects fulfil the NINCDS/ADRDA
criteria for probable AD [33]. MCI subjects included have MMSE
scores between 2430 (inclusive), a memory complaint, have
objective memory loss measured by education adjusted scores on
Wechsler Memory Scale Logical Memory II, a CDR of 0.5,
absence of significant levels of impairment in other cognitive
domains, essentially preserved activities of daily living, and an
absence of dementia [30]. Healthy subjects have MMSE scores
between 2430 (inclusive), a CDR of 0, are non-depressed, non
MCI, and nondemented. A more detailed description of the ADNI
study is given in Appendix S1.
All 834 ADNI subjects (231 HC, 238 S-MCI, 167 P-MCI, 198
AD) for which a 1.5T T1-weighted MRI scan at baseline was
available were included in this study. 167 subjects in the MCI
group converted to AD as of July 2011. We therefore
independently analysed progressive MCI (P-MCI) subjects and
subjects with a stable diagnosis of MCI (S-MCI).
Table 1 shows the demographics for the 834 study subjects.
Statistically significant differences in the demographics and clinical
variables between the study groups were assessed using Students
unpaired t-test. In this work, the difference was considered
statistically significant if pv0.05 if not stated otherwise. There
were more men than women in all other groups besides the AD
group. MMSE scores were significantly different in the pairwise
comparisons between all study groups. CDR scores of the HC and
AD groups are significantly different to the ones of the two MCI
groups. Healthy subjects had a significantly lower Geriatric
Depression Scale (GDS) compared to all other groups. Compared
to all other groups, AD subjects had significantly shorter education.
MRI Acquisition
Standard 1.5T screening/baseline T1-weighted images obtained using volumetric 3D MPRAGE protocol with resolutions
Table 1. Subjects.
Group
HC
S-MCI
P-MCI
AD
231
238
167
198
Men
52%
66%
62%
52%
Age
76.02 (5.0)
74.85 (7.8)
74.6 (7.0)
75.68 (7.7)
MMSE
29.1* (1.0)
27.3* (1.8)
26.6* (1.7)
23.3* (2.0)
CDR
0 (0)
0.49 (0.05)
0.50 (0)
0.75 (0.25)
GDS
0.83* (1.14)
1.60 (1.42)
1.53 (1.30)
1.67 (1.42)
Education
16.0 (2.8)
15.6 (3.1)
15.7 (2.9)
14.7* (3.1)
APOE4 status
(e3e4/e4e4)
23%/2%
31%/8%
50%/16%
42%/18%
Months to
conversion
18.2 (10.1)
31
Feature extraction
32
Study design
Table 2 presents an overview on the features calculated for all
834 available ADNI baseline images. All feature values were
corrected for age and gender using a linear regression model
where control subjects were used as the training set, i.e., the
normal, not disease-related, age and gender related differences in
the classification features were removed. Feature selection was
then carried out on the corrected feature sets using stepwise
regression [51].
We used two subsets to perform classification:
I.
II.
Figure 1. 2D manifold embedding of a set of images acquired from healthy controls (red) and subjects with AD (blue).
doi:10.1371/journal.pone.0025446.g001
Classification methods
We used two different widely used methods to perform
classification based on individual features and their combination:
Figure 2. Results for voxelwise t-tests for statistically significant group differences with features extracted from TBM.
doi:10.1371/journal.pone.0025446.g002
33
Figure 3. Results for t-tests for statistically significant group differences based on cortical thickness measurements.
doi:10.1371/journal.pone.0025446.g003
Discussion
Results
Method
No of features
Description
Cortical thickness
9 (HC vs AD)
(CTH)
7 (HC vs P-MCI)
84
average Jacobian of atrophic voxels within a ROI, weighted based on voxel-wise p-values
20
8 (S-MCI vs P-MCI)
34
doi:10.1371/journal.pone.0025446.t002
Feature
LDA
SVM
Feature
SEN
SPE
SEN
SPE
MBL
87
83
85 [64 100]
87
83
HV
81
79
84
77
CTH
89
71
90
TBM
90
84
87 [71 100]
89
All
89 [71 100]
93
85
86 [71 100]
94
SEN
SPE
MBL
64
66
77
48
HV
63
67
62 [36 86]
83
33
73
CTH
63
45
59 [36 79]
96
03
84
TBM
65
62
77
44
78
All
68 [43 93]
67
69
60 [36 86]
92
14
SEN
SPE
LDA
SVM
Feature
LDA
SVM
SEN
SPE
SEN
SPE
MBL
81
75
84
69
HV
77
76
78 [54 92]
83
71
CTH
85
65
77 [54 100]
89
62
TBM
82
76
85
74
All
84 [62 100]
86
82
82 [62 100]
93
67
Feature
SEN
HC vs AD
SPE
HC vs P-MCI
S-MCI vs P-MCI
SEN
SEN
SPE
SPE
MBL
90
74
84
92
55
76
HV
80
69
75
76
63
70
CTH
85
75
86
59
72
35
TBM
93
76
90
84
63
59
All
94
76
94
89
69
54
doi:10.1371/journal.pone.0025446.t006
35
Table 7. Classification results of healthy control (HC), mild cognitive impairment (MCI) and Alzheimers disease subjects reported
in the recent literature.
Study
Features
HC vs AD
HC vs P-MCI
S-MCI vs P-MCI
CCR
SEN
SPE
CCR
SEN
SPE
CCR
SEN
SPE
333
Cortical volumes
91
92
90
70
Hippocampus shape
94
96
92
605
Hippocampus volume
76
75
77
64
60
65
382
Cortical thickness
85
73
75
68
312
Amygdala/caudate volumes
69
76
68
356
SPARE-AD index
56
95
38
509
Various
81
95
73
85
62
69
159
81
60
92
14
351
Various volumes
82
73
398
89
83
93
380
STAND score
86
86
36
Supporting Information
Appendix S1
(DOCX)
Appendix S2 ROI-wise features for CTH and TBM.
(DOCX)
Author Contributions
Acknowledgments
Data used in the preparation of this article were obtained from the ADNI
database (www.loni.ucla.edu/ADNI). As such, the investigators within the
References
24. Liu Y, Paajanen T, Zhang Y, Westman E, Wahlund LO, et al. (2009)
Combination analysis of neuropsychological tests and structural MRI measures
in differentiating AD, MCI and control groupsThe AddNeuroMed study.
Neurobiology of Aging In Press, Corrected Proof: -.
25. Liu Y, Paajanen T, Zhang Y, Westman E, Wahlund LO, et al. (2010) Analysis of
regional MRI volumes and thicknesses as predictors of conversion from mild
cognitive impairment to Alzheimers disease. Neurobiology of Aging 31:
13751385.
26. McEvoy L, Fennema-Notestine C, JC R, Hagler D, Holland D, et al. (2009)
Alzheimer disease: quantitative structural neuroimaging for detection and
prediction of clinical and structural changes in mild cognitive impairment.
Radiology 251: 1950205.
27. Querbes O, Aubry F, Pariente J, Lotterie JA, Dmonet JF, et al. (2009) Early
diagnosis of alzheimers disease using cortical thickness: impact of cognitive
reserve. Brain 132: 20362047.
28. Mueller SG, Weiner MW, Thal LJ, Petersen RC, Jack C, et al. (2005) The
Alzheimers Disease Neuroimaging Initiative. Neuroimaging Clinics of North
America 15: 869877.
29. Cuingnet R, Gerardin E, Tessieras J, Auzias G, Lehricy S, et al. (2010)
Automatic classification of patients with Alzheimers disease from structural
MRI: A comparison of ten methods using the ADNI database. NeuroImage In
Press, Corrected Proof: -.
30. Petersen RC, Aisen PS, Beckett LA, Donohue MC, Gamst AC, et al. (2010)
Alzheimers Disease Neuroimaging Initiative (ADNI): clinical characterization.
Neurology 74: 201209.
31. Folstein MF, Folstein SE, McHugh PR (1975) Mini-mental state: A practical
method for grading the cognitive state of patients for the clinician. Journal of
Psychiatric Research 12(3): 189198.
32. Morris J (1993) The Clinical Dementia Rating (CDR): current version and
scoring rules. Neurology 43: 24122414.
33. McKhann G, Drachman D, Folstein M, Katzman R, Price D, et al. (1984)
Clinical diagnosis of alzheimers disease: report of the nincds-adrda work group
under the auspices of department of health and human services task force on
alzheimers disease. Neurology 34: 939944.
34. Jack CR, Jr., Bernstein MA, Fox NC, Thompson P, Alexander G, et al. (2008)
The Alzheimers disease neuroimaging initiative (ADNI): MRI methods. Journal
of Magnetic Resonance Imaging 27: 685691.
35. Lotjonen JM, Wolz R, Koikkalainen JR, Thurfjell L, Waldemar G, et al. (2010)
Fast and robust multi-atlas segmentation of brain magnetic resonance images.
NeuroImage 49: 23522365.
36. Lotjonen J, Wolz R, Koikkalainen J, Julkunen V, Thurfjell L, et al. (2011) Fast
and robust extraction of hippocampus from mr images for diagnostics of
alzheimers disease. NeuroImage 56: 185196.
37. Van Leemput K, Maes F, Vandermeulen D, Suetens P (1999) Automated
model-based tissue classification of MR images of the brain. IEEE Transactions
on Medical Imaging 18: 897908.
38. Lerch J, Evans A (2005) Cortical thickness analysis examined through power
analysis and a population simulation. NeuroImage 24: 163173.
39. Mazziotta J, Toga A, Evans A, Fox P, Lancaster J, et al. (2001) A probabilistic
atlas and reference system for the human brain: International Consortium for
Brain Mapping (ICBM). Philos Trans R Soc Lond B Biol Sci 356: 12931322.
40. Sled JG, Zijdenbos AP, Evans AC (1998) A nonparametric method for
automatic correction of intensity nonuniformity in MRI data. IEEE Transactions on Medical Imaging 17: 8797.
41. Smith SM (2002) Fast robust automated brain extraction. Hum Brain Mapp 17:
143155.
42. Zijdenbos AP, Forghani R, Evans AC (1998) Automatic quantification of MS
lesions in 3D MRI brain data sets: Validation of INSECT. In: III WMW,
Colchester ACF, Delp SL, eds. MICCAI Springer, volume 1496 of Lecture
Notes in Computer Science. pp 439448.
43. Tohka J, Zijdenbos A, Evans A (2004) Fast and robust parameter estimation for
statistical partial volume models in brain mri. NeuroImage 23: 8497.
44. Kim JS, Singh V, Lee JK, Lerch J, Ad-Dabbagh Y, et al. (2005) Automated 3-d
extraction and evaluation of the inner and outer cortical surfaces using a
laplacian map and partial volume effect classification. NeuroImage 27: 210221.
45. Chung M, Taylor J (2004) Diffusion smoothing on brain surface via finite
element method. In: ISBI IEEE. pp 432435.
37
46. Brun CC, Lepor N, Pennec X, Lee AD, Barysheva M, et al. (2009) Mapping the
regional inuence of genetics on brain structure variability a tensor-based
morphometry study. NeuroImage 48: 3749.
47. Koikkalainen J, Lotjonen J, Thurfjell L, Rueckert D, Waldemar G, et al. (2011)
Multi-template tensor-based morphometry: Application to analysis of alzheimers disease. NeuroImage 56: 11341144.
48. Belkin M, Niyogi P (2003) Laplacian eigenmaps for dimensionality reduction
and data representation. Neural Computation 15: 13731396.
49. Rueckert D, Sonoda LI, Hayes C, Hill DLG, Leach MO, et al. (1999) Nonrigid
registration using free-form deformations: Application to breast MR images.
IEEE Transactions on Medical Imaging 18: 712721.
50. Wolz R, Heckemann RA, Aljabar P, Hajnal JV, Hammers A, et al. (2010)
Measurement of hippocampal atrophy using 4D graph-cut segmentation:
Application to ADNI. NeuroImage 52: 10091018.
51. Draper NR, Smith H (1998) Applied Regression Analysis (Wiley Series in
Probability and Statistics). Wiley.
52. Krzanowski WJ (1988) Principles of Multivariate Analysis: A Users Perspective
Oxford University Press.
53. Cortes C, Vapnik V (1995) Support-vector networks. Machine Learning 20:
273297.
54. Vemuri P, Gunter JL, Senjem ML, Whitwell JL, Kantarci K, et al. (2008)
Alzheimers disease diagnosis in individual subjects using structural mr images:
Validation studies. NeuroImage 39: 11861197.
55. Westman E, Simmons A, Zhang Y, Muehlboeck JS, Tunnard C, et al. (2011)
Multivariate analysis of mri data for alzheimers disease, mild cognitive
impairment and healthy controls. NeuroImage 54: 11781187.
38
PLoS ONE | www.plosone.org
Abstract
Cerebral spinal fluid (CSF) and structural imaging markers are suggested as biomarkers amended to existing diagnostic
criteria of mild cognitive impairment (MCI) and Alzheimers disease (AD). But there is no clear instruction on which markers
should be used at which stage of dementia. This study aimed to first investigate associations of the CSF markers as well as
volumes and shapes of the hippocampus and lateral ventricles with MCI and AD at the baseline and secondly apply these
baseline markers to predict MCI conversion in a two-year time using the Alzheimers Disease Neuroimaging Initiative (ADNI)
cohort. Our results suggested that the CSF markers, including Ab42, t-tau, and p-tau, distinguished MCI or AD from NC,
while the Ab42 CSF marker contributed to the differentiation between MCI and AD. The hippocampal shapes performed
better than the hippocampal volumes in classifying NC and MCI, NC and AD, as well as MCI and AD. Interestingly, the
ventricular volumes were better than the ventricular shapes to distinguish MCI or AD from NC, while the ventricular shapes
showed better accuracy than the ventricular volumes in classifying MCI and AD. As the CSF markers and the structural
markers are complementary, the combination of them showed great improvements in the classification accuracies of MCI
and AD. Moreover, the combination of these markers showed high sensitivity but low specificity for predicting conversion
from MCI to AD in two years. Hence, it is feasible to employ a cross-sectional sample to investigate dynamic associations of
the CSF and imaging markers with MCI and AD and to predict future MCI conversion. In particular, the volumetric
information may be good for the early stage of AD, while morphological shapes should be considered as markers in the
prediction of MCI conversion to AD together with the CSF markers.
Citation: Yang X, Tan MZ, Qiu A (2012) CSF and Brain Structural Imaging Markers of the Alzheimers Pathological Cascade. PLoS ONE 7(12): e47406. doi:10.1371/
journal.pone.0047406
Editor: Yong Fan, Institution of Automation, CAS, China
Received April 13, 2012; Accepted September 13, 2012; Published December 19, 2012
Copyright: 2012 Yang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The work was supported by grants: a centre grant from the National Medical Research Council (NMRC/CG/NUHS/2010), the Young Investigator Award
at the National University of Singapore (NUSYIA FY10 P07), and the National University of Singapore MOE AcRF Tier 1. The funders had no role in study design,
data collection and analysis, decision to publish, or preparation of the manuscript.
Competing Interests: The authors have declared that no competing interests exist.
* E-mail: [email protected]
Introduction
39
Methods
The ADNI was launched in 2003 by the National Institute on
Aging (NIA), the National Institute of Biomedical Imaging and
Bioengineering (NIBIB), the Food and Drug Administration [21],
private pharmaceutical companies and non-profit organizations,
as a $60 million, 5-year publicprivate partnership. The primary
goal of ADNI has been to test whether serial MRI, PET, other
biological markers, and clinical and neuropsychological assessments can be combined to measure the progression of MCI and
early AD. Determination of sensitive and specific markers of very
early AD progression is intended to aid researchers and clinicians
to develop new treatments and monitor their effectiveness, as well
as lessen the time and cost of clinical trials.
ADNI is the result of efforts of many coinvestigators from a
broad range of academic institutions and private corporations, and
subjects have been recruited from over 50 sites across the U.S. and
Canada. The initial goal of ADNI was to recruit 800 adults, ages
55 to 90, to participate in the research approximately 200
cognitively normal older individuals to be followed for 3 years, 400
people with MCI to be followed for 3 years, and 200 people with
early AD to be followed for 2years (see www.adni-info.org for upto-date information). The data were analyzed anonymously, using
publicly available secondary data from the ADNI study, therefore
no ethics statement is required for this work.
Subjects
The ADNI general eligibility criteria are described at www.
adni- info.org. Briefly, subjects are between 5590 years of age,
having a study partner able to provide an independent evaluation
of functioning. Specific psychoactive medications will be excluded.
General inclusion/exclusion criteria are as follows: 1) healthy
subjects: Mini- Mental State Examination (MMSE) scores between
2430, a Clinical Dementia Rating (CDR) of 0, non-depressed,
non-MCI, and nondemented; 2) MCI subjects: MMSE scores
between 2430, a memory complaint, having objective memory
loss measured by education adjusted scores on Wechsler Memory
Scale Logical Memory II, a CDR of 0.5, absence of significant
levels of impairment in other cognitive domains, essentially
preserved activities of daily living, and an absence of dementia;
and 3) mild AD: MMSE scores between 2026, CDR of 0.5 or 1.0,
and meets the National Institute of Neurological and Communicative Disorders and Stroke and the Alzheimers Disease and
Related Disorders Association (NINCDS/ADRDA) criteria for
probable AD.
In this study, 383 subjects were chosen from our previous study
[22]. Within this group, 218 have both MRI and CSF baseline
data (age: 74.467.2 years), with 72 normal controls (NC), 35 AD
subjects, and 111 MCI patients. Amongst these 383 subjects, 25
subjects with MCI converted to AD within 24 months.
Structural MR scans were collected across a variety of scanners
with protocols individualized for each scanner, as defined at www.
loni.ucla.edu/ADNI/Research/Cores/index.shtml. The CSF
Ab42, t-tau and p-tau data were downloaded from the ADNI
web site (www.loni.ucla.edu/ADNI).
MRI analysis
40
Statistical Analysis
Hippocampal volume and shape markers. Bilateral hippocampal volumes distinguished MCI and AD subjects from
normal controls at an accuracy of 61.9% and 65.5% respectively,
with relatively high specificity (MCI: 66.1%; AD: 73.3%) but low
sensitivity (MCI: 57.7%; AD: 57.8%) (Table 2). However, the
hippocampal volumes lost statistical power in the separation of
subjects with MCI and AD (classification accuracy: 42.3%;
sensitivity: 45.3%; specificity: 39.2%, Table 2).
In shape analysis, the first 20 ISOMAP components characterized bilateral hippocampal shapes among all 383 subjects. The 1st,
2nd, 3rd, 4th, 9th, 13th, 18th components contributed to hippocampal shape differences between NC and AD. Among these
components, most of them (the 1st, 2nd, 3rd, 7th, 13th) also
contributed to shape differences between NC and MCI (Table 3).
Interestingly, only the 4th and 9th components showed shape
differences between MCI and AD. Moreover, among these
ISOMAP components, the 1st and 2nd components were highly
correlated with the hippocampal volume (Pearson correlations:
r = 0.7, p,0.01 for the 1st component; r = 0.5, p,0.01 for the 2nd
component), which were dominant components for group
differences in hippocampal shapes between NC and MCI. This
is illustrated as relatively homogeneous shrinkage over the bilateral
hippocampi in Figure 2 (a,b). But the 4th and 9th components
were not associated with the hippocampal volume (p.0.05),
suggesting that only local hippocampal shapes contributed to the
difference between MCI and AD, as seen in Figure 2 (c,d). This
can be further supported by evidence of increased classification
accuracy rates for the classifications between NC and AD (79.2%),
NC and MCI (67.4%), and MCI and AD (57.2%) (Table 2) when
the ISOMAP embedding of hippocampal shapes were used in the
SVM classifiers.
Lateral ventricular volume and shape markers. The
volumes of bilateral lateral ventricles distinguished subjects with
MCI and AD from normal controls at the accuracy of 63.1% and
65.5% respectively, with relatively high specificity (MCI: 75.3%;
AD: 72.3%) but low sensitivity (MCI: 50.9%; AD: 58.8%)
(Table 2). However, volumes of the lateral ventricles lost statistical
Results
Demographic Information
Demographic information for different diagnosis groups at
baseline are shown in Table 1. No significant differences in age
were found among the NC, MCI, and AD groups (ANOVA,
p = 0.454). 25 out of 111 MCI subjects were diagnosed as AD at
the two-year follow up and were denoted as the MCI-c group.
Rest of the MCI subjects were placed in the MCI-s group. No
significant MMSE difference was found between the two MCI
groups at baseline (p = 0.10).
PLOS ONE | www.plosone.org
41
Table 1. Demographic information for each of the diagnosis groups (normal controls (NC), mild cognitive impairment (MCI), and
Alzheimers disease (AD)) at the baseline.
Group
Subjects n
Age (meanSD)
Gender (female/male)
MMSE (meanSD)
NC
72
75.265.2
35/37
2961
MCI-s
86
7467.7
28/58
26.761.8
MCI-c
25
73.566.9
6/19
26.161.5
AD
35
74.669.3
15/20
22.961.8
Note: SD standard deviation; MMSE mini-mental state examination; MCI-s subjects with MCI who remained as MCI at the two-year follow up; MCI-c subjects with
MCI who converted as AD at the two-year follow up.
doi:10.1371/journal.pone.0047406.t001
Table 2. The classification accuracy, sensitivity, and specificity of the support vector machine (SVM) classifiers are given for
distinguishing normal controls (NC) and subjects with Alzheimers disease (AD), NC and subjects with mild cognitive impairment
(MCI), and subjects with MCI and AD.
NC.vs. AD
NC.vs. MCI
MCI.vs. AD
42.3% (CI:32.2%,52.3%)
Sensitivity = 45.3%, Specificity = 39.2%
Hp shapes
LV shapes
60.1% (CI:57.1%,63.1%)
Sensitivity = 62%, Specificity = 58.3%
CSF, Hp volumes,
LV volumes
CSF, Hp shapes,
LV shapes
CSF
CSF markers
Combination
42
The volumes and shapes of the hippocampus (Hp) and lateral ventricles (LV) as well as cerebral spinal fluid (CSF) markers are respectively used as features in the SVM.
doi:10.1371/journal.pone.0047406.t002
Discussion
Our study demonstrated the dynamic trajectories of the CSF,
hippocampal and lateral ventricular markers in the Alzheimers
pathological cascade using a cross-sectional ADNI sample and also
showed the feasibility of predicting future MCI-to-AD conversion
using baseline CSF and imaging markers. The CSF markers,
including Ab42, t-tau, and p-tau, distinguished MCI or AD from
NC, while only the Ab42 CSF marker contributed to the
differentiation between MCI and AD. The hippocampal shapes
performed better than the hippocampal volumes in classifying NC
and MCI, NC and AD, as well as MCI and AD. Interestingly, as
compared to the ventricular shape, ventricular volume performed
better in distinguishing MCI or AD from NC. The ventricular
shape, however, showed better accuracy in the classification for
MCI and AD. As the CSF and structural markers were
complementary, their combination showed great improvement
in the classification accuracies at all the stages of AD. Moreover,
the combination of these baseline markers also showed high
sensitivity but low specificity for predicting MCI conversion to AD
during a two-year period.
Our findings supported the conclusion drawn in previous
studies [31]; [32], where abnormality of both CSF Ab42 and
neurodegenerative biomarkers, including CSF tau and MRI
markers, precedes clinical symptoms; all these markers showed
significant differences between NC and MCI groups. However,
our findings did not support the hypothesis where CSF Ab42
reaches a plateau before the appearance of MRI atrophy and
cognitive symptoms, and remain static thereafter [7]. In our study,
CSF Ab42 continued to show appreciable power discriminant
between MCI and AD. In contrast, CSF tau lost its discriminating
power in distinguishing MCI and AD patients, suggesting that
CSF Ab42 reaches its plateau after CSF tau in the Alzheimers
pathological cascade.
MCI and AD patients were not well-separated using hippocampal volume, implying that the overall tissue loss in the
hippocampus may not be a good marker for monitoring AD
progression. However, we may not conclude that the hippocampus
Table 3. ISOMAP components of the hippocampus and lateral ventricles as well as CSF markers contribute to the group
differences between normal controls (NC) and subjects with Alzheimers disease (AD), NC and subjects with mild cognitive
impairment (MCI), and subjects with MCI and AD.
NC.vs. AD
NC.vs. MCI
MCI vs AD
Imaging Markers
ISOMAP components
ISOMAP components
ISOMAP components
hippocampal shapes
1,2,3,4,9,13,18
1,2,3,7,13
4,9
1,7,8,9,17,19
1,8,12,14
9, 13,17
Ab42
CSF Markers
CSF Ab42, t-tau, p-tau
doi:10.1371/journal.pone.0047406.t003
43
5
Figure 3. Shape differences of the lateral ventricles among normal controls (NC), mild cognitive impairment (MCI), and Alzheimers
Disease (AD). Panels (a,b) respectively show group differences in the left and right lateral ventricular surface deformations between MCI and NC.
Panels (c,d) respectively show group differences in the left and right lateral ventricular surface deformations between AD and MCI. Warm color
denotes regions where structures have surface outward-deformation in the former group when compared with the latter group, while cool color
denotes regions where structures have surface inward-deformation in the former group when compared with the latter group.
doi:10.1371/journal.pone.0047406.g003
reached its abnormality peak before the late stage of AD, as its
local shape variations were significantly associated with progression from MCI to AD. Our results also showed that such local
shape markers aided the hippocampal markers in achieving
slightly better accuracy than the CSF markers in the prediction for
MCI conversion to AD. This was also supported by previous
studies, suggesting that MRI markers (e.g. cortical thickness of the
medial temporal lobe) correlate well with severity of cognitive
impairment and have greater predictive power than the CSF tau
[33]. Based on these evidences, we may conclude that the
hippocampal shape marker reaches its plateau after CSF tau.
However, the order in which the CSF Ab42 and MRI markers
reach their abnormality peaks is still unclear based on our current
study.
The volume of the lateral ventricles cannot distinguish MCI and
AD patients. Interestingly, the overall expansion of the lateral
ventricles showed better performance in identifying MCI or AD
from NC when compared with their shapes. This result agrees
with previous studies, suggesting that rates of ventricular
expansion were significantly different between AD (or MCI) and
NC groups [8]. This implies that complicated shape analysis might
not always be necessary to provide better structural morphological
Table 4. The accuracy, sensitivity, and specificity for predicting the MCI converters are listed when the volumes or shapes of the
hippocampus (Hp) or the lateral ventricles (LV), or the CSF markers, or their combination were used as features in the classification.
Markers
Sensitivity
Specificity
54%(53.1%,54.9%)
54.4%
53.6%
Hp shapes
63%(62.1%,64%)
74.9%
51.2%
55.6%(54.7%,56.5)
63.7%
47.5%
LV shapes
62.7%(61.5%,63.9%)
67.9%
57.5%
62.2%(61.3%,63.1%)
80.4%
44%
81.1%(80.5%,81.6%)
37.2%(35.6%,38.9%)
82%(81.5%,82.6%)
51.4% (49.4%,53.3%)
CSF
CSF markers
Combination
44
66.7%(65.7%,67.8%)
doi:10.1371/journal.pone.0047406.t004
Acknowledgments
Data used in preparation of this article were obtained from the Alzheimers
Disease Neuroimaging Initiative (ADNI) database (adni.loni.ucla.edu). As
such, the investigators within the ADNI contributed to the design and
implementation of ADNI and/or provided data but did not participate in
analysis or writing of this report. A complete listing of ADNI investigators
can be found at: http://adni.loni.ucla.edu/wp-content/uploads/how_to_
apply/ADNI_Acknowledgement_List.pdf.
Author Contributions
Conceived and designed the experiments: XY MZT AQ. Performed the
experiments: XY MZT AQ. Analyzed the data: XY MZT AQ.
Contributed reagents/materials/analysis tools: XY MZT AQ. Wrote the
paper: XY MZT AQ.
References
1. Jack CR, Jr., Albert MS, Knopman DS, McKhann GM, Sperling RA, et al.
(2011) Introduction to the recommendations from the National Institute on
Aging-Alzheimers Association workgroups on diagnostic guidelines for Alzheimers disease. Alzheimers Dement 7: 257262.
2. McKhann GM, Knopman DS, Chertkow H, Hyman BT, Jack CR, Jr., et al.
(2011) The diagnosis of dementia due to Alzheimers disease: recommendations
from the National Institute on Aging-Alzheimers Association workgroups on
diagnostic guidelines for Alzheimers disease. Alzheimers Dement 7: 263269.
3. Albert MS, DeKosky ST, Dickson D, Dubois B, Feldman HH, et al. (2011) The
diagnosis of mild cognitive impairment due to Alzheimers disease: recommendations from the National Institute on Aging-Alzheimers Association workgroups on diagnostic guidelines for Alzheimers disease. Alzheimers Dement 7:
270279.
4. Sperling RA, Aisen PS, Beckett LA, Bennett DA, Craft S, et al. (2011) Toward
defining the preclinical stages of Alzheimers disease: recommendations from the
National Institute on Aging-Alzheimers Association workgroups on diagnostic
guidelines for Alzheimers disease. Alzheimers Dement 7: 280292.
5. Ferrarini L, Palm WM, Olofsen H, van der Landen R, van Buchem MA, et al. (2008)
Ventricular shape biomarkers for Alzheimers disease in clinical MR images.
Magnetic resonance in medicine : official journal of the Society of Magnetic
Resonance in Medicine/Society of Magnetic Resonance in Medicine 59: 260267.
6. Ferrarini L, Frisoni GB, Pievani M, Reiber JH, Ganzola R, et al. (2009)
Morphological hippocampal markers for automated detection of Alzheimers
disease and mild cognitive impairment converters in magnetic resonance images.
Journal of Alzheimers disease : JAD 17: 643659.
7. Jack CR, Jr., Knopman DS, Jagust WJ, Shaw LM, Aisen PS, et al. (2010)
Hypothetical model of dynamic biomarkers of the Alzheimers pathological
cascade. Lancet Neurol 9: 119128.
8. Jack CR, Jr., Lowe VJ, Weigand SD, Wiste HJ, Senjem ML, et al. (2009) Serial
PIB and MRI in normal, mild cognitive impairment and Alzheimers disease:
implications for sequence of pathological events in Alzheimers disease. Brain
132: 13551365.
9. Davatzikos C, Bhatt P, Shaw LM, Batmanghelich KN, Trojanowski JQ (2011)
Prediction of MCI to AD conversion, via MRI, CSF biomarkers, and pattern
classification. Neurobiol Aging 32: 2322 e23192327.
10. Kantarci K, Weigand SD, Przybelski SA, Shiung MM, Whitwell JL, et al. (2009)
Risk of dementia in MCI: combined effect of cerebrovascular disease, volumetric
MRI, and 1H MRS. Neurology 72: 15191525.
11. Qiu A, Younes L, Miller MI, Csernansky JG (2008) Parallel transport in
diffeomorphisms distinguishes the time-dependent pattern of hippocampal
surface deformation due to healthy aging and the dementia of the Alzheimers
type. Neuroimage 40: 6876.
12. Costafreda SG, Dinov ID, Tu Z, Shi Y, Liu CY, et al. (2011) Automated
hippocampal shape analysis predicts the onset of dementia in mild cognitive
impairment. Neuroimage 56: 212219.
13. Westman E, Muehlboeck JS, Simmons A (2012) Combining MRI and CSF
measures for classification of Alzheimers disease and prediction of mild
cognitive impairment conversion. Neuroimage.
14. Cui Y, Liu B, Luo S, Zhen X, Fan M, et al. (2011) Identification of conversion
from mild cognitive impairment to Alzheimers disease using multivariate
predictors. PLoS One 6: e21896.
15. Cardenas VA, Du AT, Hardin D, Ezekiel F, Weber P, et al. (2003) Comparison
of methods for measuring longitudinal brain change in cognitive impairment
and dementia. Neurobiol Aging 24: 537544.
16. Ridha BH, Barnes J, Bartlett JW, Godbolt A, Pepple T, et al. (2006) Tracking
atrophy progression in familial Alzheimers disease: a serial MRI study. Lancet
Neurol 5: 828834.
17. Apostolova LG, Dutton RA, Dinov ID, Hayashi KM, Toga AW, et al. (2006)
Conversion of mild cognitive impairment to Alzheimer disease predicted by
hippocampal atrophy maps. Arch Neurol 63: 693699.
18. Fox NC, Warrington EK, Freeborough PA, Hartikainen P, Kennedy AM, et al.
(1996) Presymptomatic hippocampal atrophy in Alzheimers disease. A
longitudinal MRI study. Brain 119 (Pt 6): 20012007.
19. Wang L, Swank JS, Glick IE, Gado MH, Miller MI, et al. (2003) Changes in
hippocampal volume and shape across time distinguish dementia of the
Alzheimer type from healthy aging. Neuroimage 20: 667682.
20. Apostolova LG, Dinov ID, Dutton RA, Hayashi KM, Toga AW, et al. (2006) 3D
comparison of hippocampal atrophy in amnestic mild cognitive impairment and
Alzheimers disease. Brain 129: 28672873.
21. Alves Fda S, Figee M, Vamelsvoort T, Veltman D, de Haan L (2008) The revised
dopamine hypothesis of schizophrenia: evidence from pharmacological MRI studies
with atypical antipsychotic medication. Psychopharmacol Bull 41: 121132.
22. Qiu A, Fennema-Notestine C, Dale AM, Miller MI, Initiative AsDN (2009)
Regional shape abnormalities in mild cognitive impairment and Alzheimers
disease. Neuroimage 45: 656661.
23. Fischl B, Salat DH, Busa E, Albert M, Dieterich M, et al. (2002) Whole Brain
Segmentation: Automated Labeling of Neuroanatomical Structures in the
Human Brain. Neuron 33: 341355.
24. Miller MI, Qiu A (2009) The emerging discipline of Computational Functional
Anatomy. Neuroimage 45: S1639.
25. Qiu A, Brown T, Fischl B, Ma J, Miller MI (2010) Atlas Generation for
Subcortical and Ventricular Structures with Its Applications in Shape Analysis.
IEEE Transactions on Image Processing 19: 15391547.
26. Qiu A, Miller MI (2008) Multi-structure network shape analysis via normal
surface momentum maps. Neuroimage 42: 14301438.
27. Qiu A, Brown T, Fischl B, Ma J, Miller MI (2010) Atlas generation for
subcortical and ventricular structures with its applications in shape analysis.
IEEE Trans Image Process 19: 15391547.
28. Tenenbaum JB, de Silva V, Langford JC (2000) A global geometric framework
for nonlinear dimensionality reduction. Science 290: 23192323.
29. Yang X, Goh A, Qiu A (2011) Locally Linear Diffeomorphic Metric Embedding
(LLDME) for surface-based anatomical shape modeling. Neuroimage 56: 149
161.
30. Vapnik VN (1995) The Nature of Statistical Learning Theory: Springer.
31. Chetelat G, Desgranges B, De La Sayette V, Viader F, Eustache F, et al. (2002)
Mapping gray matter loss with voxel-based morphometry in mild cognitive
impairment. Neuroreport 13: 19391943.
32. Shaw LM, Vanderstichele H, Knapik-Czajka M, Clark CM, Aisen PS, et al.
(2009) Cerebrospinal fluid biomarker signature in Alzheimers disease neuroimaging initiative subjects. Ann Neurol 65: 403413.
33. Vemuri P, Wiste HJ, Weigand SD, Shaw LM, Trojanowski JQ, et al. (2009)
MRI and CSF biomarkers in normal, MCI, and AD subjects: predicting future
clinical change. Neurology 73: 294301.
34. Vemuri P, Wiste HJ, Weigand SD, Knopman DS, Trojanowski JQ, et al. (2010)
Serial MRI and CSF biomarkers in normal aging, MCI, and AD. Neurology 75:
143151.
35. Ewers M, Walsh C, Trojanowski JQ, Shaw LM, Petersen RC, et al. (2010)
Prediction of conversion from mild cognitive impairment to Alzheimers disease
dementia based upon biomarkers and neuropsychological test performance.
Neurobiol Aging.
36. Zhang D, Wang Y, Zhou L, Yuan H, Shen D (2011) Multimodal classification of
Alzheimers disease and mild cognitive impairment. Neuroimage 55: 856867.
45
Abstract
Background: Visual assessment rating scales for medial temporal lobe (MTL) atrophy have been used by neuroradiologists
in clinical practice to aid the diagnosis of Alzheimers disease (AD). Recently multivariate classification methods for magnetic
resonance imaging (MRI) data have been suggested as alternative tools. If computerized methods are to be implemented in
clinical practice they need to be as good as, or better than experienced neuroradiologists and carefully validated. The aims
of this study were: (1) To compare the ability of MTL atrophy visual assessment rating scales, a multivariate MRI classification
method and manually measured hippocampal volumes to distinguish between subjects with AD and healthy elderly
controls (CTL). (2) To assess how well the three techniques perform when predicting future conversion from mild cognitive
impairment (MCI) to AD.
Methods: High resolution sagittal 3D T1w MP-RAGE datasets were acquired from 75 AD patients, 101 subjects with MCI and
81 CTL from the multi-centre AddNeuroMed study. An automated analysis method was used to generate regional volume
and regional cortical thickness measures, providing 57 variables for multivariate analysis (orthogonal partial least squares to
latent structures using seven-fold cross-validation). Manual hippocampal measurements were also determined for each
subject. Visual rating assessment of MTL atrophy was performed by an experienced neuroradiologist according to the
approach of Scheltens et al.
Results: We found prediction accuracies for distinguishing between AD and CTL of 83% for multivariate classification, 81%
for the visual rating assessments and 89% for manual measurements of total hippocampal volume. The three different
techniques showed similar accuracy in predicting conversion from MCI to AD at one year follow-up.
Conclusion: Visual rating assessment of the MTL gave similar prediction accuracy to multivariate classification and manual
hippocampal volumes. This suggests a potential future role for computerized methods as a complement to clinical
assessment of AD.
Citation: Westman E, Cavallin L, Muehlboeck J-S, Zhang Y, Mecocci P, et al. (2011) Sensitivity and Specificity of Medial Temporal Lobe Visual Ratings and
Multivariate Regional MRI Classification in Alzheimers Disease. PLoS ONE 6(7): e22506. doi:10.1371/journal.pone.0022506
Editor: Jerson Laks, Federal University of Rio de Janeiro, Brazil
Received April 14, 2011; Accepted June 22, 2011; Published July 21, 2011
Copyright: 2011 Westman et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This study was supported by InnoMed (Innovative Medicines in Europe) an Integrated Project funded by the European Union of the Sixth Framework
program priority FP6-2004-LIFESCIHEALTH-5, Life Sciences, Genomics and Biotechnology for Health. AS and SL were supported by funds from the NIHR
Biomedical Research Centre for Mental Health at the South London (www.nihr.ac.uk), Maudsley NHS Foundation Trust and Institute of Psychiatry, Kings College,
London. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing Interests: The authors have the following competing interest: this study was supported by InnoMed (Innovative Medicines in Europe). There are no
patents, products in development or marketed products to declare. This does not alter the authors adherence to all the PLoS ONE policies on sharing data and
materials, as detailed online in the guide for authors.
* E-mail: [email protected]
. These authors contributed equally to this work.
Introduction
Dementia is the third most common cause of death in society
today, exceeded only by cancer and cardiovascular disorders and
Alzheimers disease (AD) is the most common form of dementia.
46
AddNeuroMed
MCI
AD
Number
81
101
75
Female/Male
45/36
52/49
50/25
Age
73.666.3
74.065.8
74.266.0
Education
11.064.8
8.764.3
8.364.2
MMSE
29.061.2
27.261.6
21.364.6
CDR
0.5
1.160.5
ADAS1
3.461.5
5.361.2
6.661.5
CTL
47
MRI
Data acquisition for the AddNeuroMed study was designed to be
compatible with the Alzheimer Disease Neuroimaging Initiative
(ADNI) [16]. The imaging protocol for both studies included a high
resolution sagittal 3D T1-weighted MPRAGE volume (voxel size
1.161.161.2 mm3) and axial proton density/T2-weighted fast spin
echo images. The MPRAGE volume was acquired using a custom
pulse sequence specifically designed for the ADNI study to ensure
compatibility across scanners [16]. Full brain and skull coverage was
required for both of the latter datasets and detailed quality control
carried out on all MR images from both studies according to the
AddNeuroMed quality control procedure [13,14].
48
Figure 1. Representations of ROIs included as candidate input variables in the multivariate OPLS model. (A) Regional volumes. (B)
Regional cortical thickness measures.
doi:10.1371/journal.pone.0022506.g001
Scale
Width of
Choroid fissure
Width of
temporal horn
Hippocampal
thickness
qqq
qqq
qqq
QQQ
49
Results
Subject cohort
252 subjects were included in this study: 75 AD patients, 101
MCI patients and 81 control subjects. Using age as an x-variable
in the OPLS models did not have any effect on the predictive
power of the models separating the groups when all image
variables were included. Therefore, age was excluded from further
analysis. All MRI volumetric measures were normalised by
dividing by each subjects intracranial volume. As expected,
performance on the MMSE, CDR and ADAS1 was poorest
among AD patients and best among controls (Table 1). The MCI
group had scores between the AD and the control groups (Table 1).
50
Figure 3. OPLS cross validated score plots and MRI measures of importance for the separation between AD and CTL. (A) The scatter
plot visualises group separation and the predictability of the AD vs. CTL model. Each black square represents an AD subject and each gray circle a
control subject. Control subjects to the left of zero and AD subjects to the right of zero are falsely predicted. Q2(Y).0.05 (statistically significant
model). (B) Measures above zero have a larger value in controls compared to AD and measures below zero have a lower value in controls compared
to AD. A measure with a high covariance is more likely to have an impact on group separation than a measure with a low covariance. Measures with
jack knifed confidence intervals that include zero have low reliability.
doi:10.1371/journal.pone.0022506.g003
Discussion
Automated computerized MRI methods to aid in the diagnosis
of AD will only be implemented in clinical practice if they are
carefully investigated and validated. The aim of this study was to
further validate the OPLS technique with fully automated MRI
measures as input and compare it to the Scheltens scale for visual
assessment of medial temporal lobe atrophy and to that of total
hippocampal volume. To our knowledge this is the first time that
6
51
Sensitivity
Specificity
Accuracy
LR+
LR2
Visual assessment
78.7 (68.186.4)
82.7 (73.189.4)
80.8 (73.986.2)
4.6 (2.87.4)
0.26 (0.160.40)
77.3(66.785.3)
87.7 (78.793.2)
82.7 (75.687.8)
6.3 (3.511.3)
0.26 (0.170.40)
Manual outlining
93.3 (85.397.1)
85.2(75.991.3)
89.1 (83.293.1)
6.3(3.710.7)
0.08 (0.030.18)
Confidence intervals within parentheses, LR+ = positive likelihood ratio and LR2 = negative likelihood ratio.
doi:10.1371/journal.pone.0022506.t003
52
Method
Number
AD-like
CTL-like
19
68% (13)
32% (6)
19
74% (14)
26% (5)
19
79% (15)
21% (4)
82
32% (26)
68% (56)
82
30% (24)
70% (58)
82
46% (38)
54% (44)
Number
MCI-c
as AD
MCI-s
as CTL
Visual assessment
101
68%
68%
101
79%
68%
101
68%
68%
70%. Although the OPLS method gave the best results, subjects
diagnosed at baseline with MCI who are still classified as MCI at
one year follow up may subsequently convert to AD at a later stage
and thus assessing whether subjects will eventually convert to AD
may require longer follow up times. Further large studies with
longer follow up times are warranted to investigate this issue. To
make the comparison of the performance of the three methods
easier we calculated the sensitivity (MCI-c predicted as AD) at a
fixed specificity (MCI-s predicted as CTL). This yielded the best
sensitivity for the OPLS model with automated regional MRI
measures as input. This model predicted 79% of the MCI-c
subjects who converted to AD at one year follow-up as more ADlike, compared to 68% of the other two methods. Although the
number of converters is relatively small, a likely explanation for
these results is that combining multiple regions across the brain
aids in the prediction of MCI conversion.
Conclusion
Visual rating assessment of the medial temporal lobe gave
similar prediction accuracy to computerized multivariate classification and both accuracies are comparable to that of manual
hippocampal volume measurements. While manual hippocampal
volumes are a valuable research tool and have found utility in
clinical trials, they are not however practical for use in routine
clinical work. Our results demonstrate that computerized
multivariate classification is as good as expert radiological review
of MRI using a validated and widely used visual rating scale,
which suggests a potential future role for computerized methods as
a complement to clinical assessment. Improving classification
models by adding other biomarkers will hopefully increase the
predictive power of the models and improve our understanding of
the etiology of the disease. Two limitations of the current study are
however that the data is not neuropathologically confirmed and
that the follow-up of the MCI subjects was for one year.
In conclusion we believe that this study and previous work
[6,7,9] has shown that the OPLS model with automated MRI
measures as input has the potential to serve as a complement to
clinical assessment of AD, and to target appropriate populations
for clinical trials.
Acknowledgments
The authors thank the Gamla Tjanarinnor Foundation, the Swedish
Alzheimers Association, Swedish Brain Power, Health Research Council
of Academy of Finland and Stockholm Medical Image Laboratory and
Education (SMILE).
Author Contributions
Conceived and designed the experiments: EW LC AS L-OW HS CS J-SM
YZ PM BV MT IK SL. Performed the experiments: AS HS PM BV MT
IK EW J-SM. Analyzed the data: EW LC YZ AS. Contributed reagents/
materials/analysis tools: EW LC AS L-OW HS CS J-SM YZ PM BV MT
IK SL. Wrote the paper: EW LC AS L-OW.
References
1. Kantarci K (2005) Magnetic resonance markers for early diagnosis and
progression of Alzheimers disease. Expert Rev Neurother 5: 663670.
2. Dubois B, Feldman HH, Jacova C, DeKosky ST, Barberger-Gateau P, et al.
(2007) Research criteria for the diagnosis of Alzheimers disease: revising the
NINCDS-ADRDA criteria. The Lancet Neurology 6: 734746.
3. McKhann GM, Knopman DS, Chertkow H, Hyman BT, Jack CR, Jr., et al.
(2011) The diagnosis of dementia due to Alzheimers disease: Recommendations
from the National Institute on Aging-Alzheimers Association workgroups on
diagnostic guidelines for Alzheimers disease. Alzheimers Dement 7: 263269.
53
PLoS ONE | www.plosone.org
28. Fischl B, van der Kouwe A, Destrieux C, Halgren E, Segonne F, et al. (2004)
Automatically parcellating the human cerebral cortex. Cereb Cortex 14: 1122.
29. Liu Y, Paajanen T, Zhang Y, Westman E, Wahlund L-O, et al. (2011)
Combination analysis of neuropsychological tests and structural MRI measures
in differentiating AD, MCI and control groups-The AddNeuroMed study.
Neurobiol Aging 32: 11981206.
30. Liu Y, Paajanen T, Zhang Y, Westman E, Wahlund L-O, et al. (2010) Analysis
of regional MRI volumes and thicknesses as predictors of conversion from mild
cognitive impairment to Alzheimers disease. Neurobiology of Aging 31:
13751385.
31. Tunnard C, Whitehead D, Hurt C, Wahlund L-O, Mecocci P, et al. (2011)
Apathy and cortical atrophy in Alzheimers disease. Int J Geriatr Psychiatry 26:
741748.
32. Liu Y, Paajanen T, Westman E, Wahlund L-O, Simmons A, et al. (2010) Effect
of APOE epsilon4 allele on cortical thicknesses and volumes: the AddNeuroMed
study. J Alzheimers Dis 21: 947966.
33. Liu Y, Paajanen T, Westman E, Zhang Y, Wahlund L-O, et al. (2010) APOE
epsilon2 allele is associated with larger regional cortical thicknesses and volumes.
Dement Geriatr Cogn Disord 30: 229237.
34. Thambisetty M, Simmons A, Velayudhan L, Hye A, Campbell J, et al. (2010)
Association of plasma clusterin concentration with severity, pathology, and
progression in Alzheimer disease. Arch Gen Psychiatry 67: 739748.
35. Pantel J, OLeary DS, Cretsinger K, Bockholt HJ, Keefe H, et al. (2000) A new
method for the in vivo volumetric measurement of the human hippocampus with
high neuroanatomical accuracy. Hippocampus 10: 752758.
36. Bylesjo M, Eriksson D, Kusano M, Moritz T, Trygg J (2007) Data integration in
plant biology: the O2PLS method for combined modeling of transcript and
metabolite data. Plant J 52: 11811191.
37. Johan Trygg SW (2002) Orthogonal projections to latent structures (O-PLS).
Journal of Chemometrics 16: 119128.
38. Rantalainen M, Cloarec O, Beckonert O, Wilson ID, Jackson D, et al. (2006)
Statistically integrated metabonomic-proteomic studies on a human prostate
cancer xenograft model in mice. J Proteome Res 5: 26422655.
39. Wiklund S, Johansson E, Sjostrom L, Mellerowicz EJ, Edlund U, et al. (2008)
Visualization of GC/TOF-MS-based metabolomics data for identification of
biochemically interesting compounds using OPLS class models. Anal Chem 80:
115122.
40. Levine B, Kovacevic N, Nica EI, Cheung G, Gao F, et al. (2008) The Toronto
traumatic brain injury study: Injury severity and quantified MRI. Neurology 70:
771778.
41. McIntosh AR, Lobaugh NJ (2004) Partial least squares analysis of neuroimaging
data: applications and advances. Neuroimage 23: S250S263.
42. Oberg J, Spenger C, Wang FH, Andersson A, Westman E, et al. (2007) Age
related changes in brain metabolites observed by (1)H MRS in APP/PS1 mice.
Neurobiol Aging 29(9): 142333.
43. Westman E, Spenger C, Oberg J, Reyer H, Pahnke J, et al. (2009) In vivo 1Hmagnetic resonance spectroscopy can detect metabolic changes in APP/PS1
mice after donepezil treatment. BMC Neurosci 10: 33.
44. Westman E, Spenger C, Wahlund L-O, Lavebratt C (2007) Carbamazepine
treatment recovered low N-acetylaspartate+N-acetylaspartylglutamate (tNAA)
levels in the megencephaly mouse BALB/cByJ-Kv1.1mceph/mceph. Neurobiology of Disease 26: 221228.
45. Eriksson L, Johansson E, Kettaneh-Wold N, Trygg J, Wiksrom C, et al. (2006)
Multi- and Megavariate Data Analysis (Part I -Basics and Principals and
Applications). Umea: Umetrics AB.
46. Qizilbash N, Schneider L, Chui H, Tariot P, Brodaty H, et al. (2002)
Evindenced-based Dementia Pratice. Oxford, UK: Blackwell Publishing. pp
2023.
47. Blennow K, de Leon MJ, Zetterberg H (2006) Alzheimers disease. The Lancet
368: 387403.
54
PLoS ONE | www.plosone.org
Abstract
The aim of the study was to evaluate the value of assessing white matter integrity using diffusion tensor imaging (DTI) for
classification of mild cognitive impairment (MCI) and prediction of cognitive impairments in comparison to brain atrophy
measurements using structural MRI. Fifty-one patients with MCI and 66 cognitive normal controls (CN) underwent DTI and
T1-weighted structural MRI. DTI measures included fractional anisotropy (FA) and radial diffusivity (DR) from 20
predetermined regions-of-interest (ROIs) in the commissural, limbic and association tracts, which are thought to be involved
in Alzheimers disease; measures of regional gray matter (GM) volume included 21 ROIs in medial temporal lobe, parietal
cortex, and subcortical regions. Significant group differences between MCI and CN were detected by each MRI modality: In
particular, reduced FA was found in splenium, left isthmus cingulum and fornix; increased DR was found in splenium, left
isthmus cingulum and bilateral uncinate fasciculi; reduced GM volume was found in bilateral hippocampi, left entorhinal
cortex, right amygdala and bilateral thalamus; and thinner cortex was found in the left entorhinal cortex. Group
classifications based on FA or DR was significant and better than classifications based on GM volume. Using either DR or FA
together with GM volume improved classification accuracy. Furthermore, all three measures, FA, DR and GM volume were
similarly accurate in predicting cognitive performance in MCI patients. Taken together, the results imply that DTI measures
are as accurate as measures of GM volume in detecting brain alterations that are associated with cognitive impairment.
Furthermore, a combination of DTI and structural MRI measurements improves classification accuracy.
Citation: Zhang Y, Schuff N, Camacho M, Chao LL, Fletcher TP, et al. (2013) MRI Markers for Mild Cognitive Impairment: Comparisons between White Matter
Integrity and Gray Matter Volume Measurements. PLoS ONE 8(6): e66367. doi:10.1371/journal.pone.0066367
Editor: Yong Fan, Institution of Automation, CAS, China
Received October 10, 2012; Accepted May 7, 2013; Published June 6, 2013
Copyright: 2013 Zhang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This research was supported in part by National Institutes of Health/National Institute on Aging grants P01AG19724, P01AG12435 and grant from NIH/
National Center for Research Resources P41RR023953, which were administered by the Northern California Institute for Research and Education, and with
resources of the Veterans Affairs Medical Center, San Francisco, California. The funders had no role in study design, data collection and analysis, decision to
publish, or preparation of the manuscript.
Competing Interests: The authors have declared that no competing interests exist.
* E-mail: [email protected]
Introduction
attracted considerable interest, in part because of the low risk, noninvasive examinations and broad availability nowadays.
Aside from MRI measurements of brain atrophy [614],
diffusion tensor imaging (DTI), a variant of MRI that is sensitive
to microstructural tissue changes, has also shown promise in
characterizing MCI [1524]. DTI studies in MCI generally rely
on diffusion summary measures, such as fractional anisotropy (FA),
an index of the spatial directionality of tissue water diffusion and
mean diffusivity (MD), a measure of the diffusion magnitude.
These measures uniquely detect abnormalities in white matter
(WM) regions. In particular, FA alterations in MCI though often
subtle [15,17,2528] have been found in limbic circuits
(including posterior and parahippocampal tracts and fornix)
[19,2932] and posterior callosal areas [33,34]. To gain further
insight into WM alterations in MCI, DTI studies have computed
55
subjects had predominant memory impairments) and nonamnestic MCI (naMCI, if the MCI subjects had predominant
non-memory impairments, such as executive, language, visual,
behavioral, motor, or other deficits). Sixty-six cognitively normal
subjects, whose neuropsychological scores were no worse than one
standard deviation below the mean on standardized tests, but also
not better than 2 standard deviations above the mean were
included as controls. Exclusion criteria for all participants included
any poorly controlled illness, use of medication or recreational
drugs that could affect brain function, a history of brain trauma,
brain surgery, ischemic events, or skull defects.
MRI Acquisitions
All scans were preformed on a 4 Tesla (Bruker/Siemens) MRI
system with a single housing birdcage transmit and 8-channel
receiver head coil. T1-weighted images were obtained using a 3D
volumetric magnetization prepared rapid gradient echo
(MPRAGE) sequence with TR/TE/TI = 2300/3/950 ms, 7degree flip angle, 1.061.061.0 mm3 resolution, 157 continuous
sagittal slices. In addition, FLAIR (fluid attenuated inversion
recovery) images (TR/TE/TI = 5000/355/1900 ms) and T2weighted images (TR/TE = 4000/30 ms) were acquired for
clinical reading. DTI was acquired based on a dual-refocused
spin-echo EPI sequence supplemented with twofold parallel
imaging acceleration (GRAPPA) [43] to reduce susceptibility
distortions. Other DTI parameters were TR/TE = 6000/77 ms,
field of view 2566224 cm, 1286112 matrix size, yielding
262 mm2 in-plane resolution, and 40 contiguous slices each
3 mm thick. One reference image (b = 0) and six diffusionweighted images (b = 800 s/mm2, along 6 noncollinear directions)
were acquired.
56
Statistics
Differences between patients and controls were tested using a ttest for each ROI and MRI modality after effects of age and
gender were linearly regressed out. For the rest of the analysis, FA,
DR, and GM volume measures were centered and normalized to
their respective standard deviations (i.e. transformed to z-scores) to
eliminate scaling differences across the ROIs and MRI modalities.
In addition, a non-parametric test using permutations of the
Figure 1. GM and WM parcellations. A. Automated parcellation of 21 cortical and subcortical ROIs for GM measurement, performed by Freesurfer
software. B. Automated parcellation of 20 deep WM ROIs for DTI measurement, performed by SPM8.
doi:10.1371/journal.pone.0066367.g001
57
Control
MCI subtypes
MCI
naMCI
Subject Number
66
51
Age (years)
67.2610.0
72.868.7
Sex
31 M : 35 F
Years of Education*
16.762.4
MMSE*
29.461.0
27.861.9
28.264.1
19.166.8
7.1261.5
4.5262.7
7.4561.3
4.6562.4
14.665.6
13.466.0
21.666.1
15.766.2
aMCI
17
31
71.9610.9
73.167.5
30 M : 21 F
9M:8F
20 M : 11 F
18.6612.4
22.1621.5
17.562.4
28.561.7
27.561.9
21.666.0
17.166.8
5.7661.9
3.9362.8
5.9461.9
3.8762.5
`
`
`
`
ICV (cm )
10766125
10916147
Total GM/ICV
0.4060.04
0.3760.04
Total WM/ICV
0.4160.05
0.3960.05
14.365.9
12.966.1
18.967.0
13.965.3
10466111
11086164
0.3860.03
0.3660.04
0.3960.05
0.3960.05
*6 subjects years of education was missing. 4 subjects MMSE was missing. 28 subjects verbal fluency and semantic fluency were missing. Note, smaller scores of
neurocognitive measures indicate greater impairment.
Significance of group differences between paired groups (MCI vs. Control, aMCI vs. naMCI):
{
0.05,p#0.001,
`
p,0.001.
doi:10.1371/journal.pone.0066367.t001
diagnostic labels was applied to determine whether the heterogeneity of MCI has to be taken into account for findings of
differences between uniform MCI and CN. The concept of false
discovery rate (FDR) was used to account for the multiple
comparisons problem [51]. An adjusted level of p,0.05 was
selected as threshold of significance. The classification accuracy of
FA, DR, and GM volume was tested using regularized logistic
regressions with all respective ROIs included simultaneously as
factors. An L2-norm regularization term (Ridge regression) was
used to control for overfitting [52]. The classification accuracy of
FA, DR, and GM were then compared pairwise by testing
differences in the area under the curve of an operator characteristics analysis (ROC) based on bootstraping (using the roc.test
function in the pROC package of R, URL: http://www.r-project.
org/index.html). Similarly, to test the accuracy of FA, DR, and
GM volume in predicting neurocognitive performance, regularized linear regressions were used with cognitive scores as outcomes
and with all respective ROIs included simultaneously as factors.
Accuracy was expressed as the root mean square error (RMSE)
between predicted and clinical values, where a smaller RMSE
indicating a more accurate prediction. Differences in RMSE
distributions of FA, DR, and GM volume where compared using ttests. Lastly, classifications and predictions were performed with 5fold cross-validation and bootstrapping to estimate reliability in
terms of 95% confidence intervals (CI).
Results
Demographic and Clinical Characteristics
58
59
PLOS ONE | www.plosone.org
Figure 2. Mean differences of the DTI and GM measures between MCI and control. A. Mean differences and standard errors of regional DTI
and GM measures, expressed as Z-scores, between MCI patients and controls for each ROI. Abbreviations:CC = corpus callosum; Post. CG = posterior
cingulum; Isth. CG = isthmus cingulum; FX-ST = fornix (cres) and stria terminalis; IFO = inferior fronto-occipital fasciculus; SLF = superior
longitudinal fasciculus; ILF = inferior longitudinal fasciculus. B. Mean differences and standard errors of regional DTI and MRI measures between
aMCI, naMCI group and controls for each ROI. DTI and GM measures in ROIs with significant heterogeneities between aMCI and naMCI groups were
labeled as #.
doi:10.1371/journal.pone.0066367.g002
Discussion
Group classifications
Classification accuracy of MCI and control subjects based on
either FA, DR, GM volume alone or in combinations is
summarized in Table 2. The corresponding receiver operating
characteristic curves of each classification are illustrated in
Figure 4. The table indicates that classifications based on FA or
DR alone were significant (lower bound of the confidence interval
.50%) in contrast to classifications based on GM volumes that
were not better than chance. However, using either DR or FA
together with GM volume further improved accuracy.
60
Figure 3. The correlations between DTI and GM volumes in MCI patients. The p value maps of the Pearsons correlation between DTI values
and GM volumes. The green and warmer colors indicate significant correlations.
doi:10.1371/journal.pone.0066367.g003
Table 2. Group classifications based on either DTI or GM volume measures separately or used together.
Measure
median
(95% CI Upper)
FA
70.6
77.2
74.4
82.8
51.5
70.8
90.0
DR
66.7
78.8
73.5
83.5
58.0
78.5
95.0
GM volume
64.7
83.3
75.2
79.1
49.5
67.3
86.2
FA + GM volume
88.2
90.9
89.7
94.9
63.9
82.3a
96.2
DR + GM volume
82.3
86.4
84.6
92.6
61.9
81.2b
96.2
AUC area under a receiver operating characteristic curve. AUCs were tested using bootstrap (2000 boots):
Differences between modalities (i.e. FA + GM volume vs. GM volume) were significant at 95% confidence interval.
b
Differences between modalities (i.e. DR + GM volume vs. GM volume) were significant at 90% confidence interval.
doi:10.1371/journal.pone.0066367.t002
a
61
Figure 5. Accuracies of predicting cognitive scores based on DTI and GM volume measures. Distribution of root mean square errors in
predicting MMSE and CVLT based on FA, DR or GM volume measures.
doi:10.1371/journal.pone.0066367.g005
62
Acknowledgments
The authors thank all the participants in this study. We also thank Mr.
Matthew Nordstrom, Mr. Frank Ezekiel for the assistance with image
processing, Mr. Derek Flenniken for clinical data organizing, and Mr.
Philip Insel for the help with statistics.
Author Contributions
Conceived and designed the experiments: YZ NS MWW. Performed the
experiments: YZ NS MC LLC. Analyzed the data: YZ NS MC LLC.
Contributed reagents/materials/analysis tools: NS LLC TPF KY SCW
CM HJR BLM MWW. Wrote the paper: YZ NS. Designed the software
used in analysis: TPF. Results interpretation and manuscript editing: YZ
NS LLC KY SCW CM HJR BLM MWW.
References
1. Albert MS, DeKosky ST, Dickson D, Dubois B, Feldman HH, et al. (2011) The
diagnosis of mild cognitive impairment due to Alzheimers disease: recommendations from the National Institute on Aging-Alzheimers Association workgroups on diagnostic guidelines for Alzheimers disease. Alzheimers Dement 7:
270279.
2. Ramani A, Jensen JH, Helpern JA (2006) Quantitative MR imaging in
Alzheimer disease. Radiology 241: 2644.
3. Frisoni GB, Padovani A, Wahlund LO (2003) The diagnosis of Alzheimer
disease before it is Alzheimer dementia. Arch Neurol 60: 1023; author reply
10231024.
4. Vemuri P, Wiste HJ, Weigand SD, Shaw LM, Trojanowski JQ, et al. (2009)
MRI and CSF biomarkers in normal, MCI, and AD subjects: diagnostic
discrimination and cognitive correlations. Neurology 73: 287293.
5. Sluimer JD, Bouwman FH, Vrenken H, Blankenstein MA, Barkhof F, et al.
(2010) Whole-brain atrophy rate and CSF biomarker levels in MCI and AD: a
longitudinal study. Neurobiol Aging 31: 758764.
6. Jack CR, Jr., Petersen RC, Xu YC, O9Brien PC, Smith GE, et al. (1999)
Prediction of AD with MRI-based hippocampal volume in mild cognitive
impairment. Neurology 52: 13971403.
7. Du AT, Schuff N, Amend D, Laakso MP, Hsu YY, et al. (2001) Magnetic
resonance imaging of the entorhinal cortex and hippocampus in mild cognitive
impairment and Alzheimers disease. J Neurol Neurosurg Psychiatry 71: 441
447.
8. Kantarci K, Xu Y, Shiung MM, O9Brien PC, Cha RH, et al. (2002)
Comparative diagnostic utility of different MR modalities in mild cognitive
impairment and Alzheimers disease. Dement Geriatr Cogn Disord 14: 198
207.
9. Grundman M, Jack CR, Jr., Petersen RC, Kim HT, Taylor C, et al. (2003)
Hippocampal volume is associated with memory but not monmemory cognitive
performance in patients with mild cognitive impairment. J Mol Neurosci 20:
241248.
10. Pennanen C, Kivipelto M, Tuomainen S, Hartikainen P, Hanninen T, et al.
(2004) Hippocampus and entorhinal cortex in mild cognitive impairment and
early AD. Neurobiol Aging 25: 303310.
11. Whitwell JL, Petersen RC, Negash S, Weigand SD, Kantarci K, et al. (2007)
Patterns of atrophy differ among specific subtypes of mild cognitive impairment.
Arch Neurol 64: 11301138.
12. Morra JH, Tu Z, Apostolova LG, Green AE, Avedissian C, et al. (2009)
Automated 3D mapping of hippocampal atrophy and its clinical correlates in
400 subjects with Alzheimers disease, mild cognitive impairment, and elderly
controls. Hum Brain Mapp 30: 27662788.
13. Petersen RC (2009) Early diagnosis of Alzheimers disease: is MCI too late? Curr
Alzheimer Res 6: 324330.
14. Frisoni GB, Fox NC, Jack CR, Jr., Scheltens P, Thompson PM (2010) The
clinical use of structural MRI in Alzheimer disease. Nat Rev Neurol 6: 6777.
15. Rogalski EJ, Murphy CM, deToledo-Morrell L, Shah RC, Moseley ME, et al.
(2009) Changes in parahippocampal white matter integrity in amnestic mild
cognitive impairment: a diffusion tensor imaging study. Behav Neurol 21: 5161.
16. Kiuchi K, Morikawa M, Taoka T, Nagashima T, Yamauchi T, et al. (2009)
Abnormalities of the uncinate fasciculus and posterior cingulate fasciculus in
mild cognitive impairment and early Alzheimers disease: a diffusion tensor
tractography study. Brain Res 1287: 184191.
17. Mielke MM, Kozauer NA, Chan KC, George M, Toroney J, et al. (2009)
Regionally-specific diffusion tensor imaging in mild cognitive impairment and
Alzheimers disease. Neuroimage 46: 4755.
18. Stahl R, Dietrich O, Teipel SJ, Hampel H, Reiser MF, et al. (2007) White
matter damage in Alzheimer disease and mild cognitive impairment: assessment
19.
20.
21.
22.
23.
24.
25.
26.
27.
28.
29.
30.
31.
32.
33.
34.
63
35. Song SK, Yoshino J, Le TQ, Lin SJ, Sun SW, et al. (2005) Demyelination
increases radial diffusivity in corpus callosum of mouse brain. Neuroimage 26:
132140.
36. Acosta-Cabronero J, Williams GB, Pengas G, Nestor PJ (2010) Absolute
diffusivities define the landscape of white matter degeneration in Alzheimers
disease. Brain 133: 529539.
37. Wang L, Goldstein FC, Veledar E, Levey AI, Lah JJ, et al. (2009) Alterations in
cortical thickness and white matter integrity in mild cognitive impairment
measured by whole-brain cortical thickness mapping and diffusion tensor
imaging. AJNR Am J Neuroradiol 30: 893899.
38. Cui Y, Wen W, Lipnicki DM, Beg MF, Jin JS, et al. (2011) Automated detection
of amnestic mild cognitive impairment in community-dwelling elderly adults: a
combined spatial atrophy and white matter alteration approach. Neuroimage
59: 12091217.
39. Folstein MF, Folstein SE, McHugh PR (1975) Mini-mental state. A practical
method for grading the cognitive state of patients for the clinician. J Psychiatr
Res 12: 189198.
40. Morris JC (1993) The Clinical Dementia Rating (CDR): current version and
scoring rules. Neurology 43: 24122414.
41. Delis DC, Kaplan E, Kramer JH (2001) Delis-Kaplan Executive Function
System. San Antonio: The Psychological Corporation.
42. Delis DC, Kramer JH, Kaplan E, Ober BA (2000) California Verbal Learning
Test, 2nd Edition. San Antonio, TX: The Psychological Corporation.
43. Griswold MA, Jakob PM, Heidemann RM, Nittka M, Jellus V, et al. (2002)
Generalized autocalibrating partially parallel acquisitions (GRAPPA). Magn
Reson Med 47: 12021210.
44. Fischl B, Salat DH, Busa E, Albert M, Dieterich M, et al. (2002) Whole brain
segmentation: automated labeling of neuroanatomical structures in the human
brain. Neuron 33: 341355.
45. Fischl B, van der Kouwe A, Destrieux C, Halgren E, Segonne F, et al. (2004)
Automatically parcellating the human cerebral cortex. Cereb Cortex 14: 1122.
46. Callen DJ, Black SE, Gao F, Caldwell CB, Szalai JP (2001) Beyond the
hippocampus: MRI volumetry confirms widespread limbic atrophy in AD.
Neurology 57: 16691674.
47. Jones BF, Barnes J, Uylings HB, Fox NC, Frost C, et al. (2006) Differential
regional atrophy of the cingulate gyrus in Alzheimer disease: a volumetric MRI
study. Cereb Cortex 16: 17011708.
48. Jacobs HI, Van Boxtel MP, Jolles J, Verhey FR, Uylings HB (2011) Parietal
cortex matters in Alzheimers disease: an overview of structural, functional and
metabolic findings. Neurosci Biobehav Rev 36: 297309.
49. Tao R, Fletcher PT, Gerber S, Whitaker RT (2009) A variational image-based
approach to the correction of susceptibility artifacts in the alignment of diffusion
weighted and structural MRI. Inf Process Med Imaging 21: 664675.
50. Chua TC, Wen W, Slavin MJ, Sachdev PS (2008) Diffusion tensor imaging in
mild cognitive impairment and Alzheimers disease: a review. Curr Opin Neurol
21: 8392.
51. Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I (2001) Controlling the false
discovery rate in behavior genetics research. Behavioural Brain Research 125:
279284.
52. Hastie T, Tibshirani RJ, Friedman J (2009) The Elements of Statistical
Learning: Data Mining, Inference, and Prediction.: Springer.
53. Clerx L, Visser PJ, Verhey F, Aalten P (2012) New MRI markers for Alzheimers
disease: a meta-analysis of diffusion tensor imaging and a comparison with
medial temporal lobe measurements. J Alzheimers Dis 29: 405429.
54. Cho H, Yang DW, Shon YM, Kim BS, Kim YI, et al. (2008) Abnormal integrity
of corticocortical tracts in mild cognitive impairment: a diffusion tensor imaging
study. J Korean Med Sci 23: 477483.
55. Ringman JM, O9Neill J, Geschwind D, Medina L, Apostolova LG, et al. (2007)
Diffusion tensor imaging in preclinical and presymptomatic carriers of familial
Alzheimers disease mutations. Brain 130: 17671776.
56. Chen TF, Lin CC, Chen YF, Liu HM, Hua MS, et al. (2009) Diffusion tensor
changes in patients with amnesic mild cognitive impairment and various
dementias. Psychiatry Res 173: 1521.
57. Bai F, Zhang Z, Watson DR, Yu H, Shi Y, et al. (2009) Abnormal integrity of
association fiber tracts in amnestic mild cognitive impairment. J Neurol Sci 278:
102106.
58. Rose SE, McMahon KL, Janke AL, O9Dowd B, de Zubicaray G, et al. (2006)
Diffusion indices on magnetic resonance imaging and neuropsychological
59.
60.
61.
62.
63.
64.
65.
66.
67.
68.
69.
70.
71.
72.
73.
74.
75.
76.
77.
78.
79.
64
PLOS ONE | www.plosone.org
10
Neuroimaging
biomarkers:
PET
65
Abstract
The aim of this study is to identify mild cognitive impairment (MCI) due to Alzheimers disease (AD) using amyloid imaging
of beta amyloid (Ab) deposition and FDG imaging of reflecting neuronal dysfunction as PET biomarkers. Sixty-eight MCI
patients underwent cognitive testing, [11C]-PIB PET and [18F]-FDG PET at baseline and follow-up. Regions of interest were
defined on co-registered MRI. PIB distribution volume ratio (DVR) was calculated using Logan graphical analysis, and the
standardized uptake value ratio (SUVR) on the same regions was used as quantitative analysis for [18F]-FDG. Thirty (44.1%)
of all 68 MCI patients converted to AD over 19.267.1 months. The annual rate of MCI conversion was 23.4%. A positive Ab
PET biomarker significantly identified MCI due to AD in individual MCI subjects with a sensitivity (SS) of 96.6% and specificity
(SP) of 42.1%. The positive predictive value (PPV) was 56.8%. A positive Ab biomarker in APOE e4/4 carriers distinguished
with a SS of 100%. In individual MCI subjects who had a prominent impairment in episodic memory and aged older than 75
years, an Ab biomarker identified MCI due to AD with a greater SS of 100%, SP of 66.6% and PPV of 80%, compared to FDG
biomarker alone or both PET biomarkers combined. In contrast, when assessed in precuneus, both Ab and FDG biomarkers
had the greatest level of certainty for MCI due to AD with a PPV of 87.8%. The Ab PET biomarker primarily defines MCI due
to AD in individual MCI subjects. Furthermore, combined FDG biomarker in a cortical region of precuneus provides an
added diagnostic value in predicting AD over a short period.
Citation: Hatashita S, Yamasaki H (2013) Diagnosed Mild Cognitive Impairment Due to Alzheimers Disease with PET Biomarkers of Beta Amyloid and Neuronal
Dysfunction. PLoS ONE 8(6): e66877. doi:10.1371/journal.pone.0066877
Editor: Stephen D Ginsberg, Nathan Kline Institute and New York University School of Medicine, United States of America
Received February 8, 2013; Accepted May 10, 2013; Published June 14, 2013
Copyright: 2013 Hatashita and Yamasaki. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which
permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The authors have no support or funding to report.
Competing Interests: GE Healthcare Ltd. (Little Chalfont, Buckinghamshire, HP7 9NA, United Kingdom) entered into the sublicense agreement with the
Shonan-Atsugi Hospital to allow the use of the [11C]-PIB PET imaging. There are no further patents, products in development or marketed products to declare.
This does not alter the authors adherence to all the PLOS ONE policies on sharing data and materials, as detailed online in the guide for authors.
* E-mail: [email protected]
Introduction
Mild cognitive impairment (MCI) has a multitude of causes,
including predementia of Alzheimers disease (AD) and non-AD as
well as depression and various physical disorders. A recent study
has shown that subjects with MCI are at increased risk of
developing AD and that their overall rate of progression to AD is
typically 10%215% per year [1]. Identifying subjects with MCI
who are most likely to decline in cognition over time is a major
focus in AD research. Revised research criteria for the diagnosis of
AD have been proposed, and a framework has been developed to
capture the earliest stage of AD [2]. Recently, the National
Institute on Aging (NIA) - Alzheimers Association working group
has proposed the diagnostic criteria for the symptomatic predementia phase of AD, referred to a MCI due to AD [3]. The
MCI due to AD could be used to identify individuals with AD
pathophysiological processes as the primary cause of their
progressive cognitive dysfunction.
The clinical research diagnosis of MCI due to AD is based on
the Core Clinical Criteria used in all clinical settings without the
need of highly specialized tests and the Clinical Research Criteria
with the use of biomarkers. Some biomarker measurements
directly reflect the pathology of AD, including positron emission
tomography (PET) amyloid imaging and cerebrospinal fluid (CSF)
66
Image analysis
All subjects underwent T1-weighed MRI (1.5 T) for coregistration with the PET images. MRI-based correction of PET
data for partial volume effects was carried out using the PMOD
software (PMOD Technologies Ltd., Adliswil, Switzerland).
Regions of interest were manually drawn on the co-registered
MR image, including the following cortical regions: lateral
temporal cortex, medial temporal cortex, frontal cortex, occipital
cortex, parietal cortex, sensory motor cortex, anterior cingulate
gyrus, posterior cingulate gyrus, precuneus cortex and cerebellar
cortex.
Data management
Levels of PIB retention were determined by the distribution
volume ratio (DVR) with Logan graphical analysis for 35 to 60
minutes with cerebellar gray matter as reference [14]. Regional
DVR values of each cortical region and the cortical DVR value of
the whole cortical regions were calculated. DVR PET images were
created for visual inspection with a rainbow color scale. Scans of
all subjects were visually assessed with 2 readers as amyloidpositive or amyloid-negative for cortical PIB retention. Amyloidpositive images showed greatly increased PIB retention in the
cortex compared to white matter and a typical PIB pattern of
cortical amyloid deposition, whereas amyloid-negative images
showed no PIB retention in any of the cortical regions.
The same co-registration method was applied for quantification
of [18F]-FDG. A standardized uptake value (SUV) of the same
region was obtained and subsequently normalized to the cerebellar
cortex as reference. Glucose metabolism was referred to as the
SUV ratio (SUVR).
Statistical analysis
Group differences were evaluated with analysis of variance
(ANOVA), followed by Bonferroni post hoc tests to assess the
significance. Analysis correlation between DVR values, FDG
SUVR values, age, MMSE scores, CDR SB scores and WMS-R
recall scores yielded Pearsons product moment correlation
coefficient (r). Categorical variables were examined with Fishers
exact test. Paired t-tests were used to study changes between
baseline and follow-up data. Results were considered significant at
p,0.05. Data are presented as means 6 standard deviations (SD).
Statistical analyses were performed with Statcel 3 software (OMS
Inc. Japan).
Results
PET imaging
All PET scans were performed on the same day as the cognitive
testing, using a Siemens ECAT ACCEL scanner in 3-dimensional
scanning mode, providing 63 contiguous 2.46-mm slices with a
5.6-mm transaxial and a 5.4-mm axial resolution. All imaging data
were reconstructed into a 1286128 matrix. The amyloid PET
imaging was accomplished with the radiotracer [11C]-PIB, as
previously described in detail [9]. Briefly, [11C]-PIB was produced
in our PET center using the one-step [11C] methyl triflate
approach. [11C]-PIB was injected intravenously as a bolus with a
mean dose of 561.5611.2 MBq (n = 68). Dynamic PET scanning
was performed for 60 minutes according to a predetermined
protocol. Sixty minutes following the completion of [11C]-PIB
image, the subjects were injected intravenously with 249.9628.8
MBq (n = 68) of [18F]-FDG and remained in a dark, quiet room.
PLOS ONE | www.plosone.org
67
Ab deposition
Fifty-one (75%) of the 68 MCI patients had marked PIB
retention in the frontal, parietal and lateral temporal cortical
regions as well as the cingulate gyrus and precuneus, showing the
typical AD pattern at baseline (Figure 1). Of 51 MCI patients who
had Ab deposition, 29 (56.8%) patients converted to AD during
the follow-up period while the remaining 22 patients did not
convert despite similar Ab deposition. In contrast, 16 of 17 MCI
patients who had no PIB retention in any cortical region were
stable and only 1 MCI patient converted to AD.
Cortical PIB DVR values for MCI converters and stable
patients are presented in Figure 2. At baseline, the mean value of
cortical DVR in amyloid-positive converters (2.0260.30, n = 29,
p,0.01) was significantly higher than that in amyloid-negative
stable patients (1.2660.16, n = 16). In the amyloid-positive stable
patients, in addition, the mean cortical DVR value (2.0560.37,
n = 22) was similar to amyloid-positive converters. Also, there were
no significant differences in the regional DVR values of any
cortical region between amyloid-positive converters and stable
patients. The mean value of cortical DVR in amyloid-positive
converters at follow-up (2.1160.37, n = 29) was not significantly
different from baseline (Figure 2). There were no significant
changes in cortical DVR values between baseline and follow-up in
amyloid-positive or amyloidnegative stable patients.
Glucose metabolism
Fifty-seven (83.8%) of 68 MCI patients had reduced glucose
metabolism (FDG SUVR #0.99) in whole cortical regions at
baseline, 28 patients (49.1%) of whom progressed to AD during
the follow-up period. A mean FDG SUVR value for the whole
cortical regions in converters was 0.9360.06 (n = 30, p,0.01) at
baseline, significantly different from stable patients (0.9860.06,
n = 38). In cortical regions of the lateral temporal cortex, parietal
cortex and precuneus, in particular, MCI converters had
significantly lower regional FDG SUVR values at baseline
compared to stable patients (Figure 3). At follow-up, MCI
converters did not have significant decreases in mean values of
FDG SUVR in whole cortical regions or each cortical region.
Baseline
68
Follow up
Characteristic
Converter
Stable
Converter
Stable
study popu
30
38
30
38
Education (yr)
12.262.4
12.062.0
12.262.4
12.062.0
MMSE
26.561.5
27.361.6
22.862.5*
27.361.6
global CDR
0.5
0.5
0.660.2
0.5
CDR SB
0.860.2
0.860.3
2.561.0*
0.760.4
Immediate Rec.
4.962.7
6.963.1
4.062.8
7.263.9
Delayed Rec.
1.762.0*
4.663.4
1.562.3
4.964.3
14 (46.6%)
14 (36.8%)
14 (46.6%)
14 (36.8%)
Biomarker
SS
SP
PPV
NPV
aged 5089
Ab PET
96.6
42.1
56.8
94.1
FDG PET
93.3
23.6
49.1
81.8
96.5
21.7
60.8
83.3
100
66.6
80.0
100
FDG PET
91.6
44.4
68.7
80.0
100
60.0
84.6
100
APOE e4 carriers
Ab PET
100
28.5
58.3
100
FDG PET
100
21.4
56.0
100
100
18.1
60.8
100
SS: sensitivity, SP: specificity, PPV: positive predictive value, NPV: negative
predictive value, Ab: b-amyloid, MCI: mild cognitive impairment, APOE:
apolipoprotein E, Data are presented as percentages.
doi:10.1371/journal.pone.0066877.t002
69
Discussion
We demonstrated that 44.1% of 68 MCI patients converted to
AD during a follow-up of 19.267.1 months and the rate of MCI
conversion was 23.4% per year. A recent study has reported that
48% of 31 subjects with MCI converted to AD during a follow-up
period of 2.9 years [15]. The Alzheimers Disease Neuroimaging
Initiative has reported that 32.9% of 85 single-domain or multipledomain amnestic MCI patients converted to AD during a followup period of 1.960.4 years (an annual rate of 17.2%) [1]. Thus,
conversion rates vary greatly depending on the criteria applied, the
nature of the subject population and the period of observation.
In the present study, 50% of APOE e4 carriers with MCI
converted to AD while 40% of APOE e4 non-carriers with MCI
converted. This supports that the progressive value of genetic
assessment alone is comparatively low although a high percentage
of progressive MCI patients carry the APOE e4 allele [18]. In
contrast, this study showed that 83.3% of APOE e4/4 carriers
with MCI converted to AD. These findings indicate that the
presence of the APOE e4/4 in MCI patients is strongly associated
with conversion to AD. The homozygous carriers of the APOE e4
allele could have an increased risk of MCI progression to AD.
The present study demonstrated, using [11C]-PIB PET
imaging, that 56.8% of 51 MCI patients with Ab deposition
converted to AD over 2 years, compared to 5.8% of 17 MCI
patients without Ab deposition. In addition, when MCI patients
with Ab deposition, aged 75289 years, had an impaired episodic
memory, the rate of MCI conversion was 80%. It has previously
been reported that 82% of PIB-positive amnestic subjects with
MCI converted to AD over a 3-year follow-up period [15], and
67% of MCI subjects with high PIB retention converted over 2
years [16]. These findings indicate that a high proportion of MCI
patients with Ab deposition convert to AD within a short period.
In vivo detection of Ab deposition with PET imaging is a useful
prognostic tool for identifying MCI progression to AD.
We found that MCI converters had the regional hypometabolism in cortical regions of lateral temporal cortex, parietal cortex
and precuneus at baseline compared to stable patients. Our results
are consistent with the previous study that regional metabolic
reduction, in temporoparietal or posterior cingulate cortices, was
detected at baseline in MCI patients who converted to AD within
1 year, compared with non-converters [17]. The AD pathophysiological process has been demonstrated to be temporoparietal
and/or precuneus hypometabolism, which is a topographic
pattern characteristic of AD. Therefore, if a regional hypometabolism in the temporal, parietal and/or precuneus cortices is
detected on FDG PET imaging, the progression to AD may occur
within a few years.
An evaluation in predictive accuracy by combining different
biomarkers and neuropsychological variables for the prediction of
AD in MCI has been reported [19]. Inconsistencies in recent
findings are likely due to a variability in neuroimaging processing
techniques, CSF protein immunoassay techniques, setting cutoff
Acknowledgments
GE Healthcare Ltd. (Little Chalfont, Buckinghamshire, HP7 9NA, United
Kingdom) entered into the sublicense agreement with the Shonan-Atsugi
Hospital to allow the use of the [11C]-PIB PET imaging. The authors
thank all the patients and their relatives who have participated in this study,
and the staff at Shonan-Atsugi Hospital.
Author Contributions
Conceived and designed the experiments: SH HY. Performed the
experiments: SH HY. Analyzed the data: SH HY. Contributed
reagents/materials/analysis tools: SH HY. Wrote the paper: SH HY.
References
70
1. Petersen RC, Aise PS, Beckett LA, Donohue MC, Gamst AC, et al. (2010)
Alzheimers Disease Neuroimaging Initiative (ADNI). Neurology 74: 201209.
2. Dubois B, Feldman HH, Jacova C, Dekosky ST, Barberger-Gateauet P, et al.
(2007) Research criteria for the diagnosis of Alzheimers disease: revising the
NINCDS-ADRDA criteria. Lancet Neurol 6: 734746.
3. Albert MS, DeKosky ST, Diskson D, Dubois B, Feldman HH, et al. (2011) The
diagnosis of mild cognitive impairment due to Alzheimers disease: Recommen-
5. Jack CR, Lowe VJ, Weigand SD, Wiste HJ, Senjem ML, et al. (2009) Serial PIB
and MCI in normal, mild cognitive impairment and Alzheimers disease:
implications for sequence of pathologic events in Alzheimers disease. Brain 132:
13551365.
6. Landau SM, Harvey D, Madison CM, Reiman EM, Foster NL, et al. (2010)
Comparing predictors of conversion and decline in mild cognitive impairment.
Neurology 75: 230238.
7. Mathis CA, Wang Y, Holt DP, Huang GF, Debnath ML, et al. (2003) Synthesis
and evaluation of [11C]-labeled 6-substituted 2-arylbenzothiazoles as amyloid
imaging agents. J Med Chem 46: 27402754.
8. Klunk WE, Engler H, Nordberg A, Wang Y, Blomqvist G, et al. (2004) Imaging
brain amyloid in Alzheimer disease with Pittsburgh Compound-B. Ann Neurol
55: 306319.
9. Hatashita S, Yamasaki H (2010) Clinically different stages of Alzheimers disease
associated by amyloid deposition with [11C]-PIB PET imaging. J Alzheimers
Dis 21: 9951003.
10. Folstein MF, Folstein SE, McHugh PR (1975) Mini-mental state A practical
method for grading the cognitive state of patients for the clinician. J Psychiatr
Res 12: 189198.
11. Morris JC (1993) The clinical dementia rating (CDR): current version and
scoring rules. Neurology 43: 24122414.
12. Wechsler DA (1987) Wechsler Memory Scale-Revised. New York: Psychological
Corporation. New York.
13. Mckhann G, Drachman DA, Folstein M, KatzmanR, Price D, et al. (1984)
Clinical diagnosis of Alzheimers disease-report of the NINCDS-ADRDA work
14.
15.
16.
17.
18.
19.
20.
group under the auspices of Department of Health and Human Services Task
Force on Alzheimers disease. Neurology 34: 939944.
Lopresti BJ, Klunk WE, Mathis CA, Hoge JA, Ziolko SK, et al. (2005) Simplified
quantification of Pittsburgh Compound B amyloid imaging PET studies: a
comparative analysis. J Nucl Med 46: 19591972.
Okello A, Koivumen J, Edison P, Acher HA, Turkheimer FE, et al. (2009)
Conversion of amyloid positive and negative MCI to AD over 3 years: An
[11C]-PIB PET study. Neurology 73: 754760.
Villemagne VL, Pike KE, Chetelat G, Ellis KA, Mulligan RS, et al. (2011)
Longitudinal Assessment of AB and cognition in Aging and Alzheimer Disease.
Ann Neurol 69: 181192.
Drzezga A, Lautenschlager N, Siebner H, Riemenschneider M, WiIloch F, et al.
(2003) Cerebral metabolic change accompanying conversion of mild cognitive
impairment into Alzheimers disease: a PET follow-up study. Eur J Nucl Med
Mol Imaging 30: 11041113.
Drzezga A, Grimer T, Riemenschneider M, Lautenschlager N, Siebner H, et al.
(2005) Prediction of individual clinical outcome in MCI by means of genetic
assessment and [18F]-FDG PET. J Nucl Med 46: 16251632.
Ewers M, Walsh C, Trojanowski JQ, Shaw LM, Peterson RC, et al. (2012)
Prediction of conversion from mil cognitive impairment to Alzheimers disease
dementia based upon biomarkers and neuropsychological test performance.
Neurobiol Aging 33: 12031214.
Zhang S, Han D, Tan X, Feng J, Guo Y, et al. (2012) Diagostic accuracy of
[18F]-FDG and [11C]-PIB-PET for prediction of short-term conversion to
Alzheimers disease in subjects with mild cognitive impairment. Int J Clin Pract
66: 185198.
71
PLOS ONE | www.plosone.org
Abstract
Background: To compare the neocortical amyloid loads among cognitively normal (CN), amnestic mild cognitive
impairment (aMCI), and Alzheimers disease (AD) subjects with [18F]AV-45 positron emission tomography (PET).
Materials and Methods: [18F]AV-45 PET was performed in 11 CN, 13 aMCI, and 12 AD subjects to compare the cerebral
cortex-to-whole cerebellum standard uptake value ratios (SUVRs) of global and individual volumes of interest (VOIs) cerebral
cortex. The correlation between global cortical [18F]AV-45 SUVRs and Mini-Mental State Examination (MMSE) scores was
analyzed.
Results: The global cortical [18F]AV-45 SUVRs were significantly different among the CN (1.0860.08), aMCI (1.2760.06), and
AD groups (1.3460.13) (p = 0.0003) with amyloidosis positivity rates of 9%, 62%, and 92% in the three groups respectively.
Compared to CN subjects, AD subjects had higher SUVRs in the global cortical, precuneus, frontal, parietal, occipital,
temporal, and posterior cingulate areas; while aMCI subjects had higher values in the global cortical, precuneus, frontal,
occipital and posterior cingulate areas. There were negative correlations of MMSE scores with SUVRs in the global cortical,
precuneus, frontal, parietal, occipital, temporal, posterior cingulate and anterior cingulate areas on a combined subject pool
of the three groups after age and education attainment adjustment.
Conclusions: Amyloid deposition occurs relatively early in precuneus, frontal and posterior cingulate in aMCI subjects.
Higher [18F]AV-45 accumulation is present in parietal, occipital and temporal gyri in AD subjects compared to the aMCI
group. Significant correlation between MMSE scores and [18F]AV-45 SUVRs can be observed among CN, aMCI and AD
subjects.
Citation: Huang K-L, Lin K-J, Hsiao I-T, Kuo H-C, Hsu W-C, et al. (2013) Regional Amyloid Deposition in Amnestic Mild Cognitive Impairment and Alzheimers
Disease Evaluated by [18F]AV-45 Positron Emission Tomography in Chinese Population. PLoS ONE 8(3): e58974. doi:10.1371/journal.pone.0058974
Editor: Karl Herholz, University of Manchester, United Kingdom
Received November 28, 2012; Accepted February 8, 2013; Published March 14, 2013
Copyright: 2013 Huang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This study was carried out with financial support from the National Science Council, Taiwan (Grant NSC-97-2314-B-182A-050, NSC-98-2314-B-182-056,
NSC-100-2314-B-182-003, NSC-100-2314-B-182A-091-MY3, and NSC-101-2314-B-182A-061-MY2), and Chang Gung Memorial Hospital (Grant CMRPG300071,
CMRPG390793, and CMRPD1A0312). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing Interests: The authors have declared that no competing interests exist.
* E-mail: [email protected] (CCH); [email protected] (TCY)
. These authors contributed equally to this work.
Introduction
72
Methods
Participants
Thirty-six subjects were studied at the Center for Dementia of
the Linkou Chang Gung Memorial Hospital (CGMH), Taiwan,
including 11 CN, 13 aMCI, and 12 AD subjects. All subjects
underwent a neurological examination, neurocognitive evaluation,
and routine blood analysis. Study subjects were categorized on the
basis of the consensus of panels composed of neurologists,
neuropsychologists, neuroradiologists and experts in nuclear
medicine. The CN subjects were volunteers with normal cognitive
performance; the methods for CN subject recruitment have been
described previously [8]. The diagnostic criteria for aMCI were
based on those proposed by Petersen et al [13]: (1) subjective
memory complaints by the patient or an informant, (2) relatively
normal performance in other cognitive domains and (3) normal
activities of daily living, (4) objective memory impairment on at
least one neurocognitive test of memory performance, and (5) no
dementia according to DSM-IV criteria [14]. No rigid cutoff score
was applied to determine objective memory impairment, but
aMCI was generally determined when memory measures fell 1.0
1.5 standard deviations below the means for age-matched norms
in Taiwan. The diagnosis of AD was made when subjects fulfilled
the National Institute of Neurological and Communicative
Disorders and StrokeAlzheimers Disease and Related Disorders
Association (NINCDS-ADRDA) criteria for probable AD [15]. A
comprehensive battery of neurocognitive tests was administered to
obtain objective evidence of cognitive impairment, including the
Mini-Mental State Examination (MMSE) [16], Clinical Dementia
Rating (CDR), logical memory in the Wechsler memory scalerevised (LM), visual-association memory test (VAMT), category
verbal fluency test (CVFT), trail-making A test (TMAT), and
clock-drawing test (CDT) [17]. Prior to enrollment, the study
objectives and protocol were well informed to all subjects and also
PLOS ONE | www.plosone.org
73
2
Both PET and MRI images from each subject were coregistered first, and then MRI images were spatially normalized
into the MNI T1 template as in Statistical Parametric Mapping 5
(SPM5) (Wellcome Department of Cognitive Neurology, University College London, London, UK). The transformation parameters of MR normalization were then applied to corresponding
PET images. All the image data were processed in the PMOD
software (version 3.2; PMOD Technologies Ltd., Zurich, Switzerland). For all subjects, a binary gray matter mask was first created
by taking the non-zero pixels from an averaged gray matter
probability map generated from the segmentation of the spatially
normalized MR images from normal subjects. The final VOI
template for VOI analysis was then obtained by taking the
intersection of both the binary gray matter mask and the
automated anatomic labeling (AAL) atlas (VOI template) for all
subjects [21,22]; this step was performed to minimize the inclusion
of both cerebrospinal fluid and white matter (and thus non-specific
white matter [18F]AV-45 retention) in the statistical measures of all
VOIs [23]. The original AAL atlas which contains 116 regions was
combined into the following 8 labeled VOIs for further analysis:
the precuneus, frontal, parietal, occipital, temporal, posterior
cingulate and anterior cingulate areas, and the whole cerebellum.
The standard uptake value ratio (SUVR) was calculated by
determining the ratios of the integrated activities between the
target VOI (cortical regions) and whole cerebellum reference VOI.
Voxel-based statistical analysis was conducted with SPM5 [24].
The spatial normalization of individual [18F]AV-45 images to
MNI MRI T1-template space was performed in the same way as
that for the VOI analysis described above. All voxels in the
normalized [18F]AV-45 images were divided by the mean
[18F]AV-45 uptake of the whole cerebellum VOI in each subject
to form SUVR images. Moreover, all PET images in SPM analysis
were smoothed by 8 mm FWHM Gaussian filtering. Voxel-wise
[18F]AV-45 uptake differences among the CN, aMCI, and AD
subjects were assessed in SPM5. Statistical maps displaying group
differences were shown at a significance value of p,0.01 and
extent voxels of 100. Due to a relatively small sample size of our
study, and to increase the sensitivity in the SPM analysis, the lower
threshold of two-tailed p,0.01 under uncorrected statistics was
applied here [25].
Results
The study recruited 11, 13, and 12 subjects in the CN, aMCI,
and AD groups, respectively. The demographic data were shown
in Table 1. None of the subjects had adverse events directly related
to [18F]AV-45 administration.
VOI Analysis
The VOI analysis of [18F]AV-45 images revealed significant
differences in the global cortical SUVR among subjects with CN,
aMCI and AD (Figure 1). Regarding to the differences in regional
SUVRs, the precuneus, frontal, parietal, occipital, temporal, and
posterior and anterior cingulate cortical SUVRs were also
different among the three groups. In the post hoc analyses, all of
these areas except the anterior cingulate had a higher SUVR of
[18F]AV-45 in AD subjects than CN subjects. Furthermore, the
SUVRs of the global cortical, precuneus, frontal, occipital and
posterior cingulate areas were higher in aMCI subjects than CN
subjects. There was no difference in [18F]AV-45 SUVRs of these
predefined areas between aMCI and AD subjects. With the SUVR
cutoff point of 1.178, the positive rate of cerebral cortical
amyloidosis was 8%, 62%, and 92% for CN, aMCI and AD
groups, respectively. Interestingly, there were varying degrees of
[18F]AV-45 retention in CN, aMCI negative for cerebral
amyloidosis (aMCI (2)), aMCI positive for cerebral amyloidosis
(aMCI (+)), and AD subjects, whose patterns closely followed the
amyloid deposition distribution observed in post mortem brain
tissue (Figure 2) [26].
SPM Analysis
The SPM analysis showed that aMCI subjects had significantly
higher [18F]AV-45 uptake than CN subjects in the frontal,
temporal, parietal, occipital areas, as well as in the precuneus,
which showed the most prominent differences (Figure 3A). These
results were supported by the automated VOI analysis in which
aMCI patients showed 15% increased [18F]AV-45 uptake in the
Statistical Analyses
74
CN (n = 11)
aMCI (n = 13)
AD (n = 12)
Gender (M/F)
6/5
5/8
3/9
0.35
Age, years
69.3 (6.6)
70.0 (11.2)
75.8 (5.7)
0.15
CDR
0.0 (0.0)
0.5 (0.0)*
1.1 (0.3)*
,0.0001
MMSE
27.9 (1.6)
19.7 (2.4)*
12.9 (3.7)*
,0.0001
3.0 (1.2)
5.9 (2.2)
,0.001
6.2 (4.1)
7.2 (4.3)
0.19
Figure 1. Scatter plots showing general and individual region-to-cerebellum ratios between different groups. Comparisons of [18F]AV45 SUVR in the predefined VOIs of the global cortical (A), precuneus (B), frontal (C), parietal (D), occipital (E), temporal (F), posterior cingulate (G) and
anterior cingulate (H) areas among the cognitively normal (CN), amnestic mild cognitive impairment (aMCI), and Alzheimers disease (AD) subjects. (*
p,0.05, ** p,0.005, *** p,0.0005, when comparing SUVR among CN, aMCI and AD subjects with the KruskalWallis ANOVA test. { p,0.05 versus CN
subjects with Dunns post hoc analysis).
doi:10.1371/journal.pone.0058974.g001
Figure 2. Representative sagittal [18F]AV-45 PET images. Varying degrees of [18F]AV-45 PET uptake in CN, aMCI negative for cerebral
amyloidosis (aMCI (-)), aMCI positive for cerebral amyloidosis (aMCI (+)), and AD subjects (bottom row) were closely corresponding to the pathological
amyloid deposition distribution adapted from Braak and Braak (top row) [26].
doi:10.1371/journal.pone.0058974.g002
75
Figure 3. Statistical parametric mapping analysis: Localization of increased [18F]AV-45 retention between CN, aMCI and AD
subjects. Comparisons of [18F]AV-45 SUVRs between cognitively normal (CN) and amnestic mild cognitively impairment (aMCI) subjects (A), and
between aMCI and Alzheimers disease (AD) subjects (B) (Uncorrected for multiple comparisons and the color bar values indicate the value of the Tstatistic in each display). Surface rendering was used to illustrate the cortical areas where [18F]AV-45 SUVRs were increased in aMCI than CN subjects
(red) and increased in AD than aMCI subjects (green) (C).
doi:10.1371/journal.pone.0058974.g003
fluorine-18-labeled radiotracers, including [18F]AV-45 (florbetapir), [18F]AV-1 (florbetaben) and [18F]GE067 (flutemetamol), have
been extensively employed for amyloid plaque detection both in
vivo and ex vivo [29]. Although [11C] PIB PET has been applied
to compare the amyloid plaque distribution among CN, MCI and
AD subjects in the Asian population [30], there are limited reports
about fluorine-18-labeled amyloid PET studies in non-Caucasians
[20,31,32]. In a recent report by Barthel et al, the Caucasian
elderly normal subjects had higher SUVRs in different neocortical
Discussion
Although the cause of AD is still unknown, the presence of
amyloid plaques plays a pivotal role in the pathogenesis of the
disease. Studies show that the presence of amyloid plaques in
patients with mild cognitive impairment may increase the potential
of conversion to AD [27]. Therefore, the development of imaging
techniques to detect amyloid plaques in vivo would be helpful for
early diagnosis and monitoring disease progress [28]. Several
76
Figure 4. Relationship between MMSE scores and global and regional cortical [18F]AV-45 retention. Results of the negative regressions
between MMSE scores and global and regional cortical [18F]AV-45 SUVRs among the cognitive normal, amnestic mild cognitive impairment, and
Alzheimers disease subjects with age and years of education as covariates.
doi:10.1371/journal.pone.0058974.g004
77
6
Conclusions
The findings from our study show that [18F]AV-45 PET can be
applied in vivo to differentiate the cerebral amyloid burden among
CN, aMCI, and AD subjects using continuous or binary
measurements, and [18F]AV-45 SUVR is associated with MMSE
scores. The clinical diagnosis-based threshold for cerebral
amyloidosis derived from AD patients is similar to the pathology-based threshold. Therefore, [18F]AV-45 PET is a potential
biomarker for AD in multicenter monitoring and treatment trials.
Acknowledgments
Supporting Information
We would like to thank Ko-Ting Chao and Yu-Chen Hsieh for their
technical assistance. We also thank Avid Radiopharmaceuticals, Inc.
(Philadelphia, PA, USA) for providing the precursor for the preparation of
[18F]AV-45.
Author Contributions
Conceived and designed the experiments: KLH MPK TCY CCH.
Performed the experiments: KLH KJL HCK WCH WLC TCY CCH.
Analyzed the data: KLH KJL ITH MPK SPW CJH YYW TCY CCH.
Contributed reagents/materials/analysis tools: KLH KJH ITH TCY
CCH. Wrote the paper: KLH KJL MPK YYW TCY CCH.
Comparing [18F]AV-45 uptake between amnestic mild cognitive impairment (aMCI) patients and
cognitively normal (CN) subjects. The locations and values
Table S2
References
78
18. Yao CH, Lin KJ, Weng CC, Hsiao IT, Ting YS, et al. (2010) GMP-compliant
automated synthesis of [(18)F]AV-45 (Florbetapir F 18) for imaging betaamyloid plaques in human brain. Appl Radiat Isot 68: 22932297.
19. Wong DF, Rosenberg PB, Zhou Y, Kumar A, Raymont V, et al. (2010) In vivo
imaging of amyloid deposition in Alzheimer disease using the radioligand 18FAV-45 (Florbetapir F 18). J Nucl Med 51: 913920.
20. Hsiao IT, Huang CC, Hsieh CJ, Hsu WC, Wey SP, et al. (2012) Correlation of
early-phase 18F-florbetapir (AV-45/Amyvid) PET images to FDG images:
preliminary studies. Eur J Nucl Med Mol Imaging 39: 613620.
21. Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, et al.
(2002) Automated anatomical labeling of activations in SPM using a
macroscopic anatomical parcellation of the MNI MRI single-subject brain.
NeuroImage 15: 273289.
22. Huang KL, Lin KJ, Ho MY, Chang YJ, Chang CH, et al. (2012) Amyloid
deposition after cerebral hypoperfusion: evidenced on [18F]AV-45 positron
emission tomography. J Neurol Sci 319: 124129.
23. Hsieh CJ, Fan L, Lin KJ, Hsiao IT, Huang CC, et al. (2011) Precise F-18
florbetapir (AV-45) PET image spatial normalization using early-phase image
template. J Nucl Med 52: 549.
24. Ashburner J, Friston KJ (2005) Unified segmentation. NeuroImage 26: 839851.
25. Juh R, Kim J, Moon D, Choe B, Suh T (2004) Different metabolic patterns
analysis of Parkinsonism on the 18F-FDG PET. European Journal of Radiology
51: 223233.
26. Braak H, Braak E (1997) Frequency of Stages of Alzheimer-Related Lesions in
Different Age Categories. Neurobiol Aging 18: 351357.
27. Okello A, Koivunen J, Edison P, Archer HA, Turkheimer FE, et al. (2009)
Conversion of amyloid positive and negative MCI to AD over 3 years: an 11CPIB PET study. Neurology 73: 754760.
28. Small GW (2006) Diagnostic issues in dementia: neuroimaging as a surrogate
marker of disease. J Geriatr Psychiatry Neurol 19: 180185.
29. Herholz K, Ebmeier K (2011) Clinical amyloid imaging in Alzheimers disease.
The Lancet Neurology 10: 667670.
30. Hatashita S, Yamasaki H (2010) Clinically different stages of Alzheimers disease
associated by amyloid deposition with [11C]-PIB PET imaging. J Alzheimers
Dis 21: 9951003.
31. Barthel H, Senda M, Yamane T, Patt M, Sattler B, et al. (2010) Phase 1 trials on
florbetaben brain PET in elderly normal controls. J Nucl Med 51: 385.
32. Senda M, Yamane T, Sasaki M, Shimizu K, Barthel H, et al. (2010) Ethnic
comparability of safety and pharmacokinetics of florbetaben (BAY 94-9172) in
normal Caucasian and Asian subjects. J Nucl Med 51: 1754.
33. Fleisher AS, Chen K, Liu X, Roontiva A, Thiyyagura P, et al. (2011) Using
positron emission tomography and florbetapir F 18 to image cortical amyloid in
patients with mild cognitive impairment or dementia due to Alzheimer disease.
Arch Neurol 68: 14041411.
34. Jobst KA, Barnetson LPD, Shepstone BJ (1998) Accurate prediction of
histologically confirmed Alzheimers disease and the differential diagnosis of
dementia: the use of NINCDS-ADRDA and DSM-III-R criteria, SPECT, XRay CT, and Apo E4 in medial temporal lobe dementias. Int Psychogeriatr 10:
271302.
35. Ranginwala NA, Hynan LS, Weiner MF, White CLI (2008) Clinical criteria for
the diagnosis of Alzheimer disease: still good after all these years. Am J Geriatr
Psychiatry 16: 384-388.
36. LeVine III H, Walker LC (2010) Molecular polymorphism of Ab in Alzheimers
disease. Neurobiol Aging 31: 542548.
37. Lockhart A, Lamb JR, Osredkar T, Sue LI, Joyce JN, et al. (2007) PIB is a nonspecific imaging marker of amyloid-beta (Ab) peptide-related cerebral amyloidosis. Brain 130: 26072615.
38. Kemppainen NM, Aalto S, Wilson IA, Nagren K, Helin S, et al. (2007) PET
amyloid ligand [11C]PIB uptake is increased in mild cognitive impairment.
Neurology 68: 16031606.
39. Forsberg A, Engler H, Almkvist O, Blomquist G, Hagman G, et al. (2008) PET
imaging of amyloid deposition in patients with mild cognitive impairment.
Neurobiol Aging 29: 14561465.
40. Rowe CC, Ng S, Ackermann U, Gong SJ, Pike K, et al. (2007) Imaging betaamyloid burden in aging and dementia. Neurology 68: 17181725.
41. Aizenstein HJ, Nebes RD, Saxton JA, Price JC, Mathis CA, et al. (2008)
Frequent amyloid deposition without significant cognitive impairment among
the elderly. Arch Neurol 65: 15091517.
42. Rowe CC, Ellis KA, Rimajova M, Bourgeat P, Pike KE, et al. (2010) Amyloid
imaging results from the Australian Imaging, Biomarkers and Lifestyle (AIBL)
study of aging. Neurobiol Aging 31: 12751283.
43. Chetelat G, Villemagne VL, Bourgeat P, Pike KE, Jones G, et al. (2010)
Relationship between atrophy and beta-amyloid deposition in Alzheimer
disease. Ann Neurol 67: 317324.
44. Jack CR Jr, Lowe VJ, Weigand SD, Wiste HJ, Senjem ML, et al. (2009) Serial
PIB and MRI in normal, mild cognitive impairment and Alzheimers disease:
implications for sequence of pathological events in Alzheimers disease. Brain
132: 13551365.
45. Teipel SJ, Sabri O, Grothe M, Barthel H, Prvulovic D, et al. (2013) Perspectives
for Multimodal Neurochemical and Imaging Biomarkers in Alzheimers Disease.
J Alzheimers Dis 33: S329S347.
46. Roe CM, Mintun MA, DAngelo G, Xiong C, Grant EA, et al. (2008) Alzheimer
Disease and Cognitive Reserve: Variation of Education Effect With Carbon 11Labeled Pittsburgh Compound B Uptake. Arch Neurol 65: 14671471.
47. Barthel H, Luthardt J, Becker G, Patt M, Hammerstein E, et al. (2011)
Individualized quantification of brain b-amyloid burden: results of a proof of
mechanism phase 0 florbetaben PET trial in patients with Alzheimers disease
and healthy controls. Eur J Nucl Med Mol Imaging 38: 17021714.
48. Jack Jr CR, Knopman DS, Jagust WJ, Shaw LM, Aisen PS, et al. (2010)
Hypothetical model of dynamic biomarkers of the Alzheimers pathological
cascade. The Lancet Neurology 9: 119128.
79
PLOS ONE | www.plosone.org
Abstract
The role of structural brain magnetic resonance imaging (MRI) is becoming more and more emphasized in the early
diagnostics of Alzheimers disease (AD). This study aimed to assess the improvement in classification accuracy that can be
achieved by combining features from different structural MRI analysis techniques. Automatically estimated MR features
used are hippocampal volume, tensor-based morphometry, cortical thickness and a novel technique based on manifold
learning. Baseline MRIs acquired from all 834 subjects (231 healthy controls (HC), 238 stable mild cognitive impairment (SMCI), 167 MCI to AD progressors (P-MCI), 198 AD) from the Alzheimers Disease Neuroimaging Initiative (ADNI) database
were used for evaluation. We compared the classification accuracy achieved with linear discriminant analysis (LDA) and
support vector machines (SVM). The best results achieved with individual features are 90% sensitivity and 84% specificity
(HC/AD classification), 64%/66% (S-MCI/P-MCI) and 82%/76% (HC/P-MCI) with the LDA classifier. The combination of all
features improved these results to 93% sensitivity and 85% specificity (HC/AD), 67%/69% (S-MCI/P-MCI) and 86%/82% (HC/
P-MCI). Compared with previously published results in the ADNI database using individual MR-based features, the presented
results show that a comprehensive analysis of MRI images combining multiple features improves classification accuracy and
predictive power in detecting early AD. The most stable and reliable classification was achieved when combining all
available features.
Citation: Wolz R, Julkunen V, Koikkalainen J, Niskanen E, Zhang DP, et al. (2011) Multi-Method Analysis of MRI Images in Early Diagnostics of Alzheimers
Disease. PLoS ONE 6(10): e25446. doi:10.1371/journal.pone.0025446
Editor: Celia Oreja-Guevara, University Hospital La Paz, Spain
Received July 12, 2011; Accepted September 5, 2011; Published October 13, 2011
Copyright: 2011 Wolz et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This project is partially funded under the 7th Framework Programme by the European Commission (http://cordis.europa.eu/ist/). Data collection and
sharing for this project was funded by the Alzheimers Disease Neuroimaging Initiative (ADNI; Principal Investigator: Michael Weiner; NIH grant U01 AG024904).
ADNI is funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering (NIBIB), and through generous contributions
from the following: Pzer Inc., Wyeth Research, Bristol-Myers Squibb, Eli Lilly and Company, GlaxoSmithKline, Merck & Co. Inc., AstraZeneca AB, Novartis
Pharmaceuticals Corporation, Alzheimers Association, Eisai Global Clinical Development, Elan Corporation, plc, Forest Laboratories, and the Institute for the Study
of Aging, with participation from the U.S. Food and Drug Administration. Industry partnerships are coordinated through the Foundation for the National Institutes
of Health. The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimers Disease
Cooperative Study at the University of California San Diego. ADNI data are disseminated by the Laboratory of Neuro Imaging at the University of California Los
Angeles. Data used in the preparation of this article were obtained from the ADNI database (www.loni.ucla.edu/ADNI). As such, the investigators within the ADNI
contributed to the design and implementation of ADNI and/or provided data but did not participate in analysis or writing of this report. ADNI investigators
include (complete listing available at http://adni.loni.ucla.edu/wp-content/uploads/how_to_apply/ADNI_Authorship_List.pdf).
Competing Interests: The authors have declared that no competing interests exist.
* E-mail: [email protected]
. These authors contributed equally to this work.
" For information on the Alzheimers Disease Neuroimaging Initiative please see the Acknowledgments section.
Introduction
80
Subjects
In the ADNI study, brain MR images were acquired at regular
intervals after an initial baseline scan from approximately 200
cognitively normal older subjects (HC), 400 subjects with mild
cognitive impairment (MCI), and 200 subjects with early AD.
Detailled inclusion/exclusion criteria used for the different subject
groups in ADNI are defined in [30]. The AD group has scores
between 2026 (inclusive) on the Mini-Mental State Examination
(MMSE) [31], and a Clinical Dementia Rating (CDR) [32] of 0.5
or 1.0. Furthermore, these subjects fulfil the NINCDS/ADRDA
criteria for probable AD [33]. MCI subjects included have MMSE
scores between 2430 (inclusive), a memory complaint, have
objective memory loss measured by education adjusted scores on
Wechsler Memory Scale Logical Memory II, a CDR of 0.5,
absence of significant levels of impairment in other cognitive
domains, essentially preserved activities of daily living, and an
absence of dementia [30]. Healthy subjects have MMSE scores
between 2430 (inclusive), a CDR of 0, are non-depressed, non
MCI, and nondemented. A more detailed description of the ADNI
study is given in Appendix S1.
All 834 ADNI subjects (231 HC, 238 S-MCI, 167 P-MCI, 198
AD) for which a 1.5T T1-weighted MRI scan at baseline was
available were included in this study. 167 subjects in the MCI
group converted to AD as of July 2011. We therefore
independently analysed progressive MCI (P-MCI) subjects and
subjects with a stable diagnosis of MCI (S-MCI).
Table 1 shows the demographics for the 834 study subjects.
Statistically significant differences in the demographics and clinical
variables between the study groups were assessed using Students
unpaired t-test. In this work, the difference was considered
statistically significant if pv0.05 if not stated otherwise. There
were more men than women in all other groups besides the AD
group. MMSE scores were significantly different in the pairwise
comparisons between all study groups. CDR scores of the HC and
AD groups are significantly different to the ones of the two MCI
groups. Healthy subjects had a significantly lower Geriatric
Depression Scale (GDS) compared to all other groups. Compared
to all other groups, AD subjects had significantly shorter education.
MRI Acquisition
Standard 1.5T screening/baseline T1-weighted images obtained using volumetric 3D MPRAGE protocol with resolutions
Table 1. Subjects.
Group
HC
S-MCI
P-MCI
AD
231
238
167
198
Men
52%
66%
62%
52%
Age
76.02 (5.0)
74.85 (7.8)
74.6 (7.0)
75.68 (7.7)
MMSE
29.1* (1.0)
27.3* (1.8)
26.6* (1.7)
23.3* (2.0)
CDR
0 (0)
0.49 (0.05)
0.50 (0)
0.75 (0.25)
GDS
0.83* (1.14)
1.60 (1.42)
1.53 (1.30)
1.67 (1.42)
Education
16.0 (2.8)
15.6 (3.1)
15.7 (2.9)
14.7* (3.1)
APOE4 status
(e3e4/e4e4)
23%/2%
31%/8%
50%/16%
42%/18%
Months to
conversion
18.2 (10.1)
81
Feature extraction
82
Study design
Table 2 presents an overview on the features calculated for all
834 available ADNI baseline images. All feature values were
corrected for age and gender using a linear regression model
where control subjects were used as the training set, i.e., the
normal, not disease-related, age and gender related differences in
the classification features were removed. Feature selection was
then carried out on the corrected feature sets using stepwise
regression [51].
We used two subsets to perform classification:
I.
II.
Figure 1. 2D manifold embedding of a set of images acquired from healthy controls (red) and subjects with AD (blue).
doi:10.1371/journal.pone.0025446.g001
Classification methods
We used two different widely used methods to perform
classification based on individual features and their combination:
Figure 2. Results for voxelwise t-tests for statistically significant group differences with features extracted from TBM.
doi:10.1371/journal.pone.0025446.g002
83
Figure 3. Results for t-tests for statistically significant group differences based on cortical thickness measurements.
doi:10.1371/journal.pone.0025446.g003
Discussion
Results
Method
No of features
Description
Cortical thickness
9 (HC vs AD)
(CTH)
7 (HC vs P-MCI)
84
average Jacobian of atrophic voxels within a ROI, weighted based on voxel-wise p-values
20
8 (S-MCI vs P-MCI)
84
doi:10.1371/journal.pone.0025446.t002
Feature
LDA
SVM
Feature
SEN
SPE
SEN
SPE
MBL
87
83
85 [64 100]
87
83
HV
81
79
84
77
CTH
89
71
90
TBM
90
84
87 [71 100]
89
All
89 [71 100]
93
85
86 [71 100]
94
SEN
SPE
MBL
64
66
77
48
HV
63
67
62 [36 86]
83
33
73
CTH
63
45
59 [36 79]
96
03
84
TBM
65
62
77
44
78
All
68 [43 93]
67
69
60 [36 86]
92
14
SEN
SPE
LDA
SVM
Feature
LDA
SVM
SEN
SPE
SEN
SPE
MBL
81
75
84
69
HV
77
76
78 [54 92]
83
71
CTH
85
65
77 [54 100]
89
62
TBM
82
76
85
74
All
84 [62 100]
86
82
82 [62 100]
93
67
Feature
SEN
HC vs AD
SPE
HC vs P-MCI
S-MCI vs P-MCI
SEN
SEN
SPE
SPE
MBL
90
74
84
92
55
76
HV
80
69
75
76
63
70
CTH
85
75
86
59
72
35
TBM
93
76
90
84
63
59
All
94
76
94
89
69
54
doi:10.1371/journal.pone.0025446.t006
85
Table 7. Classification results of healthy control (HC), mild cognitive impairment (MCI) and Alzheimers disease subjects reported
in the recent literature.
Study
Features
HC vs AD
HC vs P-MCI
S-MCI vs P-MCI
CCR
SEN
SPE
CCR
SEN
SPE
CCR
SEN
SPE
333
Cortical volumes
91
92
90
70
Hippocampus shape
94
96
92
605
Hippocampus volume
76
75
77
64
60
65
382
Cortical thickness
85
73
75
68
312
Amygdala/caudate volumes
69
76
68
356
SPARE-AD index
56
95
38
509
Various
81
95
73
85
62
69
159
81
60
92
14
351
Various volumes
82
73
398
89
83
93
380
STAND score
86
86
86
Supporting Information
Appendix S1
(DOCX)
Appendix S2 ROI-wise features for CTH and TBM.
(DOCX)
Author Contributions
Acknowledgments
Data used in the preparation of this article were obtained from the ADNI
database (www.loni.ucla.edu/ADNI). As such, the investigators within the
References
24. Liu Y, Paajanen T, Zhang Y, Westman E, Wahlund LO, et al. (2009)
Combination analysis of neuropsychological tests and structural MRI measures
in differentiating AD, MCI and control groupsThe AddNeuroMed study.
Neurobiology of Aging In Press, Corrected Proof: -.
25. Liu Y, Paajanen T, Zhang Y, Westman E, Wahlund LO, et al. (2010) Analysis of
regional MRI volumes and thicknesses as predictors of conversion from mild
cognitive impairment to Alzheimers disease. Neurobiology of Aging 31:
13751385.
26. McEvoy L, Fennema-Notestine C, JC R, Hagler D, Holland D, et al. (2009)
Alzheimer disease: quantitative structural neuroimaging for detection and
prediction of clinical and structural changes in mild cognitive impairment.
Radiology 251: 1950205.
27. Querbes O, Aubry F, Pariente J, Lotterie JA, Dmonet JF, et al. (2009) Early
diagnosis of alzheimers disease using cortical thickness: impact of cognitive
reserve. Brain 132: 20362047.
28. Mueller SG, Weiner MW, Thal LJ, Petersen RC, Jack C, et al. (2005) The
Alzheimers Disease Neuroimaging Initiative. Neuroimaging Clinics of North
America 15: 869877.
29. Cuingnet R, Gerardin E, Tessieras J, Auzias G, Lehricy S, et al. (2010)
Automatic classification of patients with Alzheimers disease from structural
MRI: A comparison of ten methods using the ADNI database. NeuroImage In
Press, Corrected Proof: -.
30. Petersen RC, Aisen PS, Beckett LA, Donohue MC, Gamst AC, et al. (2010)
Alzheimers Disease Neuroimaging Initiative (ADNI): clinical characterization.
Neurology 74: 201209.
31. Folstein MF, Folstein SE, McHugh PR (1975) Mini-mental state: A practical
method for grading the cognitive state of patients for the clinician. Journal of
Psychiatric Research 12(3): 189198.
32. Morris J (1993) The Clinical Dementia Rating (CDR): current version and
scoring rules. Neurology 43: 24122414.
33. McKhann G, Drachman D, Folstein M, Katzman R, Price D, et al. (1984)
Clinical diagnosis of alzheimers disease: report of the nincds-adrda work group
under the auspices of department of health and human services task force on
alzheimers disease. Neurology 34: 939944.
34. Jack CR, Jr., Bernstein MA, Fox NC, Thompson P, Alexander G, et al. (2008)
The Alzheimers disease neuroimaging initiative (ADNI): MRI methods. Journal
of Magnetic Resonance Imaging 27: 685691.
35. Lotjonen JM, Wolz R, Koikkalainen JR, Thurfjell L, Waldemar G, et al. (2010)
Fast and robust multi-atlas segmentation of brain magnetic resonance images.
NeuroImage 49: 23522365.
36. Lotjonen J, Wolz R, Koikkalainen J, Julkunen V, Thurfjell L, et al. (2011) Fast
and robust extraction of hippocampus from mr images for diagnostics of
alzheimers disease. NeuroImage 56: 185196.
37. Van Leemput K, Maes F, Vandermeulen D, Suetens P (1999) Automated
model-based tissue classification of MR images of the brain. IEEE Transactions
on Medical Imaging 18: 897908.
38. Lerch J, Evans A (2005) Cortical thickness analysis examined through power
analysis and a population simulation. NeuroImage 24: 163173.
39. Mazziotta J, Toga A, Evans A, Fox P, Lancaster J, et al. (2001) A probabilistic
atlas and reference system for the human brain: International Consortium for
Brain Mapping (ICBM). Philos Trans R Soc Lond B Biol Sci 356: 12931322.
40. Sled JG, Zijdenbos AP, Evans AC (1998) A nonparametric method for
automatic correction of intensity nonuniformity in MRI data. IEEE Transactions on Medical Imaging 17: 8797.
41. Smith SM (2002) Fast robust automated brain extraction. Hum Brain Mapp 17:
143155.
42. Zijdenbos AP, Forghani R, Evans AC (1998) Automatic quantification of MS
lesions in 3D MRI brain data sets: Validation of INSECT. In: III WMW,
Colchester ACF, Delp SL, eds. MICCAI Springer, volume 1496 of Lecture
Notes in Computer Science. pp 439448.
43. Tohka J, Zijdenbos A, Evans A (2004) Fast and robust parameter estimation for
statistical partial volume models in brain mri. NeuroImage 23: 8497.
44. Kim JS, Singh V, Lee JK, Lerch J, Ad-Dabbagh Y, et al. (2005) Automated 3-d
extraction and evaluation of the inner and outer cortical surfaces using a
laplacian map and partial volume effect classification. NeuroImage 27: 210221.
45. Chung M, Taylor J (2004) Diffusion smoothing on brain surface via finite
element method. In: ISBI IEEE. pp 432435.
87
46. Brun CC, Lepor N, Pennec X, Lee AD, Barysheva M, et al. (2009) Mapping the
regional inuence of genetics on brain structure variability a tensor-based
morphometry study. NeuroImage 48: 3749.
47. Koikkalainen J, Lotjonen J, Thurfjell L, Rueckert D, Waldemar G, et al. (2011)
Multi-template tensor-based morphometry: Application to analysis of alzheimers disease. NeuroImage 56: 11341144.
48. Belkin M, Niyogi P (2003) Laplacian eigenmaps for dimensionality reduction
and data representation. Neural Computation 15: 13731396.
49. Rueckert D, Sonoda LI, Hayes C, Hill DLG, Leach MO, et al. (1999) Nonrigid
registration using free-form deformations: Application to breast MR images.
IEEE Transactions on Medical Imaging 18: 712721.
50. Wolz R, Heckemann RA, Aljabar P, Hajnal JV, Hammers A, et al. (2010)
Measurement of hippocampal atrophy using 4D graph-cut segmentation:
Application to ADNI. NeuroImage 52: 10091018.
51. Draper NR, Smith H (1998) Applied Regression Analysis (Wiley Series in
Probability and Statistics). Wiley.
52. Krzanowski WJ (1988) Principles of Multivariate Analysis: A Users Perspective
Oxford University Press.
53. Cortes C, Vapnik V (1995) Support-vector networks. Machine Learning 20:
273297.
54. Vemuri P, Gunter JL, Senjem ML, Whitwell JL, Kantarci K, et al. (2008)
Alzheimers disease diagnosis in individual subjects using structural mr images:
Validation studies. NeuroImage 39: 11861197.
55. Westman E, Simmons A, Zhang Y, Muehlboeck JS, Tunnard C, et al. (2011)
Multivariate analysis of mri data for alzheimers disease, mild cognitive
impairment and healthy controls. NeuroImage 54: 11781187.
88
PLoS ONE | www.plosone.org
Abstract
Introduction: Various biomarkers have been reported in recent literature regarding imaging abnormalities in different types
of dementia. These biomarkers have helped to significantly improve early detection and also differentiation of various
dementia syndromes. In this study, we systematically applied whole-brain and region-of-interest (ROI) based support vector
machine classification separately and on combined information from different imaging modalities to improve the detection
and differentiation of different types of dementia.
Methods: Patients with clinically diagnosed Alzheimers disease (AD: n = 21), with frontotemporal lobar degeneration (FTLD:
n = 14) and control subjects (n = 13) underwent both [F18]fluorodeoxyglucose positron emission tomography (FDG-PET)
scanning and magnetic resonance imaging (MRI), together with clinical and behavioral assessment. FDG-PET and MRI data
were commonly processed to get a precise overlap of all regions in both modalities. Support vector machine classification
was applied with varying parameters separately for both modalities and to combined information obtained from MR and
FDG-PET images. ROIs were extracted from comprehensive systematic and quantitative meta-analyses investigating both
disorders.
Results: Using single-modality whole-brain and ROI information FDG-PET provided highest accuracy rates for both,
detection and differentiation of AD and FTLD compared to structural information from MRI. The ROI-based multimodal
classification, combining FDG-PET and MRI information, was highly superior to the unimodal approach and to the wholebrain pattern classification. With this method, accuracy rate of up to 92% for the differentiation of the three groups and an
accuracy of 94% for the differentiation of AD and FTLD patients was obtained.
Conclusion: Accuracy rate obtained using combined information from both imaging modalities is the highest reported up
to now for differentiation of both types of dementia. Our results indicate a substantial gain in accuracy using combined
FDG-PET and MRI information and suggest the incorporation of such approaches to clinical diagnosis and to differential
diagnostic procedures of neurodegenerative disorders.
Citation: Dukart J, Mueller K, Horstmann A, Barthel H, Moller HE, et al. (2011) Combined Evaluation of FDG-PET and MRI Improves Detection and Differentiation of
Dementia. PLoS ONE 6(3): e18111. doi:10.1371/journal.pone.0018111
Editor: Wang Zhan, University of California, San Francisco, United States of America
Received August 30, 2010; Accepted February 25, 2011; Published March 23, 2011
Copyright: 2011 Dukart et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The work was funded by the Max Planck Society. The funders had no role in study design, data collection and analysis, decision to publish, or
preparation of the manuscript.
Competing Interests: The authors have declared that no competing interests exist.
* E-mail: [email protected]
Introduction
In recent research, various biomarkers have been reported to
differentiate between early stages of dementia and healthy control
subjects or between different types of neurodegenerative disorders,
suggesting an integration of these would improve diagnostic
accuracy of dementia [1][9].
For the detection of dementia, accuracy rates significantly above
90% have recently been reported using univariate and multivariate statistical approaches in magnetic resonance imaging (MRI)
and [F18]fluorodeoxyglucose positron emission tomography
(FDG-PET) [1],[10][14]. However, the differentiation of the
two most common types of dementia, namely Alzheimers disease
(AD) and frontotemporal lobar degeneration (FTLD), is still
89
Controls
ANOVA
(df,F,P)
13
21
14
Male/Female
7/6
9/12
7/7
Age (years)
53.966.0
61.166.7
60.866.4
2, 5.76, 0.006
CDR (score)
0.2360.26
0.7160.25
0.8260.42
2, 13.93, 0.000
MMSE (score)
n.a.
23.263.9
24.464.2
10.763.1
11.663.8
2,1.02,0.368
Data acquisition
MRI data. For each subject, a high-resolution T1-weighted
MRI scan was obtained, consisting of 128 sagittal slices adjusted to
AC-PC line and a with slice thickness of 1.5 mm and pixel size of
161 mm2. MRI was performed on two different 3T scanners
(MedSpec 30/100, Bruker Biospin, Ettlingen Germany and
Magnetom Trio, Siemens, Erlangen, Germany) using two
different T1-weighted sequences (MDEFT or MP-RAGE with
TR = 1300 ms, TI = 650 ms, TE = 3.93 ms or TE = 10 ms; FOV
25625 cm2; matrix = 2566256 voxels). On the MedSpec
scanner, only the MDEFT-sequence and on the Magnetom Trio
scanner, either MDEFT or MP-RAGE sequences were used. The
distribution of scanner types and sequences used to obtain the
MRI data was random across subjects and did not differ
significantly in its distribution between the groups nor for
scanner type nor for sequence.
PET data. Each subject also underwent FDG-PET imaging
either a few a weeks before or after the MRI scan. All PET data
were acquired on a Siemens ECAT EXACT HR+ scanner (CTI/
Siemens, Knoxville, TN, USA) under a standard resting condition
in 2-dimensional (2D) mode. The 2D acquisition mode was used
because it allows a better quantification of the PET data due to
lower scatter radiation. Sixty-three slices were simultaneously
collected with an axial resolution of 5 mm full width at half
maximum (FWHM) and in-plane resolution of 4.6 mm. After
correction for attenuation, scatter, decay and scanner-specific
dead-time, images were reconstructed by filtered back-projection
using a Hann-filter of 4.9 mm FWHM. The 63 transaxial slices
obtained had a matrix of 1286128 voxels with an edge length of
2.45 mm.
Ethics Statement
The research protocol was approved by the Ethics Committee
of the University of Leipzig, and was in accordance with the latest
version of the Declaration of Helsinki. Informed consent was
obtained from all subjects.
Subjects
We analyzed FDG-PET and T1-weighted MRI data of 21
patients (Table 1) with an early stage of probable AD, 14 patients
with an early stage of FTLD and 13 control subjects. Patients were
recruited from the Day Clinic of Cognitive Neurology at the
University of Leipzig. Probable AD was diagnosed according to
NINCDS-ADRDA criteria [28]. Although all AD subjects also
fulfilled the revised NINCDS-ADRDA criteria suggested by
Dubois et al. [17] the fulfillment of the original McKhann criteria
was sufficient for the inclusion into the study. Diagnosis of FTLD
was based on criteria suggested by Neary et al. [29]. The control
group included subjects who visited the Day Clinic with subjective
cognitive complaints, which were not objectively confirmed by a
comprehensive neuropsychological and clinical evaluation. FDGPET and MRI for these subjects was conducted for diagnostic
reasons within the clinical assessment. This control group was
chosen because, in clinical practice, it is crucial to discriminate
between these subjects showing a normal age-related decrease in
cognitive performance and patients with an early stage of
dementia. Patients were excluded if structural imaging revealed
PLoS ONE | www.plosone.org
FTLD
Number
Methods
90
AD
Figure 1. Schematic representation of the procedure for FDG-PET and MRI data handling and processing steps. FDG-PET
[F18]fluorodeoxyglucose positron emission tomography, MRI magnetic resonance imaging, MNI Montreal Neurological Institute, PVE partial volume
effect correction.
doi:10.1371/journal.pone.0018111.g001
91
92
BA
Lat x
24/32
25
220
Medial thalamus
22
219 6
34 21
23
Anterior insula
15/16
247
10 27
10
54 0
25
22 225
Amygdala
MRI
BA
Lat x
25
49 25
24/32
25
34 21
44
254
11 20
220
247
10 27
15/16
Subcallosal/septal area
25
23
22
10 210
Amygdala
25
22 225
226
25 225
Temporal pole
38
229
11 243
Temporal pole
38
247
244
BA
Lat x
AD vs. Controls
FDG-PET
Angular gyrus
39
238
268 37
Angular gyrus
39
43
268 33
21/22
R
31 31
236 27
9/10
32
Inferior precuneus
31
23
21/22
251
261 23
20/21
59
231 223
MRI
BA
Lat x
Posterior insula
13
238
225 15
Medial thalamus
25
213 3
Hippocampal body/tail
31
238 26
21/22
263
221 25
25
28 218
226
28 218
Results
Clinical characteristics
The chi-square test for independent samples did not reveal any
statistical differences in sex between the groups [x2(2) = 0.42;
p = 0.809]. The three groups did not differ in education (Table 1).
CDR scores differed significantly in the three groups. The post-hoc
test revealed no differences in the mean CDR scores between both
groups of dementia patients [t(33) = 0.94;p = 0.977]. As expected,
both early AD [t(32) = 5.36;p,0.001] and early FTLD
[t(25) = 4.35;p,0.001] had significantly higher CDR scores compared to the control subjects. MMSE scores also did not differ
between both groups of dementia patients indicating a similar severity
of dementia syndrome [t(29) = 281;p = 0.95]. The ANOVA also
revealed a significant group difference in age. The two groups of
dementia patients did not differ significantly in age
[t(33) = 0.16;p = 1.0]. There was a minor but significant difference
between AD patients and controls [t(32) = 3.18;p = 008] and FTLD
patients and controls [t(25) = 2.86; p = 024].
93
Table 3. Accuracy rates for whole-brain and ROI-based SVM classification for FDG-PET and MRI.
AD vs FTLD
AD vs Controls
FTLD vs Controls
GM
whole-brain
72.9%
80.0%
88.2%
77.8%
WM
whole-brain
66.7%
74.3%
79.4%
77.8%
FDG-PET
whole-brain
81.3%
82.9%
94.1%
92.6%
GM/ FDG-PET
whole-brain
79.2%
82.9%
94.1%
88.9%
GM/WM/FDG-PET
whole-brain
77.1%
82.9%
91.2%
85.2%
GM + FDG-PET
whole-brain
81.3%
88.6%
91.2%
88.9%
GM
ROIs
56.3%
60.0%
82.4%
85.2%
FDG-PET
ROIs
75.0%
80.0%
94.1%
85.2%
GM/FDG-PET
ROIs
91.7%
94.3%
100.0%
92.6%
Accuracy represents the percentage of subjects correctly assigned to the correct condition. AD Alzheimers disease, FDG-PET [18F]fluorodeoxyglucose positron emission
tomography, FTLD frontotemporal lobar degeneration, GM gray matter, MRI magnetic resonance imaging, ROI region-of-interest, SVM support vector machine, WM
white matter.
doi:10.1371/journal.pone.0018111.t003
Discussion
In this study we performed a multimodal comparison and
discrimination of dementia patients using FDG-PET and MRI. To
enable a more accurate coevaluation of both imaging modalities,
we developed a new preprocessing algorithm. This algorithm was
designed to enable an accurate anatomical registration of both
modalities. All processing steps were performed as far as possible
simultaneously by applying the same deformations and preprocessing parameters to both modalities of the same subject. This
procedure resulted in an accurate anatomical overlap of both
imaging modalities and in an accurate between-subject registration, with both images having the same voxel size and
approximately the same effective smoothness.
SVM
94
Accuracy
Sensitivity
Specificity
92.9%
AD vs FTLD
94.3%
95.2%*
AD vs Controls
100.0%
100.0%
100.0%
FTLD vs Controls
92.6%
85.7%
100.0%
95
Author Contributions
Conceived and designed the experiments: JD KM AH HB HEM AV OS
MLS. Performed the experiments: JD KM AH HB HEM AV OS MLS.
Analyzed the data: JD KM AH HB HEM AV OS MLS. Contributed
reagents/materials/analysis tools: JD KM AH HB HEM AV OS MLS.
Wrote the paper: JD KM AH HB HEM AV OS MLS.
References
1. Hoffman JM, Welsh-Bohmer KA, Hanson M, Crain B, Hulette C, et al. (2000)
FDG PET imaging in patients with pathologically verified dementia. J Nucl Med
41: 19201928.
2. Rosen HJ, Gorno-Tempini ML, Goldman WP, Perry RJ, Schuff N, et al. (2002)
Patterns of brain atrophy in frontotemporal dementia and semantic dementia.
Neurology 58: 198208.
3. Diehl J, Grimmer T, Drzezga A, Riemenschneider M, Forstl H, et al. (2004)
Cerebral metabolic patterns at early stages of frontotemporal dementia and
semantic dementia. A PET study. Neurobiol Aging 25: 10511056.
4. Jeong Y, Cho SS, Park JM, Kang SJ, Lee JS, et al. (2005) 18F-FDG PET
findings in frontotemporal dementia: an SPM analysis of 29 patients. J Nucl Med
46: 233239.
5. Diehl-Schmid J, Grimmer T, Drzezga A, Bornschein S, Riemenschneider M,
et al. (2007) Decline of cerebral glucose metabolism in frontotemporal dementia:
a longitudinal 18F-FDG-PET-study. Neurobiol Aging 28: 4250.
6. Fung G, Stoeckel J (2007) SVM feature selection for classification of SPECT
images of Alzheimers disease using spatial information. Knowledge and
Information Systems 11: 243258.
7. Sabri O, Kendziorra K, Wolf H, Gertz HJ, Brust P (2008) Acetylcholine
receptors in dementia and mild cognitive impairment. Eur J Nucl Med Mol
Imaging 35(Suppl 1): S3045.
8. Schroeter ML, Raczka K, Neumann J, Yves von Cramon D (2007) Towards a
nosology for frontotemporal lobar degenerations-a meta-analysis involving 267
subjects. Neuroimage 36: 497510.
9. Schroeter ML, Stein T, Maslowski N, Neumann J (2009) Neural correlates of
Alzheimers disease and mild cognitive impairment: a systematic and
quantitative meta-analysis involving 1351 patients. Neuroimage 47: 11961206.
10. Kloppel S, Stonnington CM, Chu C, Draganski B, Scahill RI, et al. (2008)
Automatic classification of MR scans in Alzheimers disease. Brain 131:
681689.
11. Matsunari I, Samuraki M, Chen WP, Yanase D, Takeda N, et al. (2007)
Comparison of 18F-FDG PET and optimized voxel-based morphometry for
detection of Alzheimers disease: aging effect on diagnostic performance. J Nucl
Med 48: 19611970.
12. Davatzikos C, Resnick SM, Wu X, Parmpi P, Clark CM (2008) Individual
patient diagnosis of AD and FTD via high-dimensional pattern classification of
MRI. Neuroimage 41: 12201227.
13. Fan Y, Resnick SM, Wu X, Davatzikos C (2008) Structural and functional
biomarkers of prodromal Alzheimers disease: a high-dimensional pattern
classification study. Neuroimage 41: 277285.
14. Sadeghi N, Foster NL, Wang AY, Minoshima S, Lieberman AP, et al. (2008)
Automatic classification of Alzheimers disease vs. Frontotemporal dementia: A
spatial decision tree approach with FDG-PET. 2008 Ieee International
Symposium on Biomedical Imaging: From Nano to Macro Vols 1-4: 408411.
15. Knopman DS, Boeve BF, Parisi JE, Dickson DW, Smith GE, et al. (2005)
Antemortem diagnosis of frontotemporal lobar degeneration. Ann Neurol 57:
480488.
16. Rascovsky K, Salmon DP, Ho GJ, Galasko D, Peavy GM, et al. (2002) Cognitive
profiles differ in autopsy-confirmed frontotemporal dementia and AD.
Neurology 58: 18011808.
17. Dubois B, Feldman HH, Jacova C, Dekosky ST, Barberger-Gateau P, et al.
(2007) Research criteria for the diagnosis of Alzheimers disease: revising the
NINCDS-ADRDA criteria. Lancet Neurol 6: 734746.
18. Kipps CM, Hodges JR, Fryer TD, Nestor PJ (2009) Combined magnetic
resonance imaging and positron emission tomography brain imaging in
behavioural variant frontotemporal degeneration: refining the clinical phenotype. Brain 132: 25662578.
19. Ishii K, Willoch F, Minoshima S, Drzezga A, Ficaro EP, et al. (2001) Statistical
brain mapping of 18F-FDG PET in Alzheimers disease: validation of anatomic
standardization for atrophied brains. J Nucl Med 42: 548557.
20. Ishii K, Sasaki H, Kono AK, Miyamoto N, Fukuda T, et al. (2005) Comparison
of gray matter and metabolic reduction in mild Alzheimers disease using FDGPET and voxel-based morphometric MR studies. Eur J Nucl Med Mol Imaging
32: 959963.
21. Kanda T, Ishii K, Uemura T, Miyamoto N, Yoshikawa T, et al. (2008)
Comparison of grey matter and metabolic reductions in frontotemporal
22.
23.
24.
25.
26.
27.
28.
29.
30.
31.
32.
33.
34.
35.
36.
37.
38.
39.
40.
41.
42.
43.
96
PLoS ONE | www.plosone.org
Patrizia Mecocci9, Simon Lovestone1,2, Stephen Newhouse1,2., Richard Dobson1,2.*, for the Alzheimers
Disease Neuroimaging Initiative
1 National Institute of Health Research Biomedical Research Centre for Mental Health, South London and Maudsley National Health Service Foundation Trust, London,
United Kingdom, 2 Kings College London, Institute of Psychiatry, London, United Kingdom, 3 Laboratory of Behavioral Neuroscience, National Institute on Aging,
Baltimore, Maryland, United States of America, 4 Proteome Sciences plc, Cobham, Surrey, United Kingdom, 5 School of Neurology, University of Eastern Finland and
University Hospital of Kuopio, Kuopio, Finland, 6 Department of Old Age Psychiatry and Psychotic Disorders, Medical University of Lodz, Lodz, Poland, 7 3rd Department of
Neurology, G. Papanicolaou Hospital, Aristotle University of Thessaloniki, Thessaloniki, Greece, 8 Department of Geriatric Medicine, Grontople de Toulouse, Toulouse
University Hospital, Toulouse, France, 9 Institute of Gerontology and Geriatrics, University of Perugia, Perugia, Italy
Abstract
Changes in brain amyloid burden have been shown to relate to Alzheimers disease pathology, and are believed to precede
the development of cognitive decline. There is thus a need for inexpensive and non-invasive screening methods that are
able to accurately estimate brain amyloid burden as a marker of Alzheimers disease. One potential method would involve
using demographic information and measurements on plasma samples to establish biomarkers of brain amyloid burden; in
this study data from the Alzheimers Disease Neuroimaging Initiative was used to explore this possibility. Sixteen of the
analytes on the Rules Based Medicine Human Discovery Multi-Analyte Profile 1.0 panel were found to associate with [11 C]PiB PET measurements. Some of these markers of brain amyloid burden were also found to associate with other AD related
phenotypes. Thirteen of these markers of brain amyloid burden c-peptide, fibrinogen, alpha-1-antitrypsin, pancreatic
polypeptide, complement C3, vitronectin, cortisol, AXL receptor kinase, interleukin-3, interleukin-13, matrix metalloproteinase-9 total, apolipoprotein E and immunoglobulin E were used along with co-variates in multiple linear regression, and
were shown by cross-validation to explain w30% of the variance of brain amyloid burden. When a threshold was used to
classify subjects as PiB positive, the regression model was found to predict actual PiB positive individuals with a sensitivity of
0.918 and a specificity of 0.545. The number of APOE E 4 alleles and plasma apolipoprotein E level were found to contribute
most to this model, and the relationship between these variables and brain amyloid burden was explored.
Citation: Kiddle SJ, Thambisetty M, Simmons A, Riddoch-Contreras J, Hye A, et al. (2012) Plasma Based Markers of [11C] PiB-PET Brain Amyloid Burden. PLoS
ONE 7(9): e44260. doi:10.1371/journal.pone.0044260
Editor: Ashley I. Bush, Mental Health Research Institute of Victoria, Australia
Received April 13, 2012; Accepted July 31, 2012; Published September 24, 2012
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for
any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Funding: Alzheimers Disease Neuroimaging Initiative (ADNI) data collection and sharing for this project was funded by the ADNI (National Institutes of Health Grant
U01 AG024904). ADNI is funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering, and through generous
contributions from the following: Abbott; Alzheimers Association; Alzheimers Drug Discovery Foundation; Amorfix Life Sciences Ltd.; AstraZeneca; Bayer HealthCare;
BioClinica, Inc.; Biogen Idec Inc.; Bristol-Myers Squibb Company; Eisai Inc.; Elan Pharmaceuticals Inc.; Eli Lilly and Company; F. Hoffmann-La Roche Ltd and its affiliated
company Genentech, Inc.; GE Healthcare; Innogenetics, N.V.; Janssen Alzheimer Immunotherapy Research & Development, LLC.; Johnson & Johnson Pharmaceutical
Research & Development LLC.; Medpace, Inc.; Merck & Co., Inc.; Meso Scale Diagnostics, LLC.; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Servier; Synarc Inc.; and
Takeda Pharmaceutical Company. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. Private sector contributions are
facilitated by the Foundation for the National Institutes of Health (www.fnih.org). The grantee organization is the Northern California Institute for Research and
Education, and the study is coordinated by the Alzheimers Disease Cooperative Study at the University of California, San Diego. ADNI data are disseminated by the
Laboratory for Neuro Imaging at the University of California, Los Angeles. This research was also supported by National Institutes of Health grants P30 AG010129 and K01
AG030514. This study was supported by InnoMed (Innovative Medicines in Europe www.imi.europa.eu/) an Integrated Project funded by the European Union of the
Sixth Framework program priority FP6-2004-LIFESCIHEALTH-5, Life Sciences, Genomics and Biotechnology for Health. This study was also supported by funds from the
National Institutes for Health Research Biomedical Research Centre for Mental Health at the South London and Maudsley National Health Service Foundation Trust and
Institute of Psychiatry, Kings College London. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing Interests: Intellectual property has been registered on the use of plasma proteins for use as biomarkers for AD by Kings College London and Proteome
Sciences, with SL and M. Thambisetty named as inventors. IP and MW were full-time employees of Proteome Sciences, London, United Kingdom, at the time of their
contribution to the work described in this manuscript. This does not alter the authors adherence to all the PLOS ONE policies on sharing data and materials. Patent Title
Methods and Compositions Relating to Alzheimers Disease Subject Covers utility of around 30 proteins, specifically listing 16 in the Dependent claims for diagnosis of
Alzheimers disease Filing Proteome Sciences with Kings Business United Kingdom Priority GB0421639.6 dated 29/09/2004 PCT Application PCT/GB2005/003756 dated
29/09/2005 Application in Europe, Japan, United States, Australia, and Canada, dated 15/03/2007 to 16/10/2007. In addition Alzheimers Disease Neuroimaging
Initiative recieved funding from: Abbott; Alzheimers Association; Alzheimers Drug Discovery Foundation; Amorfix Life Sciences Ltd.; AstraZeneca; Bayer HealthCare;
BioClinica, Inc.; Biogen Idec Inc.; Bristol-Myers Squibb Company; Eisai Inc.; Elan Pharmaceuticals Inc.; Eli Lilly and Company; F. Hoffmann-La Roche Ltd and its affiliated
company Genentech, Inc.; GE Healthcare; Innogenetics, N.V.; Janssen Alzheimer Immunotherapy Research & Development, LLC.; Johnson & Johnson Pharmaceutical
Research & Development LLC.; Medpace, Inc.; Merck & Co., Inc.; Meso Scale Diagnostics, LLC.; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Servier; Synarc Inc.; and
Takeda Pharmaceutical Company. This does not alter the authors adherence to all the PLOS ONE policies on sharing data and materials.
* E-mail: [email protected] (SJK); [email protected] (RD)
. These authors contributed equally to this work.
97
Introduction
The failure of several clinical trials targeting brain amyloid
deposition in patients with Alzheimers disease (AD) has led to the
suggestion that these agents may be useful if targeted at older
individuals in pre-symptomatic stages of the disease [1,2].
Screening methods that accurately identify at-risk non-demented
older individuals who are most likely to benefit from such
treatments will therefore represent a major advance in our ability
to effectively test these disease-modifying treatments [3]. If clinical
trials of amyloid lowering interventions were successful in the presymptomatic stages of AD, then there would be a desire to identify
non-demented elderly individuals with elevated brain amyloid
burden who could potentially benefit from early intervention.
However identifying suitable individuals poses a great challenge in
terms of feasibility and cost. To date, the two methods that are
most likely to be useful in estimating levels of brain amyloid
burden are in vivo imaging with positron emission tomography
(PET) using radioligands binding to fibrillar amyloid beta (Ab),
such as [11 C] Pittsburgh B compound (PiB), and assays of Ab levels
in cerebrospinal fluid (CSF) [4,5]. However, both these methods
have inherent drawbacks that limit their utility as screening tools,
especially in resource-poor settings. While PET scanning is
expensive and limited to specialised centres, lumbar puncture to
obtain CSF is associated with some patient discomfort and is
unlikely to be used in primary health care centres to routinely
screen large numbers of elderly patients. An inexpensive, noninvasive screening method that accurately estimates brain amyloid
burden would therefore fulfil a critical unmet need in the care of
the elderly.
The identification of blood-based biomarkers associated with
AD diagnosis [69] or distinct endophenotypes of AD pathology
such as brain atrophy [1012], hippocampal metabolite abnormalities [13] and amyloid burden [14], have previously been
reported. In these studies, proteomic analyses were combined with
neuroimaging methods to identify plasma signals associated with
measures of AD pathology. In this study, a different strategy was
used by examining the association between brain amyloid burden
and a panel of 146 plasma analytes proteins, complexes and
metabolites measured by Rules Based Medicine, Inc. (RBM)
(Austin, TX) using the Human Discovery Multi-Analyte Profile
(MAP) 1.0 panel and a Luminex 100 platform. Some of the
analytes on this panel, such as apolipoprotein E (APOE) and
Results
RBM analytes associate with Ab levels in the brain
Levels of analytes measured by the RBM Human Discovery
MAP 1.0 from ADNI plasma samples were compared to fibrillar
amyloid in the RBM-PiB PET cohort (N = 71). Characteristics of
this subcohort are summarised in Table 1 where it can be seen that
brain amyloid burden was almost significantly different at the 0.05
level between diagnostic groups (Kruskal-Wallis (KW) x2 test pvalue 0.055). The distribution of brain amyloid burden in the
RBM-PiB PET cohort is shown in Figure S1. In the slightly larger
ADNI-PiB PET cohort (i.e. all ADNI subjects with [11 C] PiB-PET
scans performed at baseline), whose sample characteristics are
shown in Table S1, brain amyloid burden was found to be
significantly different across diagnostic groups (KW p-value 0.022).
The analytes most associated with brain amyloid burden in the
RBM-PiB PET cohort, after taking into account co-variates (age,
gender, years of education, number of APOE E 4 alleles and the
number of days between [11 C]-PiB PET scan and plasma sample),
are shown in Table 2.
Control (3)
MCI (52)
AD (16)
P-value
77.4 [5.6]
75.4 [11.1]
72.3 [8.2]
0.290
(Median [IQR])
Sex (Male/Female)
1/2
37/15
10/6
0.263
13.0 [3.0]
16.0 [5.0]
16.0 [5.3]
0.398
2/1/0
25/22/5
7/7/2
0.977
Days between [ C]-PiB PET scan and plasma sample (Median [IQR])
5.0 [15.5]
23.5 [60.5]
21.5 [42.3]
0.288
1.31 [0.108]
1.98 [0.723]
1.90 [0.438]
0.055
11
(Median [IQR])
98
Characteristics of the ADNI RBM-PiB PET subcohort by diagnostic group. P-values were calculated for differences across diagnostic groups using a Kruskal-Wallis x2 test
for continuous characteristics and simulated contingency p-values for categorical characteristics.
doi:10.1371/journal.pone.0044260.t001
RBM analyte
Gene name
Uniprot ID
C-peptide
INS
P01308
20.310
0.010
0.351
FG(A/B/G)
P02671 P02675
P02679
20.307
0.010
0.351
Alpha-1-antitrypsin
SERPINA1
P01009
20.302
0.012
0.351
Pancreatic polypeptide
PPY
P01298
20.296
0.014
0.351
Complement C3
C3
P01024
20.296
0.014
0.351
Vitronectin
VTN
P04004
20.295
0.014
0.351
VWF
P04275
20.287
0.017
0.363
Cortisol
(NA)
(NA)
0.271
0.025
0.412
APCS
P02743
20.268
0.027
0.412
AXL
P30530
0.266
0.028
0.412
Interleukin-3
IL3
P08700
0.261
0.032
0.412
Interleukin-13
IL13
P35225
0.252
0.038
0.412
MMP9
P14780
20.250
0.040
0.412
APOE
APOE
P02649
20.248
0.042
0.412
Leptin
LEP
P41159
20.248
0.042
0.412
Immunoglobulin E (IgE)
(NA)
(NA)
20.243
0.046
0.424
Analytes with a partial SRC p-value of v0.05 are shown. Benjamini-Hochberg corrected p-values were calculated to take into account the comparisons against all RBM
analytes.
doi:10.1371/journal.pone.0044260.t002
99
Figure 1. Analyte correlations. Heatmaps of partial SRC between RBM analytes significantly associated (p-value v0.05) with brain amyloid
burden, taking into account: age, gender, years of education, number of APOE E 4 alleles and the difference, in days, between plasma sampling and
[11 C]-PiB PET scan date. RBM analytes have been ordered by hierarchical clustering, the final cut-off is shown in purple. Variables that have been
grouped are shown in blue boxes, and the representative for each group is shown in a red box.
doi:10.1371/journal.pone.0044260.g001
100
PLOS ONE | www.plosone.org
Discussion
In this study, fibrillar amyloid beta levels, in ADNI subjects,
have been compared and related to the level of analytes on the
RBM panel in plasma. Brain amyloid burden appears to be
distributed bimodally in the RBM-PiB PET cohort, as has been
previously reported for a larger ADNI subcohort by Ewers et al.
[21]. Associations of C-peptide, fibrinogen, alpha-1-antitrypsin,
pancreatic polypeptide, complement C3, vitronectin, von willebrand factor, cortisol, serum amyloid p-component, AXL
receptor tyrosine kinase, interleukin-3, interleukin-13, matrix
metallproteinase-9, APOE, leptin and immunoglobulin E with
brain amyloid burden have been found in this study. Some of these
markers of brain amyloid burden were also found to associate with
other AD related phenotypes, such as CSF Ab1{42 , MRI features,
cognitive tests and diagnostic groups. In regression models it was
found that models including both RBM analytes and co-variates
performed better than those using only co-variate information,
suggesting that the RBM panel of analytes can be used as markers
of brain amyloid burden. Combining highly correlated variables
was found to reduce overfitting and led to a set of 13 RBM
analytes that together with co-variates could explain w30% of the
variance of brain amyloid burden and predict PiB positive
individuals with a high sensitivity. This result, and the increased
predictive accuracy of the 13 RBM model in comparison to using 101
n=0
ADNI-RBM
1.79 [0.200]
ANM + KHPDCR
1.91 [9.75610
n=1
1.66 [0.203]
22
Kruskal-Wallis x2 Pvalue
n=2
,2.20610216
1.52 [0.233]
1.88 [9.00610
22
1.82 [7.00610
22
,2.20610216
ANM + KHPDCR
controls
1.9 [0.100]
1.88 [8.0061022]
1.83 [3.0061022]
1.5761024
102
ADNI Data
Data used in the preparation of this article were obtained from
the ADNI database (adni.loni.ucla.edu). The ADNI was launched
in 2003 by the National Institute on Aging (NIA), the National
Institute of Biomedical Imaging and Bioengineering (NIBIB), the
Food and Drug Administration (FDA), private pharmaceutical
companies and non-profit organisations, as a $60 million, 5- year
public-private partnership. The primary goal of ADNI has been to
test whether serial magnetic resonance imaging (MRI), positron
emission tomography (PET), other biological markers, and clinical
and neuropsychological assessment can be combined to measure
the progression of mild cognitive impairment (MCI) and early
Alzheimers disease (AD). Determination of sensitive and specific
markers of very early AD progression is intended to aid
researchers and clinicians to develop new treatments and monitor
their effectiveness, as well as lessen the time and cost of clinical
trials.
The Principal Investigator of this initiative is Michael W.
Weiner, MD, VA Medical Center and University of California
San Francisco. ADNI is the result of efforts of many coinvestigators from a broad range of academic institutions and
private corporations, and subjects have been recruited from over
50 sites across the U.S. and Canada. The initial goal of ADNI was
to recruit 800 adults, ages 55 to 90, to participate in the research,
approximately 200 cognitively normal older individuals to be
followed for 3 years, 400 people with MCI to be followed for 3
years and 200 people with early AD to be followed for 2 years. For
up-to-date information, see www.adni-info.org.
Demographic (age, gender, years of education), genetic (number
of APOE E 4 alleles), diagnosis (control, MCI or AD at a given date)
and analyte (metabolite/protein/complex) levels in plasma were
compared with brain amyloid burden (or other markers of AD
pathology such as: CSF Ab1{42 level, MRI features, cognitive
scores or diagnostic groups). Diagnoses were recorded for each
subject at each visit. Plasma and CSF were collected from fasted
subjects using the procedures described previously [29,40]. Levels
of 190 analytes were measured from subject plasma using the
Rules Based Medicine (RBM, rulesbasedmedicine.com, Austin,
TX) Human Discovery Multi-Analyte Profile (MAP) 1.0 panel and
a Luminex 100 platform [41]. Measurements of 44 analytes were
excluded on the basis of quality control, leaving 146 analytes in the
subsequent analysis. The levels of all 146 analytes (except:
apolipoprotein H, complement factor H, E selectin, epidermal
growth factor, fibrinogen, interleukin-12 subunit p40, placenta
growth factor, serum glutamic oxaloacetic transaminase and 103
7
Linear regression
Linear regression was performed using the lm function in the R
stats package. This was appropriate because although many
variables were not distributed normally, the residuals of the
regression models used were approximately. Before regression,
measurements of each RBM analyte were transformed to a
standard deviation of one to allow each analyte to have equal
influence on the model (but not transformed to a mean of zero, to
make the analysis more comparable with that used in Thambisetty
et al. (2010) [14]). LOO CV was performed by fitting linear
regression models to the data, leaving out one subject at a time,
and using the model to predict brain amyloid burden in that
subject based on the fitted model. LOO CV R2 was calculated as
the square of the Pearsons correlation coefficient, calculated using
cor.test, between the predicted and observed brain amyloid
burden.
A cut-off of 1.5 was used to dichotomise brain amyloid burden
in the RBM-PiB PET cohort as PiB positive (49 subjects) or
negative (22 subjects), as previously suggested by Jagust et al.
(2010) [16]. Predicted values were similarly dichotomised.
Sensitivity and specificity of this prediction was calculated in R
using epi.test in the epiR package [50].
Permutation tests were performed by permuting brain amyloid
burden across subjects in the data, and re-calculating the LOO
CV R2 that was achieved by fitting the regression model to the
resulting data. 100,000 permutations were used, and the number
of cases in which the LOO CV R2 exceeded that achieved with
the original dataset was recorded. Relative importance was
calculated using the LMG score [51], using the R package
relaimpo [52].
Correlation analysis
104
Supporting Information
(PDF)
Acknowledgments
We are grateful to peer-reviewers, whose suggestions helped us to refine
this manuscript.
Data used in preparation of this article were obtained from the
Alzheimers Disease Neuroimaging Initiative (ADNI) database (adni.loni.ucla.edu). As such, the investigators within the ADNI contributed to the
design and implementation of ADNI and/or provided data but did not
participate in analysis or writing of this report. A complete listing of ADNI
investigators can be found at: http://adni.loni.ucla.edu/wp-content/
uploads/how_to_apply/ADNI_Acknowledgement_List.pdf.
by diagnostic group. P-values were calculated when appropriate for differences across diagnostic groups, using a Kruskal-Wallis
x2 test for continuous characteristics and simulated contingency
table p-values for discrete characteristics.
(PDF)
Author Contributions
Conceived and designed the experiments: M. Thambisetty AS SL SN RD.
Analyzed the data: SJK. Wrote the paper: SJK M. Thambisetty AS SN
RD. APOE measurements in AddNeuroMed and DCR: JRC AH IP MW
MKL KL. Input on image analysis: EW. Provided cluster computing
support: CJ. AddNeuroMed clinical centre lead: HS IK M. Tsolaki BV
PM.
References
1. Golde TE, Schneider LS, Koo EH (2008) Anti-Ab therapeutics in Alzheimers
disease: the need for a paradigm shift. Neuron 69: 203213.
2. Reiman EM, Langbaum JB, Fleisher AS, Caselli RJ, Chen K, et al. (2011)
Alzheimers Prevention Initiative: a plan to accelerate the evaluation of
presymptomatic treatments. J Alzheimers Dis 26: 321329.
3. Aisen PS, Andrieu S, Sampaio C, Carrillo M, Khachaturian ZS, et al. (2011)
Report of the task force on designing clinical trials in early (predementia) AD.
Neurology 76: 280286.
4. Grimmer T, Riemenschneider M, Forstl H, Henriksen G, Klunk WE, et al.
(2009) Beta amyloid in Alzheimers disease: increased deposition in brain is
reected in reduced concentration in cerebrospinal fluid. Biol Psychiatry 65: 927
934.
5. Koivunen J, Pirttila T, Kemppainen N, Aalto S, Herukka SK, et al. (2008) PET
amyloid ligand [11C] PIB uptake and cerebrospinal uid beta-amyloid in mild
cognitive impairment. Dement Geriatr Cogn Disord 26: 378383.
6. Zhang R, Barker L, Pinchev D, Marshall J, Rasamoelisolo M, et al. (2004)
Mining biomarkers in human sera using proteomic tools. Proteomics 4: 244
256.
7. Hye A, Lynham S, Thambisetty M, Causevic M, Campbell J, et al. (2006)
Proteome-based plasma biomarkers for Alzheimers disease. Brain 129: 3042
3050.
8. Soares HD, Potter WZ, Pickering E, Kuhn M, Immermann FW, et al. (Ahead of
print) Plasma biomarkers associated with the apolipoprotein e genotype and
alzheimer disease. Arch Neurol.
9. Doecke JD, Laws SM, Faux NG, Wilson W, Burnham SC, et al. (Ahead of print)
Blood-based protein biomarkers for diagnosis of Alzheimer disease. Arch
Neurol.
10. Thambisetty M, Simmons A, Velayudhan L, Hye A, Campbell J, et al. (2010)
Association of plasma clusterin concentration with severity, pathology, and
progression in Alzheimer disease. Arch Gen Psychiatry 67: 739748.
11. Thambisetty M, Simmons A, Hye A, Campbell J, Westman E, et al. (2011)
Plasma biomarkers of brain atrophy in Alzheimers disease. PloS one 6: e28527.
12. Thambisetty M, An Y, Kinsey A, Koka D, Saleem M, et al. (2012) Plasma
clusterin concentration is associated with longitudinal brain atrophy in mild
cognitive impairment. Neuroimage 59: 212217.
13. Thambisetty M, Hye A, Foy C, Daly E, Glover A, et al. (2008) Proteome-based
identification of plasma proteins associated with hippocampal metabolism in
early Alzheimers disease. J Neurol 255: 17121720.
14. Thambisetty M, Tripaldi R, Riddoch-Contreras J, Hye A, An Y, et al. (2010)
Proteome-based plasma markers of brain amyloid-b deposition in non-demented
older individuals. J Alzheimers Dis 22: 10991109.
15. Gupta VB, Laws SM, Villemagne VL, Ames D, Bush AI, et al. (2011) Plasma
apolipoprotein E and Alzheimer disease risk: the AIBL study of aging.
Neurology 76: 10911098.
16. Jagust WJ, Bandy D, Chen K, Foster NL, Landau SM, et al. (2010) The
Alzheimers Disease Neuroimaging Initiative positron emission tomography
core. Alzheimers Dement 6: 221229.
17. Slooter AJ, de Knijff P, Hofman A, Cruts M, Breteler MM, et al. (1998) Serum
apolipoprotein E level is not increased in Alzheimers disease: the Rotterdam
study. Neurosci Lett 248: 2124.
18. Siest G, Bertrand P, Qin B, Herbeth B, Serot JM, et al. (2000) Apolipoprotein E
polymorphism and serum concentration in Alzheimers disease in nine European
centres: the ApoEurope study. Clin Chem Lab Med 38: 721730.
19. Evans AE, Zhang W, Moreel JF, Bard JM, Ricard S, et al. (1993)
Polymorphisms of the apolipoprotein B and E genes and their relationship to
plasma lipid variables in healthy Chinese men. Hum Genet 92: 191197.
20. Johnstone D, Milward EA, Berretta R, Moscato P, ADNI (2012) Multivariate
protein signatures of pre-clinical Alzheimers disease in the Alzheimers Disease
Neuroimaging Initiative (ADNI) plasma proteome dataset. PLoS one 7: e34341.
21. Ewers M, Insel P, Jagust WJ, Shaw L, Trojanowski JQ, et al. (2012) CSF
Biomarker and PIBPETDerived Beta-Amyloid Signature Predicts Metabolic,
Gray Matter, and Cognitive Changes in Nondemented Subjects. Cerebral
Cortex 22: 19932004.
22. Lunnon K, Ibrahim Z, Proitsi P, Lourdusamy A, Newhouse S, et al. (2012)
Mitochondrial Dysfunction and Immune Activation are Detectable in Early
Alzheimers Disease Blood. J Alz Dis 30: 685710.
23. Maier M, Peng Y, Jiang L, Seabrook TJ, Carroll MC, et al. (2008) Complement
C3 deficiency leads to accelerated amyloid b plaque deposition and
neurodegeneration and modulation of the microglia/macrophage phenotype
in amyloid precursor protein transgenic mice. J Neurosci 28: 63336341.
24. Wetterling T, Tegtmeyer KF (1994) Serum alpha 1-antitrypsin and alpha 2macroglobulin in Alzheimers and Binswangers disease. Clin Investig 72: 196
199.
25. Bemelmans KJ, Noort A, Rijk RD, Middelkoop HAM, Kempen GMJV, et al.
(2007) Plasma cortisol and norepinephrine in Alzheimers disease: opposite
relations with recall performance and stage of progression. Acta Neuropsychiatr
19: 231237.
26. Lee BK, Glass TA, Wand GS, McAtee MJ, Bandeen-Roche K, et al. (2008)
Apolipoprotein E genotype, cortisol, and cognitive function in communitydwelling older adults. Am J Psychiat 165: 14561464.
27. Ray S, Britschgi M, Herbert C, Takeda-Uchimura Y, Boxer A, et al. (2007)
Classification and prediction of clinical Alzheimers diagnosis based on plasma
signaling proteins. Nat Med 13: 13591362.
28. Lorenzl S, Albers DS, Relkin N, Ngyuen T, Hilgenberg SL, et al. (2003)
Increased plasma levels of matrix metalloproteinase-9 in patients with
Alzheimers disease. Neurochem Int 43: 191196.
29. OBryant SE, Xiao G, Barber R, Huebinger R, Wilhelmsen K, et al. (2011) A
Blood-Based Screening Tool for Alzheimers Disease That Spans Serum and
Plasma: Findings from TARC and ADNI. PLoS one 6: e28092.
30. Nishiyama E, Iwamoto N, Kimura M, Arai H (1996) Serum amyloid P
component level in Alzheimers disease. Dementia 7: 256259.
31. OBryant SE, Xiao G, Barber R, Reisch J, Doody R, et al. (2010) A serum
protein-based algorithm for the detection of Alzheimer disease. Arch Neurol 67:
10771081.
105
PLOS ONE | www.plosone.org
32. Mold M, Shrive A (2012) Serum amyloid P component accelerates the formation
and enhances the stability of amyloid fibrils in a physiologically significant
under-saturated solution of amyloid-Ab42. J Alzheimers Dis 29: 875881.
33. de la Monte SM, Wands JR (2005) Review of insulin and insulin-like growth
factor expression, signaling, and malfunction in the central nervous system:
relevance to Alzheimers disease. J Alzheimers Dis 7: 4561.
34. Zambrano A, Otth C, Mujica L, Concha II, Maccioni RB (2007) Interleukin-3
prevents neuronal death induced by amyloid peptide. BMC Neurosci 8: 82.
35. Greco SJ, Sarkar S, Johnston JM, Tezapsidis N (2009) Leptin regulates tau
phosphorylation and amyloid through AMPK in neuronal cells. Biochem Bioph
Res Co 380: 98104.
36. Szczepanik AM, Funes S, Petko W, Ringheim GE (2001) IL-4, IL-10 and IL-13
modulate Ab(1-42)-induced cytokine and chemokine production in primary
murine microglia and a human monocyte cell line. J Neuroimmunol 113: 4962.
37. Bell RD, Winkler EA, Singh I, Sagare AP, Deane R, et al. (2012) Apolipoprotein
E controls cerebrovascular integrity via cyclophilin A. Nature 485: 512516.
38. Kohavi R (1995) A study of cross-validation and bootstrap for accuracy
estimation and model selection. Proceedings of IJCAI95 2: 11371143.
39. Alzheimers Association (2004) Research consent for cognitively impaired adults:
recommendations for institutional review boards and investigators. Alzheimer
Dis Assoc Disord 18: 171175.
40. Vemuri P, Wiste HJ, Weigand SD, Shaw LM, Trojanowski JQ, et al. (2009)
MRI and CSF biomarkers in normal, MCI, and AD subjects: diagnostic
discrimination and cognitive correlations. Neurology 73: 287297.
41. Trojanowski JQ, Vandeerstichele H, Korecka M, Clark CM, Aisen PS, et al.
(2010) Update on the biomarker core of the Alzheimers Disease Neuroimaging
Initiative subjects. Alzheimers Dement 6: 230238.
42. Jagust WJ, Landau SM, Shaw LM, Trojanowski JQ, Koeppe RA, et al. (2009)
Relationships between biomarkers in aging and dementia. Neurology 73: 1193.
43. R Development Core Team (2011) R: A Language and Environment for
Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria.
URL http://www.R-project.org. ISBN 3-900051-07-0.
106
PLOS ONE | www.plosone.org
10
Abstract
The failure of current strategies to provide an explanation for controversial findings on the pattern of pathophysiological
changes in Alzheimers Disease (AD) motivates the necessity to develop new integrative approaches based on multi-modal
neuroimaging data that captures various aspects of disease pathology. Previous studies using [18F]fluorodeoxyglucose
positron emission tomography (FDG-PET) and structural magnetic resonance imaging (sMRI) report controversial results
about time-line, spatial extent and magnitude of glucose hypometabolism and atrophy in AD that depend on clinical and
demographic characteristics of the studied populations. Here, we provide and validate at a group level a generative
anatomical model of glucose hypo-metabolism and atrophy progression in AD based on FDG-PET and sMRI data of 80
patients and 79 healthy controls to describe expected age and symptom severity related changes in AD relative to a
baseline provided by healthy aging. We demonstrate a high level of anatomical accuracy for both modalities yielding
strongly age- and symptom-severity- dependant glucose hypometabolism in temporal, parietal and precuneal regions and a
more extensive network of atrophy in hippocampal, temporal, parietal, occipital and posterior caudate regions. The model
suggests greater and more consistent changes in FDG-PET compared to sMRI at earlier and the inversion of this pattern at
more advanced AD stages. Our model describes, integrates and predicts characteristic patterns of AD related pathology,
uncontaminated by normal age effects, derived from multi-modal data. It further provides an integrative explanation for
findings suggesting a dissociation between early- and late-onset AD. The generative model offers a basis for further
development of individualized biomarkers allowing accurate early diagnosis and treatment evaluation.
Citation: Dukart J, Kherif F, Mueller K, Adaszewski S, Schroeter ML, et al. (2013) Generative FDG-PET and MRI Model of Aging and Disease Progression in
Alzheimers Disease. PLoS Comput Biol 9(4): e1002987. doi:10.1371/journal.pcbi.1002987
Editor: Olaf Sporns, Indiana University, United States of America
Received December 11, 2012; Accepted January 27, 2013; Published April 4, 2013
Copyright: 2013 Dukart et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: JD is supported by the Synapsy NCCR of the Swiss National Science Foundation (http://www.nccr-synapsy.ch/). BD is supported by the Swiss National
Science Foundation (grant 320030_135679 and NCCR Synapsy), Foundation Parkinson Switzerland, Foundation Synapsys, Novartis Foundation for medical
biological research and Deutsche Forschungsgemeinschaft (Kfo 247). FK is supported by the Velux Stiftung. MLS is supported by the German Federal Ministry of
Education and Research (BMBF; German FTLD consortium, www.ftld.de/) and LIFE - Leipzig Research Center for Civilization Diseases at the University of Leipzig
(http://life.uni-leipzig.de/). LIFE is funded by means of the European Union, by the European Regional Development Fund (ERFD) and by means of the Free State
of Saxony within the framework of the excellence initiative. Data collection and sharing for this project was funded by the Alzheimers Disease Neuroimaging
Initiative (ADNI, adni.loni.ucla.edu) (National Institutes of Health Grant U01 AG024904). ADNI is funded by the National Institute on Aging, the National Institute of
Biomedical Imaging and Bioengineering, and through generous contributions from the following: Abbott, AstraZeneca AB, Bayer Schering Pharma AG, BristolMyers Squibb, Eisai Global Clinical Development, Elan Corporation, Genentech, GE Healthcare, GlaxoSmithKline, Innogenetics, Johnson and Johnson, Eli Lilly and
Co., Medpace, Inc., Merck and Co., Inc., Novartis AG, Pfizer Inc, F. Hoffman-La Roche, Schering-Plough, Synarc, Inc., as well as non-profit partners the Alzheimers
Association and Alzheimers Drug Discovery Foundation, with participation from the U.S. Food and Drug Administration. Private sector contributions to ADNI are
facilitated by the Foundation for the National Institutes of Health (www.fnih.org). The grantee organization is the Northern California Institute for Research and
Education, and the study is coordinated by the Alzheimers Disease Cooperative Study at the University of California, San Diego. ADNI data are disseminated by
the Laboratory for Neuro Imaging at the University of California, Los Angeles. This research was also supported by NIH grants P30 AG010129, K01 AG030514, and
the Dana Foundation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing Interests: The authors have declared that no competing interests exist.
* E-mail: [email protected]
" Please see the acknowledgments for further details of the Alzheimers Disease Neuroimaging Initiative.
Introduction
107
Author Summary
Establishing an accurate diagnosis of Alzheimers disease
has been a major challenge in the past decades. With an
increasing amount of studies aiming at detection and
validation of imaging biomarkers for this disease, many
apparently controversial findings have been reported over
the time. The failure of current strategies to provide a
consistent explanation for these differential findings
motivates the necessity to develop new integrative
approaches based on multi-modal data that capture
various aspects of disease pathology. In our study we
propose such a generative model providing a comprehensive approach towards integration of previously
published differential findings in early- and late-onset
AD. We believe that our analytical strategy not only
provides the link between imaging biomarkers and clinical
phenotype considering the effects of aging, but could also
lead to new areas of research in terms of creation of new,
individualized biomarkers for a more accurate diagnosis of
Alzheimers disease.
108
Methods
Subjects
To derive a generative model of age and symptom severity
related changes, we extracted from the Alzheimers disease
Neuroimaging Initiative (ADNI) database (www.adni-info.org)
sMRI and FDG-PET data from multiple centres of 80 patients
with a clinical diagnosis of AD and 79 healthy controls (Table 1). A
full list of subject and scan IDs used in this study is provided at the
following location: http://www.unil.ch/webdav/site/lren/shared/
Juergen/Overview_patandcon_MRIandPETdate.xlsx. For all subjects, follow-up evaluations were available for up to 5 years after
initial examination. Control subjects and AD patient groups were
matched for gender and age. A diagnosis of AD was based on
NINCDS/ARDRA criteria [23]. Exclusion criteria for the ADNI
data included the presence of any significant neurological disease
other than AD, history of head trauma followed by persistent
2
FDG-PET data
Image pre-processing
All data processing steps were carried out using the SPM5
software package (Statistical Parametric Mapping software:
http://www.fil.ion.ucl.ac.uk/spm/) implemented in Matlab 7.7
(MathWorks Inc., Sherborn, MA). The same pre-processing
algorithm was used for sMRI and FDG-PET data, as described
elsewhere [11]. This procedure includes co-registration and
interpolation of both FDG-PET and sMR images to an isotropic
resolution of 16161 mm3, bias correction for inhomogeneity
artefacts for sMRI data, segmentation of sMRI data into different
tissue classes (only the grey matter tissue class is used for further
analyses), and masking of non-GM voxels in FDG-PET data. PVE
correction using the modified Muller-Gartner method [21,24] was
in the PVElab software package [25] that is compatible with
SPM5 only. This procedure uses the segmented sMR images to
account for PVE and for potential atrophy effects in FDG-PET.
DARTEL (Diffeomorphic Anatomical Registration using Exponentiated Lie algebra) based on grey matter tissue probability
maps was used for spatial normalization of data to an average size
template created from all study participants [26]. Structural MR
images were additionally modulated to preserve the total amount
of signal from each region. The same deformation matrices used to
normalise sMRI scans to a template were used to co-register the
FDG-PET images. After spatial normalization anatomical regions
of all subjects were located at same location in the images.
Smoothing with a Gaussian kernel of 12 mm FWHM (full width at
half maximum) accounted for minor misalignment errors.. FDGPET data were intensity normalized to cerebellar mean [27] and
masked to avoid big edge effects. The cerebellar region was chosen
for intensity normalization of FDG-PET as it has been shown to
be a region of choice for intensity normalization which is
unaffected in healthy aging and early stages of AD when
correcting for PVE caused by atrophy [2729].
sMRI data
The sMRI dataset included standard T1-weighted images
obtained with different scanner types using a 3D MP-RAGE
(magnetization-prepared 180 degrees radio-frequency pulses and
rapid gradient-echo) sequence varying in TR and TE (repetition
and echo time) with an in-plane resolution of 1.2561.25 mm and
1.2 mm slice thickness acquired at 1.5T magnetic field strength.
All raw data were pre-processed to correct for distortion and B1
non-uniformity as described on the ADNI webpage (http://www.
loni.ucla.edu/ADNI/Data/ADNI_Data.shtml).
AD
T-test (df,t,p)
Number
79
80
Male/Female
41/38
40/40
Age (years)
75.864.9
75.767.0
157,0.1,.887
6287
5588
MMSE (score)
28.761.6
23.662.2
157,16.6,,.001
All statistical analyses were also carried out using the SPM5
software package and Matlab 7.7. The effect of aging in healthy
control subjects was estimated separately for FDG-PET and sMRI
with voxel-wise linear regressions. To obtain the healthy aging
component of our generative model we used the beta coefficients
of aging in healthy controls to simulate voxel-wise changes in both
imaging modalities for the age range 50 to 80 years (Figure 1a).
The estimated values at 50 years were used as a 100% baseline.
Estimated age-related changes for the whole age range for both
FDG-PET and sMRI were expressed as percent decreases relative
to this baseline.
109
Figure 1. Schematic representation of voxel-wise age- and symptom severity- related models. a) Schematic representation of agerelated changes in one voxel (in %) considering GM volume at the age of 50 years as baseline. b) Schematic representation of changes related to
healthy aging (black line) and age-related differences in AD (red line) in one voxel. Intersection age (dotted line) represents the age at which healthy
aging in this voxel becomes similar to changes observed in AD. The hinge in the red line (aging in AD) at the intersection point indicates that after the
intersection age, according to our assumption of the additive impact of AD related processes to healthy aging, the healthy aging model would apply
in AD patients as no pathological processes in terms of atrophy or glucose hypometabolism are longer observable after this time point. c) Decrease
(in %) in GM volume observed in an exemplary voxel in AD depending on the constellation of age and symptom severity (MMSE) relative to the
baseline provided by healthy aging (violet line). AD Alzheimers disease, GM grey matter, MMSE Mini Mental State Examination.
doi:10.1371/journal.pcbi.1002987.g001
110
cCon {cAD
:
bAD {bCon
Results
Subject demographic characteristics
AD patients and control subjects did not differ in age
[t(157) = 0.1;p = .887]. As expected MMSE differed significantly
between the groups [t(157) = 16.6;p,.001]. The comparison of
AD patients and control subjects in relation to sex showed no
statistical differences [x2(1) = 0.06;p = .811].
112
PLOS Computational Biology | www.ploscompbiol.org
Figure 2. Healthy aging related changes observed in MRI (a) and FDG-PET (b) considering the expression of GM volume and
glucose metabolism at the age of 50 years as baseline. FDG-PET [18F]fluorodeoxyglucose positron emission tomography, GM grey matter, MRI
structural magnetic resonance imaging.
doi:10.1371/journal.pcbi.1002987.g002
Discussion
In this study we demonstrate that a generative model captures
the anatomical and metabolic features associated with AD and
Figure 3. A linear model of age and MMSE related changes observed in AD in FDG-PET considering the healthy control group as
baseline. AD Alzheimers disease, FDG-PET [18F]fluorodeoxyglucose positron emission tomography, MMSE Mini Mental State Examination.
doi:10.1371/journal.pcbi.1002987.g003
113
Figure 4. A linear model of age and MMSE related changes observed in AD in MRI considering the healthy control group as
baseline. AD Alzheimers disease, MMSE Mini Mental State Examination, MRI structural magnetic resonance imaging.
doi:10.1371/journal.pcbi.1002987.g004
114
Figure 5. Voxel-wise intersections of healthy aging and changes observed in AD in MRI (a) and FDG-PET (b). Colour bars represent
the intersection age. AD Alzheimers disease, FDG-PET [18F]fluorodeoxyglucose positron emission tomography, MRI structural magnetic resonance
imaging.
doi:10.1371/journal.pcbi.1002987.g005
115
Figure 6. Positive linear (a) or quadratic (b) relationship observed between age and MMSE in MRI (blue) and FDG-PET (red) in AD
patients after removing the variance explained by heathy aging. Only clusters are shown exceeding a significance threshold of p = 0.001
uncorrected on voxel level and p = 0.05 FWE-corrected on cluster level. AD Alzheimers disease, FDG-PET [18F]fluorodeoxyglucose positron emission
tomography, MRI structural magnetic resonance imaging.
doi:10.1371/journal.pcbi.1002987.g006
Supporting Information
Text S1 Acknowledgment list for ADNI publications.
(PDF)
Acknowledgments
We thank Jessica Peters for preparing parts of the data used in this project.
Data used in preparation of this article were obtained from the
Alzheimers Disease Neuroimaging Initiative (ADNI) database (adni.loni.ucla.edu). As such, the investigators within the ADNI contributed to the
design and implementation of ADNI and/or provided data but did not
participate in analysis or writing of this report. A complete listing of ADNI
investigators can be found in Text S1.
Author Contributions
Conceived and designed the experiments: JD MLS RSJF BD. Performed
the experiments: JD. Analyzed the data: JD. Contributed reagents/
materials/analysis tools: JD FK KM SA. Wrote the paper: JD FK KM SA
MLS RSJF BD.
116
PLOS Computational Biology | www.ploscompbiol.org
10
References
1. Mielke R, Pietrzyk U, Jacobs A, Fink GR, Ichimiya A, et al. (1994) HMPAO
SPET and FDG PET in Alzheimers disease and vascular dementia: comparison
of perfusion and metabolic pattern. Eur J Nucl Med 21: 10521060.
2. Laakso MP, Soininen H, Partanen K, Helkala EL, Hartikainen P, et al. (1995)
Volumes of hippocampus, amygdala and frontal lobes in the MRI-based
diagnosis of early Alzheimers disease: correlation with memory functions.
J Neural Transm Park Dis Dement Sect 9: 7386.
3. Silverman DH, Gambhir SS, Huang HW, Schwimmer J, Kim S, et al. (2002)
Evaluating early dementia with and without assessment of regional cerebral
metabolism by PET: a comparison of predicted costs and benefits. J Nucl Med
43: 253266.
4. Davatzikos C, Resnick SM, Wu X, Parmpi P, Clark CM (2008) Individual
patient diagnosis of AD and FTD via high-dimensional pattern classification of
MRI. Neuroimage 41: 12201227.
5. Habeck C, Foster NL, Perneczky R, Kurz A, Alexopoulos P, et al. (2008)
Multivariate and univariate neuroimaging biomarkers of Alzheimers disease.
Neuroimage 40: 15031515.
6. Kloppel S, Stonnington CM, Chu C, Draganski B, Scahill RI, et al. (2008)
Automatic classification of MR scans in Alzheimers disease. Brain 131: 681
689.
7. Sadeghi N, Foster NL, Wang AY, Minoshima S, Lieberman AP, et al. (2008)
Automatic classification of Alzheimers disease vs. Frontotemporal dementia: A
spatial decision tree approach with FDG-PET. In: 5th IEEE International
Symposium on Biomedical Imaging: From Nano to Macro; 1417 May 2008;
Paris, France. Vols 14: 408411.
8. Haense C, Herholz K, Jagust WJ, Heiss WD (2009) Performance of FDG PET
for detection of Alzheimers disease in two independent multicentre samples
(NEST-DD and ADNI). Dement Geriatr Cogn Disord 28: 259266.
9. Schroeter ML, Stein T, Maslowski N, Neumann J (2009) Neural correlates of
Alzheimers disease and mild cognitive impairment: a systematic and
quantitative meta-analysis involving 1351 patients. Neuroimage 47: 11961206.
10. Morinaga A, Ono K, Ikeda T, Ikeda Y, Shima K, et al. (2010) A comparison of
the diagnostic sensitivity of MRI, CBF-SPECT, FDG-PET and cerebrospinal
fluid biomarkers for detecting Alzheimers disease in a memory clinic. Dement
Geriatr Cogn Disord 30: 285292.
11. Dukart J, Mueller K, Horstmann A, Barthel H, Moller HE, et al. (2011)
Combined Evaluation of FDG-PET and MRI Improves Detection and
Differentiation of Dementia. PLoS One 6: e18111.
12. Sakamoto S, Ishii K, Sasaki M, Hosaka K, Mori T, et al. (2002) Differences in
cerebral metabolic impairment between early and late onset types of Alzheimers
disease. J Neurol Sci 200: 2732.
13. Ishii K, Kawachi T, Sasaki H, Kono AK, Fukuda T, et al. (2005) Voxel-based
morphometric comparison between early- and late-onset mild Alzheimers
disease and assessment of diagnostic performance of z score images. AJNR
Am J Neuroradiol 26: 333340.
14. Frisoni GB, Pievani M, Testa C, Sabattoli F, Bresciani L, et al. (2007) The
topography of grey matter involvement in early and late onset Alzheimers
disease. Brain 130: 720730.
15. Franke K, Ziegler G, Kloppel S, Gaser C (2010) Estimating the age of healthy
subjects from T1-weighted MRI scans using kernel methods: exploring the
influence of various parameters. Neuroimage 50: 883892.
16. Good CD, Johnsrude IS, Ashburner J, Henson RN, Friston KJ, et al. (2001) A
voxel-based morphometric study of ageing in 465 normal adult human brains.
Neuroimage 14: 2136.
17. Baron JC, Chetelat G, Desgranges B, Perchey G, Landeau B, et al. (2001) In
vivo mapping of gray matter loss with voxel-based morphometry in mild
Alzheimers disease. Neuroimage 14: 298309.
18. Teipel SJ, Wegrzyn M, Meindl T, Frisoni G, Bokde AL, et al. (2012) Anatomical
MRI and DTI in the Diagnosis of Alzheimers Disease: A European Multicenter
Study. J Alzheimers Dis 31: S3347.
19. Dukart J, Schroeter ML, Mueller K, Initiative TAsDN (2011) Age Correction in
Dementia - Matching to a Healthy Brain. PLoS One 6: e22193.
20. Karow DS, McEvoy LK, Fennema-Notestine C, Hagler DJ, Jennings RG, et al.
(2010) Relative capability of MR imaging and FDG PET to depict changes
associated with prodromal and early Alzheimer disease. Radiology 256: 932
942.
21. Muller-Gartner HW, Links JM, Prince JL, Bryan RN, McVeigh E, et al. (1992)
Measurement of radiotracer concentration in brain gray matter using positron
22.
23.
24.
25.
26.
27.
28.
29.
30.
31.
32.
33.
34.
35.
36.
37.
38.
39.
40.
41.
42.
43.
117
PLOS Computational Biology | www.ploscompbiol.org
11
Abstract
Alzheimers disease (AD) is characterized by neurofibrillary tangle and neuropil thread deposition, which ultimately results in
neuronal loss. A large number of magnetic resonance imaging studies have reported a smaller hippocampus in AD patients
as compared to healthy elderlies. Even though this difference is often interpreted as atrophy, it is only an indirect
measurement. A more direct way of measuring the atrophy is to use repeated MRIs within the same individual. Even though
several groups have used this appropriate approach, the pattern of hippocampal atrophy still remains unclear and difficult
to relate to underlying pathophysiology. Here, in this longitudinal study, we aimed to map hippocampal atrophy rates in
patients with AD, mild cognitive impairment (MCI) and elderly controls. Data consisted of two MRI scans for each subject.
The symmetric deformation field between the first and the second MRI was computed and mapped onto the threedimensional hippocampal surface. The pattern of atrophy rate was similar in all three groups, but the rate was significantly
higher in patients with AD than in control subjects. We also found higher atrophy rates in progressive MCI patients as
compared to stable MCI, particularly in the antero-lateral portion of the right hippocampus. Importantly, the regions
showing the highest atrophy rate correspond to those that were described to have the highest burden of tau deposition.
Our results show that local hippocampal atrophy rate is a reliable biomarker of disease stage and progression and could also
be considered as a method to objectively evaluate treatment effects.
Citation: Franko E, Joly O, for the Alzheimers Disease Neuroimaging Initiative (2013) Evaluating Alzheimers Disease Progression Using Rate of Regional
Hippocampal Atrophy. PLoS ONE 8(8): e71354. doi:10.1371/journal.pone.0071354
Editor: Karl Herholz, University of Manchester, United Kingdom
Received March 13, 2013; Accepted June 28, 2013; Published August 12, 2013
Copyright: 2013 Franko et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The authors of this study have only been funded by the University of Caen and INSERM. Data collection and sharing for this project was funded by the
Alzheimers Disease Neuroimaging Initiative (ADNI) (National Institutes of Health Grant U01 AG024904). ADNI was funded by the National Institute on Aging, the
National Institute of Biomedical Imaging and Bioengineering, and through generous contributions from the following: Abbott, AstraZeneca AB, Bayer Schering
Pharma AG, Bristol-Myers Squibb, Eisai Global Clinical Development, Elan Corporation, Genentech, GE Healthcare, GlaxoSmithKilne, Innogenetics, Johnson and
Johnson, Eli Lilly and Co., Medpace, Inc., Merck and Co., Inc., Novartis AD, Pfizer Inc, F. Hoffman-La Roche, Schering-Plough, Synarc, Inc., as well as non-profit
partners the Alzheimers Association and Alzheimers Drug Discovery Foundation, with participation from the U.S. Food and Drug Administration. Private sector
contributions to ADNI are facilitated by the Foundation for the National Institutes of Health (www.fnih.org). The grantee organization is the Northern California
Institute for Research and Education, and the study is coordinated by the Alzheimers Disease Cooperative Study at the University of California, San Diego. ADNI
data are disseminated by the Laboratory for Neuro Imaging at the University of California, Los Angeles. This research was also supported by NIH grants P30
AG010129, K01 AG030514, and the Dana Foundation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of
the manuscript.
Competing Interests: The authors have declared that no competing interests exist.
* E-mail: [email protected]
Introduction
118
3.9 ms, flip angle 8, and slice thickness 1.2 mm. In-plane
resolution differed slightly among subjects: 0.9460.94 mm,
1.2561.25 mm, or 1.361.3 mm.In this study, we used two scans
from each subject. MRI_acq1 was the baseline scan (first scan of
the subject), and MRI_acq2 was the last available scan (through
December 2011). All subjects included in the analyses had a
duration of at least 200 days between the first and last scans.
Hippocampal Segmentation
A template hippocampus (right hemisphere) was manually
segmented from the single-subject MNI-T1 template available in
SPM (www.fil.ion.ucl.ac.uk/spm). The segmentation followed the
description of Franko et al. [23] and Insausti and Amaral [24] and
contains the hippocampus proper (cornu ammonis (CA) 13
fields), the dentate gyrus and the subiculum. For guided extraction
of the subjects hippocampus in MRI_acq1, two sets of seven
points were defined in the template and in the MRI_acq1,
respectively. The locations of the points are similar to those used
by Pluta et al. [25]; three points were placed on coronal slices and
four points on sagittal slices. The first point on the coronal view
marks the appearance of the hippocampus inferior to the
amygdala, the second marks the most medial point of the
hippocampus at the level of the hippocampal-amygdaloid
transitional area, and the third point indicates the most posterior
part of the hippocampus [23]. On the sagittal view, we placed the
first point on the most lateral slice of the hippocampus, the second
and the third were placed 23 mm more medial on the anterior
and posterior borders of hippocampus, respectively, and the last
point was placed at the end of the intralimbic gyrus lateral to the
hippocampal fissure. Finally, to extract the 3D hippocampus from
MRI_acq1, the Iterative Closest Point (ICP) algorithm (www.vtk.
org) was used to compute the affine transformation between the
template and the subjects hippocampus (Fig. 1A).
Methods
Subjects
Data used in the preparation of this article were obtained from
the Alzheimers Disease Neuroimaging Initiative (ADNI) database
(adni.loni.ucla.edu), for up-to-date information, see www.adniinfo.org. The data were analysed anonymously, using publicly
available secondary data from the ADNI study, therefore no ethics
statement is required for this work.
From the ADNI database, a preliminary dataset containing AD
and normal subjects were downloaded first while developing the
method [22]. From this first dataset, we selected the subjects who
had at least 2 scans with at least 200 days apart. This resulted in
90 AD patients and 54 control subjects. Later we expanded the
dataset to have similar number of controls as patients and added
the MCI group. From the downloaded subjects, we included the
first about 90 subjects who fulfilled the inclusion criteria, namely
having at least two MRIs with at least 200 days apart. Few control
subjects, who turned to MCI or AD during the follow-up period,
were then removed from the statistical tests. The final number of
subjects included in this study were 85 healthy controls, 102
patients diagnosed with MCI and 90 patients with AD. An
overview of the subject groups is provided in Table 1, including
the total number of subjects, their age, sex, mini-mental state
examination (MMSE) and clinical dementia rating (CDR) scores.
The interscan intervals are also shown for each group. We further
divided the MCI group into progressive MCI (pMCI, n = 39),
containing subjects who converted to AD between the baseline
and the last scan (CDR.0.5), and stable MCI (sMCI, n = 44),
containing subjects with CDR scores of 0.5 at the time of the last
scan. The subjects, whose CDR score was not available at the time
of the last scan, were not included into this analysis. Almost all
patients with AD (88 out of 90) received medication (cholinesterase
inhibitors and NMDA-receptor antagonists), as did 78% of the
patients with MCI (80 out of 102) whereas none of the controls.
Measurements of Atrophy
Following the recommendation of Yushkevich et al. [26], we
performed an unbiased symmetric deformation field estimation to
measure the hippocampal atrophy rate. ANTS software (www.
picsl.upenn.edu/ANTS) was used to compute the symmetric
deformation field (SyN method) between MRI_acq1 and
MRI_acq2 using a normalised cross-correlation metric. The
displacement of each voxel in MRI_acq1 was written in the
WarpX, WarpY, and WarpZ images of the deformation fields.
The parameters of our deformation-based morphometry included
isotropic 2 mm (FWHM Gaussian kernel) image smoothing.
Finally, for each vertex of the hippocampal surface (MRI_acq1),
the dot product between the normal and the deformation field
defined the signed displacement of the vertex (in mm) (Fig. 1B).
The resulting values (divided by the duration between acquisition
MRI_acq1 and MRI_acq2) represent the atrophy rate in mm/
year. The texture was displayed on the 3D hippocampus in the
radiological convention with visualisation software (anatomist,
brainvisa.info/). The 3D mesh together with statistical maps is also
available: brainsenses.x10host.com/hc.htm. The hippocampal
volume was computed from the volumetric meshes. The technique
described here was partly published in abstract form with
preliminary results [22].
MRI Acquisitions
Statistical Analyses
Group
Age in years
MMSE
CDR
76 (5)
29.2 (1)
0 (0)
1217 (303)
75 (7)
26.8 (1.8)
0.5 (0)
1030 (393)
75 (7)
23.3 (2.2)
0.83 (0.35)
652 (188)
Mean values are followed by standard deviations (SD). MMSE: mini-mental state examination; CDR: clinical dementia rating; F: female; M:male.
doi:10.1371/journal.pone.0071354.t001
in controls in both hippocampi (right: W = 5390, p-value = 1.501e06; left: W = 4690, p-value = 0.004927). The volume loss was
significantly higher in patients with MCI than in controls in the
right hemisphere but not the left (right: W = 5273, p-value = 0.005484; left: W = 4719, p-value = 0.149 n.s.). Similarly,
the AD group showed significantly greater volume loss than the
MCI group in the right but not in the left hippocampus (right:
W = 5507,
p-value = 0.008535;
left:
W = 5190,
p-value = 0.05936 n.s.). No significant difference in volume change
was found between progressive and stable MCI (right: W = 1025,
p-value = 0.06436 n.s.; left: W = 996, p-value = 0.1048 n.s.). Collectively, the volumetric measurements revealed significant loss in
all groups and a larger loss in the AD and the MCI groups than in
controls, but they failed to show significant difference between
pMCI and sMCI.
signed-rank test, which does not require the normality of the data.
For both type of tests, significance level (p,0.05) was corrected for
multiple comparison (number of vertices) using Bonferroni
correction.
Hippocampal Subfields
Even though the hippocampal subfields are determined by
histology and their boundaries are not visible on MRI, the
subfields are often indicated on in vivo MRIs based solely on
visual cues derived from histological images. However, a more
reliable and unbiased subfield delineation can be achieved by
using an independent MRI atlas. Hence, to help localise significant
atrophy rates on the hippocampus, we projected the hippocampal
subfields as defined from the high-resolution atlas of the human
hippocampus [27] computed from post-mortem MRI at 9.4 Tesla
(www.nitrc.org/projects/pennhippoatlas) onto our template using
ICP algorithm. The borders of the subfields including CA1, CA2
3, the dentate gyrus, and the subiculum are illustrated as outlines
on the 3D surface of the right hippocampus in Fig. 2A.
Group Differences
The significantly higher rate of atrophy in patients with AD
than in controls (Fig. 4A and online 3D mesh, see Methods) was
found mainly in the medial part of the head and body and along
the lateral side of the hippocampi (Wilcoxon test, pv0.05 corr.). A
higher atrophy rate in the MCI group as compared to controls
(Fig. 4B) was found in much smaller regions, mainly on the lateral
side of the hippocampal head. Finally, significant differences
between the AD and MCI groups (Fig. 4C) were found in a few
vertices along the lateral side.
Results
120
Figure 1. Illustration of the methods. (A) On the left side of this panel, we illustrate the triangulated surface of the template hippocampus and
the manually defined seven points (in red). On the right side, the seven points drawn on the subjects hippocampus are shown in blue, and the ICP
algorithm is illustrated which minimises the sum of distances (di) between the red and blue points. (B) Illustration of the deformation field mapping
onto the triangulated hippocampal surface resulting in different colours for inward (blue) and outward (orange) deformation.
doi:10.1371/journal.pone.0071354.g001
Discussion
In the present study, we used two MRI scans to map the rate of
hippocampal atrophy in patients with AD and MCI and in healthy
elderly controls. We also measured the total hippocampal volume
121
Figure 2. Mean atrophy rate. Mean atrophy rate in mm/year for the left and right hippocampi in patients with AD (A), MCI (B) and controls (C). On
the right hippocampus, projected borders of subfields CA1, CA23, dentate gyrus (DG) and subiculum (Sub) are indicated with grey outlines.
doi:10.1371/journal.pone.0071354.g002
122
PLOS ONE | www.plosone.org
Group
2*Control (n = 85)
2402 (476)/
29 (30)/
2059 (404)
26 (22)
2181 (468)/
220 (31)/
1946 (431)
210 (30)
2129 (446)/
228 (41)/
1805 (397)
214 (34)
2*MCI (n = 102)
2*AD (n = 90)
Average hippocampal volume in mm3 (SD) at the first scan, average volume loss in mm3/year between the two scans for right and left (R/L) hippocampi (SD), and the
statistical tests on volume loss (Wilcoxon test) in each group.
doi:10.1371/journal.pone.0071354.t002
and compared its decrease among the three groups. We found the
lowest baseline volume in the AD group, and the greatest in the
group of healthy elderlies. The average hippocampal volume of
patients with MCI was in between the other two groups. Looking
at the volume loss over time, we found the fastest loss in the AD
group. This was significantly higher than the loss in the normal
group both in the left and the right hippocampi. The difference in
volume loss between the AD and MCI groups was significant only
in the right side. Similarly, only the right hippocampus showed
greater volume loss in the MCI group than in controls.
Comparison of the atrophy rates among the three subject groups
revealed similar patterns. The regions with the highest rate of
atrophy were located on the medial side of the hippocampal head
and body and along the lateral side. These regions also showed
significantly higher atrophy rates in patients with AD compared
with controls. The difference in hippocampal atrophy rate
between MCI and controls was located at the lateral side of the
right hippocampal head, whereas the higher atrophy rate in AD
compared to MCI was mainly found on the lateral side of the right
hippocampal body and tail. To relate these regions to the
histological subparts of the hippocampus, we used the atlas of the
human hippocampus segmented from high-resolution MRIs [27].
Based on this atlas, the fastest atrophy occurred in the CA1 zone.
Similarly, the strongest differences between the groups were
localised in the CA1 zone and in the subiculum.
Our volumetric measurements are in accordance with previous
studies measuring the hippocampal volume in different populations, namely sequential decrease in volume among healthy
elderlies, MCI and AD patients [29]. We also found smaller
volume in the left side compared to the right, as did many groups
previously [3032]. Our results are also in line with studies
examining the volume loss over time [11,31,3337]. They
reported the greatest annual percent change in the hippocampal
volume in patients with AD. This was significantly higher than
that in patients with MCI and in normal subjects. Leung et al. [37]
showed that the volume loss was even accelerating in the MCI
group over time.
Although the difference in hippocampal volume had the same
trend in many studies previously, the volume itself differed
substantially. Groups, using the ADNI database of patients with
AD and MCI and healthy elderly controls, reported hippocampal
volumes in AD varying between 1600 [35] and 3000 mm3 [20]. In
the present study, the hippocampal volume of AD patients was
about 2000 mm3. The remarkable difference among studies might
derive from the segmentation of the template hippocampus used to
extract the structure in each individual. There are different
recommendations for hippocampal segmentation which include
different amount of the subiculum and the tail of the hippocampus
PLOS ONE | www.plosone.org
Figure 3. Significant atrophy rates. Statistical maps showing the -log(p-values) for significant (pv0.05 corr.) atrophy rates for the left and right
hippocampi in patients with AD (A), MCI (B) and controls (C). Significance level (pv0.05 corr.) is indicated in the colour bar.
doi:10.1371/journal.pone.0071354.g003
124
PLOS ONE | www.plosone.org
Figure 4. Statistical maps of atrophy rate differences among groups. Statistical maps showing significantly higher atrophy rate in AD versus
control (A), MCI versus control (B), and AD versus MCI (C). Significance level (pv0.05 corr.) is indicated in the colour bar.
doi:10.1371/journal.pone.0071354.g004
125
Figure 5. ROI analysis of atrophy rates in progressive and stable MCI. Bar plots of atrophy rates in the left and right hippocampal ROIs
(shown in Figure 4) for pMCI (n = 39) and sMCI (n = 44). *: pv0.05, **: pv0.001 Bonferroni adjusted p-values.
doi:10.1371/journal.pone.0071354.g005
126
along the middle of the superior and inferior surfaces. Wang et al.
[32,39] compared patients with MCI (CDR = 0.5) to healthy
controls. They found more inward deformation in patients on the
superior side of the head and the inferior and lateral parts of the
hippocampal body and tail. They also reported the inward
deformation of these regions to predict conversion from healthy
state to MCI [49]. Apostolova et al. [42] compared subjects who
remained cognitively normal with those who converted to MCI or
AD and showed significant differences in hippocampal atrophy
rate in all parts of the hippocampus without subfield specificity
(CA13 subfields and subiculum).
Although Gerardin et al. [40] used a single scan to examine the
hippocampal shape, the topography of atrophy found in the
present study is the most like their results. They found
morphological differences between the hippocampi of patients
with AD and those of controls, predominantly on the medial part
of the head and the lateral part of the body and tail. However, the
use of a single scan prevented the investigation of atrophy
evolution. Our data suggest that the morphological differences
among groups in their study may be the consequence of different
atrophy rates. Costafreda et al. [43] also used a single scan to
predict conversion from MCI to AD using a baseline atrophic
phenotype. Similarly to our findings, the atrophy of the anterior
part of the CA1 field best predicted the conversion, with 80%
accuracy. We also found the lateral side of the hippocampal head
(corresponding to CA1) to be the best predictor of the progression
to AD from MCI based on the atrophy rate. Although the
accuracy is slightly lower, our method can be used to identify
patients at high risk of progression. This can help identify patients
most likely to be helped by treatments that are mainly effective in
the early phase of AD.
We found similar patterns of atrophy rate among the three
subject groups which might suggest that certain hippocampal
subregions are more prone to atrophy caused by normal or
pathological ageing. However, the current method cannot
Acknowledgments
The authors are indebted to T.E Cope for comments on earlier versions of
the manuscript.
Data used in preparation of this article were obtained from the
Alzheimers Disease Neuroimaging Initiative (ADNI) database (adni.loni.ucla.edu). As such, the investigators within the ADNI contributed to the
design and implementation of ADNI and/or provided data but did not
participate in analysis or writing of this report. A complete listing of ADNI
investigators can be found at: http://adni.loni.ucla.edu/wp-content/
uploads/how_to_ apply/ADNI_Acknowledgement_List.pdf.
Author Contributions
Conceived and designed the experiments: EF OJ. Performed the
experiments: EF OJ. Analyzed the data: EF OJ. Contributed reagents/
materials/analysis tools: EF OJ. Wrote the paper: EF OJ.
References
20. Morra JH, Tu Z, Apostolova LG, Green AE, Avedissian C, et al. (2009)
Automated 3d mapping of hippocampal atrophy and its clinical correlates in 400
subjects with alzheimers disease, mild cognitive impairment, and elderly
controls. Hum Brain Mapp 30: 276688.
21. Morra JH, Tu Z, Apostolova LG, Green AE, Avedissian C, et al. (2009)
Automated mapping of hippocampal atrophy in 1-year repeat mri data from 490
subjects with alzheimers disease, mild cognitive impairment, and elderly
controls. Neuroimage 45: S315.
22. Franko E, Joly O (2010) Mapping changes in the hippocampal shape over time
in alzheimers disease. Program No 75620 2010 Neuroscience Meeting Planner
San Diego, CA: Society for Neuroscience, 2010 Online.
23. Franko E, Insausti AM, Artacho-Perula E, Insausti R, Chavoix C (2012)
Identification of the human medial temporal lobe regions on magnetic resonance
images. Hum Brain Mapp : In Press, DOI: 10.1002/hbm.22170.
24. Insausti R, Amaral D (2004) Hippocampal formation. The human nervous
system: 871914.
25. Pluta J, Avants BB, Glynn S, Awate S, Gee JC, et al. (2009) Appearance and
incomplete label matching for diffeomorphic template based hippocampus
segmentation. Hippocampus 19: 56571.
26. Yushkevich PA, Avants BB, Das SR, Pluta J, Altinay M, et al. (2010) Bias in
estimation of hippocampal atrophy using deformation-based morphometry
arises from asymmetric global normalization: an illustration in adni 3 t mri data.
Neuroimage 50: 43445.
27. Yushkevich PA, Avants BB, Pluta J, Das S, Minkoff D, et al. (2009) A highresolution computational atlas of the human hippocampus from postmortem
magnetic resonance imaging at 9.4 t. Neuroimage 44: 38598.
28. Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, et al. (2011) proc: an opensource package for r and s+ to analyze and compare roc curves. BMC
Bioinformatics 12: 77.
29. Convit A, Leon MJD, Tarshish C, Santi SD, Tsui W, et al. (1997) Specific
hippocampal volume reductions in individuals at risk for alzheimers disease.
Neurobiol Aging 18: 131138.
30. Jack CR, Twomey CK, Zinsmeister AR, Sharbrough FW, Petersen RC, et al.
(1989) Anterior temporal lobes and hippocampal formations: normative
volumetric measurements from mr images in young adults. Radiology 172:
549554.
31. Laakso MP, Lehtovirta M, Partanen K, Riekkinen PJ, Soininen H (2000)
Hippocampus in alzheimers disease: a 3-year follow-up mri study. Biol
Psychiatry 47: 557561.
32. Wang L, Swank JS, Glick IE, Gado MH, Miller MI, et al. (2003) Changes in
hippocampal volume and shape across time distinguish dementia of the
alzheimer type from healthy aging. Neuroimage 20: 66782.
33. Jack CR, Shiung MM, Gunter JL, OBrien PC, Weigand SD, et al. (2004)
Comparison of different mri brain atrophy rate measures with clinical disease
progression in ad. Neurology 62: 591600.
34. Fox NC, Black RS, Gilman S, Rossor MN, Griffth SG, et al. (2005) Effects of
abeta immunization (an1792) on mri measures of cerebral volume in alzheimer
disease. Neurology 64: 15631572.
35. Schuff N,Woerner N, Boreta L, Kornfield T, Shaw LM, et al. (2009) Mri of
hippocampal volume loss in early alzheimers disease in relation to apoe
genotype and biomarkers. Brain 132: 10671077.
36. Macdonald KE, Bartlett JW, Leung KK, Ourselin S, Barnes J, et al. (2012) The
value of hippocampal and temporal horn volumes and rates of change in
predicting future conversion to ad. Alzheimer Dis Assoc Disord.
37. Leung KK, Bartlett JW, Barnes J, Manning EN, Ourselin S, et al. (2013)
Cerebral atrophy in mild cognitive impairment and alzheimer disease: Rates and
acceleration. Neurology 80: 648654.
10
127
38. Konrad C, Ukas T, Nebel C, Arolt V, Toga AW, et al. (2009) Defining the
human hippocampus in cerebral magnetic resonance imagesan overview of
current segmentation protocols. Neuroimage 47: 118595.
39. Wang L, Miller JP, Gado MH, McKeel DW, Rothermich M, et al. (2006)
Abnormalities of hippocampal surface structure in very mild dementia of the
alzheimer type. Neuroimage 30: 5260.
40. Gerardin E, Chetelat G, Chupin M, Cuingnet R, Desgranges B, et al. (2009)
Multidimensional classification of hippocampal shape features discriminates
alzheimers disease and mild cognitive impairment from normal aging.
Neuroimage 47: 147686.
41. Wang L, Harms MP, Staggs JM, Xiong C, Morris JC, et al. (2010) Donepezil
treatment and changes in hippocampal structure in very mild alzheimer disease.
Arch Neurol 67: 99106.
42. Apostolova LG, Mosconi L, Thompson PM, Green AE, Hwang KS, et al. (2010)
Subregional hippocampal atrophy predicts alzheimers dementia in the
cognitively normal. Neurobiol Aging 31: 107788.
43. Costafreda SG, Dinov ID, Tu Z, Shi Y, Liu CY, et al. (2011) Automated
hippocampal shape analysis predicts the onset of dementia in mild cognitive
impairment. Neuroimage 56: 2129.
44. Gutman B, Wang Y, Morra J, Toga AW, Thompson PM (2009) Disease
classification with hippocampal shape invariants. Hippocampus 19: 572578.
45. Frisoni GB, Sabattoli F, Lee AD, Dutton RA, Toga AW, et al. (2006) In vivo
neuropathology of the hippocampal formation in ad: a radial mapping mr-based
study. Neuroimage 32: 104110.
46. Frisoni GB, Ganzola R, Canu E, Rb U, Pizzini FB, et al. (2008) Mapping local
hippocampal changes in alzheimers disease and normal ageing with mri at 3
tesla. Brain 131: 32663276.
47. Pluta J, Yushkevich P, Das S, Wolk D (2012) In vivo analysis of hippocampal
subfield atrophy in mild cognitive impairment via semi-automatic segmentation
of t2-weighted mri. J Alzheimers Dis 31: 8599.
48. Lorenzi M, Ayache N, Frisoni GB, Pennec X, for the Alzheimers Disease
Neuroimaging Initiative (ADNI) (2013) Lcc-demons: A robust and accurate
diffeomorphic registration algorithm. Neuroimage.
49. Csernansky JG, Wang L, Swank J, Miller JP, Gado M, et al. (2005) Preclinical
detection of alzheimers disease: hippocampal shape and volume predict
dementia onset in the elderly. Neuroimage 25: 783792.
50. Leemput KV, Bakkour A, Benner T, Wiggins G, Wald LL, et al. (2009)
Automated segmentation of hippocampal subfields from ultra-high resolution
in vivo mri. Hippocampus 19: 549557.
51. Chupin M, Mukuna-Bantumbakulu AR, Hasboun D, Bardinet E, Baillet S, et al.
(2007) Anatomically constrained region deformation for the automated
segmentation of the hippocampus and the amygdala: Method and validation
on controls and patients with alzheimers disease. Neuroimage 34: 9961019.
52. Chupin M, Grardin E, Cuingnet R, Boutet C, Lemieux L, et al. (2009) Fully
automatic hippocampus segmentation and classification in alzheimers disease
and mild cognitive impairment applied on data from adni. Hippocampus 19:
579587.
128
PLOS ONE | www.plosone.org
11
130