Original Research ARTICLE
Tau Pathology Distribution in Alzheimer's disease Corresponds Differentially to Cognition-Relevant Functional Brain Networks
- 1Clinical Memory Research Unit, Department of Clinical Sciences, Lund University, Malmö, Sweden
- 2Memory Clinic, Skåne University Hospital, Malmö, Sweden
- 3German Center for Neurodegenerative Diseases, Rostock, Germany
- 4Department of Radiation Physics, Skåne University Hospital, Lund, Sweden
- 5Department of Clinical Neurophysiology, Skåne University Hospital, Lund, Sweden
- 6Department of Clinical Physiology and Nuclear Medicine, Skåne University Hospital, Lund, Sweden
- 7Department of Neurology, Skåne University Hospital, Lund, Sweden
- 8MedTech West and the Department of Psychiatry and Neurochemistry, University of Gothenburg, Gothenburg, Sweden
Neuropathological studies have shown that the typical neurofibrillary pathology of hyperphosphorylated tau protein in Alzheimer's disease (AD) preferentially affects specific brain regions whereas others remain relatively spared. It has been suggested that the distinct regional distribution profile of tau pathology in AD may be a consequence of the intrinsic network structure of the human brain. The spatially distributed brain regions that are most affected by the spread of tau pathology may hence reflect an interconnected neuronal system. Here, we characterized the brain-wide regional distribution profile of tau pathology in AD using 18F-AV 1451 tau-sensitive positron emission tomography (PET) imaging, and studied this pattern in relation to the functional network organization of the human brain. Specifically, we quantified the spatial correspondence of the regional distribution pattern of PET-evidenced tau pathology in AD with functional brain networks characterized by large-scale resting state functional magnetic resonance imaging (rs-fMRI) data in healthy subjects. Regional distribution patterns of increased PET-evidenced tau pathology in AD compared to controls were characterized in two independent samples of prodromal and manifest AD cases (the Swedish BioFINDER study, n = 44; the ADNI study, n = 35). In the BioFINDER study we found that the typical AD tau pattern involved predominantly inferior, medial, and lateral temporal cortical areas, as well as the precuneus/posterior cingulate, and lateral parts of the parietal and occipital cortex. This pattern overlapped primarily with the dorsal attention, and to some extent with higher visual, limbic and parts of the default-mode network. PET-evidenced tau pathology in the ADNI replication sample, which represented a more prodromal group of AD cases, was less pronounced but showed a highly similar spatial distribution profile, suggesting an earlier-stage snapshot of a consistently progressing regional pattern. In conclusion, the present study indicates that the regional deposition of tau aggregates in AD predominantly affects higher-order cognitive over primary sensory-motor networks, but does not appear to be specific for the default-mode or related limbic networks.
The brain can be subdivided into a number of intrinsic connectivity networks (ICNs), which have been established mainly through interpretation of large resting-state functional magnetic resonance imaging (rs-fMRI) datasets (Power et al., 2011; Yeo et al., 2011). These networks have been found to be disrupted in a distinct manner in Alzheimer's disease (AD), directly affecting cognitive performance (Pievani et al., 2011). The accumulation of fibrillar amyloid-β (Aβ), one of the major traits of AD pathogenesis, has consistently been demonstrated to follow spatial patterns (Braak and Braak, 1991; Villeneuve et al., 2015) that overlap strongly with the so-called default-mode network (DMN; Buckner et al., 2009; Sperling et al., 2009), a network of brain regions generally activated in “task-negative” states but also during autobiographical and introspective tasks (Greicius et al., 2003, 2004). However, a recent study generally confirming these findings, also observed considerable spatial overlap of Aβ accumulating regions with other ICNs, such as a frontoparietal-control network (FPN) and the dorsal-attention network (DAN; Grothe and Teipel, 2016). This study further found substantial dissociation between correspondence of imaging-derived measures of Aβ pathology, measures of neurodegeneration proxied by glucose hypometabolism and gray matter atrophy, and ICNs (Grothe and Teipel, 2016).
Hyperphosphorylated tau, another misfolded protein to accumulate in AD, is known to be closer related to neurodegeneration and cognitive impairment than Aβ (Serrano-Pozo et al., 2011; Nelson et al., 2012). The recent introduction of tau-sensitive positron emission tomography (PET) ligands binding to paired-helical filaments of tau protein has enabled the in vivo assessment of hyperphosphorylated tau pathology (Xia et al., 2013; Marquié et al., 2015). Neuropathological studies have established a rather consistent spatial pattern of progressing tau pathology that has been categorized into so-called Braak stages (Braak and Braak, 1991; Braak et al., 2006). These nested stages of progressing tau pathology could recently be replicated in vivo using tau PET (Schwarz et al., 2016; Schöll et al., 2016).
Studies in cell and transgenic mouse models further suggest that tau in its harmful hyperphosphorylated form likely spreads in a trans-synaptic manner along neuronal networks, eventually causing synaptic damage, functional disruption, and neurodegeneration of the affected networks (Fox et al., 2011; Liu et al., 2012; Spires-Jones and Hyman, 2014; Menkes-Caspi et al., 2015).
In the present study, we aimed at exploring the spatial overlap of pre-defined functional brain networks and the extent of neurofibrillary tau pathology assessed with 18F-AV-1451 PET in prodromal and dementia stages of AD, also answering the question whether tau pathology preferrably affects the DMN over other ICNs.
We studied samples from two independent cohorts. The first sample was recruited from the prospective and longitudinal Swedish BioFINDER study (further information available at: www.BioFINDER.se). We included data from 17 cognitively healthy elderly control (HC) participants and 27 AD patients with prodromal (MCI due to AD, n = 11) or clinically manifest disease (AD dementia, n = 16). Inclusion criteria for HC were: (1) aged ≥ 60 years old, (2) scored 28–30 points on the Mini-Mental State Examination (MMSE) at the screening visit, (3) absence of cognitive symptoms as evaluated by a physician, (4) fluent in Swedish, (5) did not fulfill the criteria of MCI or any dementia; with the following exclusion criteria: (1) significant systemic illness making it difficult to participate, (2) presence of significant neurologic or psychiatric disease (e.g., stroke, Parkinson's disease, multiple sclerosis, major depression), (3) significant substance abuse. Inclusion criteria for MCI due to AD were: (1) age 40–100 years, (2) referred to the memory clinics due to cognitive symptoms experienced by the patient and/or informant, (3) MMSE score of 24–30, (4) objective impairment according to neuropsychological testing, (5) no fulfillment of the criteria for any dementia disorder (major neurocognitive disorder) according to DSM-V, (6) abnormal cerebrospinal fluid (CSF) Aβ42 indicative of prodromal AD and (7) fluency in Swedish; and the exclusion criteria were: (1) cognitive impairment that without doubt could be explained by another condition (other than prodromal AD) and (2) severe somatic disease. AD dementia patients met the DSM-IIIR criteria for dementia (American Psychiatric Association and American Psychiatric Association Work Group to Revise DSM-III., 1987) as well as the NINCDS-ADRDA criteria for AD (McKhann et al., 1984). Exclusion criteria were significant systemic illness and significant alcohol abuse. AD diagnoses were confirmed by physicians who were blinded to any PET and CSF data. We only recruited AD patients with a late-onset amnestic clinical presentation to ensure a typical pattern of tau ligand uptake as differences in ligand retention between different AD variants had previously been reported (Ossenkoppele et al., 2016). All participants provided written informed consent to participate in the study according to the Declaration of Helsinki, ethical approval was given by the Ethics Committee of Lund University, Lund, Sweden, and all methods were carried out in accordance with the approved guidelines. Approval for PET imaging was obtained from the Swedish Medicines and Products Agency and the local Radiation Safety Committee at Skåne University Hospital, Sweden.
The second sample was recruited from the ADNI2 cohort, an extension of the original ADNI1 and ADNI-GO studies. ADNI (Alzheimer's Disease Neuroimaging Initiative) is a multisite longitudinal biomarker study that has enrolled over 1,500 cognitively normal older individuals, people with early or late amnestic MCI, and people with early AD (www.adni-info.org). Detailed inclusion and exclusion criteria for the diagnostic categories can be found on the ADNI website (http://adni.loni.usc.edu/methods/). Briefly, healthy subjects had MMSE scores between 24 and 30 (inclusive), a clinical dementia rating (CDR) = 0, were non-depressed, non-MCI, and non-demented. MCI subjects had MMSE scores between 24 and 30 (inclusive), a subjective memory concern reported by subject, informant, or clinician, objective memory loss measured by education adjusted scores on delayed recall, a CDR = 0.5, absence of significant levels of impairment in other cognitive domains, essentially preserved activities of daily living, and an absence of dementia. Subjects with AD dementia had initial MMSE scores between 20 and 26 (inclusive), a CDR = 0.5 or 1.0, and fulfill NINCDS-ADRDA criteria for clinically probable Alzheimer's disease (McKhann et al., 1984). For the present study, we recruited data from 19 HC, 26 prodromal, and four patients with manifest AD.
Diagnostic groups were dichotomized into β-amyloid-positive (Aβ+) and -negative (Aβ−) subgroups, based on 18F-flutemetamol (BioFINDER) or 18F-AV45 (ADNI) PET evidence of global Aβ pathology indicative of AD. For ADNI subjects, cortex-to-whole cerebellum AV45 standard uptake value ratios (SUVR) had been calculated and made available on the ADNI server by the Jagust Lab, UC Berkeley. Aβ-positivity was established using the recommended threshold for cross-sectional analyses of SUVR ≥ 1.11 (Landau et al., 2013). For BioFINDER subjects, global 18F-flutemetamol SUVR had been calculated using a composite whole cerebellum, the pons/brainstem region, and eroded cortical white matter reference region; a cutoff of SUVR ≥ 0.79 was derived using a mixture modeling analysis in a large BioFINDER cohort (n = 406) to describe Aβ-positivity.
Given that we aimed at examining a typical late-onset amnestic AD type, and with all recruited BioFINDER MCI and AD subjects being Aβ-positive, Aβ-negative subjects were omitted from the ADNI2 sample, which resulted in a final ADNI sample size of 19 HC and 16 prodromal dementia and clinically manifest AD patients (12 Aβ + MCI, 4 Aβ + AD; Albert et al., 2011; McKhann et al., 2011; Sperling et al., 2011).
All BioFINDER subjects underwent structural magnetic resonance imaging (MRI) on a Siemens Tim Trio 3T scanner (Siemens Medical Solutions, Erlangen, Germany). High resolution T1-weighted anatomical magnetization-prepared rapid gradient echo (MPRAGE) images (TR = 1950 ms TE = 3.4 ms, 1 mm isotropic voxels and 176 slices) were acquired for PET image co-registration, processing, and template normalization.
18F-AV-1451 PET scans were performed on a GE Discovery 690 PET scanner (General Electric Medical Systems) as dynamic scans using LIST-mode 80–120 min after a bolus injection of 370 MBq of 18F-AV-1451. Low-dose CT-scans for attenuation correction were performed immediately prior to the PET scans. PET data were reconstructed into 5 min frames using an iterative Vue Point HD algorithm with six subsets, 18 iterations with 3 mm filter, and no time-of-flight correction. The dynamic scans were motion corrected using AFNI's 3dvolreg (Cox, 1996), time-averaged, and rigidly co-registered to the skull-stripped MRI scan.
For ADNI participants, T1-weighted MPRAGE MR images were acquired on multiple 3T scanners. 18F-AV-1451 PET scans were also acquired on several different scanners at multiple sites. Subjects were examined for 30 min (6 × 5 min frames) starting at 75 min post-injection of a 370 MBq/kg bolus. Standardized image pre-processing is applied to all original ADNI scans (see http://adni.loni.usc.edu/methods for details).
Preprocessed ADNI image data was downloaded from the ADNI image database (https://ida.loni.usc.edu) and further processed in the same manner as BioFINDER image data (see below).
All image data was thus processed at Lund University employing an in-house developed pipeline. The MR scans were normalized to a common MNI152 space (Montreal Neurological Institute) with a diffeomorphic transform using the Advanced Normalization Tools (ANTs) toolbox (Avants et al., 2014) for further use in the PET processing pipeline. Cortical reconstruction and volumetric segmentation and parcellation were performed with Freesurfer v5.3 (http://surfer.nmr.mgh.harvard.edu). Reconstructed data sets were visually inspected for inaccuracies, and major segmentation errors were manually corrected. The Freesurfer parcellation in the MR space of the anatomical scan was then applied to the processed, coregistered, and time-averaged PET image to extract reference regional uptake values.
We created BioFINDER 18F-AV-1451 standardized uptake value (SUV) PET images based on mean uptake over 80–120 min postinjection normalized to uptake in a gray matter masked cerebellum reference region to create voxelwise SUV ratio (SUVR) images in each participant's MRI native space. ADNI 18F-AV-1451 SUVR images were created based on mean uptake over 80–100 min post-injection and intensity normalized using the same cerebellar reference region as for the BioFINDER data. All PET images were then spatially normalized to MNI152 space, employing the ANTs spatial transformation parameters derived from the co-registered MR scans, and smoothed with a 8 mm FWHM Gaussian filter. Using FSL (v5.0.6, http://fsl.fmrib.ox.ac.uk), we finally created a mean and a standard deviation image based on the respective control group's 18F-AV-1451 images in MNI152 standard space and calculated individual Z-score maps for each patient (Z-score = (individual value—control mean)/control standard deviation; Chételat et al., 2008; Grothe and Teipel, 2016).
Analysis of the network-specificity of the AD-typical in vivo distribution of tau deposits as imaged by 18F-AV-1451 PET closely followed the approach of a previously published study assessing the network-specificity of AD-related amyloid deposition, hypometabolism, and gray matter atrophy (Grothe and Teipel, 2016). Briefly, this approach quantitatively assesses the correspondence of the cortex-wide pathologic imaging pattern with functional networks in the human brain as defined by standardized maps of ICN derived from large-scale resting-state rs-fMRI data of a healthy adult population (Yeo et al., 2011). Quantitative metrics include the mean Z-score within each ICN template, reflecting the extent of tau deposition in relation to control values, as well as a goodness-of-fit (GOF) score for each ICN template quantifying the spatial correspondence of the cortex-wide tau deposition pattern with the functional network topography. GOF scores are calculated as the difference between the mean Z-score value of voxels falling within a given ICN template and the mean Z-score value of cortical voxels outside the ICN template (Greicius et al., 2004; Lehmann et al., 2013). Thus, a positive GOF-score indicates a relative preference of cortex-wide tau deposition to occur within the respective ICN, whereas a uniform distribution of tau deposition across the cortex would result in GOF-scores close to zero for all ICNs.
One sample t-tests were used to assess the significance of ICN-specific increases in tau deposition in the patient groups compared to the respective control groups (mean Z-scores per ICN), as well as the significance of the spatial correspondence of the AD-typical tau deposition with a given ICN (positive GOF-scores). Statistical significance was set at p < 0.05 (two-tailed), Bonferroni-corrected for the number of networks assessed.
In our primary analysis we used an ICN definition based on a recently published functional parcellation scheme of the human brain into seven major ICNs, including the 1, DMN; 2, frontoparietal-control network (FPN); 3, dorsal attention network (DAN); 4, ventral attention network (VAN); 5, limbic network (LIM); 6, visual network (VIS) and 7, somatomotor network (SMN; Figure 1, Supplementary Figure 1; Yeo et al., 2011; https://sites.google.com/site/yeoyeo02/software). While these principal large-scale ICNs are consistently reproduced in rs-fMRI based functional parcellations across several independent cohorts and using diverse parcellation methods (Greicius et al., 2003; Damoiseaux et al., 2006; Fox et al., 2006; Cohen et al., 2008; Kahn et al., 2008; Vincent et al., 2008; Smith et al., 2009; Bellec et al., 2010; Power et al., 2011; Jones et al., 2012; Das et al., 2015; Pascual et al., 2015), there is currently no established way of unambiguously defining the most appropriate number of separate connectivity modules within the brain's functional connectivity architecture, and at higher parcellation resolutions, the principal large-scale ICNs typically fractionate into further submodules (Andrews-Hanna et al., 2010; Power et al., 2011). Thus, in a secondary analysis we assessed the network-specificity of the AD-typical tau deposition pattern using a more fine-grained parcellation scheme into 17 functional subnetworks (Yeo et al., 2011): 1, Higher visual; 2, Primary visual; 3, Dorsal SMN; 4, Ventral SMN-auditory; 5, Posterior DAN; 6, Frontal eye field-DAN; 7, Posterior VAN; 8, Anterior VAN; 9, Temporal pole-Anterior medial temporal lobe (MTL) LIM; 10, Orbitofrontal LIM; 11, Precuneus non-DMN; 12, FPN component 1; 13, FPN component 2; 14, Lateral-temporal DMN/language; 5, Posterior MTL-retrosplenial DMN; 16, Midline DMN; 17, Anterior DMN (these codes are consequently used in graphs and tables).
Figure 1. Brain-wide patterns of increased tau deposition in the BioFINDER and ADNI cohorts and illustration of the standardized intrinsic connectivity networks (for details see Supplementary Figure 1). The top and bottom row images represent mean Z-score images of the BioFINDER and ADNI patient samples. Note the different scales used for color-coded visualization of mean Z-score values in the BioFINDER and ADNI patient samples.
The spatial correspondence between regional tau distribution profiles in the independent BioFINDER and ADNI samples was assessed using Pearson's correlation across all cortical voxels of the mean Z-score maps from the respective cohorts (Buckner et al., 2009).
Demographics for all participants are displayed in Table 1. There was no statistically significant difference in age between the BioFINDER groups, but the ADNI HC were significantly younger than the ADNI patient group. Patient groups in both BioFINDER and ADNI performed significantly worse on the MMSE when compared to the respective HC.
Consistent Brain-Wide Pattern of Increased 18F-AV-1451 Uptake in AD
Figure 1 shows mean Z-score images of the AD patients from the BioFINDER and the ADNI cohort, respectively. The BioFINDER sample with its higher proportion of manifest AD dementia cases displayed clear-cut 18F-AV-1451 uptake in the lateral and inferior temporal lobes with partial involvement of the lateral parietal and occipital cortices, as well as in the precuneus and posterior cingulate. The ADNI sample with its greater proportion of prodromal AD cases exhibited less pronounced ligand uptake, albeit demonstrating a very similar spatial distribution profile across all cortical voxels (r = 0.79, p < 0.001; Figures 1, 2).
Figure 2. Consistency of tau deposition pattern across cohorts. Mean voxelwise cortical Z-scores for the BioFINDER and ADNI sample were highly correlated.
The Brain-Wide Pattern of Increased 18F-AV-1451 Uptake in AD Overlaps with Several Functional Brain Networks Related to Higher Cognition
Increased AV-1451 uptake was observed in all seven main ICNs at a corrected alpha level = 0.05/7 = 0.007 (Z-scores different from 0, Table 2, Figure 3A), but AV-1451 retention was not homogeneously distributed across these networks. The gradient of severity of tau accumulation followed the order: DAN > VIS > DMN > LIM > FPN > VAN > SMN. Accordingly, positive GOF-scores, indicating a preferential tau accumulation in relation to brain-wide uptake, were observed for the DAN (p < 0.001), with trends (p < 0.05) for the VIS (p = 0.03), and the DMN (p = 0.04; Table 2, Figure 3B).
Table 2. Z- and Goodness-of-fit scores for AV-1451 uptake in seven large-scale networks in the BioFINDER sample.
Figure 3. Mean Z- and Goodness-of-fit scores for each of the seven main networks in the BioFINDER (A,B) and ADNI (C,D) samples. Error bars represent 1 standard error.
Fifteen of the 17 subnetworks showed significantly increased AV-1451 uptake at a corrected alpha level = 0.05/17 = 0.003 (Z-scores different from 0, Table 3, Figure 4A). Most clear AV-1451 uptake was observed in nodes of the posterior DAN (5), a lateral temporal/language subnetwork of the DMN (14), a medial parietal/precuneal component of the FPN (11), a higher visual network (1), and a posterior MTL/retrosplenial subnetwork of the DMN (15) (all z-score >5). Least pronounced uptake was seen in the dorsal and ventral SMN (3, 4) and the orbitofrontal node of the limbic network (10) (all z-scores ≤2.5).
Table 3. Z- and Goodness-of-fit scores for AV-1451 uptake in 17 sub-networks in the BioFINDER sample.
Figure 4. Mean Z- and Goodness-of-fit scores for each of the 17 sub-networks in the BioFINDER (A,B) and ADNI (C,D) samples. Error bars represent 1 standard error.
Statistically significant positive GOF scores (p < 0.003) were observed for a total of six subnetworks: posterior DAN (5), lateral-temporal language (14), precuneus non-DMN (11), higher visual (1), posterior MTL-retrosplenial (15), and medial parietal and lateral frontoparietal components (12; Table 3, Figure 4B).
At a corrected alpha-level of p < 0.007, we observed no significantly increased AV-1451 uptake compared to the control group within any of the seven large-scale brain networks. However, in accordance with the very high similarity in the voxel-wise regional distribution profile (Figure 2), the overall rank order of network involvement was very similar to the pattern in the BioFINDER sample (DAN > VIS > LIM > DMN > FPN > VAN > SMN), and trends for increased AV-1451 uptake compared to the control group were noted in the DAN (p = 0.045), LIM (p = 0.057), and DMN (p = 0.062; Table 4, Figure 3C). The only network that showed a trend for a significantly positive GOF score was the DAN (p = 0.03; Table 4, Figure 3D).
Table 4. Z- and Goodness-of-fit scores for AV-1451 uptake in seven large-scale networks in the ADNI sample.
At a corrected alpha-level of p < 0.003, none of the 17 subnetworks showed significantly increased tau uptake compared to the control group. However, similar to the findings in the BioFINDER sample, trends (p < 0.05) were observed for the posterior DAN (5), a medial parietal/precuneal component of the FPN (11), as well as the temporal pole-anterior MTL limbic subnetwork (9), and posterior MTL/retrosplenial (15) and midline (16) subnetworks of the DMN. Further, similar to the findings in the BioFINDER sample, the weakest uptake was seen in the dorsal and ventral somatomotor networks (3, 4), the posterior VAN (7), and the orbitofrontal node of the limbic network (10) (Table 5, Figure 4C). A significant positive GOF score was only observed for the posterior DAN (p = 0.003; Table 5, Figure 4D).
In the present study, we examined the spatial distribution of 18F-AV-1451 retention, likely representing the presence of hyperphosporylated tau pathology, in the brains of prodromal and clinically manifest AD patients in relation to the spatial extent of predefined templates of functional brain networks. In the past years, emerging evidence for a consistent regional deposition pattern of in vivo tau PET ligands in AD supports a typical involvement of the inferior and lateral temporal lobes, precuneus and posterior cingulate, as well as occipital and lateral parietal lobes (Brier et al., 2016; Ishiki et al., 2015; Johnson et al., 2016; Ossenkoppele et al., 2016; Schwarz et al., 2016). Our voxel-wise (average Z-score) maps of regional AV-1451 distribution reflect these reports of regional tau PET-ligand retention in independent AD cohorts. Accordingly, the regional tau distribution profiles were very similar between the independent BioFINDER and ADNI cohorts of our study, despite obvious differences in recruitment criteria, PET scanning platforms, disease severity, and age. The ADNI sample featuring a greater proportion of prodromal AD cases consequently demonstrated a slightly alleviated AV-1451 signal when compared to the BioFINDER sample.
Previous research has suggested that AD pathology and related brain atrophy predominantly target the DMN (Seeley et al., 2009), a finding we could not confirm for the in vivo measures of AD-related tau pathology in our study. This is in line with a previous investigation of network-specific atrophy measures in AD (Grothe and Teipel, 2016), indicating little specificity for the entire DMN, but rather a distinct spatial overlap with an anterior limbic subnetwork. Here, we found increased tau deposition in limbic and DMN network components, but especially posterior cortical networks (dorsal attention and higher visual networks). A recent publication reported a significant correlation between regional AV-1451 ligand uptake and FDG-PET measures of neurodegeneration across AD patients of differing clinical phenotypes (Ossenkoppele et al., 2015, 2016). However, spatial correspondence of the two imaging modalities was far from being complete (r = −0.49 to −0.60), and future studies will have to further explore similarities and differences in the in vivo spatial distribution profiles of AD-related tau pathology and neurodegeneration, as well as their respective relations to the spatial topography of intrinsic brain networks (Sepulcre et al., 2016).
The posterior cortical predominance of in vivo tau deposition, also corroborated by the above-mentioned reports, contrasts partly with neuropathological assessments of regional tau deposition severity in post mortem tissue (Braak and Braak, 1991; Braak et al., 2006) where no such posterior predominance has been reported. It should be taken into account, however, that neuropathological sampling techniques commonly examine selected slices from different brain regions and are not necessarily comparable to an averaged quantification of tau PET ligand uptake within a specific brain region. Moreover, the distinct distribution profile of AV-1451 may be determined by a specific binding to certain types of tau aggregates, such as indicated by a recent study showing a relatively selective binding of AV-1451 to neurofibrillary tangles as compared to the paired helical tau filaments found in neuropil threads and neuritic plaques (Braak et al., 1986; Ono et al., 2017).
Nonetheless, tau deposition in AD might exhibit regional non-linearity, occurring early in allocortical regions of the MTL, but in greater absolute amounts at later disease stages in isocortical areas. This is in line with the here observed relatively more pronounced involvement of MTL/limbic subnetworks (e.g., 9, Temporal pole-Anterior MTL LIM and 15, Posterior MTL-retrosplenial DMN) in the more prodromal ADNI sample compared to BioFINDER sample (see Figure 4). Moreover, in our study tau deposition was not quantified absolutely but relative to cognitively healthy controls, a group that has previously been reported to exhibit significant tau deposition in MTL and related limbic circuits (primary age-related tauopathy, PART; Crary et al., 2014; Schöll et al., 2016), which may result in an underestimation of absolute MTL/limbic AV-1451 retention in patient groups as represented in our mean Z-maps.
On a methodological note, tau-sensitive PET ligands might exhibit regionally differing saturation rates of ligand uptake, leading to over- or underestimation of actual regional tau pathology by interpretation of tau PET data; however, this has hitherto only been reported for another tau PET ligand, not AV-1451 (Lemoine et al., 2015; Marquié et al., 2015; Smith et al., 2016). In addition, several nuisance factors could have influenced our results and have not been taken into account when interpreting our results. Both study samples were relatively small and were not matched for age, education level, or disease severity. We used plain Z-scores to examine spatial AV-1451 retention patterns, not adjusting for any covariates which might have affected group differences.
According to our present findings, AD-typical tau deposition predominantly targets higher-order cognitive networks over primary sensory-motor networks, but is not specific for the DMN or any other single large-scale functional brain network as a whole. This finding has implications for popular models of network specific spread of tau pathology derived from observations of prion-like mechanisms of transsynaptic tau transmission in animal models. According to these models, tau deposition within a given seed region would primarily spread within the interconnected network of this brain region, before it spreads to separate networks (Fox et al., 2011; Liu et al., 2012; Spires-Jones and Hyman, 2014; Menkes-Caspi et al., 2015). Our finding of a disproportionate affection of specific submodules within large-scale functional brain networks could be explained by a combination of (i) transsynaptic spread process and (ii) inherent regional differences in the vulnerability to tau accumulation. Tau may spread to all intrinsic network connections of one specific seed region, but only accumulate in those areas that exhibit an inherent susceptibility to tau aggregation, possibly determined by its specific functional role and molecular architecture (Zhou et al., 2012; Freer et al., 2016). Furthermore, our observation of a fairly parallel affection of several distinct (higher-order) brain networks in AD may point to the existence of several parallel seed regions of tau progression, rather than a single disease “epicenter” (Seeley et al., 2009; Zhou et al., 2012).
In conclusion, the present study indicates that the regional deposition of hyperphosphorylated tau aggregates in AD does not specifically target the default-mode or related limbic networks, but more generally affects higher-order cognitive over primary sensory-motor networks.
MS, MG, and OH designed the study; TS and MS conducted data preprocessing; MG and MS performed data analyses; OH, TO, DH, JJ, and RS contributed with data acquisition; MS and MG drafted the manuscript; all authors revised and approved the final manuscript.
Conflict of Interest Statement
Dr. Hansson has served at advisory boards for Eli Lilly and received research support from GE Healthcare and Hoffmann LaRoche.
The other 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 authors would like to thank the collaborators of this study and the entire BioFINDER Study group (www.BioFINDER.se). Work in the authors' laboratory was supported by the European Research Council, the Swedish Research Council, the Strategic Research Area MultiPark (Multidisciplinary Research in Parkinson's disease) at Lund University, the Swedish Brain Foundation, the Skåne University Hospital Foundation, the Parkinson foundation of Sweden, the Swedish Alzheimer Association, and the Swedish federal government under the ALF agreement. The BioFINDER team wishes to thank Avid Radiopharmaceuticals for providing the precursor for AV-1451, and GE Healthcare for providing 18F-flutemetamol. The authors would further like to thank the ADNI investigators, and all study participants for their valuable contributions. 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 (www.fnih.org). 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.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fnins.2017.00167/full#supplementary-material
Supplementary Figure 1. Intrinsic functional brain network segmentation based on rs-MRI data (Yeo et al., 2011).
Albert, M. S., Dekosky, S. T., Dickson, D., Dubois, B., Feldman, H. H., Fox, N. C., et al. (2011). The diagnosis of mild cognitive impairment due to Alzheimer's disease: recommendations from the national institute on aging-alzheimer's association workgroups on diagnostic guidelines for Alzheimer's disease. Alzheimers. Dement. 7, 270–279. doi: 10.1016/j.jalz.2011.03.008
Andrews-Hanna, J. R., Reidler, J. S., Sepulcre, J., Poulin, R., and Buckner, R. L. (2010). Functional-anatomic fractionation of the brain's default network. Neuron 65, 550–562. doi: 10.1016/j.neuron.2010.02.005
Avants, B. B., Libon, D. J., Rascovsky, K., Boller, A., McMillan, C. T., Massimo, L., et al. (2014). Sparse canonical correlation analysis relates network-level atrophy to multivariate cognitive measures in a neurodegenerative population. Neuroimage 84, 698–711. doi: 10.1016/j.neuroimage.2013.09.048
Bellec, P., Rosa-Neto, P., Lyttelton, O. C., Benali, H., and Evans, A. C. (2010). Multi-level bootstrap analysis of stable clusters in resting-state FMRI. Neuroimage 51, 1126–1139. doi: 10.1016/j.neuroimage.2010.02.082
Braak, H., Alafuzoff, I., Arzberger, T., Kretzschmar, H., and Del Tredici, K. (2006). Staging of Alzheimer disease-associated neurofibrillary pathology using paraffin sections and immunocytochemistry. Acta Neuropathol. 112, 389–404. doi: 10.1007/s00401-006-0127-z
Braak, H., Braak, E., Grundke-Iqbal, I., and Iqbal, K. (1986). Occurrence of neuropil threads in the senile human brain and in Alzheimer's disease: a third location of paired helical filaments outside of neurofibrillary tangles and neuritic plaques. Neurosci. Lett. 65, 351–355. doi: 10.1016/0304-3940(86)90288-0
Brier, M. R., Gordon, B., Friedrichsen, K., McCarthy, J., Stern, A., Christensen, J., et al. (2016). Tau and abeta imaging, CSF measures, and cognition in Alzheimer's disease. Sci. Transl. Med. 8:338ra66. doi: 10.1126/scitranslmed.aaf2362
Buckner, R. L., Sepulcre, J., Talukdar, T., Krienen, F. M., Liu, H., Hedden, T., et al. (2009). Cortical hubs revealed by intrinsic functional connectivity: mapping, assessment of stability, and relation to Alzheimer's disease. J. Neurosci. 29, 1860–1873. doi: 10.1523/JNEUROSCI.5062-08.2009
Chételat, G., Desgranges, B., Landeau, B., Mézenge, F., Poline, J. B., de La Sayette, V., et al. (2008). Direct voxel-based comparison between grey matter hypometabolism and atrophy in Alzheimer's disease. Brain 131, 60–71. doi: 10.1093/brain/awm288
Cohen, A. L., Fair, D. A., Dosenbach, N. U., Miezin, F. M., Dierker, D., Van Essen, D. C., et al. (2008). Defining functional areas in individual human brains using resting functional connectivity MRI. Neuroimage 41, 45–57. doi: 10.1016/j.neuroimage.2008.01.066
Crary, J. F., Trojanowski, J. Q., Schneider, J. A., Abisambra, J. F., Abner, E. L., Alafuzoff, I., et al. (2014). Primary age-related tauopathy (part): a common pathology associated with human aging. Acta Neuropathol. 128, 755–766. doi: 10.1007/s00401-014-1349-0
Damoiseaux, J. S., Rombouts, S. A., Barkhof, F., Scheltens, P., Stam, C. J., Smith, S. M., et al. (2006). Consistent resting-state networks across healthy subjects. Proc. Natl. Acad. Sci. U.S.A. 103, 13848–13853. doi: 10.1073/pnas.0601417103
Das, S. R., Pluta, J., Mancuso, L., Kliot, D., Yushkevich, P. A., and Wolk, D. A. (2015). Anterior and posterior MTL networks in aging and MCI. Neurobiol Aging 36(Suppl. 1) S141-50, S150 ES1. doi: 10.1016/j.neurobiolaging.2014.03.041
Fox, L. M., William, C. M., Adamowicz, D. H., Pitstick, R., Carlson, G. A., Spires-Jones, T. L., et al. (2011). Soluble tau species, not neurofibrillary aggregates, disrupt neural system integration in a tau transgenic model. J. Neuropathol. Exp. Neurol. 70, 588–595. doi: 10.1097/NEN.0b013e318220a658
Fox, M. D., Corbetta, M., Snyder, A. Z., Vincent, J. L., and Raichle, M. E. (2006). Spontaneous neuronal activity distinguishes human dorsal and ventral attention systems. Proc. Natl. Acad. Sci. U.S.A. 103, 10046–10051. doi: 10.1073/pnas.0604187103
Freer, R., Sormanni, P., Vecchi, G., Ciryam, P., Dobson, C. M., and Vendruscolo, M. (2016). A protein homeostasis signature in healthy brains recapitulates tissue vulnerability to Alzheimer's disease. Sci. Adv. 2:E1600947. doi: 10.1126/sciadv.1600947
Greicius, M. D., Krasnow, B., Reiss, A. L., and Menon, V. (2003). Functional connectivity in the resting brain: a network analysis of the default mode hypothesis. Proc. Natl. Acad. Sci. U.S.A. 100, 253–258. doi: 10.1073/pnas.0135058100
Greicius, M. D., Srivastava, G., Reiss, A. L., and Menon, V. (2004). Default-mode network activity distinguishes Alzheimer's disease from healthy aging: evidence from functional MRI. Proc. Natl. Acad. Sci. U.S.A. 101, 4637–4642. doi: 10.1073/pnas.0308627101
Grothe, M. J., and Teipel, S. J. (2016). Spatial patterns of atrophy, hypometabolism, and amyloid deposition in Alzheimer's disease correspond to dissociable functional brain networks. Hum. Brain Mapp. 37, 35–53. doi: 10.1002/hbm.23018
Ishiki, A., Okamura, N., Furukawa, K., Furumoto, S., Harada, R., Tomita, N., et al. (2015). Longitudinal assessment of tau pathology in patients with Alzheimer's disease using [18f]Thk-5117 positron emission tomography. PLoS ONE 10:E0140311. doi: 10.1371/journal.pone.0140311
Johnson, K. A., Schultz, A., Betensky, R. A., Becker, J. A., Sepulcre, J., Rentz, D., et al. (2016). Tau positron emission tomographic imaging in aging and early Alzheimer disease. Ann. Neurol. 79, 110–119. doi: 10.1002/ana.24546
Jones, D. T., Vemuri, P., Murphy, M. C., Gunter, J. L., Senjem, M. L., Machulda, M. M. Jr., et al. (2012). Non-stationarity in the “Resting brain's” modular architecture. PLoS ONE 7:E39731. doi: 10.1371/journal.pone.0039731
Kahn, I., Andrews-Hanna, J. R., Vincent, J. L., Snyder, A. Z., and Buckner, R. L. (2008). Distinct cortical anatomy linked to subregions of the medial temporal lobe revealed by intrinsic functional connectivity. J. Neurophysiol. 100, 129–139. doi: 10.1152/jn.00077.2008
Landau, S. M., Breault, C., Joshi, A. D., Pontecorvo, M., Mathis, C. A., Jagust, W. J., et al. (2013). Amyloid-β imaging with pittsburgh compound B and florbetapir: comparing radiotracers and quantification methods. J. Nucl. Med. 54, 70–77. doi: 10.2967/jnumed.112.109009
Lehmann, M., Madison, C. M., Ghosh, P. M., Seeley, W. W., Mormino, E., Greicius, M. D., et al. (2013). Intrinsic connectivity networks in healthy subjects explain clinical variability in Alzheimer's disease. Proc. Natl. Acad. Sci. U.S.A. 110, 11606–11611. doi: 10.1073/pnas.1221536110
Lemoine, L., Saint-Aubert, L., Marutle, A., Antoni, G., Eriksson, J. P., Ghetti, B., et al. (2015). Visualization of regional tau deposits using (3)H-Thk5117 in Alzheimer brain tissue. Acta Neuropathol. Commun. 3:40. doi: 10.1186/s40478-015-0220-4
Marquié, M., Normandin, M. D., Vanderburg, C. R., Costantino, I. M., Bien, E. A., Rycyna, L. G., et al. (2015). Validating novel tau positron emission tomography tracer [F-18]-Av-1451 (T807) on postmortem brain tissue. Ann. Neurol. 78, 787–800. doi: 10.1002/ana.24517
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
McKhann, G. M., Knopman, D. S., Chertkow, H., Hyman, B. T., and Jack, C. R. Jr. (2011). The diagnosis of dementia due to Alzheimer's disease: recommendations from The National Institute on Aging-Alzheimer's association workgroups on diagnostic guidelines for Alzheimer's disease. Alzheimers Dement. 7, 263–269. doi: 10.1016/j.jalz.2011.03.005
Menkes-Caspi, N., Yamin, H. G., Kellner, V., Spires-Jones, T. L., Cohen, D., and Stern, E. A. (2015). Pathological tau disrupts ongoing network activity. Neuron 85, 959–966. doi: 10.1016/j.neuron.2015.01.025
Nelson, P. T., Alafuzoff, I., Bigio, E. H., Bouras, C., Braak, H., Cairns, N. J., et al. (2012). Correlation of Alzheimer disease neuropathologic changes with cognitive status: a review of the literature. J. Neuropathol. Exp. Neurol. 71, 362–381. doi: 10.1097/NEN.0b013e31825018f7
Ono, M., Sahara, N., Kumata, K., Ji, B., Ni, R., Koga, S., et al. (2017). Distinct binding of PET ligands PBB3 and AV-1451 to tau fibril strains in neurodegenerative tauopathies. Brain. doi: 10.1093/brain/aww339. [Epub ahead of print].
Ossenkoppele, R., Schonhaut, D. R., Baker, S. L., O'neil, J. P., Janabi, M., Ghosh, P. M., et al. (2015). Tau, Amyloid, and hypometabolism in a patient with posterior cortical atrophy. Ann. Neurol. 77, 338–342. doi: 10.1002/ana.24321
Ossenkoppele, R., Schonhaut, D. R., Schöll, M., Lockhart, S. N., Ayakta, N., Baker, S. L., et al. (2016). Tau pet patterns mirror clinical and neuroanatomical variability in Alzheimer's disease. Brain 139, 1551–1567. doi: 10.1093/brain/aww027
Pascual, B., Masdeu, J. C., Hollenbeck, M., Makris, N., Insausti, R., Ding, S. L., et al. (2015). Large-scale brain networks of the human left temporal pole: a functional connectivity MRI study. Cereb. Cortex 25, 680–702. doi: 10.1093/cercor/bht260
Power, J. D., Cohen, A. L., Nelson, S. M., Wig, G. S., Barnes, K. A., Church, J. A., et al. (2011). Functional network organization of the human brain. Neuron 72, 665–678. doi: 10.1016/j.neuron.2011.09.006
Schöll, M., Lockhart, S. N., Schonhaut, D. R., O'neil, J. P., Janabi, M., Ossenkoppele, R., et al. (2016). PET imaging of tau deposition in the aging human brain. Neuron 89, 971–982. doi: 10.1016/j.neuron.2016.01.028
Schwarz, A. J., Yu, P., Miller, B. B., Shcherbinin, S., Dickson, J., Navitsky, M., et al. (2016). Regional profiles of the candidate tau pet ligand 18f-Av-1451 recapitulate key features of braak histopathological stages. Brain 139, 1539–1550. doi: 10.1093/brain/aww023
Seeley, W. W., Crawford, R. K., Zhou, J., Miller, B. L., and Greicius, M. D. (2009). Neurodegenerative diseases target large-scale human brain networks. Neuron 62, 42–52. doi: 10.1016/j.neuron.2009.03.024
Sepulcre, J., Schultz, A. P., Sabuncu, M., Gomez-Isla, T., Chhatwal, J., Becker, A., et al. (2016). In vivo tau, amyloid, and gray matter profiles in the aging brain. J. Neurosci. 36, 7364–7374. doi: 10.1523/JNEUROSCI.0639-16.2016
Smith, R., Puschmann, A., Schöll, M., Ohlsson, T., Van Swieten, J., Honer, M., et al. (2016). 18f-Av-1451 Tau PET imaging correlates strongly with tau neuropathology in MAPT mutation carriers. Brain 139, 2372–2379. doi: 10.1093/brain/aww163
Smith, S. M., Fox, P. T., Miller, K. L., Glahn, D. C., Fox, P. M., Mackay, C. E., et al. (2009). Correspondence of the brain's functional architecture during activation and rest. Proc. Natl. Acad. Sci. U.S.A. 106, 13040–13045. doi: 10.1073/pnas.0905267106
Sperling, R. A., Aisen, P. S., Beckett, L. A., Bennett, D. A., Craft, S., Fagan, A. M., et al. (2011). Toward defining the preclinical stages of Alzheimer's disease: recommendations from The National Institute on Aging-Alzheimer's Association Workgroups on diagnostic guidelines for Alzheimer's disease. Alzheimers Dement. 7, 280–292. doi: 10.1016/j.jalz.2011.03.003
Sperling, R. A., Laviolette, P. S., O'keefe, K., O'brien, J., Rentz, D. M., Pihlajamaki, M., et al. (2009). Amyloid deposition is associated with impaired default network function in older persons without dementia. Neuron 63, 178–188. doi: 10.1016/j.neuron.2009.07.003
Villeneuve, S., Rabinovici, G. D., Cohn-Sheehy, B. I., Madison, C., Ayakta, N., Ghosh, P. M., et al. (2015). Existing pittsburgh compound-B positron emission tomography thresholds are too high: statistical and pathological evaluation. Brain 138, 2020–2033. doi: 10.1093/brain/awv112
Vincent, J. L., Kahn, I., Snyder, A. Z., Raichle, M. E., and Buckner, R. L. (2008). Evidence for a frontoparietal control system revealed by intrinsic functional connectivity. J. Neurophysiol. 100, 3328–3342. doi: 10.1152/jn.90355.2008
Xia, C. F., Arteaga, J., Chen, G., Gangadharmath, U., Gomez, L. F., Kasi, D., et al. (2013). [(18)F]T807, a novel tau positron emission tomography imaging agent for Alzheimer's disease. Alzheimers Dement. 9, 666–676. doi: 10.1016/j.jalz.2012.11.008
Yeo, B. T., Krienen, F. M., Sepulcre, J., Sabuncu, M. R., Lashkari, D., Hollinshead, M., et al. (2011). The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J. Neurophysiol. 106, 1125–1165. doi: 10.1152/jn.00338.2011
Keywords: Alzheimer's disease, functional brain networks, tau, positron emission tomography
Citation: Hansson O, Grothe MJ, Strandberg TO, Ohlsson T, Hägerström D, Jögi J, Smith R and Schöll M (2017) Tau Pathology Distribution in Alzheimer's disease Corresponds Differentially to Cognition-Relevant Functional Brain Networks. Front. Neurosci. 11:167. doi: 10.3389/fnins.2017.00167
Received: 24 December 2016; Accepted: 13 March 2017;
Published: 31 March 2017.
Edited by:Betty M. Tijms, VU University Medical Center, Netherlands
Reviewed by:Kim A. Staats, University of Southern California, USA
Efthimios M. C. Skoulakis, Biomedical Sciences Research Centre “Alexander Fleming”, Greece
Copyright © 2017 Hansson, Grothe, Strandberg, Ohlsson, Hägerström, Jögi, Smith and Schöll. 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) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Michael Schöll, email@example.com
†These authors have contributed equally to this work.
‡Data used in preparation of this article were obtained from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database (http://adni.loni.usc.edu/). As such, the investigators within the ADNI contributed to the design and implementation of ADNI and/or provided data but did not participate in analysis or writing of this report. A complete listing of ADNI investigators can be found at: http://adni.loni.usc.edu/wp-content/uploads/how_to_apply/ADNI_Acknowledgement_List.pdf.