Abstract
Keywords
INTRODUCTION
Damage to the cerebral white matter (WM), known by several terminologies including leukoaraiosis, is a pathology commonly found in older adults. Magnetic resonance imaging (MRI) is sensitive to the detection of WM damage, which we refer to here as white matter signal abnormalities (WMSA) [1, 2].Although histopathological studies have indicated a mix of heterogeneous findings that correlate with WMSA including demyelination and gliosis [3, 4], they are often denoted ‘white matter hyperintensities of presumed vascular origin’ [5, 6], and a loss of surrounding vascular integrity is typically indicated as the upstream pathogenic mechanism responsible for their manifestation in both cognitively normal aging as well as in Alzheimer’s disease (AD) [4, 8]. The clinical significance of WMSA has been unclear and in some cases this damage appears to be a benign comorbidity of typical aging with a prevalence as high as 80–95% in older populations [9, 10]. However, it is also known that WMSA are associated with some degree of cognitive decline in otherwise healthy older adult [11–13], exhibit alterations in imaging properties around the time of a dementia diagnosis [2], and are increased in volume in AD [14, 15], suggesting a more direct link between this damage and the clinical expression of dementia [9, 16–19].
Although prior work has described the ‘periventricular’ and ‘deep’ spatial distributions of white matter lesions [20, 21], no juxtacortical localization has been performed to date, and it is unknown whether such cortically proximal regional distributions differ in typical aging compared to AD. This information is critical in that it may provide new mechanisms for identifying individuals at risk for the development of AD, as well as yield a clearer picture of the underlying pathologies of the disease, which is becoming increasingly recognized as a heterogeneous process with potential vascular contributions [22, 23]. In addition, associating WMSA with proximal cortex may be important for understanding the functional consequences of damages white matter as well as for designing and testing potential therapeutic interventions. The current understanding is that WMSA first appear in periventricular WM where tissue is most susceptible to hypoperfusion and likely progresses outwards toward the juxtacortical WM over time in ‘normal’ aging [9, 24], although regionally isolated WMSA exist as well. Several studies have demonstrated WMSA localized to parietal [25, 26] and periventricular regions [18, 27–30] as being indicators of AD, yet other studies have found limited or no spatial differences between AD and controls in terms of distribution of WMSA [15, 31].
A potential reason for the discrepancies between studies may be the methodologies used to assess WMSA. In general, WMSA are rated using semi-quantitative measures such as the Fazekas scale [32] or they are automatically quantified with computational segmentation techniques. In either case, their spatial patterns are typically described as belonging to ‘deep’ versus ‘periventricular’ WM, sometimes further subdivided into the four main lobes of the brain [15, 33–35]. This has provided some information about the spatial distribution of tissue damage, yet the large volume of such descriptions likely results in regional heterogeneity within each lobar measurement, and additionally may miss critical areas in which WMSA may reflect specific focal functional consequence. In addition to the lack of understanding regarding the precise spatial accumulation of WMSA, it is unclear if certain regions accumulate more WMSA than other regions as the total amount of WMSA in the brain increases, and whether regional differences may differentiate cognitively healthy older adults form individuals with AD.
In the present study, we demonstrate a novel mapping procedure for the spatial assessment of WMSA based on proximity to 35 anatomically-based cortical regions [36, 37]. We use this technique to determine whether the global increase in WMSA reported for individuals with AD is due to preferential regional accumulation as well as whether spatial patterns in WMSA volume differ in AD compared to cognitively healthy older adults when accounting for global WMSA volume (i.e., does accumulation in certain regions scale disproportionately with global WMSA volume in AD compared to what is seen is typical aging). The overall findings in this study demonstrate spatial differences in WMSA accumulation in AD that are related to the conversion from mild cognitive impairment (MCI) to AD providing information linking specific regional accumulation of WMSA to the development of dementia.
METHODS
Data
Data used in the preparation of this article were obtained from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database (http://adni.loni.usc.edu). The ADNI was launched in 2003 as a public-private partnership, led by Principal Investigator Michael W. Weiner, MD. The primary goal of ADNI has been to test whether serial MRI, positron emission tomography, other biological markers, and clinical and neuropsychological assessment can be combined to measure the progression of MCI and early AD. For up-to-date information, see http://www.adni-info.org.
The present study used data from the ADNI-1 database that were processed locally at the time of analysis and were selected based on availability of T1-weighted, T2-weighted, and proton density (PD)-weighted images for WMSA processing. Three hundred and eight individuals were used in the following analyses; 127 of these had a diagnosis of AD, 107 were cognitively healthy age-matched older controls (OC), and 74 had a diagnosis of amnestic MCI (Table 1). ADNI-1 follows individuals every 6 months and so MCI individuals were stratified by ADNI-1 into three groups based on their distance in time from AD conversion: those who were 6 months away from converting (MCI-6; n = 16), those 12 months away from converting (MCI-12; n = 35), and those 24 months away from converting (MCI-24; n = 23). Data were selected based on the having sufficient quality to be processed by FreeSurfer and limited to datasets with available T1-weighted, T2-weighted, and proton density (PD)-weighted images for WMSA processing. Demographic data such as age, sex, years of education, history of hypertension, history of endocrine-metabolic disorder, and Mini-Mental State Examination (MMSE) scores were additionally acquired from the ADNI database (Table 1).
MRI acquisition
All data were acquired on a 1.5-T scanner at rigorously validated sites, which all followed a previously described standardized protocol [38]. The protocol included a high-resolution, T1- weighted sagittal volumetric magnetization prepared rapid gradient echo sequence and axial PD and/or T2-weighted fast spin echo sequence. The ADNI MRI core optimized the acquisition parameters of these sequences for each make and model of scanner included in the study. All scanner sites were required to pass a strict scanner validation test before being allowed to scan ADNI participants. Additionally, each scan of ADNI participants included a scan of the phantom, which was required to pass additional strict validation tests. For the purposes of this study, it is important to note that although T2/PD data were acquired at a lower resolution than the T1 data, the T2/PD sequences having a 5 mm slice thickness, identification of juxtacortical white matter was done solely with the 1 mm3 T1-weighted data through FreeSurfer’s main recon-all processing stream. T2/PD data were then registered to the T1 data and partial volume correction was performed [39].
MRI preprocessing
Each individual’s T1-weighted MRI was processed using FreeSurfer’s main recon-all processing stream with a recently described extension for the segmentation of WMSAs [2]. Cortical reconstruction and volumetric segmentation was performed using FreeSurfer (http://surfer.nmr.mgh.harvard.edu/, version 5.1). The technical details of these procedures are described in prior publications [36, 40–43]. Importantly, the FreeSurfer atlas used in the recon-all processing stream is based on manual labelings of 40 subjects distributed in age and pathology, with 10 each of young, middle-aged, and elderly subjects and an additional 10 with an in vivo diagnosis of probable AD. The explicit inclusion of pathologicalexamples in the atlas makes it significantly more robust to neurodegeneration in morphometric studies. All processed data were manually inspected for accurate volumetric segmentation and quality of cortical boundary delineation.
WMSA segmentation
Images were processed with an automatic WMSA segmentation stream that has been previously described [2]. Briefly, this procedure registers T2- and PD-weighted images to an individual’s T1-weighted image after it has been processed through FreeSurfer’s recon-all stream. It then performs intensity normalization of all three modalities using a multimodal atlas and segments WMSA from normal-appearing white matter (NAWM) using a spatial array of multimodal Gaussian classifiers as well as individual-based heuristics. The principles implemented are similar to those used in the FreeSurfer automatic segmentation stream for normal brain structures, but it is important to note that this method is distinct from the ‘WMHypointensity’ segmentation method performed in the FreeSurfer recon-all processing stream which uses T1-weighted images alone. The use of three different imaging modalities with this newer method (more thoroughly described in Lindemer et al. [2]) provides more information on WM integrity and therefore allows for more accurate WMSA labeling. Given the use of multiple image types (T1/T2/PD), the signal may be ‘hyper’ or ‘hypo’ intense depending on the image modality, motivating our use of the term ‘WMSA’. Specifically, on T1-weighted images gray matter (GM) and cerebrospinal fluid (CSF) appear darker than WM, but WMSA appear darker than the surrounding NAWM. On T2-weighted imaging, CSF is the brightest followed by GM and then WM, but WMSA appear brighter than its surrounding NAWM. On PD-weighted imaging, GM is brighter than WM but WMSA appear brighter than surroundingNAWM.
The WM was divided into 70 regions of interest (ROIs) that encompassed the entire subcortical WM as well as the periventricular WM as described previously [44]. Briefly, the cortical parcellations defined by the FreeSurfer recon-all processing stream are used to subsequently assign a label to the underlying WM by the construction of a Voronoi diagram in the WM voxels based on distance to the nearest cortical parcellation label. Using a distance constraint of 5 mm, each Voronoi polygon inherited the label of the cortical parcellation unit, yielding a complete labeling of the cerebral WM. As WMSA has been shown to be a generally symmetrical pathology [45] and we expected and then verified regional WMSA symmetry in this work, we therefore combined across hemispheres for a total of 35 final ROIs. Regional WMSA differences between diagnostic groups were first assessed using the raw WMSA volume as a metric. Due to the known decrease in WM volume with increasing age [44], the WMSA volume as a percent of the individual’s native-space ROI total WM was also calculated for each ROI. In this way, a clearer picture of ‘functional’ impairment may be disentangled, as it indicates how much healthy tissue remains in the region.
Statistical analyses
A general linear model (GLM) was implemented to test the difference in WMSA burden in each ROI between AD (n = 127) and OC (n = 107) groups while controlling for age. This was performed using both raw WMSA volume and WMSA as a percent of each ROI. To insure that regional increases were not simply due to greater global WMSA volume, these GLMs were also conducted using regional WMSA burden normalized to the individual’s total WMSA burden. In total, each ROI was tested for group differences four times using two different WMSA metrics (raw WMSA volume and WMSA% of ROI) with and without normalization for total WMSA. Results are reported with and without correction for multiple comparisons using the false discovery rate (FDR) method [46], correcting for 140 tests. Since raw WMSA volume and WMSA as a percent of each ROI showed similar results when comparing AD to OC, the WMSA percent of each ROI was calculated for each MCI (n = 74) individual. A linear regression was performed on WMSA burdens in each ROI across the three MCI subgroups, and a linear hypothesis test was conducted to determine if there was a significant increase in WMSA burden in MCI individuals who were closer in time to ADconversion.
Correction for global WMSA reduced the number of regions exhibiting group differences in WMSA volume between AD and OC. We therefore conducted secondary analyses to determine whether regional group by global WMSA burden interactions could be detected which would indicate differential scaling of regional WMSA with global WMSA between the groups suggesting differential regional damage propensity.
RESULTS
Demographics
Demographic information for the individuals in each of the diagnostic groups is listed in Table 1. As determined by a one-way ANOVA with post-hoc pairwise comparisons, groups did not differ significantly in age or years of education, but the AD group had significantly higher total WMSA (p < 0.001). As expected due to the use of MMSE score in diagnoses, AD individuals had lower scores over OC and all MCI subgroups, and OC were significantly higher than all MCI subgroups. Chi-squared goodness-of-fit tests demonstrated that there was no significant difference in group distributions of sex, history of hypertension, hypercholesteremia, or endocrine-metabolic disorders. In the OC group, age was significantly correlated with total WMSA burden (R = 0.3299, p < 0.001) but not in AD (R = 0.1369, p = 0.1248), nor in the combined group of MCI individuals (R = 0.1512, p = 0.1985). There was no significant group by global WMSA interaction with age. Due to these findings, all further group comparisons with WMSA were controlled for age.
Spatial differences in raw WMSA volume
Out of the 35 ROIs, 26 demonstrated significantly greater raw WMSA volume in AD compared to OC (p < 0.05 uncorrected; Fig. 1, Top Left). All of these survived FDR correction except for the banks of the superior temporal sulcus, lateral orbitofrontal, precuneus, and superior temporal WM (p = 0.02). After normalizing the WMSA volume in each ROI by each individual’s total WMSA burden, nine regions remained significant at p < 0.05 (Fig. 1, bottom left). Of these, fusiform, inferior parietal, inferior temporal, pars orbitalis, rostral middle frontal, and supramarginal WM survived FDR correction. Insular and periventricular WMSA showed greater normalized WMSA volumes in OC compared to AD, but only insular WM survived FDR correction.
There was a significant group by global WMSA interaction for WMSA within several regions including cuneus, fusiform, inferior parietal, inferior temporal, lateral occipital, middle temporal, periventricular and frontal pole WM (p < 0.05). This indicates that increasing global WMSA is linked to a greater accumulation of local WMSA in these regions compared to the accumulation in OC demonstrating the differential lesion properties in these groups.
Follow-up analyses assessed whether spatial differences were seen when examining only individuals with hypertension as well as only individuals without hypertension in each diagnostic group. In AD individuals with hypertension, there was a significantly higher WMSA burden in caudal anterior cingulate WM than in those without hypertension. In OC individuals with hypertension, there was a significantly higher WMSA burden in the banks of the superior temporal sulcus, lateral orbitofrontal, middle temporal, precuneus, superior temporal, supramarginal, and periventricular WM over those without hypertension.
Spatial Differences in WMSA as a percent of total ROI WM
Measures used in the prior results were based on raw WMSA within a region. However, WMSA volumes corrected for regional total tissue volume is likely a more relevant marker of the degree of regional damage and compromise. We therefore repeated the ROI comparisons using a corrected measure of regional WMSA as a % of total WM volume in the region. Results remained essentially unchanged except in precuneus and superior temporal regions, which no longer reached significance (Fig. 1, top right). Of the 24 significant ROIs, only lateral orbitofrontal, posterior cingulate, and superior parietal did not survive FDR correction. We next performed a secondary normalization of the WMSA percent measure to each individual’s total WMSA (normalization for global WMSA) and found 8 ROIs to have significantly higher WMSA in AD, and insular WM to have higher WMSA in OC (Fig. 1, bottom right). Of these, only inferior parietal and frontal pole regions did not survive FDR correction.
When comparing individuals with and without hypertension in each diagnostic group, only precuneus WM showed a greater WMSA/total WMSA ratio in the OC group with hypertension over the OC group without hypertension.
Spatial WMSA and MCI time to conversion
Of the 35 ROIs tested, all but four exhibited the highest mean WMSA in the MCI-6 group (the MCI closest to conversion) relative to the other groups. A linear regression demonstrated that there was a significant relationship between regional WMSA burden and time-to-AD-conversion in fusiform, inferior temporal, pericalcarine, rostral middle frontal, and superior frontal WM (p < 0.05 uncorrected) (Fig. 2). After normalizing for total WMSA, only inferior temporal and rostral middle frontal regions exhibited a significant relationship with time-to-AD-conversion.
DISCUSSION
We demonstrate here for the first time a set of juxtacortical WM regions where a greater WMSA volume is found in individuals with AD, where the scaling of WMSA with global WMSA is different in patients with AD compared to OC, and where larger volumes of WMSA are indicative of individuals with MCI who are temporally close to conversion. These results provide unique insight into the pathological basis of AD dementia and demonstrate the utility of the regional WMSA mapping procedure based on proximity to specific cortical regions. Future work will investigate if this regionally-specific WMSA manifestation is linked to distinct functional and clinical consequences. In addition to the regional analyses, we explored the relationship between regional and global WMSA volume and also measure WMSA as a percent of the ROI’s total WM to account for WM tissue atrophy that is known to occur in aging and dementia [44]. While prior studies have demonstrated regional differences in periventricular [27, 28] and parietal [25, 26] WMSA between the two diagnostic groups, we provide evidence for more widespread increases in WMSA in AD that include temporal and frontal WM, and that some of these regions also show increases in WMSA as individuals with MCI are closer to AD conversion. It has been previously established that individuals with AD have greater global WMSA burdens [2, 48], and to ensure that this factor is not solely responsible for any regional increases, our technique in which we normalize for global WMSA describes what we believe to be true regional increases in cuneus, inferior temporal, middle frontal, and inferior parietal WM. Based on our findings, we hypothesize that the spatial distribution of WMSA may be important to the understanding of the pathological processes that occur separately from normal aging during the development of AD. This is supported by the relationship between increasing regional WMSA burden and decreasing time-to-AD-conversion in MCI individuals. While the histopathological profile of WMSA presents as complex and heterogeneous and includes gliosis and demyelination, WMSA are generally regarded as a vascular pathology [5, 6]. Furthermore, studies that compare the histopathological profiles of WMSA between AD and OC individuals have demonstrated that they are similarly comprised of lacunes, arteriosclerosis, cerebral amyloid angiopathy, and microinfarcts [49]. Individuals with AD, however, demonstrate more marked deep WM demyelination as well as nearly doubled thickness of the adventitia of deep WM arteries over OC individuals [4] and WM atherosclerosis as well as WMSA burden are associated with higher Braak [50] scores in aging populations [8]. Thus, our findings along with previous findings provide evidence that vascular pathologies may be contributing factors to AD development, acting either independently from more classically known AD pathologies or due to a common upstream process. While some studies have found WMSA to be associated with classic AD pathologies such as cerebral amyloid [51], others have found WMSA and cerebral amyloid to be independent predictors of AD [14] and that vascular risk factors are associated with WMSA independent of CSF amyloid levels [52]. The findings of increased regional WMSA volume independent of global volume suggests that there is something unique about AD that results in this regional WMSA accumulation, and that WMSA do not simply represent a comorbid vascular condition that is found to be heightened in individuals with AD.
As we demonstrate in our analyses in which we normalize for total WMSA, regional differences between groups are strongly influenced by total WMSA burden. Prior to normalization for global WMSA burden, there were widespread regions where AD had greater WMSA volume compared to OC including medial frontal, postcentral, and precentral WM. This finding perhaps demonstrates a vascular comorbidity of AD that is expressed over cognitively healthy older controls but which is not a unique component of the AD disease process itself. After normalization for global WMSA however, 9 out of 35 ROIs remain significantly greater in WMSA volume in AD including inferior temporal, middle frontal, inferior parietal, and cuneus WM. We believe that this is a critical concept in understanding the contribution of WMSA to AD in that it isolates regional differences independent of global WMSA burden. To our knowledge, with the exception of parietal WMSA [25, 26], these regions have not been reported to demonstrate higher WMSA in AD and we attribute these findings to our methodology which normalizes for global WMSA burden. WMSA in temporal WM fit into the classic picture of AD that demonstrates marked cortical changes in relation to memory loss symptoms, but whether these are two separate processes or part of the same pathological pathway is currently unknown. The increase in frontal WMSA in AD may explain some of the executive dysfunction seen in AD [53, 54] which is highly dependent on frontal lobe function. Furthermore, others have shown that increased WMSA anywhere in the brain is associated with greater executive dysfunction [55]. Interestingly, we find that normalizing for total WMSA brings out a significant increase in periventricular and insular WMSA in cognitively healthy older controls over individuals with AD, which to our knowledge has not been previously demonstrated. In fact, several studies have reported greater WMSA in periventricular regions in individuals with AD [18, 28], but these studies did not control for global WMSA and therefore may have been detecting global WMSA increases that were most obvious in periventricular WM as this is where the majority of WMSA occur in the course of aging. We interpret this finding to mean that WMSA preferentially accumulate in deep regions where perfusion is low during the course of normal aging such as periventricular and insular WM, but that the WMSA that accumulate elsewhere in AD may be due to a selective type of vascular dysfunction that is unique to AD and not found in general vascularconditions.
Our results indicate that as MCI individuals are closer in time to AD conversion, WMSA burden increases in regions that demonstrate higher WMSA in AD over OC. In our linear regression analyses, we demonstrate that the regions showing a statistically significant relationship between WMSA and time-to-AD-conversion are limited to the inferior temporal and the frontal WM. We note that this does not include parietal WM, which did show significantly increased WMSA in AD individuals over the OC group, and has also been shown by other groups to be predictive of AD conversion [26]. We interpret this finding to suggest that perhaps parietal WMSA follow a different time course than temporal and frontal WMSA, as we see a significant relationship in these regions within the two years prior to AD conversion, but the duration of MCI can be much longer depending on the individual. Additionally, with greater sample sizes within each of the MCI subgroups, significant relationships may be revealed as trend-level relationships are seen in parietal WMSA but were not statistically significant, and we acknowledge the small MCI sample sizes as a limitation of the present study which likely contributes to the high variance in MCI subgroups. Taken together, however, these findings indicate that regionally-specific WMSA may in fact be a critical component of AD development as their progression can be viewed prior to AD onset and is directly related to time-to-conversion.
Interestingly, the regions demonstrating higher WMSA in AD as well as a relationship in the MCI individuals regardless of the WMSA metric used and after normalizing for total WMSA align to known cerebrovascular boundary zones (also known as watershed regions; Supplementary Figure 1), which are areas that are at the border of blood supply zones from two different major cerebral arteries and are thus most susceptible to hypoperfusion [56]. Specifically, these include superior frontal, supramarginal, inferior temporal, and cuneus regions. Furthermore, these regions are in close proximity to cortical regions in which cerebral blood flow (CBF) is known to decrease with advancing age [57]. The relationship between reductions in WM integrity and reductions in cortical CBF has also been previously demonstrated [58, 59], supporting the idea that the effects seen in these regions may be linked to neighboring cortical perfusion changes. We therefore interpret our findings to suggest that there may be vascular changes in AD that are at least regionally accelerated relative to global changes in comparison to those seen in normal aging. Our initial findings with WMSA, hypertension, and AD diagnosis suggest that this relationship is complex and future work will examine whether or not there are interactive effects between vascular risk factors and a clinical diagnosis of AD on the manifestation of WMSA, and how these effects are spatially represented.
Overall, the current data demonstrate a regionally stereotyped increase in WMSA burden with increasing overall burden in cognitively healthy older adults as well as adults with AD. There is a widespread increase in WMSA in individuals with AD, although the spatial extent of these increases is reduced when regional WMSA is normalized for total WMSA burden and are limited to inferior temporal, inferior parietal, cuneus, and middle frontal WM. Our results indicate that regionally-specific WMSA may be an important pathological component of AD development and we hypothesize that these WMSA may be linked to a vascular component of diseaseprogression.
Footnotes
ACKNOWLEDGMENTS
Support for this work was provided by the NIHNeuroimaging Training Program Grant #T32 EB001680 from the National Institute of Biomedical Imaging and Bioengineering (NIBIB), a training grant from the NIH Blueprint for Neuroscience Research (T90DA022759/R90DA023427), and National Institute for Nursing Research R01NR010827. Additional support was provided in part by the National Institute for Biomedical Imaging and Bioengineering (P41EB015896, R01EB006758, R21EB018907, R01EB019956), the National Institute on Aging (5R01AG008122, R01AG016495), the National Institute of Diabetes and Digestive and Kidney Diseases (1-R21-DK-108277-01), the National Institute for Neurological Disorders and Stroke (R01NS0525851, R21NS072652, R01NS070963, R01NS083534, 5U01NS086625), and was made possible by the resources provided by Shared Instrumentation Grants 1S10RR023401, 1S10RR019307, and 1S10RR023043. Additional support was provided by the NIH Blueprint for Neuroscience Research (5U01-MH093765), part of the multi-institutional Human Connectome Project.
Data collection and sharing for this project was funded by the Alzheimer’s Disease Neuroimaging Initiative (ADNI) (National Institutes of Health Grant U01 AG024904) and DOD ADNI (Department of Defense award number W81XWH-12-2-0012). ADNI is funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering, and through generous contributions from the following: AbbVie, Alzheimer’s Association; Alzheimer’s Drug Discovery Foundation; Araclon Biotech; BioClinica, Inc.; Biogen; Bristol-Myers Squibb Company; CereSpir, Inc.; Cogstate; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; EuroImmun; F. Hoffmann-La Roche Ltd and its affiliated company Genentech, Inc.; Fujirebio; GE Healthcare; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research & Development, LLC.; Johnson & Johnson Pharmaceutical Research & Development LLC.; Lumosity; Lundbeck; Merck & Co., Inc.; Meso Scale Diagnostics, LLC.; NeuroRx Research; Neurotrack Technologies; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Piramal Imaging; Servier; Takeda Pharmaceutical Company; and Transition Therapeutics. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health (
). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer’s Therapeutic Research Institute at the University of Southern California. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of Southern California.
