Original Research ARTICLE
Association of Subcortical Structural Shapes With Tau, Amyloid, and Cortical Atrophy in Early-Onset and Late-Onset Alzheimer’s Disease
- 1School of Biomedical Engineering, Korea University, Seoul, South Korea
- 2Department of Psychiatry, Gil Medical Center, Gachon University College of Medicine, Incheon, South Korea
- 3Department of Neuroscience, College of Medicine, Gachon University, 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
- 6Department of Neurology, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul, South Korea
- 7Neuroscience Center, Samsung Medical Center, Seoul, South Korea
- 8Department of Health Science and Technology, GAIHST, Gachon University, Incheon, South Korea
- 9Department of Artificial Intelligence, Korea University, Seoul, South Korea
The objectives of this study were to compare the topographical subcortical shape and to investigate the effects of tau or amyloid burden on atrophic patterns in early onset Alzheimer’s disease (EOAD) and late-onset Alzheimer’s disease (LOAD). One hundred and sixty-one participants (53 EOAD, 44 LOAD, 33 young controls, and 31 older controls) underwent [18F]THK5351 positron emission tomography (PET), [18F]flutemetamol (FLUTE) PET, and 3T MRI scans. We used surface-based analysis to evaluate subcortical structural shape, permutation-based statistics for group comparisons, and Spearman’s correlations to determine associations with THK, FLUTE, cortical thickness, and neuropsychological test results. When compared to their age-matched controls, EOAD patients exhibited shape reduction in the bilateral amygdala, hippocampus, caudate, and putamen, while in LOAD patients, the bilateral amygdala and hippocampus showed decreased shapes. In EOAD, widespread subcortical shrinkage, with less association of the hippocampus, correlated with THK retention and cortical thinning, while in LOAD patients, subcortical structures were limited which had significant correlation with THK or mean cortical thickness. Subcortical structural shape showed less correlation with FLUTE global retention in both EOAD and LOAD. Multiple cognitive domains, except memory function, correlated with the bilateral amygdala, caudate, and putamen in EOAD patients, while more restricted regions in the subcortical structures were correlated with neuropsychological test results in LOAD patients. Subcortical structures were associated with AD hallmarks in EOAD. However, the correlation was limited in LOAD. Moreover, relationship between subcortical structural atrophy and cognitive decline were quite different between EOAD and LOAD. These findings suggest that the effects of Alzheimer’s pathologies on subcortical structural changes in EOAD and LOAD and they may have different courses of pathomechanism.
Alzheimer’s disease (AD) is a neurodegenerative disease associated with cognitive decline. Accumulations of neurofibrillary tangles and beta-amyloid (Aβ) are the two pathologic hallmarks of AD, as well as structural changes in the brain cortices (Jack et al., 2016). Early onset AD (EOAD) is characterized by onset before 65 years of age and accounts for 10–15% of all cases (Panegyres and Chen, 2014; Wattmo and Wallin, 2017). Studies have found that EOAD patients have typical neuropsychological features with more progressive cognitive decline (Seltzer and Sherwin, 1983; Smits et al., 2012) as well as heterogenous neuroimaging findings (Frisoni et al., 2007; Migliaccio et al., 2015).
Compared with LOAD patients, EOAD patients manifest with more widespread tau PET retention and heavier amyloid PET uptake in the brain cortices, and more severe gray matter loss suggesting different etiologies and predisposing factors for EOAD (Frisoni et al., 2007; Choo et al., 2011; Scholl et al., 2017; Wattmo and Wallin, 2017). Like gray matter change, more rapid volumetric declines in several subcortical structures including the caudate, putamen, and thalamus have been highlighted as distinct to LOAD patients (Cho et al., 2013; Pievani et al., 2013). It has also been reported that characteristic symptoms of EOAD such as severe extrapyramidal signs (Chui et al., 1985) and non-memory cognitive dysfunction such as executive function, visuospatial functioning, and attention (Smits et al., 2012) are related to deterioration of the basal ganglia and the thalamus, both of which play major roles in movement symptoms, brain connectivity, and memory function (Hahn et al., 2016; Dipasquale and Cercignani, 2017).
Although there have been previous investigations into the association between tau, Aβ, and atrophy in brain cortices, little is known of the relationship between the well-known AD biomarkers and subcortical structures. Unraveling the effect of tau, Aβ, and atrophy in the subcortical structures may contribute to a better understanding of the pathological mechanisms of AD and the different developmental paths of EOAD and LOAD. We sought to compare topographical changes in the subcortical structures in EOAD and LOAD patients and investigated their association with cortical tau and amyloid global retention, as well as cognitive functions.
Materials and Methods
Ninety-one participants who had been clinically diagnosed with AD dementia (EOAD; n = 55, LOAD; n = 46) and 66 cognitively normal (CN) participants were prospectively recruited. All participants underwent [18F]THK5351 PET scans, [18F]flutemetamol (FLUTE) PET scans, and 3.0-Tesla MRI scans at Gachon University Gil Medical Center, from October 2015 to June 2017. Of the 167 participants, 2 patients diagnosed with AD were excluded due to head motion issues during the [18F]THK5351 PET scan acquisition. An additional four participants (two AD and two CN) were excluded due to errors with the FreeSurfer software. Thus, data from 161 participants (53 EOAD, 44 LOAD, 33 YC, and 31 OC) were used in the analyses.
All patients with AD dementia met the probable AD criteria as proposed by the National Institute of Neurological and Communicative Disorders and Stroke and the AD and Related Disorders Association (McKhann et al., 1984). Patients with familial AD with autosomal dominant inheritance were excluded. Participants with brain lesions on brain MRI such as intracranial hemorrhage, traumatic brain injury, hydrocephalus, territorial infarction, severe white matter hyperintensity (WMH) or WMH associated with radiation, multiple sclerosis or vasculitis were excluded. Severe WMH was defined as both periventricular WMH as cap or band ≥10 mm and deep WMH ≥25 mm using modified Fazekas visual rating scale (Fazekas et al., 1987; Noh et al., 2014). Secondary causes of dementia were ruled out through serum laboratory tests including for thyroid function, metabolic profile, vitamin B12, folate, complete blood counts, and syphilis serology. APOE genotyping was also performed. Clinical interviews with a neurologist or psychiatrist and standardized comprehensive neuropsychological tests were undertaken with all participants.
CN participants were either volunteers from the community or spouses of the patients (age range 44–92, female 45.31%). They had no history of neurologic diseases, psychiatric disorders, abnormalities on neurologic examination, structural lesions on brain MRI such as cerebral infarction, intracranial hemorrhage, traumatic brain injury, hydrocephalus, or severe WMH. Control participants had a clinical dementia rating (CDR) score of 0, and normal cognitive function defined as within 1.5 standard deviations of the age- and education-corrected normative mean on detailed neuropsychological tests. Written informed consent was obtained from each participant, and the Institutional Review Boards of Gachon University Gil Medical Center approved this study.
Image Data Acquisition and Parcellation
3D T1-MPRAGE (Repetition time: 1,900 ms, echo time: 2.93 ms, flip angle: 8°, pixel bandwidth: 170 Hz/pixel, matrix size: 256 × 208, field of view: 256 mm2 × 208 mm2, slice thickness: 1.0 mm, voxel size: 1.0 mm3 × 1.0 mm3 × 1.0 mm3, NEX: 1, total acquisition time: 4 min 9 s, GRAPPA acceleration factor: 2 along phase-encoding direction) was acquired with a 3.0T MRI (Verio, Siemens with a Siemens matrix coil). Images were analyzed using FreeSurfer 5.11 to define regions-of-interest (ROIs) in native-space for each participant and to support a correction of gray matter atrophy and white matter spillover in the PET data. In addition, subcortical region parcellation including the left and right of the amygdala, caudate, hippocampus, pallidum, putamen, and thalamus was also conducted using Desikan-Killiany atlas in MNI 152 space. After subcortical segmentation, each subcortical region on the template surface was sampled as a mesh surface with 2,562 vertices and then transformed into its own original surface. This procedure was conducted for every 12 subcortical regions as detailed in a previous study (Chung et al., 2017; Koo et al., 2017).
The pretreatment process for the PET images proceeded according to the process outlined in a previous study (Kang et al., 2017). [18F]THK5351 was synthesized and radiolabeled at Gachon University Neuroscience Research Institute. All participants underwent a 20-min emission scan beginning 50 min after 185 MBq of [18F]THK5351 was injected intravenously (50–70 min), and a 20-min emission scan 90 min after the intravenous injection of 185 MBq of [18F]FLUTE (90–110 min).
The quantification process for the [18F]THK and [18F]FLUTE PET images was conducted using the same methodology outlined in a previous study (Kang et al., 2017). After co-registration of PET images with the corresponding T1 image, region-based partial volume correction was performed. The regional mean values were then calculated and weight-averaged for pre-defined ROIs (Greve et al., 2014, 2016; Kang et al., 2017). Regional standardized uptake value ratios (SUVRs) were calculated using cerebellar gray matter as the reference region for THK images and the pons for FLUTE images (Okamura et al., 2014; Thurfjell et al., 2014; Lockhart et al., 2016). THK5351 global SUVR was calculated based on the AD-related regions including the frontal, parietal, precuneus, occipital, temporal, anterior cingulate, and posterior cingulate, and the cortical composite FLUTE SUVR was calculated based on AD-related regions including the frontal, parietal, lateral temporal, anterior and posterior cingulate cortices (Thurfjell et al., 2014).
Assessment of Local Shape of Subcortical Structures
To measure the degree of deformation of the subcortical structures, a subcortical mesh was constructed similarly to that described in previous studies (Cho et al., 2011; Chung et al., 2017; Koo et al., 2017). The local shape at every location was then calculated using the surface-based method proposed by Shapira et al. (2008). Weighted sum of the depths at every location was measured in various directions at every location, using the weight as an angle with the normal vector (Chung et al., 2017; Koo et al., 2017). This depth value obtained at 2,562 locations for a single structure reflects local contraction or local shrinkage at that location (i.e., the smaller the value, the more shrinkage caused by atrophy occurred at that position, which expands as the value becomes larger). Therefore, these values were set as the subcortical local shape. This method was undertaken for each of the subcortical regions (amygdala, hippocampus, caudate, pallidum, putamen, and left and right areas of the thalamus).
All participants underwent the Korean version of the mini-mental state examination (MMSE), CDR, and a detailed neuropsychological function test battery (Kang et al., 2003) including attention, praxis, language, verbal and visual memory, visuospatial function, frontal/executive function, and elements of Gerstmann syndrome. Details of the comprehensive test battery are described in our previous study (Lee et al., 2018) and all scores were calculated as age- and level of education-matched z-scores.
A pairwise group comparisons comparing EOAD and LOAD, EOAD and YC, LOAD and OC, YC and OC for subcortical structure shape were performed by permutation tests with 5,000 iterations with ANCOVA. For the comparison between EOAD and LOAD, or between YC and OC, the local shape of the subcortical structures were adjusted for sex, years of education, and intracranial volume (ICV). Additionally, for the comparison between EOAD and YC, or between LOAD and OC, the values for the subcortical local shape were adjusted for age, sex, years of education, and ICV. After the permutation test, the Benjamini-Hochberg false discovery rate was used for region-wise multiple comparison (Benjamini and Hochberg, 1995). A p-value threshold (0.05) was applied to extract the regions showing significant differences between two groups. If the size of the statistically significant region was smaller than the specific size, it was regarded as noise and removed.
We calculated Spearman’s rank correlation coefficient to identify the relationship between the deformity at every location in the subcortical structures and the three types of neuroimaging findings such as THK5351 SUVR, FLUTE SUVR, and cortical thickness (Han et al., 2013). The vertices with a stronger negative correlation of −0.3 or a stronger positive correlation of 0.3 were then selected. The significance level of each vertex was corrected using the Cluster-Based statistics (Han et al., 2013). Through this process we examined the correlation between THK5351 global SUVR, FLUTE global SUVR, mean cortical thickness and the subcortex structures with age, sex, years of education and ICV as covariates separately within EOAD and LOAD groups. In this process, a p-value threshold (0.05) was applied to extract the regions with a significant correlation between the subcortical structure and neuroimaging findings. In cases where the size of the statistically significant region was smaller than the specific size, it was regarded as noise and removed.
To compare the effects of tau, amyloid, and mean cortical thickness on subcortical shape deformity, the relative importance for each factor was calculated using general linear model. For the general linear model, age, sex, education and ICV were used as covariates, and the degree of tau global retention, amyloid global retention, and mean cortical thickness were used as predictors, and local shape was encoded as dependent variable. The r-squared value of the model was calculated from the predictor combinations in all cases in which a specific predictor was included in the general linear model, and the accumulated value was called relative importance and is shown in the Supplementary Figure 1.
Comparisons of the demographic and clinical information were conducted using independent t-test and chi square test followed by Yate’s continuity correction for nominal variables. For volumetric analyses, group comparisons were conducted using independent t-test and analysis of covariance. Correlations between volumes of each subcortical structure and other variables such as THK, FLUTE, cortical thickness, and cognitive function were evaluated using Pearson’s correlation analysis and general linear model. Volume data analysis was conducted with PASW Statistics 23 software (SPSS Inc., Chicago, IL, United States) with significance set at p < 0.05 (two-way).
Demographics and Clinical Characteristics
Table 1 shows the demographic information and clinical characteristics of the participants. There were no significant differences in MMSE scores or CDR-SOB between EOAD and LOAD patients. Global retention of THK5351 PET and mean cortical thickness were not different between EOAD and LOAD, however, the EOAD group showed significantly higher amyloid global accumulation than the LOAD group (p < 0.001).
Comparison of Subcortical Structures
Figure 1 shows the results of group comparisons of the subcortical structural shapes. The EOAD group showed widespread subcortical shape difference in the bilateral amygdala, caudate, hippocampus, putamen, thalamus and the right pallidum relative to YC. The LOAD group showed significant regional shape difference in the bilateral amygdala, hippocampus, and putamen compared to OC. Comparisons of the subcortical structures between the EOAD and LOAD groups showed shrinkage in the bilateral thalamus, left amygdala and left hippocampus in LOAD group.
Figure 1. Group comparison of shape of subcortical structures. (A,B) Group-wise comparison of the subcortical local shape adjusted with age, sex, education, ICV; (C,D) sex, education, ICV was adjusted. Corrected p-values after FDR multiple correction for each 12 structures (left and right amygdala, hippocampus, caudate, pallidum, putamen, and thalamus) were mapped. Abbreviations: EOAD, early-onset Alzheimer’s disease; LOAD, late-onset Alzheimer’s disease; YC, young control; OC, old control.
Volume analysis of the subcortical structures between the groups are presented in Supplementary Table 1. EOAD patients showed reduced volume in the bilateral amygdala, hippocampus, caudate, putamen, and the left thalamus compared to YC. LOAD patients showed reduced volumes in the bilateral amygdala and hippocampus compared to OC.
Correlation Between THK5351, FLUTE, Cortical Thickness, and Subcortical Structures
Figure 2 shows the correlation between THK global retention, FLUTE global retention, mean cortical thickness, and subcortical shape in the EOAD and LOAD groups. In EOAD patients, as THK global retention increased, subcortical structure showed contraction significantly in the bilateral amygdala, caudate, putamen, hippocampus, left pallidum and the left thalamus. In LOAD patients, as THK global retention increased, subcortical shape decreased in the bilateral caudate, putamen, left pallidum and the left amygdala and hippocampus.
Figure 2. Spearman correlation between subcortical local shape and THK5351 retention, FLUTE global retention, and mean cortical thickness in EOAD and LOAD, adjusting with age, sex, years of education, and intracranial volume (ICV) in each EOAD and LOAD group. The mapped p-values are multiple corrected using cluster-based statistics. Abbreviations: EOAD, early onset Alzheimer’s disease; LOAD, late-onset Alzheimer’s disease.
As FLUTE global retention increased, subcortical shape deformity significantly increased in the bilateral hippocampus and the right amygdala and left putamen in EOAD. LOAD patients showed significant subcortical shape deformity in the bilateral putamen, the left amygdala and hippocampus, and the right caudate as FLUTE increased. There was no subcortical structure with positive correlation with THK or FLUTE global retention in EOAD or LOAD patients. Mean cortical thickness was positively associated with subcortical structure shape in the bilateral amygdala, caudate, putamen, and the left hippocampus in EOAD.
Supplementary Figure 1 shows association between THK, FLUTE global retention, cortical thickness and subcortical structures, which was analyzed using the general linear model. In EOAD group, similar to the results of our main correlation analyses, tau and mean cortical thickness were found to have great importance in the bilateral caudate and putamen, whereas amyloid showed little importance. On the other hand, in the LOAD patient group, the importance was higher in correlation with tau in both caudate. And mean cortical thickness was found to have importance than other variables in predicting the local shape in the left putamen.
Supplementary Figure 2 shows the voxel-wise analyses results to evaluate the association between subcortical structural volume and THK or FLUTE global retention. In EOAD, decreased volumes of the amygdala, caudate, putamen had association with THK retention in the nearly whole association cortices. Hippocampal volume was associated with THK retention in the bilateral medial temporal region and the right precuneus and parietal regions in EOAD. The hippocampus was the only region with significant correlation with FLUTE global retention in EOAD. In LOAD patients, no subcortical structural volume showed significant correlation with THK or FLUTE global retention in voxel-wise analyses.
Volumetric correlation between THK global retention, FLUTE global retention, mean cortical thickness, and subcortical structures are presented in Supplementary Table 2. In EOAD patients, THK, FLUTE, and cortical thickness correlated with volumes of the subcortical structures except for pallidum. In LOAD patients, no correlation was found between THK, FLUTE, and mean cortical thickness and subcortical volume.
Correlation Between Cognitive Function and Subcortical Structures
The neuropsychological test results are presented in Table 2.
Patients with EOAD had significantly worse performance in digit span forward and backward, Rey-Osterrieth complex figure test (RCFT) copy memory test other than RCFT recognition, controlled oral word association test supermarket, Stroop test color reading, trail making test type A and B compared to LOAD patients.
Figure 3 shows the association between subcortical shapes and neuropsychological test scores in the EOAD patients. MMSE scores were correlated with the bilateral amygdala, caudate, putamen, and thalamus and the left hippocampus (Figure 3A). CDR-SOB was correlated with the bilateral amygdala, hippocampus, caudate, putamen and thalamus (Figure 3B). Digit span backward was correlated with the bilateral amygdala, caudate, and putamen and the left pallidum (Figure 3C). BNT was correlated with the bilateral amygdala and putamen and the left caudate and hippocampus (Figure 3D). RCFT copy test was correlated with the bilateral caudate and putamen and right amygdala (Figure 3E). SVLT delayed recall showed correlation with the left thalamus (Figure 3F) and RCFT delayed recall showed association with the bilateral putamen (Figure 3G). Executive function measured by Stroop test color reading was correlated with the bilateral amygdala, caudate, and putamen and the left pallidum (Figure 3H).
Figure 3. Spearman correlation between neuropsychological test results and subcortical local shape in EOAD group adjusted with intracranial volume. (A) MMSE, (B) CDR-SOB, (C) digit span backward score in attention, (D) BNT score in language, (E) RCFT copy score in visuospatial, (F) SVLT delayed recall score in memory, (G) RCFT delayed recall score in memory, and (H) color reading test score in frontal/executive function. The mapped p-values are multiple corrected using cluster-based statistics. Abbreviations: EOAD, early onset Alzheimer’s disease; MMSE, mini-mental state examination; CDR-SOB, clinical dementia rating-sum of boxes.
Figure 4 shows the association between subcortical shape and neuropsychological test scores in LOAD patients. MMSE scores were correlated with the bilateral amygdala, hippocampus, and putamen and the left caudate (Figure 4A) and CDR-SOB was correlated with the bilateral amygdala, caudate, hippocampus, putamen, thalamus, and the left pallidum (Figure 4B). Digit span backward was correlated with the left caudate and hippocampus (Figure 4C). BNT was correlated with the left amygdala and hippocampus (Figure 4D). RCFT copy showed no association with subcortical structural shape (Figure 4E). SVLT delayed recall was correlated with the bilateral hippocampus and the left caudate, and the right thalamus (Figure 4F), and RCFT delayed recall was correlated with the right hippocampus (Figure 4G). Stroop test color reading was correlated with the bilateral hippocampus and the left caudate (Figure 4H). Results of the correlation analysis between volumes of the subcortical structures and the detailed neuropsychological test results were presented in Supplementary Tables 3 and 4.
Figure 4. Spearman correlation between neuropsychological test results and subcortical local shape in LOAD group adjusted with intracranial volume. (A) MMSE, (B) CDR-SOB, (C) digit span backward score in attention, (D) BNT score in language, (E) RCFT copy score in visuospatial, (F) SVLT delayed recall score in memory, (G) RCFT delayed recall score in memory, and (H) color reading test score in frontal/executive function. The mapped p-values are multiple corrected using cluster-based statistics. Abbreviations: LOAD, late onset Alzheimer’s disease; MMSE, mini-mental state examination; CDR-SOB, clinical dementia rating-sum of boxes.
In this study, we evaluated the degree of subcortical structural shape deformity in EOAD and LOAD patients, in terms of association with tau and amyloid global retention, cortical thickness, and neuropsychological test results. We observed greater subcortical shape deformity in EOAD and more limited deformation especially in the amygdala and hippocampus in LOAD patients, compared to their age-matched controls. The associations between subcortical structural shape and THK global retention, mean cortical thickness, or cognitive functions were also more prominent in EOAD compared to LOAD patients.
Subcortical shape differed according to the age of onset. Compared to each age-matched control group, EOAD patients exhibited shape deformation in nearly all the subcortical structures, including the bilateral amygdala, hippocampus, caudate, putamen, and thalamus. Meanwhile, LOAD patients manifested decreased shapes in the more restricted region such as the amygdala, hippocampus, and putamen. We also found that the amygdala, hippocampus and thalamus were the major regions that differed in the direct group comparison with more shrinkage in LOAD patients than in EOAD patients. This observation is consistent with previous studies that have shown greater hippocampal atrophy in LOAD compared to EOAD (Frisoni et al., 2007). It was also previously found that EOAD patients mapped with subcortical deformity in wider regions including the amygdala, hippocampus, and putamen compared to LOAD patients with subcortical deformity mapped in the amygdala and hippocampus only, which was also reflected in our findings (Cho et al., 2013). It has been previously reported that EOAD patients show more rapid declines in volume in the bilateral caudate, putamen, and thalamus compared to LOAD (Cho et al., 2013). In another study, researchers found putamen (dorsal striatum) atrophy in EOAD and atrophy in the nucleus accumbens (the ventral striatum) in LOAD patients (Pievani et al., 2013). The limbic area including the hippocampus and amygdala is well known to be affected in the initial stages and the basal ganglia and thalamus are affected in the later stages of LOAD (Braak and Braak, 1991). Our observations may be indicative of a dorsal atrophic pattern of the subcortical structures in EOAD and predominant hippocampal volume loss in LOAD (Pievani et al., 2013).
In terms of tau, amyloid, and cortical thickness, broad regions of the subcortical shape deformation in the caudate and putamen were associated with THK global retention and mean cortical thickness in EOAD patients. There was association with relatively low effect sizes between subcortical shape and FLUTE global retention in both groups. These results were expected because the mean cortical thickness and tau retention reportedly proceeds in a relatively parallel manner (Whitwell et al., 2008, 2018) compared to amyloid retention (Thal et al., 2002). A/T/N classification systems for AD biomarkers note a specific temporal ordering where Aβ (A) precedes tau (T) and neurodegeneration (N), which correlate with clinical symptoms (Jack et al., 2016). Neurofibrillary tangle retention is also known to accumulate to a higher degree in most of the cortical regions sparing the primary sensorimotor cortices in EOAD and in the confined regions including the medial temporal and lateral parietal regions in LOAD (Scholl et al., 2017) compared to age-matched controls. Our results appear to support these past findings and provide evidence that the dorsal subcortical structures such as the caudate and putamen are affected more in EOAD compared to LOAD and it is associated with tau global retention and cortical thinning more significantly than amyloid global accumulation. It was interesting result that amyloid global retentions were correlated with hippocampal shapes in EOAD. APOE e4 genotype might be associated with this correlation. Further study would be needed to investigate this issue.
Subcortical structures correlated with the neuropsychological test results, with greater association found more consistently in EOAD than LOAD patients. Specifically, attention, language, visuospatial function, frontal/executive function, and global cognition was correlated with subcortical shapes in most regions of the caudate and putamen. This corresponds the known function of the caudate and putamen (Koziol and Budding, 2009). The striatum, including the caudate nucleus and putamen, receives inputs from most of the major cortical regions including motor, sensory, perceptual and association, frontal executive, and affective and motivational regions (Koziol and Budding, 2009). The head of the caudate nucleus is known to receive inputs from the dorsolateral prefrontal cortex, while the tail of the caudate receives projections from the parietal cortex, and the putamen receives inputs from the motor cortex (Koziol and Budding, 2009). Regarding the circuitry of the striatum, our results indicate that the subcortical structures may associate with cortical region-specific cognitive impairment in EOAD patients. Contrary to EOAD patients, LOAD patients have shown more limited associations: language impairment was associated with shapes in the left amygdala and hippocampus, verbal memory impairment was associated with shapes in the bilateral hippocampus and the left caudate, and visual memory impairment was associated with shapes in the right hippocampus with lower levels of significance. A previous study found that EOAD patients showed more rapid cognitive decline as determined by the MMSE and frontal and total scores in the comprehensive neuropsychological test battery when compared to LOAD patients (Cho et al., 2013). Memory impairment is known to be a predominant phenotype of LOAD and non-memory impairment, including visuospatial dysfunction and apraxia, is predominant in EOAD (Koedam et al., 2010). The present findings are generally in accordance with these previous reports (Koedam et al., 2010; Cho et al., 2013) showing region-specific associations between subcortical structures and neuropsychological test results in the EOAD and LOAD groups.
We note several limitations to this study. The participants were recruited from a single memory clinic in a tertiary hospital, which may limit generalizability of the results. The lack of information on genetic mutations is also limiting because cognitive function and imaging biomarkers can vary according to genetic differences, especially for comparisons between EOAD and LOAD patients. Although we excluded participants with familial AD with autosomal dominant inheritance, the possibility which participants with genetic mutation could not excluded completely as we did not perform PSEN 1, APP or PSEN2 gene tests to all participants. Another major limitation is the known off-target binding of the tracer used in this study. [18F]THK5351, a first generation tau tracer, is known to bind to monoamine oxidase-B (Ng et al., 2017), and thus can be reactive to both neurofibrillary tangles and reactive astrocytes (Harada et al., 2018). Due to its limited utility as a sole biomarker for tau retention, our results should thus be interpreted with some degree of caution. However, we believe our study remains the first to evaluate subcortical structural shapes in EOAD and LOAD patients, using the three major imaging biomarkers tau PET, amyloid PET, and MRI which can provide valuable evidence of in vivo findings in AD.
We observed pattern differences in subcortical structural shape between EOAD and LOAD. Widespread subcortical atrophy was associated with tau and cortical thinning in most of the cognitive domains, except for memory, in EOAD. In contrast, smaller subcortical regions were associated only with THK retention and restricted cognitive impairment including memory in LOAD. These results may be indicative of a different underlying pathomechanism between EOAD and LOAD. In addition, the imaging results are in close agreement with the known progression of the A/T/N in vivo biomarkers and differences in cognitive function in both groups. The employment of all available imaging biomarkers and age of onset in AD studies may be helpful in more clearly elucidating the association between changes in subcortical regions and other AD biomarkers.
Data Availability Statement
The data set generated and/or analyzed using the current study is available from the corresponding author, YN on reasonable request.
The studies involving human participants were reviewed and approved by Institutional Review Boards of Gachon University Gil Medical Center. The patients/participants provided their written informed consent to participate in this study.
YN and J-KS conceptualized and designed the study. E-CL and JMK drafted the manuscript. JMK, H-ES, S-YL, KHP, DLN, and YN acquired the data. E-CL, JMK, and SS analyzed the data. J-KS took part in methodology. YN was in charge of funding acquisition and resources. YN and J-KS revised the manuscript for intellectual content. All authors reviewed and approved for publication.
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 & Welfare, Republic of 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); by Institute of Information & communications Technology Planning & Evaluation (IITP) grant funded by the Korean Government (MSIT) [No. 2019-0-00079, Department of Artificial Intelligence (Korea University)]; and by Brain Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science and ICT (No. 2020M3C7A101835721).
Conflict of Interest
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.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnagi.2020.563559/full#supplementary-material
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B Methodol. 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
Cho, Y., Seong, J.-K., Shin, S. Y., Jeong, Y., Kim, J. H., Qiu, A., et al. (2011). A multi-resolution scheme for distortion-minimizing mapping between human subcortical structures based on geodesic construction on Riemannian manifolds. Neuroimage 57, 1376–1392. doi: 10.1016/j.neuroimage.2011.05.066
Choo, I. H., Lee, D. Y., Kim, J. W., Seo, E. H., Lee, D. S., Kim, Y. K., et al. (2011). Relationship of amyloid-beta burden with age-at-onset in Alzheimer disease. Am. J. Geriatr. Psychiatry 19, 627–634. doi: 10.1097/jgp.0b013e318202bf3a
Chung, S. J., Shin, J. H., Cho, K. H., Lee, Y., Sohn, Y. H., Seong, J. K., et al. (2017). Subcortical shape analysis of progressive mild cognitive impairment in Parkinson’s disease. Mov. Disord. 32, 1447–1456. doi: 10.1002/mds.27106
Dipasquale, O., and Cercignani, M. (2017). Network functional connectivity and whole-brain functional connectomics to investigate cognitive decline in neurodegenerative conditions. Funct. Neurol. 31, 191–203.
Fazekas, F., Chawluk, J. B., Alavi, A., Hurtig, H. I., and Zimmerman, R. A. (1987). MR signal abnormalities at 1.5 T in Alzheimer’s dementia and normal aging. Am. J. Roentgenol. 149, 351–356. doi: 10.2214/ajr.149.2.351
Frisoni, G. B., Pievani, M., Testa, C., Sabattoli, F., Bresciani, L., Bonetti, M., et al. (2007). The topography of grey matter involvement in early and late onset Alzheimer’s disease. Brain 130, 720–730. doi: 10.1093/brain/awl377
Greve, D. N., Salat, D. H., Bowen, S. L., Izquierdo-Garcia, D., Schultz, A. P., Catana, C., et al. (2016). Different partial volume correction methods lead to different conclusions: an 18F-FDG-PET study of aging. Neuroimage 132, 334–343. doi: 10.1016/j.neuroimage.2016.02.042
Greve, D. N., Svarer, C., Fisher, P. M., Feng, L., Hansen, A. E., Baare, W., et al. (2014). Cortical surface-based analysis reduces bias and variance in kinetic modeling of brain PET data. Neuroimage 92, 225–236. doi: 10.1016/j.neuroimage.2013.12.021
Hahn, C., Lee, C. U., Won, W. Y., Joo, S. H., and Lim, H. K. (2016). Thalamic shape and cognitive performance in amnestic mild cognitive impairment. Psychiatry Investig. 13, 504–510. doi: 10.4306/pi.2016.13.5.504
Han, C. E., Yoo, S. W., Seo, S. W., Na, D. L., and Seong, J.-K. (2013). Cluster-based statistics for brain connectivity in correlation with behavioral measures. PLoS One 8:e72332. doi: 10.1371/journal.pone.0072332
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
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
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
Koedam, E. L., Lauffer, V., Van Der Vlies, A. E., Van Der Flier, W. M., Scheltens, P., and Pijnenburg, Y. A. (2010). Early-versus late-onset Alzheimer’s disease: more than age alone. J. Alzheimers Dis. 19, 1401–1408. doi: 10.3233/jad-2010-1337
Koo, D. L., Shin, J.-H., Lim, J.-S., Seong, J.-K., and Joo, E. Y. (2017). Changes in subcortical shape and cognitive function in patients with chronic insomnia. Sleep Med. 35, 23–26. doi: 10.1016/j.sleep.2017.04.002
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
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
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–944. doi: 10.1212/wnl.34.7.939
Migliaccio, R., Agosta, F., Possin, K. L., Canu, E., Filippi, M., Rabinovici, G. D., et al. (2015). Mapping the progression of atrophy in early-and late-onset Alzheimer’s disease. J. Alzheimers Dis. 46, 351–364. doi: 10.3233/jad-142292
Ng, K. P., Pascoal, T. A., Mathotaarachchi, S., Therriault, J., Kang, M. S., Shin, M., et al. (2017). Monoamine oxidase B inhibitor, selegiline, reduces 18 F-THK5351 uptake in the human brain. Alzheimers Res. Ther. 9:25.
Noh, Y., Lee, Y., Seo, S. W., Jeong, J. H., Choi, S. H., Back, J. H., et al. (2014). A new classification system for ischemia using a combination of deep and periventricular white matter hyperintensities. J. Stroke Cerebrovasc. Dis. 23, 636–642. doi: 10.1016/j.jstrokecerebrovasdis.2013.06.002
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
Pievani, M., Bocchetta, M., Boccardi, M., Cavedo, E., Bonetti, M., Thompson, P. M., et al. (2013). Striatal morphology in early-onset and late-onset Alzheimer’s disease: a preliminary study. Neurobiol. Aging 34, 1728–1739. doi: 10.1016/j.neurobiolaging.2013.01.016
Scholl, M., Ossenkoppele, R., Strandberg, O., Palmqvist, S., Jogi, J., Ohlsson, T., et al. (2017). Distinct 18F-AV-1451 tau PET retention patterns in early- and late-onset Alzheimer’s disease. Brain 140, 2286–2294. doi: 10.1093/brain/awx171
Seltzer, B., and Sherwin, I. (1983). A comparison of clinical features in early-and late-onset primary degenerative dementia: one entity or two? Arch. Neurol. 40, 143–146. doi: 10.1001/archneur.1983.04050030037006
Smits, L. L., Pijnenburg, Y. A., Koedam, E. L., Van Der Vlies, A. E., Reuling, I. E., Koene, T., et al. (2012). Early onset Alzheimer’s disease is associated with a distinct neuropsychological profile. J. Alzheimers Dis. 30, 101–108.
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
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
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
Keywords: Alzheimer’s disease, subcortical shape analysis, tau, amyloid, cortical thickness, positron emission tomography, magnetic resonance imaging
Citation: Lee E-C, Kang JM, Seo S, Seo H-E, Lee S-Y, Park KH, Na DL, Noh Y and Seong J-K (2020) Association of Subcortical Structural Shapes With Tau, Amyloid, and Cortical Atrophy in Early-Onset and Late-Onset Alzheimer’s Disease. Front. Aging Neurosci. 12:563559. doi: 10.3389/fnagi.2020.563559
Received: 19 May 2020; Accepted: 05 October 2020;
Published: 26 October 2020.
Edited by:Kuangyu Shi, Technical University of Munich, Germany
Reviewed by:Daniel Ferreira, Karolinska Institutet (KI), Sweden
Mitsuru Shinohara, Mayo Clinic, United States
Copyright © 2020 Lee, Kang, Seo, Seo, Lee, Park, Na, Noh and Seong. 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.
†These authors have contributed equally to this work