Abstract
Both the apolipoprotein E (APOE) ɛ4 allele and amnestic mild cognitive impairment (aMCI) are considered to be risk factors for Alzheimer’s disease (AD). The primary aim of this study was to determine whether the aMCI-related abnormality in gray matter (GM) cortical thickness and white matter (WM) tracts integrity would be modified by the APOE genotype. A total of 146 older adults, including 64 aMCI patients (28 ɛ4 carriers and 36 non-carriers) and 82 healthy controls (39 ɛ4 carriers and 43 non-carriers), underwent a standardized clinical interview, neuropsychological battery assessment, and multi-modal brain magnetic resonance imaging scans. Compared with control subjects, the patients with aMCI showed significantly reduced cortical thickness bilaterally in the parahippocampal gyrus and disrupted WM integrity in the limbic tracts (e.g., increased mean diffusivity in the right parahippocampal cingulum and bilateral uncinate fasciculus). However, no significant main effects of the APOE genotype and diagnosis-by-genotype interaction on GM thickness and WM integrity were observed. Further, diffusivity measures of the limbic WM tracts were significantly correlated with the parahippocampal atrophy in aMCI. Importantly, the parahippocampal thickness and diffusivity measures of the limbic WM tracts were significantly correlated with the cognitive performance (i.e., episodic memory Z score) in patients with aMCI. These results demonstrate that WM microstructural disruptions in the limbic tracts are present at the early stage of AD in an APOE-independent manner; and this degeneration may occur progressively, in parallel with parahippocampal atrophy, and may specifically contribute to early initial impairment in episodic memory.
Keywords
INTRODUCTION
Alzheimer’s disease (AD) is a neurodegenerative disorder characterized by a progressive impairment of cognitive function and behavior. Amnestic mild cognitive impairment (aMCI) is the most well-characterized cognitive risk factor for AD and is considered to be a transitional state between normal aging and dementia with rates of progression to AD reported to range between 16% and 40% [1, 2]. Thus, aMCI provides an important model to study the mechanisms of AD. Nonetheless, the outcome for any subjects with aMCI is highly variable because many subjects may remain stable or even revert to a normal state, while others progress to dementia, possibly due to the interplay of physiologic, environmental, and genetic factors [3, 4].
The apolipoprotein E (APOE) ɛ4 allele is the most common susceptibility gene for AD [5] and influences the course of disease by increasing the risk for developing AD and lowering the age at disease onset [6, 7]. Case-control studies have shown that the ɛ4 allele is overrepresented in aMCI, and increased frequency of the allele was determined to be the strongest predictor of clinical progression from aMCI to AD [8–10]. The role of APOE-ɛ4 allele in aMCI and conversion to AD has also been the subject of longitudinal studies. However, these studies have been inconsistent. Although some studies in this field have shown that APOE-ɛ4 allele leads to greater longitudinal metabolic decline in aMCI subjects converting to AD [11, 12], other studies have observed that cognitive decline and progression to AD is not strongly related to APOE-ɛ4 allele [13, 14]. Therefore, there remain large discrepancies in the literature concerning the role of APOE-ɛ4 allele in the progression from aMCI to AD.
It is becoming increasingly clearer that to achieve a comprehensive picture of AD evolution, investigations need to be based on both neuropsychological and anatomical information from subjects at high-risk of AD. Previous neuroimaging studies on AD/aMCI typically focused on investigating cortical structural changes in medial temporal grey matter (GM) such as volume reduction in the entorhinal cortex and hippocampus [15–17]. Dickerson et al. [18] also identified several cortical regions in which thickness is a potential biomarker for AD, referred to as “AD-signature” cortical thickness. This measure has predicted symptom severity, amyloid binding in asymptomatic older adults, and conversion to AD among aMCI subjects [19]. Another technique, diffusion tensor imaging (DTI), is capable of analyzing white matter (WM) integrity by using measures such as factional anisotropy (FA) and mean diffusivity (MD). Previous DTI studies on AD/aMCI have demonstrated WM changes, especially in the fiber tracts connecting to the medial temporal cortex, such as cingulum bundle and uncinate fasciculus [20–22]. Furthermore, there is also evidence that GM volume or WM integrity might be affected by the APOE genotype. For example, some [23–26], but not all [27–29], structural imaging studies of cognitively intact ɛ4-carrying adults have shown reduced brain volume (e.g., hippocampal volume) in the ɛ4 carriers. Similar DTI studies in healthy subjects also showed ɛ4-related WM integrity changes (e.g., decreased FA value) in the corpus callosum, fronto-occipital fasciculus, and parahippocampal cingulum [26, 30–32]. Overall, these above studies suggested the association of the presence of aMCI or APOE-ɛ4 status with the alterations of GM and WM integrity. However, the independent contributions of cognitive and genetic risk factors (i.e., aMCI and APOE-ɛ4 allele) for AD have received considerable attention, but less is known about the possible interactive effects.
The primary aim of this study was to determine whether the aMCI-related abnormality in GM cortical thickness and WM tract integrity would be modified by the APOE genotype. Second, several previous studies showed a significant association between hippocampal atrophy and damage to the cingulum bundle (e.g., decreased FA in the cingulum) [33] as well as between corpus callosum FA and atrophy of connected cortical areas in patients with AD [34]. Therefore, in an attempt to further explore the mechanisms associated with WM damage in the early phase of AD (i.e., aMCI), we sought to examine whether WM damage is related to GM atrophy (e.g., cortical thickness reduction). Third, we tried to clarify whether the GM and WM damage is able to explain the impairment of specific cognitive domains (e.g., episodic memory) in patients with aMCI. Finally, we conducted an exploratory classification analysis to explore the feasibility of integrating neuropsychological and neuroimaging measures to accurately distinguish patients with aMCI from healthy controls.
MATERIALS AND METHODS
Participants
Patients with aMCI and healthy controls (HCs) were recruited to establish a registry at the Affiliated Zhongda Hospital, Southeast University [35, 36]. All participants were recruited through a normal community health screening and newspaper advertisements, and they underwent a standardized clinical interview, neuropsychological battery assessment, genetic screening, and multi-modal brain MRI examinations (for details, see the below). The detailed inclusion and exclusion criteria are described in the Supplementary Material. The study was approved by the Research Ethics Committee of ZhongDa Hospital Affiliated to Southeast University, and written informed consent was obtained from all participants.
In the original registry, the frequencies of the ɛ2, ɛ3, and ɛ4 allele in the HCs (8.3%, 80.9%, and 10.8%) and the aMCI patients (6.9%, 70.7%, and 22.4%) were similar to results of previous in Chinese population [37, 38]. In this present study, we sought to investigate the main effects of diagnosis and APOE genotype and diagnosis-by-genotype interactions on GM cortical thickness and WM tracts integrity. Therefore, we selected comparable number of participants in each group to the MRI study. This experimental design has been widely used in the previous APOE ɛ4-related studies [36, 39]. 146 elderly participants (age range 60–80 years), including 64 patients with aMCI and 82 HCs, were Chinese Han and right-handed. Of these aMCI patients, there were 28 ɛ4 carriers (21 ɛ3/ɛ4 and 7 ɛ4/ɛ4) and 36 non-carriers (ɛ3/ɛ3); of these controls, there were 39 ɛ4 carriers (37 ɛ3/ɛ4 and 2 ɛ4/ɛ4) and 43 non-carriers (ɛ3/ɛ3). The participants with one or more ɛ2 allele(s) were excluded from this study due to the allele’s possible protective effect [40, 41]. Table 1 presents the demographic information, APOE status, and cognitive scores of the participants included in this study.
Neuropsychological assessment
For all participants, we assessed their general cognitive function using the Mini-Mental State Examination (MMSE) and Mattis Dementia Rating Scale-2 (MDRS-2), and performed a neuropsychological battery to evaluate their specific functions in episodic memory, visuospatial skills, information processing speed, and executive function, respectively. This battery consisted of the auditory verbal learning test with a 20-min delayed recall, the logical memory test with a 20-min delayed recall, the Rey-Osterrieth complex figure test with a 20-min delayed recall, the clock drawing test, the digital symbol substitution test, trail-making test A and B, the Stroop color-word test A, B, and C, the verbal fluency test, the digital span test, and the semantic similarity test.
In this study, a composite score analysis [36] of these neurocognitive measures was further conducted to increase statistical power by reducing random variability and floor and ceiling effects. For each subject, the raw scores from each test were first transformed to z scores with reference to the means and standard deviations of the test across all subjects. Then, the composite scores were calculated by averaging the z scores within the neuropsychological domains listed below: episodic memory (three tests, including the auditory verbal learning test with a 20-min delayed recall, the logical memory test with a 20-min delayed recall, and the Rey-Osterrieth complex figure test with a 20-min delayed recall), visuospatial function (two tests, including the Rey-Osterrieth complex figure test and the clock drawing test), information processing speed (four tests, including the digital symbol substitution test, the trail making test-A, and Stroop A and B), and executive function (five tests, including the verbal fluency test, the digital span test-backward, the trail making test-B, Stroop C, and the semantic similarity test). Notably, MMSE and MDRS-2 were used for descriptive and diagnostic classifications, but not for the composite measures.
APOE genotyping
Genomic DNA of each subject was extracted from 250 μl EDTA-anticoagulated blood using a DNA direct kit (Tiangen, China). A polymerase chain reaction-based restriction fragment length polymorphism (PCR-RFLP) assay was applied to detect the alleles of rs7412 and rs429358, the haplotype of which ultimately determined the APOE genotype. The specific process is described in the Supplementary Material.
Data acquisition
MR images were acquired in a 3.0 T Siemens Verio scanner (Siemens, Erlangen, Germany) with a 12-channel head coil at the ZhongDa Hospital affiliated to Southeast University. All participants lay supine with the head snugly fixed by a belt and foam pads to minimize head movement. High-resolution T1-weighted axial images covering the whole-brain were acquired using 3D magnetization prepared rapid gradient echo sequence as below: repetition time (TR) = 1900 ms; echo time (TE) = 2.48 ms; flip angle (FA) = 9°; acquired matrix = 256×256; field of view (FOV) = 250 mm×250 mm; thickness = 1.0 mm; gap = 0 mm; number of slices = 176. Diffusion-weighted imaging was acquired with single-shot echo-planar imaging sequence in alignment with the anterior-posterior commissural plane. The diffusion sensitizing gradients were applied along 30 non-collinear directions (b = 1000 s/mm2), together with an acquisition without diffusion weighting (b = 0). Seventy contiguous axial slices were acquired with a slice thickness of 2 mm and no gap. The acquisition parameters were as follows: TR = 10,000 ms; TE = 90 ms; FA = 90°; acquisition matrix = 128×128; FOV = 256 mm×256 mm; number of excitations = 2.
Data preprocessing
Cortical thickness mapping
We used the CIVET pipeline to measure thickness on the cortical surface as previously described [42]. Briefly, the native T1-weighted MR images were first linearly aligned into the stereotaxic space and corrected for nonuniformity artifacts using the N3 algorithm [43]. The registered and corrected images were then automatically segmented into GM, WM, cerebrospinal fluid, and background by using a partial volume classification algorithm [44]. Next, the inner and outer GM surfaces were automatically extracted from each hemisphere using the Constrained Laplacian Automated Segmentation using Proximities (CLASP) algorithm [45, 46]. Cortical thickness was further measured as the distance between corresponding vertices of inner and outer surfaces of GM across 40,962 vertices in each hemisphere [45]. Finally, the thickness data was blurred using a surface-based diffusion smoothing kernel of 30 mm FWHM that preserves corticaltopology [47].
DTI imaging analysis
The DTI data were preprocessed and analyzed using the FSL package (http://www.fmrib.ox.ac.uk/fsl/) [48] and the PANDA software (http://www.nitrc.org/projects/panda/) [49]. Briefly, preprocessing included correction for motion and eddy current effects in DTI images. Subsequently, FMRIB’s Diffusion Toolbox was used to fit the tensor model and to compute the fractional anisotropy (FA), mean diffusivity (MD), axial diffusivity (DA), and radial diffusivity (DR) maps. After that, the TBSS analysis [50] was utilized. First, all subjects’ FA images were skull-stripped and spatially registered to a target image (FMRIB58_FA) in standard space using a nonlinear transformation tool FNIRT. The mean FA image was thresholded at 0.2 and skeletonized to contain only the centers of major WM tracts in standard space. Then, the spatially normalized FA, MD, DA, and DR images for each subject were projected onto the skeleton using the derived registration and projection parameters. Finally, in order to investigate the diffusion changes in several specific tracts, the digital WM atlas JHU ICBM-DTI-81 (see http://cmrm.med.jhmi.edu/) was used to parcel the entire WM into 20 tracks of interest (TOIs) (Supplementary Figure 1). Then, the mean values of each DTI measures (i.e., FA, MD, DA, and DR) for each TOI were extracted for all subjects.
Statistical analysis
Demographic and neuropsychological variables
Statistical analyses of demographics and cognitive scores were performed using two-way analysis of covariance (ANCOVA) for continuous variables and using chi-square tests for categorical variables. Specifically, the main effects of diagnosis (i.e., aMCI versus HC) and APOE genotype (i.e., ɛ4 carriers versus non-carriers), and diagnosis-by-genotype interactions were assessed. For the ANCOVA analyses in the cognitive variables, age, gender and education were considered as covariate variables. These analyses were implemented in SPSS 17.0 (SPSS, Inc., Chicago, IL).
Group differences in cortical thickness and DTI indices
We first assessed the main effects of diagnosis and APOE genotype, and diagnosis-by-genotype interactions on cortical thickness. Specifically, a general linear model with “diagnosis”, “APOE genotype”, and “diagnosis-by-genotype” as predictor variables was applied across the entire cerebral cortex, wherein age, gender, and years of education were included as covariates and cortical clusters with an FWE-corrected p < 0.05 were considered as significant. This analysis was implemented using the SurfStat (http://www.math.mcgill.ca/keith/surfstat/). Then, to investigate the group-differences in DTI measures on WM TOIs, tract-specific ANCOVAs were used to examine the main effects of diagnosis and APOE genotype and diagnosis-by-genotype interactions, after adjusting for age, gender, andeducation.
Relationships of cortical thickness, DTI indices, and cognitive performances
When the bilateral parahippocampal gyrus (PHG) showed significant main effects of diagnosis on the thickness (see results), we further tested the associations between bilateral PHG thickness and the DTI indices of the identified WM TOIs (as determined from the ANCOVA, see results) using the Pearson’s correlation test among the aMCI patients. These analyses were applied to the aMCI patients to understand the relationship between the damage classically attributed in AD (e.g., hippocampal/parahippocampal gyrus atrophy), the deterioration of the WM tracts with the disease, and to show whether these two phenomena can be potential attributed to a linked disease pathway. Finally, multiple linear regression analyses were conducted to calculate the correlations between the neuroimaging measures (i.e., cortical thickness and DTI indices) and cognitive performance in the aMCI and HC groups separately, after controlling for age, gender, education, and APOE genotype. Bonferroni correction for multiple comparisons was used with the significant level considered at p < 0.0125 (p = 0.05/4 composite scores). To statistically compute the differences in magnitudes of the correlations between the neuroimaging measures and cognitive performance in the aMCI and HC groups, correlation coefficients obtained above were further converted into z values by using Fisher’s r-to-z transform. This transformation generated values that were approximately normally distributed. A Z statistic was then used to compare these transformed z values to determine the significance of the between-group differences in correlations, and Cohen’s q was used to quantify the magnitude of difference between correlations; according to |Cohen’s q| <0.20 is considered a small effect, 0.30 a moderate effect, and 0.50 a largeeffect [51].
Exploratory classification analysis
To measure the potential of neuropsychological and neuroimaging measures for future use in diagnosis, the combination of targeted cognitive performance, cortical thickness and DTI indices were tested to see if they could be used as a feature that could separate aMCI patients from normal controls. We took the cognitive measures, mean cortical thickness of the identified regions, and DTI indices of identified WM TOIs for each subject in the aMCI and HC groups as features in the stepwise discrimination analysis. To test the robustness of the results, we also validated the results by using the leave-one-out cross-validation method. This analysis was implemented in SPSS 17.0 (SPSS, Inc., Chicago, IL).
RESULTS
Demographic and neuropsychological variables
Table 1 illustrates the demographics and neuropsychological measures for aMCI and HC subjects stratified by APOE-ɛ4 status. The four subgroups did not differ in age, gender and education (all ps> 0.05). Two-way ANCOVA analyses revealed the main effects of diagnosis and APOE genotype and the diagnosis-by-genotype interactions on cognitive measures. A significant main effect of diagnosis on each cognitive domain was observed, with the aMCI group showing worse cognitive performance than the HC group. There was no significant main effect of APOE genotype on any cognitive measure. In addition, we found a significant diagnosis-by-genotype interaction only on the visuospatial function, with ɛ4 carriers showing worse performance than non-carriers in the aMCI group but no genotype difference in the HC group.
Group-based differences in cortical thickness and DTI indices
Cortical thickness
The main effects of diagnosis on GM cortical thickness are illustrated in Fig. 1. There were two significant cortical clusters (Fig. 1A), which were located around the bilateral parahippocampal cortex (FWE-corrected p < 0.05). The bar chart revealed a thinner thickness of the aMCI patients compared with the control subjects (Fig. 1B). No regions showed significant main effect of APOE genotype or diagnosis-by-genotype interaction.
DTI measures
The main effects of diagnosis on DTI indices are showed in Fig. 2. TOI analyses revealed increased MD in the WM tracts connecting to the medial temporal lobe in the aMCI patients comparted with the control subjects, including the right parahippocampal cingulum (PHC) and uncinate fasciculus (UF) tracts (p < 0.05). Meanwhile, increased DR and DA were found only in the right UF in the aMCI patients compared with the controls (p < 0.01). Further, DA of the right PHC and left UF were increased in the aMCI patients relative to the control subjects (p < 0.05). However, there were no significant tract-specific FA differences between the aMCI patients and control subjects. In addition, no significant main effect of APOE genotype or diagnosis-by-genotype interaction was observed on the mean diffusion metrics of each WM TOI (all ps> 0.05).
Relationship between WM abnormalities and GM atrophy
We further examined the correlation between DTI-measured WM changes and cortical thickness analysis of GM changes in the patients with aMCI. Significant correlations between changes in DTI measures and GM atrophy are shown in Table 2. In the identified WM TOIs, such as right PHC and bilateral UF, we observed that the mean diffusion metrics (i.e., MD, DA, and DR indices) correlated negatively with the bilateral PHG thickness in the aMCI group (all ps < 0.005).
Relationship between cognitive performance and neuroimaging measures
We finally examined the relationships between the cognitive composite Z scores (i.e., episodic memory, visuospatial function, information processing speed, and executive function) and neuroimaging measures (i.e., cortical thickness and DTI indices) in the aMCI and HC groups. In the aMCI group, the bilateral PHG thickness positively correlated with the episodic memory Z scores (Table 3, Supplementary Figure 2, all ps< 0.001). In addition, MD, DA, and DR indices were negatively correlated with the episodic memory and executive function Z scores in the right PHC and UF tracts for the patients with aMCI (Table 3, Supplementary Figure 3, all ps< 0.005). However, in the control subjects, the neuroimaging measures (i.e., cortical thickness and DTI indices) were not correlated significantly with any cognitive performance. Importantly, statistical comparison of the correlations using Z statistic further indicated that the associations were significantly different between the aMCI and HC groups (Table 3, all ps< 0.05, |Cohen’s q| >0.30).
Exploratory classification analysis
We selected cognitive composite scores and neuroimaging measures as independent variables for stepwise discrimination analysis in our model selection. Considering the relatively small number of subjects included in this study, we used the cognitive and neuroimaging measures with significant aMCI versus HC group differences as candidate variables. Therefore, four cognitive measures (episodic memory, visuospatial function, information processing speed, and executive function Z scores) and eight neuroimaging measures (bilateral PHG thickness, right PHC, MD, and DA indices, right UF MD, DA, and DR indices, and left UF DA index) entered in the model. Using this composite model, we could correctly distinguish the aMCI patients from the control subjects in 86.3% of the cases (sensitivity, 82.8%; specificity, 89.0%) (Table 4). The leave-one-out cross-validation results also showed that 85.6% of subjects could be correctly classified between the two groups (Table 4). Taking each subject’s discriminative score as a threshold, the performance in terms of the receiver operating characteristic (ROC) curve was shown in Fig. 3. The area under the ROC curve (AUC) of the proposed method was 0.95, indicating an excellent diagnostic power. Further analysis revealed that the most discriminative feature for classification included the episodic memory and information processing speed Z score, bilateral PHG thickness, and right UF MD index.
DISCUSSION
Neuropathological investigations suggest that the degenerative progress occurring in AD is initially localized in the medial temporal lobe (e.g., entorhinal and hippocampal cortex), and then extend to other neocortical regions [52]. This succession of pathophysiological events is thought to underlie the corresponding progressive accumulation of cognitive deficits [53]. In this perspective, we investigated the abnormalities of GM and WM in patients with aMCI, which represent the earliest stage of AD pathology. This study yielded four principal findings, which we will describe in turn. First, we found that patients with aMCI displayed significantly reduced GM thickness bilaterally in the PHG and abnormal WM integrity in the limbic tracts (i.e., right PHC and bilateral UF) in an APOE-independent manner. Second, significant and anatomically congruent correlations between parahippocampal atrophy and diffusivity measures of limbic WM tracts were observed in aMCI patients, suggesting the WM disruption of limbic tract is secondary to cortical neuronal damage. Third, the PHG thickness and diffusivity measures of limbic WM tracts (i.e., right PHC and UF) were significantly related to the cognitive performance (i.e., episodic memory Z score) in patients with aMCI. Finally, the combination of neuropsychological and neuroimaging measures could accurately distinguish patients with aMCI from healthy controls in 86.3% of the cases; the most discriminative features for classification involved the episodic memory and information processing speed Z score, bilateral PHG thickness, and right UF MD index.
aMCI-related decline in GM thickness and WM integrity
Cortical thickness of the medial temporal lobe was most severely reduced in patients with aMCI, and within the medial temporal lobe, the parahippocampal gyrus was most affected. This is consistent with previous MRI and histopathological studies showing that this region of the brain is affected early and profoundly in the course of AD [54–56]. In addition, the present study further demonstrated the microstructural changes of the limbic tracts (including PHC and UF) as manifested in higher MD, DA, and DR among the patients with aMCI. According to previous animal model studies, these changes are likely related to demyelination and axonal loss [57–59]. Interestingly, the limbic WM degeneration has been previously documented in AD/aMCI, indicating that decreased WM integrity in the limbic tracts is a highly replicable finding across DTI studies of AD/aMCI [22, 60–62]. More importantly, further correlation analyses revealed that PHG thickness and DTI indices in the WM limbic tracts were significantly correlated with the composite episodic memory Z scores across patients with aMCI. Therefore, these findings suggest that alterations within the medial temporal cortex and limbic tracts may serve as a potential imaging biomarker of early AD-related brain changes.
In the present study, the ANCOVA analyses revealed no significant main effect of APOE genotype and diagnosis-by-genotype interaction on GM thickness or WM integrity. Many previous neuroimaging studies have suggested an age-dependent effect of APOE ɛ4 on GM and WM integrity. For example, using region-of-interest or voxel-based morphometric analysis, some previous studies have shown an association between cerebral morphometry and APOE genotype in healthy young and middle-aged individuals (<65 years) [63–66] but not in healthy elderly subjects (≥65 years) [28, 29]. Likewise, using whole-brain DTI parameters, Heise et al. [67] compared the WM integrity for younger (mean age 28.6 years±4.20) and older (mean age 64.9±7.19) APOE ɛ4 carriers and non-carriers. These authors observed widespread decreases in FA values for younger but not older APOE ɛ4 carriers relative to non-carriers. Consistent with these GM/WM studies, we also revealed no association between APOE genotype and GM/WM integrity in the healthy elderly subjects (mean age 70.1 years±6.03). One plausible explanation for the lack of APOE effect in the later life is the different demographics of the nondemented subjects. It is widely accepted that inter-subjects anatomic variability increases with advancing age, especially in older subjects [29]. Therefore, it is possible that an ɛ4 effect on GM/WM integrity was present in our sample but was small enough to have been overwhelmed by the much greater effect of inter-subjects’ anatomic variability.
Relationship between WM abnormality and GM atrophy
The current study allowed us to directly correlate changes between GM and WM in patients with aMCI. DTI measures in the right PHC and bilateral UF were found to correlate with the PHG thickness across aMCI patients. In agreement, correlations between hippocampal volume and WM limbic tract (e.g., cingulum and UF) FA have been previously reported in AD patients [20, 68] or a combination of aMCI and AD patients [69]. Further studies in aMCI showed that the atrophy rate of the cingulum and UF was found to correlate with baseline hippocampal GM atrophy [70], suggesting that both WM and GM degenerative processes are not independent from each other over the course of AD development. Our correlation results in patients with aMCI fully support these findings. It is demonstrated that the UF connects the anterior part of medial temporal lobe (i.e., entorhinal and perirhinal cortices) with the orbital and polar frontal cortex [71, 72]; importantly, the entorhinal cortex is the first site of neuronal death in AD. Therefore damage in the UF could be tightly related with medial temporal GM atrophy and may result in memory impairment [70, 73]. Regarding the PHC, it connects the hippocampus with the posterior cingulate cortex (PCC) [74]. Hypometabolism in the PCC has been consistently reported in patients with aMCI [75]. Further studies indicated that PCC activity during memory tasks is correlated with hippocampal volume [76, 77] and that de-afferentation of the PCC may contribute to its dysfunction from the early-AD stage (e.g., aMCI) [70, 78]. Taken together, our present finding is consistent with the classic secondary degeneration model of AD [33, 68], and further suggests that the observed WM tract changes in patients with aMCI may be a consequence of neuronal loss in the medial temporal cortex (e.g., hippocampal and parahippocampal gyrus). However, a direct effect of amyloid-β on WM tissue cannot be excluded [79].
Exploratory classification analysis
Finally, in this present study, we selected cognitive composite scores and neuroimaging measures as independent variables for stepwise discrimination analysis in our model selection. Interestingly, this composite model could accurately distinguish patients with aMCI from healthy controls in 86.3% of the cases (sensitivity, 82.8%; specificity, 89.0%); and the most discriminative features for classification involved the episodic memory and information processing speed Z score, bilateral PHG thickness, and right UF MD index. Thus, this finding indicates that cortical thinning in the PHG, damage in the UF, a combination of episodic memory and information processing speed deficits may provide complementary information for the diagnostic classification. More generally, these finding reinforce the view that the detection of prodromal AD may benefit from the combination of targeted neuropsychological and neuroimaging markers [80, 81].
Limitations
There were several limitations in this study. First, it lacked amyloid imaging and histopathological data suggestive of AD pathology in patients with aMCI, thereby relying on clinical diagnosis only. Therefore, a considerable amount of clinical and biological heterogeneity existed in the present sample of MCI subjects. Notably, the present study only recruited aMCI subjects, who are considered to have a high risk to convert AD dementia. Importantly, as shown in the Results section, the aMCI group in our study showed significantly reduced cortical thickness in the bilateral parahippocampal gyrus compared with healthy controls. As the new guidelines for the AD-related diagnosis [82, 83] suggested the hippocampal or parahippocampal gyrus atrophy as a neurodegeneration biomarker, the aMCI subjects recruited in this study may likely for a relatively homogeneous group. Second, further validation on a larger sample is required to understand the present results, and longitudinal studies should be conducted to clarify the progression of brain changes over time. Finally, regarding the exploratory classification analysis, one might argue that the use of neuropsychological markers in this model is a tautology because they were used to diagnose the aMCI. However, this may not be correct because the diagnosis of aMCI already requires episodic memory impairment at a level consistent with AD itself [80]. Notably, use of cognitive markers has certain advantages: clear and significant effect on odds ratios, objectivity in scoring, comparative economy in terms of expense and time, and reliability[80, 84].
CONCLUSIONS
In summary, this present study reveals that WM integrity in the limbic tract (i.e., right PHC and bilateral UF) is disrupted from the early stage of AD (i.e., aMCI) in an APOE-independent manner. This degeneration seems to be closely related with, and possibly a consequence of, medial temporal GM atrophy. Additionally, this study further supports the view that the diagnosis of aMCI might benefit from the combined use of cognitive testing and brain imaging.
Footnotes
ACKNOWLEDGMENTS
We thank all the patients and volunteers for participating in this study. We would like to thank Dr. Alan Evans for kindly providing the CIVET software. This study was supported by the Projects of International Cooperation and Exchanges NSFC (grant number 81420108012), the National Key Basic Research Program of China (grant number 2014CB846102), the Natural Science Foundation of China (grant number 81601559), the Guangdong Provincial Natural Science Foundation of China (grant number 2016A030310233), and the Key Program for Clinical Medicine and Science and Technology: Jiangsu Province Clinical Medical Research Center (grant numbers BL2013025, BL2014077).
