Seun Jeon1†
Jae Myeong Kang2†
Seongho Seo3
Hye Jin Jeong4
Thomas Funck1
Sang-Yoon Lee3Kee Hyung Park5Yeong-Bae Lee5Byeong Kil Yeon2
Tatsuo Ido4
Nobuyuki Okamura6
Alan C. Evans1
Duk L. Na7,8*
Young Noh5,9*- 1McGill Centre for Integrative Neuroscience, Montreal Neurological Institute, McGill University, Montreal, QC, Canada
- 2Department of Psychiatry, Gil Medical Center, Gachon University College of Medicine, Incheon, South Korea
- 3Department of Neuroscience, Gachon University College of Medicine, Incheon, South Korea
- 4Neuroscience Research Institute, Gachon University, Incheon, South Korea
- 5Department of Neurology, Gil Medical Center, Gachon University College of Medicine, Incheon, South Korea
- 6Division of Pharmacology, Faculty of Medicine, Tohoku Medical and Pharmaceutical University, Sendai, Japan
- 7Department of Neurology, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, South Korea
- 8Neuroscience Center, Samsung Medical Center, Seoul, South Korea
- 9Department of Health Science and Technology, GAIHST, Gachon University, Incheon, South Korea
Alzheimer’s disease (AD) patients are known to have heterogeneous clinical presentation and pathologic patterns. We hypothesize that AD dementia can be categorized into subtypes based on multimodal imaging biomarkers such as magnetic resonance imaging (MRI), tau positron emission tomography (PET), and amyloid PET. We collected 3T MRI, 18F-THK5351 PET, and 18F-flutemetamol (FLUTE) PET data from 83 patients with AD dementia [Clinical Dementia Rating (CDR) ≤1] and 60 normal controls (NC), and applied surface-based analyses to measure cortical thickness, THK5351 standardized uptake value ratio (SUVR) and FLUTE SUVR for each participant. For the patient group, we performed an agglomerative hierarchical clustering analysis using the three multimodal imaging features on the vertices (n = 3 × 79,950). The identified AD subtypes were compared to NC using general linear models adjusting for age, sex, and years of education. We mapped the effect size within significant cortical regions reaching a corrected p-vertex <0.05 (random field theory). Our surface-based multimodal framework has revealed three distinct subtypes among AD patients: medial temporal-dominant subtype (MT, n = 44), parietal-dominant subtype (P, n = 19), and diffuse atrophy subtype (D, n = 20). The topography of cortical atrophy and THK5351 retention differentiates between the three subtypes. In the case of FLUTE, three subtypes did not show distinct topographical differences, although cortical composite retention was significantly higher in the P type than in the MT type. These three subtypes also differed in demographic and clinical features. In conclusion, AD patients may be clustered into three subtypes with distinct topographical features of cortical atrophy and tau deposition, although amyloid deposition may not differ across the subtypes in terms of topography.
Introduction
Alzheimer’s disease (AD) is a neurodegenerative disease characterized by histopathologic lesions of amyloid plaques and neurofibrillary tangles composed of the β-amyloid (Aβ) protein and paired helical filaments of hyperphosphorylated tau protein, respectively (Hyman and Trojanowski, 1997). These two neuropathologic hallmarks of AD are known to have typical spreading patterns. The Aβ accumulation has been known to progress from the neocortex to the brainstem (Thal et al., 2002) and neurofibrillary tangle pathology initially appears in the transentorhinal region and spreads to the limbic area, association cortices, and finally the primary cortices (Braak and Braak, 1991). Structural lesions including hippocampal and medial temporal cortical atrophy are also known to be specific to AD and can be used to screen and track the progression of AD (Scheltens et al., 1992; Frisoni et al., 2008). An A/T/N classification has been accepted for the description of multidomain biomarker findings for amyloid, tau, and neurodegeneration (Jack et al., 2016).
Previous studies have suggested anatomical and neuropathologic heterogeneity in AD. A postmortem study has found that neurofibrillary tangles can be a determinant of variability in AD (Murray et al., 2011). In terms of atrophic patterns, a voxel-based morphometry study has classified AD into four subgroups according to the regional atrophy (Shiino et al., 2006). It was reported that patterns of atrophy on magnetic resonance imaging (MRI) had three subtypes, which concomitantly correlated with pathological subtypes (Whitwell et al., 2012). Our group has also previously classified a large group of early AD dementia into three subtypes according to the regional cortical thickness (Noh et al., 2014). A recent tau positron emission tomography (PET) study using cluster analysis has reported variability of tau PET uptake in AD (Whitwell et al., 2018). In addition, these AD subtypes showed distinct clinical and demographic characteristics (Murray et al., 2011; Noh et al., 2014; Whitwell et al., 2018) and long-term disease progression (Na et al., 2016). Further investigation of the subtypes of AD dementia may facilitate a deeper understanding of its characteristics and progression.
With the development of PET tracers, researchers now can observe tau pathologies in vivo. Recent studies have shown that tau PET tracers such as 18F-THK5351 and 18F-AV-1451 significantly differentiate AD patients from old adults with normal cognition, reflect disease progression in AD, and correlates with neurofibrillary tangle retention (Cho et al., 2016; Schwarz et al., 2016; Schöll et al., 2016; Kang et al., 2017). Tau PET tracers also have been reported to bind to non-AD tauopathies, hyperphosphorylated 4R tau in tubular or straight filaments, in brain regions different from AD. Tau PET depositions have been found in patients with frontotemporal lobar degeneration such as progressive supranuclear palsy and corticobasal degeneration in the basal ganglia, thalamus, midbrain, and dentate nucleus (Chiotis et al., 2016; Kikuchi et al., 2016; Cho et al., 2017; Ishiki et al., 2017; Smith et al., 2017). A recent study found that 18F-AV-1451 deposition could differentiate dementia with Lewy bodies from AD (Kantarci et al., 2017). These studies show the utility of tau PET in evaluation of AD and non-AD pathologies.
However, off-target binding has been continuously reported for first-generation tau PET tracers. Binding affinity to β-sheet structures of 18F-THK5117 showed increased binding in the subcortical white matter (WM) retention (Harada et al., 2015). Studies have shown off-target bindings of 18F-AV-1451 in the choroid plexus due to the identification of the tau tangle-like structures (Ikonomovic et al., 2016; Johnson et al., 2016; Ossenkoppele et al., 2016). In particular, 18F-THK5351 has limited utility as a sole biomarker of AD-related tauopathy due to its binding to monoamine oxidase-B (MAO-B; Ng et al., 2017; Harada et al., 2018). A recent study has undertaken a cluster analysis based on the regional uptake of 18F-AV-1451 (Whitwell et al., 2018). Although a cautious interpretation is needed in tau PET studies due to these non-specific binding properties, tau PET provides valuable evidence of tau pathology in vivo.
Multidomain biomarker analyses based on neurodegeneration, tau, and amyloid together for AD subtypes may provide further insights into the subordinate characteristics of AD. AD subtypes have been previously defined in studies using cortical atrophy in MRI (Shiino et al., 2006; Noh et al., 2014; Whitwell et al., 2018), postmortem neurofibrillary tangle counts (Murray et al., 2011), and tau retention in PET scans (Whitwell et al., 2018), but it has not been evaluated with in vivo multimodal imaging scans. We sought to investigate whether AD dementia can be categorized into subgroups using the multimodal method comprising 3T MRI, tau PET, and amyloid PET, and whether clinical characteristics are associated with each subtype.
Materials and Methods
Participants
A total of 191 participants who had been clinically diagnosed with AD dementia or normal controls (NC) were prospectively recruited from March 2015 to November 2017. All participants underwent 3.0-Tesla MRI, 18F-THK5351 PET scans, and 18F-Flutemetamol (FLUTE) PET scans and completed neuropsychological tests at the Memory Clinic at Gachon University Gil Medical Center. Of the 191 participants, 37 participants with Clinical Dementia Rating (CDR) >1 were excluded from the study to avoid the effects of disease progression. Thus, 154 participants including patients with AD dementia with CDR ≤1 (n = 88) and NC (n = 66) were included in this study. Please note that total of 143 participants (AD dementia = 83 and NC = 60) were used since 11 participants were excluded in the quality control step as described in the “Materials and Methods” section.
AD dementia patients were recruited from memory disorder clinic at Samsung Medical Center or Gachon University Gil Medical Center and had been diagnosed with probable AD according to the National Institute of Neurological and Communicative Disorders and Stroke and the AD and Related Disorders Association (McKhann et al., 1984). Diagnoses were confirmed by follow-up for more than 1 year by a neurologist with more than 30 years of experience (DN) and a neurologist with more than 10 years clinical and research experience (YN). The AD patients were classified into early-onset AD (onset age <65) and late-onset AD (onset age ≥65). Patients were excluded if they had structural abnormalities in MRI such as intracranial hemorrhage, cerebral, cerebellar, or brainstem infarction, traumatic brain injury, hydrocephalus, tumors, severe WM hyperintensity, WM hyperintensity associated with radiation, multiple sclerosis, or vasculitis. Other causes of dementia were ruled out with laboratory tests such as complete blood count, folate levels, vitamin B12, thyroid function, metabolic profile, and syphilis serology. Patients with familial AD and vascular dementia were not included in the study. APOE4 genotyping for all participants was obtained.
The 66 participants in the NC group were either spouse of the patients or healthy volunteers from the community who did not have subjective memory complaints and objective cognitive decline. All of them had a CDR score of 0 and normal results on neuropsychological tests (defined as within 1.5 standard deviations of age- and education-corrected normative mean). Participants were excluded if they had structural abnormalities in MRI such as intracranial hemorrhage, cerebral, cerebellar, or brainstem infarction, traumatic brain injury, hydrocephalus, tumors, severe WM hyperintensity, WM hyperintensity associated with radiation, multiple sclerosis, or vasculitis.
Written informed consent was obtained from all participants and the study was approved by the Institutional Review Board of Gachon University Gil Medical Center.
Neuropsychological Assessment
Mini-Mental State Examination (MMSE), CDR, CDR-sum of boxes (CDR-SOB) results were obtained and detailed neuropsychological function tests including attention, praxis, frontal/executive function, visual and verbal memory, language, visuoconstructive ability, and elements of Gerstmann syndrome were evaluated in all participants. Detailed items of the comprehensive test battery (Kang and Na, 2003) have been described in our previous study (Lee et al., 2018).
Image Acquisition
All participants underwent brain MRI using a 3.0-T MRI scanner (Verio, Siemens with a Siemens matrix coil) as described in our previous study (Kang et al., 2017). Both 18F-THK5351 and 18F-FLUTE PET scans were acquired using a Siemens Biograph 6 Truepoint PET/computed tomography scanner (Siemens, Knoxville, TN, USA) with a list-mode emission acquisition. THK5351 scans were acquired for 20 min starting from 50 min after the injection of 185MBq of 18F-THK5351 intravenously (50–0 min), which was synthesized and radiolabeled in Gachon University Neuroscience Research Institute. 18F-FLUTE emission scans were acquired for 20 min starting from 90 min after the intravenous injection of 185 MBq of 18F-FLUTE (90–110 min), purchased from Carecamp Inc. The mean interval between the PET scans was 13.94 ± 14.02 days and detailed data of interval between the two PET scans are presented in Appendix 1 in Supplementary Materials. Attenuation correction was performed with a low-dose CT scan and data reconstruction was performed with a 2D ordered subset expectation maximization algorithm (eight iterations and 16 subsets).
Image Processing
Cortical Surface Reconstruction and Cortical Thickness Measurement
We followed the CIVET pipeline1 (version 2.1). Briefly, each subject’s T1-weighted image was corrected for intensity inhomogeneity and linearly registered to the Montreal Neurological Institute-152 template to bring the images into a common space (Collins et al., 1994). The images were then tissue classified into WM, gray matter (GM), or cerebrospinal fluid (CSF; Zijdenbos et al., 2002) and the inner (WM/GM boundary) and the outer (GM/CSF boundary) cortical surfaces were extracted resulting in 40,962 vertex points per hemisphere (Kim et al., 2005). To obtain vertex-correspondence between individuals, surfaces were registered to an unbiased group template by matching the sulcal folding pattern (Lyttelton et al., 2007). The registered surfaces were transformed back into each patient’s native space, and cortical thicknesses were calculated as the Euclidean distance between the linked vertices of the inner and outer surfaces (Lerch et al., 2005). The measured cortical thickness was smoothed using a 30 mm full width at half maximum (FWHM) surface-based diffusion smoothing kernel (Chung et al., 2003).
Surface-Based Measurement for 18F-THK5351 and 18F-FLUTE
We rigidly co-registered the PET scans to native T1-weighted images using mutual information as a cost function. The cortical surfaces and tissue classes were linearly registered into the PET scans by applying inverse transform matrices. We performed partial volume correction (PVC) within gray and WM regions using iterative deconvolution with a surface-based anatomically constrained filtering (idSURF) method that uses the representation of the volume between the inner and outer surfaces as a spatial constraint to the PET signal (Funck et al., 2014). The PVC images were normalized to the reference regions resulting in a standardized uptake value ratio (SUVR). We used cerebellum GM (Okamura et al., 2014; Lockhart et al., 2016) and pons (Thurfjell et al., 2014) as low receptor density reference regions in the THK5351 and FLUTE analyses, respectively. The SUVR signal intensities were sampled at 50% of the distance from the inner to the outer surface to minimize partial volume contamination. The measured signals were spatially blurred using a surface-based diffusion smoothing kernel (20 mm FWHM).
Image Quality Control
All raw images and the results produced from the pipeline were carefully verified (by two investigators blinded to participant information). We excluded five AD and six NC participants due to MRI motion artifacts and image processing errors in brain masking, tissue classification, and cortical surface extraction.
Cluster Analysis
We performed agglomerative hierarchical clustering analysis using z-scored multimodal imaging features without noncortical regions on the surface model (three features with 79,950 vertices each). The hierarchical clustering method combines pairs of clusters at each step while minimizing the sum of squared errors from the cluster mean (Ward, 1963). Each of the 83 patients with AD dementia was placed in their own cluster and then progressively clustered with others. The AD patients belonging to the same cluster had similar profiles, while those in the different cluster had different profiles. The dendrogram created by the surface-based multimodal cluster analysis is presented in Figure 1. To estimate the optimal cluster number, we used the Gap statistics package available in the R software (version 3.5.1, R Development Core Team). Gap compares changes in the total intra-cluster variation for the different number of clusters with the expected values under the null reference distribution of the data (i.e., a distribution with no obvious clustering; Tibshirani et al., 2001). The optimal cluster number was three which yielded the maximum Gap statistic (Appendix 2 in Supplementary Materials). The number of Monte Carlo bootstrap iterations for the computation was set to 2,000.
Figure 1. Dendrogram created by surface-based multimodal cluster analysis. AD, Alzheimer’s disease; MT, medial temporal-dominant; P, Parietal-dominant; D, Diffuse atrophy. The distance along the y-axis represents the measure of similarity between patients, such that the shorter the distance, the greater the similarity between patients. The green, blue and red lines represent the clustered subtypes of AD dementia. The three clusters were selected based on Gap statistics (see Appendix 2 in Supplementary Materials).
Regions of Interest-Based Measurement
Volume-weighted SUVR values of images were averaged for regions of interest (ROIs) predefined using the Desikan-Killiany-Tourville atlas (Desikan et al., 2006). The value of the FLUTE cortical composite ROI was calculated with cortical SUVRs including the frontal, parietal, lateral temporal, and anterior and posterior cingulate cortices (Thurfjell et al., 2014) and that of the THK5351 global ROI was based on AD-related regions. Detailed regions are presented in Appendix 3 in Supplementary Materials.
Statistical Analysis
To estimate topographical abnormalities in AD dementia subtypes, we applied general linear models and random field theory using the SurfStat toolbox (Worsley et al., 2009). The three subtypes were compared to the NC group using a general linear model adjusting for age, sex, and years of education. Intracranial volume was included in the cortical thickness analysis. We mapped the effect size (Cohen’s d, adjusted for the nuisance variables) within the significant region reaching a p-vertex < 0.05 (random field theory) on the standard cortical surface.
Group comparisons of demographic and clinical characteristics between AD and NC were performed using an independent t-test for continuous variables and chi-square test for categorical variables. In comparisons among the three AD subtypes, one-way analysis of variance followed by Bonferroni post hoc test was used (p < 0.05). Mean cortical thickness, THK SUVR, and FLUTE SUVR were compared between AD and NC using independent t-test and analysis of covariance (covariance: age, gender, and years of education) followed by pairwise comparisons for adjusted means (Bonferroni, p < 0.05). Region-wise multiple comparisons were corrected using the Benjamini-Hochberg false discovery rate method (Benjamini and Hochberg, 1995).
Results
Subtypes of Alzheimer’s Disease Based on Multimodal Cluster Analysis
Agglomerative hierarchical clustering analysis based on cortical thickness, THK5351 PET, and FLUTE PET yielded three subtypes in the AD patients (Figure 2): medial temporal-dominant subtype (MT, n = 44), parietal-dominant subtype (P, n = 19), and diffuse atrophy subtype (D, n = 20). In the MT subtype, cortical thinning was prominent in the medial temporal cortex. THK retention was most increased in the medial and lateral temporal cortices and slightly in the inferior parietal, precuneus-posterior cingulate (PC-PCC), and frontal cortices. FLUTE retention was identified in some regions of the dorsolateral prefrontal, medial frontal, lateral temporal, parietal cortices, and PC-PCC (Figures 2A–C). In the P subtype, cortical thinning was relatively distinct in the parietal, PC-PCC, and occipital cortices. THK retention was also dominant in the parietal cortex and PC-PCC and FLUTE retention was more pronounced in the parietal, PC-PCC, lateral temporal, frontal, and occipital cortices (Figures 2A–C). In the D subtype, cortical thinning was found in relatively diffuse cortices including the frontal, medial temporal, lateral temporal, inferior parietal, and PC-PCC. THK retention was found in very similar areas where cortical thinning was identified, and amyloid uptake was found in the diffuse brain cortices saving the primary sensorimotor cortex (Figures 2A–C).
Figure 2. Topographical maps of imaging biomarkers according to AD subtypes. AD, Alzheimer’s disease; MT, Medial temporal-dominant; P, parietal-dominant; D, diffuse atrophy. Comparison of regional (A) cortical thickness, (B) tau retention, and (C) amyloid uptake between identified subtypes and the control group. The color scale indicates the effect size (Cohen’s). Significance was determined based on general linear models controlling for age, sex, and years of education. Intracranial volume was added as a covariate to the cortical thickness analyses. The statistical maps were thresholded using a random field theory (p-vertex < 0.05).
The results of ROI-based comparisons of three imaging markers (cortical thickness, THK SUVR, and FLUTE SUVR) among the diagnostic groups are shown in Appendix 4 in Supplementary Materials. Cortical thinning in the superior parietal cortex was most prominent in the P subtype (2.66 ± 0.12 for MT, 2.29 ± 0.24 for P, and 2.56 ± 0.21, p < 0.001) but that in the mesial temporal cortex was similar across subtypes (2.80 ± 0.19 for MT, 2.79 ± 0.16 for P, and 2.82 ± 0.18, p = 0.831). THK retention in the mesial temporal cortex showed subtle differences across subtypes (2.16 ± 0.28 for MT, 2.15 ± 0.16 for P, and 2.09 ± 0.24, p = 0.594) but the P subtype showed dominant THK retention in the superior parietal cortex (1.36 ± 0.13 for MT, 1.98 ± 0.36 for P, and 1.54 ± 0.17, p < 0.001). In terms of FLUTE retention, there was no ROI that showed a significant difference between the MT and D subtypes. The P subtype showed greater FLUTE retention than the MT subtype in the inferior parietal areas only (Appendix Table A4-3). The FLUTE retention showed little topographical difference across the three subtypes while the cortical composite SUVR of the P type was significantly greater than that of the MT subtype (Table 1, Appendix Table A4-3).
Table 1. Demographic and clinical characteristics of the study population.
Meanwhile, clustering analyses based on cortical thickness and/or THK5351 showed less distinct classification than the present result (Appendix 5 in Supplementary Materials). The analyses with two or four subtypes also showed less significant group differences than the present result (Appendix 6 in Supplementary Materials).
Demographic and Clinical Characteristics
The demographics and clinical information showed distinct features among each subtype (Table 1). The MT subtype patients were older than the other subtypes (72.34 ± 9.37 for MT, 59.47 ± 7.21 for P, and 64.70 ± 8.07 for D, p < 0.001), and the percentage of females was greater [37 (84.1%) for MT, 13 (68.4%) for P, and 9 (45%) for D, p = 0.006]. The P subtype patients had the earliest onset age (68.75 ± 9.18 for MT, 56.12 ± 6.09 for P, and 61.27 ± 8.19 for D, p < 0.001) and were the youngest among the subtypes. There were no significant differences in MMSE or CDR-SOB among the three AD subtypes (p = 0.399 for MMSE and p = 0.094 for CDR-SOB; Table 1).
Neuropsychological Test Performances
The neuropsychological scores were different among each subtype (Table 2). The K-BNT and verbal memory (SVLT delayed recall and recognition) scores showed no significant difference between the subtypes. The P subtype patients showed poorer performance in attention (digit span backward), visuospatial function (RCFT copy), visual memory (RCFT immediate recall, delayed recall, and recognition), and frontal executive function (COWAT animal/supermarket/phonemic total, Stroop test color reading, and TMT-A/B) than the MT type patients. Patients in the D and P subtypes showed a similar decline in frontal executive function but patients in the D subtype showed better performance than those in the P subtype in the RCFT copy test.
Table 2. Neuropsychological test results for AD subtypes.
Discussion
In this study, we performed agglomerative hierarchical cluster analysis based on cortical thickness using 3T MRI, THK5351 PET, and FLUTE PET images in a group of patients with early stages of AD dementia and identified three subtypes: MT, P, and D subtypes. We also found heterogeneous topographical patterns of the three imaging biomarkers, demographic characteristics, and cognitive functions according to the AD subtypes.
Three subtypes in AD dementia were generated incorporating cortical atrophy, tau deposition, and amyloid deposition. The MT (53%), P (23%), and D (24%) subtypes are considered to be analogous to the three pathological subtypes, limbic-predominant (14.3%–21.6%), hippocampal-sparing (10.7%–19.3%), and typical AD (47.4%–70.6%) identified in previous pathologic (Murray et al., 2011; Janocko et al., 2012) and MRI (Whitwell et al., 2012) studies from another study group. Their recent tau PET study also identified three clusters according to AV-1451 uptake in the entorhinal and neocortices; entorhinal low/cortical low, entorhinal low/cortical high, entorhinal high/cortical high (Whitwell et al., 2018) and they largely corresponded to MT, P, and D subtypes in our study, respectively. The lower frequency of the D subtype (24%) in our study compared to the previous studies (47.4%–70.6%) may be attributable to the participant characteristics. There was a lower proportion of late-onset AD patients in our study participants compared to the general AD population, and the methods we used in this study are dependent on the study population.
Of the three subtypes, the P subtype was characterized as having a younger age and a lower tendency for APOE4 frequency compared to the MT and D subtypes (Table 1). These characteristics were similar to those of early-onset AD (Licht et al., 2007). The P subtype is in line with the hippocampal-sparing subtype in previous studies with features of younger age, atypical clinical presentation, and lower APOE4 frequency (Murray et al., 2011). In contrast, the MT subtype was characterized by older age, a higher percent of females, and a higher tendency for APOE4 frequency (Table 1). These results are similar with those found in previous studies showing higher alterations in the hippocampus of the female brain due to age-related estrogen reduction (Fester et al., 2012) and region-specific effects of APOE4 in the medial temporal lobe (Hashimoto et al., 2001).
A novel finding in our study is that the topographical heterogeneity was manifested based on the multimodal imaging biomarkers MRI, tau PET, and amyloid PET (Figure 2). Cortical atrophy and THK retention showed similar topography at each subtype. As for the MT subtype, cortical atrophy was observed in the medial temporal cortex and THK retention was also observed mainly in the medial temporal and basal forebrain regions which represent the earliest pathologic changes (Braak and Braak, 1991). In the P subtype, both cortical atrophy and tau retention were observed in the temporoparietal regions related to the subsequent state of AD (Braak and Braak, 1991). In the D subtype, cortical atrophy and tau retention were observed at comparable extent both in the medial temporal region and the diffuse neocortical regions. This parallel topography between tau PET and MRI scans according to AD subtypes is supported by earlier and recent studies. Autopsy studies have documented that cortical atrophy proceeds similarly to neurofibrillary tangle pathology (Braak and Braak, 1991; Whitwell et al., 2008) and in vivo imaging studies have also reported similar distribution between cortical atrophy and tau retention in patients with typical and atypical AD (Xia et al., 2017; Nasrallah et al., 2018; Whitwell et al., 2018). A previous study also suggested that CSF tau could not correspond to the cortical atrophic patterns because CSF results show only pooled information on tau in the whole brain (Hwang et al., 2015). Thus, our result may suggest that neurofibrillary tangle formation and cortical thinning can have distinct topographic patterns within AD and that MRI and tau PET scans show this considerably. On the other hand, amyloid PET findings showed only a higher uptake in the P subtype without distinct topographic patterns differentiating AD subtypes. This can be supported by the phases of Aβ deposition which are different from the cortical thickness (Chételat et al., 2010) and tau retention (Thal et al., 2002) in patients with AD. Aβ deposition is known to precede neurodegeneration and clinical decline (Jack et al., 2009) and does not correlate with cortical atrophy in AD, unlike tau retention (Josephs et al., 2008). Cortical composite SUVR for the P type was significantly greater than that of the MT subtype. This might be related to the characteristics of the subgroup population. The mean age of the MT subtype was older than the P subtype, and the main portion of patients in the MT subtype was late-onset AD. Among the late-onset AD patients, there have been known mixed pathologies in the brain such as TDP-43, argyrophilic grain disease (Ferrer et al., 2008; Landau et al., 2016), aging-related tau astrogliopathy, or hippocampal sclerosis (Nelson et al., 2011; Cairns et al., 2015).
Cognitive function differed according to the AD subtypes (Table 2). In the P subtype, visuoconstruction was significantly impaired compared to the MT and D subtypes and it is known that visuospatial cognition is predominantly mediated by parietal lobe function (Possin, 2010). Poorer results for attention, visual memory, and frontal executive function were also found in the P subtype compared to the MT subtype. The worst cognitive function and non-amnestic features in the P subtype may be characteristics of early-onset AD. However, the three subtypes did not differ in the language (K-BNT) and verbal memory tests (SVLT delayed recall and recognition).
We note several limitations to our study. We could not undertake pathologic confirmation because we only used imaging biomarkers. The cross-sectional design of the study without longitudinal follow-up is another limitation in that progression of each subtype is unknown. As the study participants were younger and had a greater proportion of early-onset AD compared to the general population, the distribution and characteristics of the AD subtypes could be affected. As noted, the limitation of the THK5351 tracer itself should be regarded in our study because THK5351 PET shows binding to MAO-B (Ng et al., 2017; Harada et al., 2018). THK5351 is known to trace not only neurofibrillary tangles but a combination of neurofibrillary tangles and reactive astrocytes. Although the analysis in this study did not include the subcortical structure that is most influenced by MAO-B availability, it should be carefully interpreted due to MAO-B distribution throughout the whole brain (Ng et al., 2017; Harada et al., 2018).
Conclusion
Our surface-based multimodal imaging cluster analysis framework has revealed three distinct subtypes among AD patients in terms of the distribution of cortical atrophy, THK5351 retention, and FLUTE retention. We used three crucial imaging biomarkers and identified the three subtypes of AD consistent with the previous pathologic or imaging studies and suggest that multimodal in vivo imaging biomarkers may differentiate the subtypes of AD, mainly by the tau deposition and cortical atrophic pattern. Future work will focus on the combinations of various biomarkers more specific to AD pathology and provide further evidence of the multifaceted basis of AD. Consideration for topographic heterogeneity may be important when planning future preventative and treatment strategies because the AD subtypes may have different courses of disease progression and different responses to treatment. In addition, since the cluster analysis algorithm is dependent on the characteristics of the participants, further evaluation in a large cohort is needed.
Data Availability
The data set generated and/or analyzed during the current study are available from the corresponding author, Prof. Young Noh on reasonable request.
Ethics Statement
Written informed consent was obtained from all participants and the study was approved by the Institutional Review Board of Gachon University Gil Medical Center.
Author Contributions
YN, DN, and AE conceptualized and designed the study. SJ and JK drafted the manuscript. JK, KP, Y-BL, BY, DN, and YN acquired the data. SJ, JK, SS, HJ, TF, AE, and YN analyzed the data. S-YL, TI, and NO also contributed to the PET data acquisition and analyses. YN revised the manuscript for intellectual content. All authors reviewed and approved for publication.
Funding
This study was supported by a grant from the Korea Healthcare Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health and Welfare, South Korea (grant No: HI14C1135), a grant from the Brain Research Program of the National Research Foundation (NRF) funded by the Korean government (MSIT; No. 2018M3C7A1056889) and Research Program 2016 funding from Seoul National University College of Medicine Research Foundation, and grants from the Canadian Institute of Health Research (CIHR; 201085 and 247003) and Brain Canada/FNC (238990) awarded to AE. SJ is the recipient of the following fellowships that contributed to his support: Bourse Fonds de Recherche Santé Québec (FRQS, dossier 34240 and 259605) and by Jeanne Timmins Costello Fellowship of the Montreal Neurological Institute (240522).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Footnotes
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2019.00211/full#supplementary-material
References
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Series B Stat. Methodol. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
Braak, H., and Braak, E. (1991). Neuropathological stageing of Alzheimer-related changes. Acta Neuropathol. 82, 239–259. doi: 10.1007/bf00308809
Cairns, N. J., Perrin, R. J., Franklin, E. E., Carter, D., Vincent, B., Xie, M., et al. (2015). Neuropathologic assessment of participants in two multi-center longitudinal observational studies: the Alzheimer Disease Neuroimaging Initiative (ADNI) and the Dominantly Inherited Alzheimer Network (DIAN). Neuropathology 35, 390–400. doi: 10.1111/neup.12205
Chételat, G., Villemagne, V. L., Bourgeat, P., Pike, K. E., Jones, G., Ames, D., et al. (2010). Relationship between atrophy and β-amyloid deposition in Alzheimer disease. Ann. Neurol. 67, 317–324. doi: 10.1002/ana.21955
Chiotis, K., Saint-Aubert, L., Savitcheva, I., Jelic, V., Andersen, P., Jonasson, M., et al. (2016). Imaging in-vivo tau pathology in Alzheimer’s disease with THK5317 PET in a multimodal paradigm. Eur. J. Nucl. Med. Mol. Imaging 43, 1686–1699. doi: 10.1007/s00259-016-3363-z
Cho, H., Choi, J. Y., Hwang, M. S., Kim, Y. J., Lee, H. M., Lee, H. S., et al. (2016). in vivo cortical spreading pattern of tau and amyloid in the Alzheimer disease spectrum. Ann. Neurol. 80, 247–258. doi: 10.1002/ana.24711
Cho, H., Choi, J. Y., Hwang, M. S., Lee, S. H., Ryu, Y. H., Lee, M. S., et al. (2017). Subcortical 18F-AV-1451 binding patterns in progressive supranuclear palsy. Mov. Disord. 32, 134–140. doi: 10.1002/mds.26844
Chung, M. K., Worsley, K. J., Robbins, S., Paus, T., Taylor, J., Giedd, J. N., et al. (2003). Deformation-based surface morphometry applied to gray matter deformation. Neuroimage 18, 198–213. doi: 10.1016/s1053-8119(02)00017-4
Collins, D. L., Neelin, P., Peters, T. M., and Evans, A. C. (1994). Automatic 3D intersubject registration of MR volumetric data in standardized talairach space. J. Comput. Assist. Tomogr. 18, 192–205. doi: 10.1097/00004728-199403000-00005
Desikan, R. S., Ségonne, F., Fischl, B., Quinn, B. T., Dickerson, B. C., Blacker, D., et al. (2006). An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage 31, 968–980. doi: 10.1016/j.neuroimage.2006.01.021
Ferrer, I., Santpere, G., and van Leeuwen, F. W. (2008). Argyrophilic grain disease. Brain 131, 1416–1432. doi: 10.1093/brain/awm305
Fester, L., Prange-Kiel, J., Zhou, L., Blittersdorf, B. V., Böhm, J., Jarry, H., et al. (2012). Estrogen-regulated synaptogenesis in the hippocampus: sexual dimorphism in vivo but not in vitro. J. Steroid Biochem. Mol. Biol. 131, 24–29. doi: 10.1016/j.jsbmb.2011.11.010
Frisoni, G. B., Ganzola, R., Canu, E., Rüb, U., Pizzini, F. B., Alessandrini, F., et al. (2008). Mapping local hippocampal changes in Alzheimer’s disease and normal ageing with MRI at 3 Tesla. Brain 131, 3266–3276. doi: 10.1093/brain/awn280
Funck, T., Paquette, C., Evans, A., and Thiel, A. (2014). Surface-based partial-volume correction for high-resolution PET. Neuroimage 102, 674–687. doi: 10.1016/j.neuroimage.2014.08.037
Harada, R., Ishiki, A., Kai, H., Sato, N., Furukawa, K., Furumoto, S., et al. (2018). Correlations of 18F-THK5351 PET with postmortem burden of tau and astrogliosis in Alzheimer disease. J. Nucl. Med. 59, 671–674. doi: 10.2967/jnumed.117.197426
Harada, R., Okamura, N., Furumoto, S., Furukawa, K., Ishiki, A., Tomita, N., et al. (2015). [18F]THK-5117 PET for assessing neurofibrillary pathology in Alzheimer’s disease. Eur. J. Nucl. Med. Mol. Imaging 42, 1052–1061. doi: 10.1007/s00259-015-3035-4
Hashimoto, M., Yasuda, M., Tanimukai, S., Matsui, M., Hirono, N., Kazui, H., et al. (2001). Apolipoprotein E ε4 and the pattern of regional brain atrophy in Alzheimer’s disease. Neurology 57, 1461–1466. doi: 10.1212/WNL.57.8.1461
Hwang, J., Kim, C. M., Jeon, S., Lee, J. M., Hong, Y. J., Roh, J. H., et al. (2015). Prediction of Alzheimer’s disease pathophysiology based on cortical thickness patterns. Alzheimers Dement 2, 58–67. doi: 10.1016/j.dadm.2015.11.008
Hyman, B. T., and Trojanowski, J. Q. (1997). Editorial on consensus recommendations for the postmortem diagnosis of Alzheimer disease from the national institute on aging and the reagan institute working group on diagnostic criteria for the neuropathological assessment of Alzheimer disease. J. Neuropathol. Exp. Neurol. 56, 1095–1097. doi: 10.1097/00005072-199710000-00002
Ikonomovic, M. D., Abrahamson, E. E., Price, J. C., Mathis, C. A., and Klunk, W. E. (2016). [F-18]AV-1451 positron emission tomography retention in choroid plexus: more than “off-target” binding. Ann. Neurol. 80, 307–308. doi: 10.1002/ana.24706
Ishiki, A., Harada, R., Okamura, N., Tomita, N., Rowe, C. C., Villemagne, V. L., et al. (2017). Tau imaging with [18F]THK-5351 in progressive supranuclear palsy. Eur. J. Neurol. 24, 130–136. doi: 10.1111/ene.13164
Jack, C. R., Bennett, D. A., Blennow, K., Carrillo, M. C., Feldman, H. H., Frisoni, G. B., et al. (2016). A/T/N: An unbiased descriptive classification scheme for Alzheimer disease biomarkers. Neurology 87, 539–547. doi: 10.1212/wnl.0000000000002923
Jack, C. R. Jr., Lowe, V. J., Weigand, S. D., Wiste, H. J., Senjem, M. L., Knopman, D. S., et al. (2009). Serial PIB and MRI in normal, mild cognitive impairment and Alzheimer’s disease: implications for sequence of pathological events in Alzheimer’s disease. Brain 132, 1355–1365. doi: 10.1093/brain/awp062
Janocko, N. J., Brodersen, K. A., Soto-Ortolaza, A. I., Ross, O. A., Liesinger, A. M., Duara, R., et al. (2012). Neuropathologically defined subtypes of Alzheimer’s disease differ significantly from neurofibrillary tangle-predominant dementia. Acta Neuropathol. 124, 681–692. doi: 10.1007/s00401-012-1044-y
Johnson, K. A., Schultz, A., Betensky, R. A., Becker, J. A., Sepulcre, J., Rentz, D., et al. (2016). Tau positron emission tomographic imaging in aging and early Alzheimer disease. Ann. Neurol. 79, 110–119. doi: 10.3410/f.725883031.793514740
Josephs, K. A., Whitwell, J. L., Ahmed, Z., Shiung, M. M., Weigand, S. D., Knopman, D. S., et al. (2008). β-amyloid burden is not associated with rates of brain atrophy. Ann. Neurol. 63, 204–212. doi: 10.1002/ana.21223
Kang, J. M., Lee, S. Y., Seo, S., Jeong, H. J., Woo, S. H., Lee, H., et al. (2017). Tau positron emission tomography using [18F]THK5351 and cerebral glucose hypometabolism in Alzheimer’s disease. Neurobiol. Aging 59, 210–219. doi: 10.1016/j.neurobiolaging.2017.08.008
Kang, Y., and Na, D. (2003). Seoul Neuropsychological Screening Battery. Seoul: Human Brain Research & Consulting Co.
Kantarci, K., Lowe, V. J., Boeve, B. F., Senjem, M. L., Tosakulwong, N., Lesnick, T. G., et al. (2017). AV-1451 tau and β-amyloid positron emission tomography imaging in dementia with Lewy bodies. Ann Neurol 81, 58–67. doi: 10.1002/ana.24825
Kikuchi, A., Okamura, N., Hasegawa, T., Harada, R., Watanuki, S., Funaki, Y., et al. (2016). In vivo visualization of tau deposits in corticobasal syndrome by 18F-THK5351 PET. Neurology 87, 2309–2316. doi: 10.1212/wnl.0000000000003375
Kim, J. S., Singh, V., Lee, J. K., Lerch, J., Ad-Dab’bagh, Y., MacDonald, D., 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, 210–221. doi: 10.1016/j.neuroimage.2005.03.036
Landau, S. M., Horng, A., Fero, A., Jagust, W. J., and Alzheimer’s Disease Neuroimaging. (2016). Amyloid negativity in patients with clinically diagnosed Alzheimer disease and MCI. Neurology 86, 1377–1385. doi: 10.1212/wnl.0000000000002576
Lee, H., Seo, S., Lee, S. Y., Jeong, H. J., Woo, S. H., Lee, K. M., et al. (2018). [18F]-THK5351 PET imaging in patients with semantic variant primary progressive aphasia. Alzheimer. Dis. Assoc. Disord. 32, 62–69. doi: 10.1097/WAD.0000000000000216
Lerch, J. P., Pruessner, J. C., Zijdenbos, A., Hampel, H., Teipel, S. J., and Evans, A. C. (2005). Focal decline of cortical thickness in Alzheimer’s disease identified by computational neuroanatomy. Cereb. Cortex 15, 995–1001. doi: 10.1093/cercor/bhh200
Licht, E. A., McMurtray, A. M., Saul, R. E., and Mendez, M. F. (2007). Cognitive differences between early- and late-onset Alzheimer’s disease. Am. J. Alzheimers Dis. Other Demen. 22, 218–222. doi: 10.1177/1533317506299156
Lockhart, S. N., Baker, S. L., Okamura, N., Furukawa, K., Ishiki, A., Furumoto, S., et al. (2016). Dynamic PET measures of tau accumulation in cognitively normal older adults and Alzheimer’s disease patients measured using [18F] THK-5351. PLoS One 11:e0158460. doi: 10.1371/journal.pone.0158460
Lyttelton, O., Boucher, M., Robbins, S., and Evans, A. (2007). An unbiased iterative group registration template for cortical surface analysis. Neuroimage 34, 1535–1544. doi: 10.1016/j.neuroimage.2006.10.041
McKhann, G., Drachman, D., Folstein, M., Katzman, R., Price, D., and Stadlan, E. M. (1984). Clinical diagnosis of Alzheimer’s disease: report of the NINCDS-ADRDA Work Group under the auspices of Department of Health and Human Services Task Force on Alzheimer’s disease. Neurology 34, 939–939. doi: 10.1212/wnl.34.7.939
Murray, M. E., Graff-Radford, N. R., Ross, O. A., Petersen, R. C., Duara, R., and Dickson, D. W. (2011). Neuropathologically defined subtypes of Alzheimer’s disease with distinct clinical characteristics: a retrospective study. Lancet Neurol. 10, 785–796. doi: 10.1016/s1474-4422(11)70156-9
Na, H. K., Kang, D. R., Kim, S., Seo, S. W., Heilman, K. M., Noh, Y., et al. (2016). Malignant progression in parietal-dominant atrophy subtype of Alzheimer’s disease occurs independent of onset age. Neurobiol. Aging 47, 149–156. doi: 10.1016/j.neurobiolaging.2016.08.001
Nasrallah, I. M., Chen, Y. J., Hsieh, M. K., Phillips, J. S., Ternes, K., Stockbower, G. E., et al. (2018). 18F-flortaucipir PET/MRI correlations in nonamnestic and amnestic variants of Alzheimer disease. J. Nucl. Med. 59, 299–306. doi: 10.2967/jnumed.117.194282
Nelson, P. T., Schmitt, F. A., Lin, Y., Abner, E. L., Jicha, G. A., Patel, E., et al. (2011). Hippocampal sclerosis in advanced age: clinical and pathological features. Brain 134, 1506–1518. doi: 10.1093/brain/awr053
Ng, K. P., Pascoal, T. A., Mathotaarachchi, S., Therriault, J., Kang, M. S., Shin, M., et al. (2017). Monoamine oxidase B inhibitor, selegiline, reduces 18F-THK5351 uptake in the human brain. Alzheimers Res. Ther. 9:25. doi: 10.1186/s13195-017-0253-y
Noh, Y., Jeon, S., Lee, J. M., Seo, S. W., Kim, G. H., Cho, H., et al. (2014). Anatomical heterogeneity of Alzheimer disease: based on cortical thickness on MRIs. Neurology 83, 1936–1944. doi: 10.1212/wnl.0000000000001003
Okamura, N., Furumoto, S., Fodero-Tavoletti, M. T., Mulligan, R. S., Harada, R., Yates, P., et al. (2014). Non-invasive assessment of Alzheimer’s disease neurofibrillary pathology using 18F-THK5105 PET. Brain 137, 1762–1771. doi: 10.1093/brain/awu064
Ossenkoppele, R., Schonhaut, D. R., Schöll, M., Lockhart, S. N., Ayakta, N., Baker, S. L., et al. (2016). Tau PET patterns mirror clinical and neuroanatomical variability in Alzheimer’s disease. Brain 139, 1551–1567. doi: 10.1093/brain/aww027
Possin, K. L. (2010). Visual spatial cognition in neurodegenerative disease. Neurocase 16, 466–487. doi: 10.1080/13554791003730600
Scheltens, P., Leys, D., Barkhof, F., Huglo, D., Weinstein, H., Vermersch, P., et al. (1992). Atrophy of medial temporal lobes on MRI in “probable” Alzheimer’s disease and normal ageing: diagnostic value and neuropsychological correlates. J. Neurol. Neurosurg. Psychiatry. 55, 967–972. doi: 10.1136/jnnp.55.10.967
Schöll, M., Lockhart, S. N., Schonhaut, D. R., O’Neil, J. P., Janabi, M., Ossenkoppele, R., et al. (2016). PET imaging of tau deposition in the aging human brain. Neuron 89, 971–982. doi: 10.1016/j.neuron.2016.01.028
Schwarz, A. J., Yu, P., Miller, B. B., Shcherbinin, S., Dickson, J., Navitsky, M., et al. (2016). Regional profiles of the candidate tau PET ligand 18F-AV-1451 recapitulate key features of Braak histopathological stages. Brain 139, 1539–1550. doi: 10.1093/brain/aww023
Shiino, A., Watanabe, T., Maeda, K., Kotani, E., Akiguchi, I., and Matsuda, M. (2006). Four subgroups of Alzheimer’s disease based on patterns of atrophy using VBM and a unique pattern for early onset disease. Neuroimage 33, 17–26. doi: 10.1016/j.neuroimage.2006.06.010
Smith, R., Schain, M., Nilsson, C., Strandberg, O., Olsson, T., Hagerstrom, D., et al. (2017). Increased basal ganglia binding of 18F-AV-1451 in patients with progressive supranuclear palsy. Mov. Disord. 32, 108–114. doi: 10.1002/mds.26813
Thal, D. R., Rub, U., Orantes, M., and Braak, H. (2002). Phases of A beta-deposition in the human brain and its relevance for the development of AD. Neurology 58, 1791–1800. doi: 10.1212/wnl.58.12.1791
Thurfjell, L., Lilja, J., Lundqvist, R., Buckley, C., Smith, A., Vandenberghe, R., et al. (2014). Automated quantification of 18F-flutemetamol PET activity for categorizing scans as negative or positive for brain amyloid: concordance with visual image reads. J. Nucl. Med. 55, 1623–1628. doi: 10.2967/jnumed.114.142109
Tibshirani, R., Walther, G., and Hastie, T. (2001). Estimating the number of clusters in a data set via the gap statistic. J. Royal Stat. Soc. B 63, 411–423. doi: 10.1111/1467-9868.00293
Ward, J. H. Jr. (1963). Hierarchical grouping to optimize an objective function. J. Am. Stat. Assoc. 58, 236–244. doi: 10.2307/2282967
Whitwell, J. L., Josephs, K. A., Murray, M. E., Kantarci, K., Przybelski, S. A., Weigand, S. D., et al. (2008). MRI correlates of neurofibrillary tangle pathology at autopsy: a voxel-based morphometry study. Neurology 71, 743–749. doi: 10.1212/01.wnl.0000324924.91351.7d
Whitwell, J. L., Dickson, D. W., Murray, M. E., Weigand, S. D., Tosakulwong, N., Senjem, M. L., et al. (2012). Neuroimaging correlates of pathologically defined subtypes of Alzheimer’s disease: a case-control study. Lancet Neurol. 11, 868–877. doi: 10.1016/s1474-4422(12)70200-4
Whitwell, J. L., Graff-Radford, J., Tosakulwong, N., Weigand, S. D., Machulda, M., Senjem, M. L., et al. (2018). [18F]AV-1451 clustering of entorhinal and cortical uptake in Alzheimer’s disease. Ann. Neurol. 83, 248–257. doi: 10.1002/ana.25142
Worsley, K. J., Taylor, J., Carbonell, F., Chung, M., Duerden, E., Bernhardt, B., et al. (2009). SurfStat: A Matlab toolbox for the statistical analysis of univariate and multivariate surface and volumetric data using linear mixed effects models and random field theory. NeuroImage 47:S102. doi: 10.1016/s1053-8119(09)70882-1
Xia, C., Makaretz, S. J., and Caso, C. (2017). Association of in vivo [18F]AV-1451 tau pet imaging results with cortical atrophy and symptoms in typical and atypical Alzheimer disease. JAMA Neurol. 74, 427–436. doi: 10.1001/jamaneurol.2016.5755
Keywords: Alzheimer’s disease, cluster analysis, tau, amyloid, cortical thickness, positron emission tomography, magnetic resonance imaging
Citation: Jeon S, Kang JM, Seo S, Jeong HJ, Funck T, Lee S-Y, Park KH, Lee Y-B, Yeon BK, Ido T, Okamura N, Evans AC, Na DL and Noh Y (2019) Topographical Heterogeneity of Alzheimer’s Disease Based on MR Imaging, Tau PET, and Amyloid PET. Front. Aging Neurosci. 11:211. doi: 10.3389/fnagi.2019.00211
Received: 05 March 2019; Accepted: 26 July 2019;
Published: 20 August 2019.
Edited by:
Changiz Geula, Northwestern University, United StatesCopyright © 2019 Jeon, Kang, Seo, Jeong, Funck, Lee, Park, Lee, Yeon, Ido, Okamura, Evans, Na and Noh. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Duk L. Na, ZHVrbmFAc2trdS5lZHU=; Young Noh, eW5vaEBnYWNob24uYWMua3I=
† These authors have contributed equally to this work