Introduction

The aggregation of abnormally hyperphosphorylated tau protein into paired helical filaments is a key aspect of the pathology of Alzheimer’s disease (AD).1 Both the regional distribution of tau pathology in the brains of patients with AD and the sequential staging of its progression have been extensively described in post-mortem studies.2,3,4,5 These studies indicated, for the first time, that an early and relatively long preclinical phase of tau aggregation precedes the symptomatic stages of AD.6,7 Despite this, the time course of tau pathology propagation, especially in relation to changes in the concomitant clinical and cognitive profiles of the individual patients, remains largely speculative because of the inherent limitations of post-mortem studies.

During the past 5 years, the development of tau-specific positron emission tomography (PET) tracers8 has provided a valuable addition to the neuroimaging arsenal. THK5317 [(S)-THK5117], a well characterised tau-specific tracer,9,10,11 showed high retention in patients with AD with a regional pattern matching that of the distribution of tau pathology described by post-mortem studies.12 Cross-sectionally, high load of tau pathology, as measured with THK5317 PET, was associated with hypometabolism in restricted brain regions,12 while more extensive associations were reported in patients with substantially greater cognitive impairment, as measured with another tau tracer.13,14,15 Studies with a longitudinal, multimodal design will shed light on the spreading of tau pathology in AD, and allow investigating whether the temporal trajectories of tau aggregation and hypometabolism are closely associated, or whether this association becomes closer with disease progression.

The aggregation of tau into filaments is not exclusively restricted to AD pathology but is also found in the less frequent amyloid-β-free degenerative tauopathies. Case series have revealed high regional retention of the tau tracers in patients with clinical diagnoses of corticobasal syndrome (CBS), but the regional distribution in these cases was different from that in patients with AD.12,16,17,18 The understanding of the in vivo propagation of the pathology in the different degenerative tauopathies, however, remains elusive.

The aim of this longitudinal, multimodal study was to assess the propagation of tau pathology (using [18F]THK5317 PET) in a cohort of patients with AD and cases with CBS, and to examine the relationships between tau pathology and markers of glucose metabolism and cognitive performance over time.

Materials and methods

Study sample

Eighteen patients who had previously participated in baseline investigations volunteered to participate in a longitudinal, multimodal follow-up study after a median of 17 months (interquartile range=(15:18)) from baseline. All patients had been referred for memory assessment to the Memory Clinic at the Department of Geriatric Medicine, Karolinska University Hospital, Stockholm, Sweden. The procedures for clinical assessment and patient recruitment are detailed elsewhere.12

At baseline, six patients were diagnosed with AD dementia (that is, probable AD and a positive [11C]PIB PET scan19,20), and 10 with prodromal AD (that is, amnestic multi-domain mild cognitive impairment and a positive [11C]PIB PET scan20,21). Two patients with CBS (fulfilled clinical criteria for possible corticobasal degeneration22) were also recruited; one of the two presented with unclear clinical features and was diagnosed with mild cognitive impairment at baseline. At follow-up, all 18 patients underwent [18F]THK5317 and [18F]FDG PET imaging as well as thorough clinical and neuropsychological investigations, within a 2-month period. During the follow-up interval, four of the patients with prodromal AD developed AD dementia.

The study was approved by the regional human ethics committee in Stockholm, as well as by the radiation safety committee of Uppsala University Hospital, Sweden. All participants and their caregivers provided written informed consent prior to the investigation and all procedures performed were in accordance with the ethical standards of the institutional and national research committee and with the 1964 Helsinki Declaration and its later amendments, or comparable ethical standards.

Neuropsychological assessment

All patients underwent neuropsychological assessments at baseline and follow-up. This included assessment of global cognition using the Mini-Mental State Examination (MMSE), and of episodic memory using, among others, the Rey Auditory-Verbal Learning (RAVL) encoding subtest (encoding scores summed across five trials). The latter was selected for further analyses because performance on this subtest was not affected by floor effects in individuals with severe cognitive deficits. Individual performance results for the RAVL encoding subtest are expressed as z-scores, in comparison with results from a reference group of healthy controls.23

Image acquisition

At baseline, all participants underwent [18F]THK5317, [11C]PIB and [18F]FDG PET imaging as well as a T1-MRI sequence. At follow-up, all participants underwent [18F]THK5317 and [18F]FDG PET imaging. For [18F]THK5317 PET, 22 frames were acquired over 60 min starting simultaneously with intravenous injection of 217±42 MBq. For [11C]PIB PET, 24 frames were acquired over 60 min after intravenous injection of 256±67 MBq. The [18F]FDG PET scans were acquired with a 15-minute static acquisition, 30 min after injection of 3 MBq/kg.

Regions of interest

For regional quantification, we used regions of interest (ROIs) derived from a probabilistic atlas,24 spatially warped in each patient’s native T1-MRI space, after application of an individual grey matter mask, as previously described.12 Our choice of ROIs was based on previously published cross-sectional data on tau PET imaging—12,25 parahippocampal gyrus, fusiform gyrus, middle and inferior temporal gyrus, posterior cingulate cortex and orbitofrontal cortex for patients with AD, and putamen for patients with CBS. A composite isocortical ROI was created for quantifying [11C]PIB retention.

Data analysis

Individual dynamic baseline and follow-up [18F]THK5317 PET images were co-registered separately onto the individual T1-MRI image using PMOD v.3.5 software (PMOD Technologies, Adliswil, Switzerland). Region-based kinetic modelling with the reference Logan graphical method was applied to extract the regional distribution volume ratio (DVR) values, with cerebellar grey matter as a reference, as previously described.10 Voxel-based modelling was also applied, to create parametric DVR images.

Summation images from 40- to 60-minute [11C]PIB PET scans were created and co-registered onto the individual T1-MRI images using SPM8 software, as were [18F]FDG PET images. Standard uptake value ratio (SUVR) images were created using the cerebellar grey matter as reference for [11C]PIB and the pons for [18F]FDG.

The annual rates of change of [18F]THK5317 DVR retention (Δ DVR/year) and [18F]FDG SUVR uptake (Δ SUVR/year) were calculated for every patient, using both ROI- and voxel-based methods, as follows: ((follow-up−baseline)/time interval between scans).

Statistical analysis

Comparisons between baseline and follow-up

Differences between baseline and follow-up assessments of cognitive performance were assessed using paired t-tests (P<0.05). Changes in [18F]THK5317 retention and [18F]FDG uptake over time were assessed using voxel-based paired t-tests as implemented in SPM8 software, with the time interval used as covariate. A cluster threshold of 20 voxels was applied, with no correction for multiple comparisons (P<0.001). The determination of the sample size was based on an earlier published, exploratory study with a similar tau PET tracer.26

To account for the inter-individual variability in spreading of [18F]THK5317 retention in our sample, tau spreading indices were calculated for every patient with AD (n=16; Supplementary Figure 1). In brief, binary tau pathology maps were created for every patient at baseline and follow-up, based on the 90% confidence interval of tau distribution in healthy controls.12 The subtraction of the baseline map from the follow-up map for each patient generated individual ‘tau spreading indices’. Positive indices would illustrate that tau pathology in the cortex was overall spreading/expanding in an individual patient, and negative indices would illustrate an overall decrease of tau pathology distribution. A single sample t-test was applied to assess whether there was a significant spreading/reduction of tau pathology at a group level (that is, whether the tau spreading indices were significantly different from 0).

To find the ROIs in which [18F]THK5317 retention increased over time, the proportion of patients in the AD groups (prodromal and dementia) in whom the annual rate of change was greater than the third quartile of the test–retest repeatability range for [18F]THK531712 was determined for each region of the probabilistic atlas, and an ROI-based map was created with the results. The same analysis was not applied to the [18F]FDG data due to the absence of [18F]FDG test–retest measures in our study.

Comparisons between modalities

The relationships between the regional retention of [18F]THK5317 at baseline and the rate of change in [18F]THK5317 retention and [18F]FDG uptake were assessed using the Pearson coefficient in the selected ROIs.

Linear mixed-effects models allowing for random intercepts between patients were employed to explore the relationship over time between regional [18F]THK5317 retention and [18F]FDG uptake as well as cognitive measures, with the time interval used as covariate. An additional linear mixed-effects model was applied to assess the relationship over time between [18F]THK5317 retention and [18F]FDG uptake across the four major brain lobes (frontal, temporal, parietal, occipital), similarly to the method from Ossenkoppele et al. 13 More specifically, [18F]FDG uptake was assigned as the dependent variable, with [18F]THK5317 retention, time point (baseline or follow-up) and their interaction as fixed factors, allowing for random intercepts across different patients and brain lobes (nested design). Simple slopes analyses were used to assess multiple linear regression two-way interactions.27

The cutoff point for statistical significance was P<0.05. For statistical analyses involving multiple ROIs, both uncorrected and Bonferroni-adjusted (five ROIs) P-values are reported. For analyses involving cognitive tests, additional adjustments were performed (five ROIs and two cognitive tests). The data met the appropriate assumptions for the different statistical comparisons. All analyses were carried out using R v.3.3.1 software.

Results

Clinical and neuropsychological profiles at follow-up in patients with AD

The characteristics of the study population are summarised in Table 1. MMSE scores were significantly lower (P=0.007) at follow-up than at baseline in all patients with AD (prodromal and dementia); there were no statistically significant changes in RAVL encoding scores.

Table 1 Demographic and clinical characteristics of the study population; follow-up investigations were performed after a median of 17 months (interquartile range=(15:18))

[18F]THK5317 retention and [18F]FDG uptake at follow-up in patients with AD

At baseline, across all patients with AD, high [18F]THK5317 retention was predominantly observed in temporal, lateral occipital, and lateral and medial frontal areas (Figure 1a). The extent of retention remained visually the same at follow-up (Figure 1b), and voxel-based paired t-tests also detected no statistically significant increase in [18F]THK5317 retention at follow-up, except for a focal area in the left inferior temporal gyrus (Figure 1c). Tau spreading index analyses revealed no spreading of [18F]THK5317 retention over time at a group level. In contrast, low [18F]FDG uptake was observed at baseline in temporoparietal areas, as seen in the average maps (Figure 1d), and this had clearly decreased at follow-up (Figure 1e). Indeed, statistical comparisons showed a significant decrease in [18F]FDG uptake at follow-up, mainly in temporoparietal areas (Figure 1f).

Figure 1
figure 1

Average baseline and follow-up images from patients with AD (prodromal AD and AD dementia, n=16), and projections of the results of voxel-based paired t-tests (SPM8 software) exploring the changes in retention/uptake over time for [18F]THK5317 (ac) and [18F]FDG (df). A cluster threshold of 20 voxels was applied, with no correction for multiple comparisons (P<0.001). AD, Alzheimer’s disease; DVR, distribution volume ratio; SUVR, standard uptake value ratio.

To further explore possible changes in [18F]THK5317 retention in the two AD groups separately (based on the baseline diagnosis), maps of the ROIs were created to illustrate the proportions of patients with increases in [18F]THK5317 retention (annual rates) that were greater than those in the test–retest repeatability range. The areas in which the greatest proportion of patients showed substantial increase were the medial orbitofrontal cortex in patients with prodromal AD (Figure 2a) and the lateral occipital cortex in patients with AD dementia (Figure 2b). Overall, a substantial increase in [18F]THK5317 retention in neocortical areas was observed more frequently in the AD dementia group than in the prodromal AD group (Figure 2).

Figure 2
figure 2

Regional maps illustrating the percentage (%) of patients with significant increase of annual [18F]THK5317 retention in patients with (a) prodromal AD and (b) AD dementia. Significant increase was defined as annual rates of increase that were greater than the third quartile of the test–retest repeatability range.12 AD, Alzheimer’s disease.

Illustrative cases of patients with AD

Figure 3a illustrates the [18F]THK5317 and [18F]FDG PET scans at baseline and the maps of respective annual rates of change from two patients with a baseline diagnosis of prodromal AD and two with baseline diagnosis of AD dementia. With regard to the prodromal AD patients, case 1, who had milder cognitive impairment at baseline, had a stable diagnosis over time, while case 2 had developed AD dementia by follow-up. In the imaging profiles, case 1 had less [18F]THK5317 retention at baseline but widespread increases in retention over time, while case 2 had only focal increases. With regard to [18F]FDG PET, case 1 did not exhibit a clear hypometabolic pattern at baseline, and did not show substantial declines in tracer uptake over time. In contrast, case 2 had very low [18F]FDG uptake in temporoparietal areas at baseline and further decreases in uptake at follow-up, mainly in temporal regions (Figure 3a).

Figure 3
figure 3

Imaging profiles from illustrative patients with baseline diagnoses of prodromal AD, AD dementia and CBS. [18F]THK5317 (DVR) and [18F]FDG (SUVR) PET scans at baseline and the maps of their respective annual rates of change are presented for each patient (a,b). (c) [18F]THK5317 retention in the putamen in the two patients with CBS over time, plotted against the baseline retention values (interquartile range) from patients with AD (n=16, brown ribbon) and those from young, healthy volunteers (n=5, grey ribbon).12 AD, Alzheimer’s disease; Apo E, apolipoprotein E; CBS, corticobasal syndrome; DVR, distribution volume ratio; FDG, fluoro-deoxyglucose; MMSE, Mini-Mental State Examination; PET, positron emission tomography; SUVR, standard uptake value ratio.

With regard to the AD dementia patients, case 3 had milder cognitive impairment at baseline (as measured with MMSE) that remained stable over the follow-up interval, while the cognitive impairment in case 4 was worse at baseline and showed a clear further decline at follow-up (loss of seven points on MMSE). [18F]THK5317 retention was less in case 3 than in case 4 at baseline but increased more over time compared with case 4. [18F]FDG PET in case 3 showed relatively restricted areas with low tracer uptake than in case 4. [18F]FDG uptake decreased widely in both cases with time, although this was more pronounced in case 4 (Figure 3a).

Relationships between [18F]THK5317, [18F]FDG and cognitive performance in AD

Significant negative correlations were observed between baseline [18F]THK5317 retention and the rate of change of [18F]THK5317 retention locally in the parahippocampal, fusiform and orbitofrontal cortex (Figure 4a). No significant correlations were observed between baseline local [18F]THK5317 retention and the rate of change of [18F]FDG in the selected ROIs.

Figure 4
figure 4

Relationship between [18F]THK5317 retention, [18F]FDG uptake and cognitive performance over time. More specifically, annual regional [18F]THK5317 retention rate of change in relation to (a) local [18F]THK5317 retention at baseline in the parahippocampal gyrus, fusiform gyrus and orbitofrontal cortex (Pearson’s correlation analyses). Performance on (b) the MMSE and (c) the RAVL encoding tests in relation to the [18F]THK5317 retention and [18F]FDG uptake in the middle and inferior temporal gyrus over time (linear mixed-effects model analyses). The RAVL encoding performance results are expressed as z-scores. Linear mixed-effects model analyses were employed. AD, Alzheimer’s disease; DVR, distribution volume ratio; e, estimate; FDG, fluoro-deoxyglucose; MMSE, Mini-Mental State Examination; P adj, Bonferroni-adjusted P-value; r, Pearson coefficient; R 2, coefficient of determination; RAVL, Rey Auditory-Verbal Learning; SUVR, standard uptake value ratio. *P<0.05.

There were no significant relationships over time between [18F]THK5317 retention and [18F]FDG uptake locally in the selected ROIs. Across the four major lobes, there was a significant effect of the interaction between [18F]THK5317 retention and the time point (baseline or follow-up), on the [18F]FDG uptake (F=5.505, P=0.021). Accordingly, the relationship between the two tracers showed significantly different slopes between baseline and follow-up (slope=−0.202 and −0.375, respectively) (Supplementary Figure 2).

Increased [18F]THK5317 retention in the middle and inferior temporal gyrus significantly predicted decreased RAVL encoding scores but not MMSE scores (Figures 4b and c). No relationships were found between [18F]THK5317 retention and cognitive performance (measured with MMSE or RAVL encoding tests) in the other ROIs. Decreased [18F]FDG uptake in the middle and inferior temporal gyrus significantly predicted decreased MMSE and RAVL encoding scores (Figures 4b and c), while decreased [18F]FDG uptake in the posterior cingulate gyrus predicted decreased MMSE scores (Figure 3); however, the statistical significance of the RAVL encoding score predictions did not survive correction for multiple comparisons for any of the tracers in the selected ROIs.

[18F]THK5317 retention and [18F]FDG uptake at follow-up in patients with CBS

Both patients with a clinical diagnosis of CBS had a negative [11C]PIB PET scan and high [18F]THK5317 retention at baseline, with a regional distribution pattern that was distinct from the pattern of [18F]THK5317 retention in patients with AD. High [18F]THK5317 retention was detected predominantly in the basal ganglia, the thalami and the white matter, but also in isocortical areas, in these patients. Asymmetrical [18F]FDG uptake was observed at baseline in isocortical areas, the basal ganglia and the thalami in both patients; however, the pattern of [18F]FDG uptake was more asymmetrical in case 6 than in case 5 (Figure 3b).

At follow-up, measures of global cognition (MMSE) had substantially declined in both patients, in comparison to baseline (Table 1). Both also showed higher [18F]THK5317 retention in the basal ganglia and isocortical areas than in the baseline investigations (Figures 3b and c). Case 6, who had higher MMSE scores and less white-matter retention at baseline, had more extensive retention in the same area at follow-up than case 5. Both patients had widespread areas of lower [18F]FDG uptake in the isocortex at follow-up; case 6 showed the greatest changes in [18F]FDG uptake over time (Figure 3b).

Discussion

Thus far, our knowledge of the temporal evolution of tau pathology has been largely based upon either post-mortem studies4,5 or theoretical models that have endeavoured to recapitulate the propagation of the pathology through cross-sectional tau PET imaging in mixed populations of patients and cognitively normal individuals.28,29,30 In this study, we used [18F]THK5317 PET imaging in a longitudinal, multimodal design to investigate in vivo the propagation of tau pathology in a cohort of patients with AD and cases with CBS, in relation to markers of glucose metabolism and cognitive performance.

The patients with AD who showed abnormally high [18F]THK5317 retention in areas exceeding the limbic lobe at baseline12 showed no significant increase in [18F]THK5317 retention in the 17-month follow-up period, at a group level. A comparison of patients with prodromal AD and AD dementia, however, revealed different patterns of increases in [18F]THK5317 retention. The medial orbitofrontal cortex, the area where greatest proportion of the prodromal AD patients showed increased [18F]THK5317 retention, although outside the limbic lobe, shows extensive tau pathology in AD,31 similar to the lateral temporal cortex,32 and has been suggested as one of the isocortical regions early affected5 due to its strong connections with the limbic system.33,34,35 In contrast, the greatest proportion of the patients with AD dementia showed increased [18F]THK5317 retention in the association area of the occipital cortex, an area thought to be affected later in the disease progression.5 Substantial changes in retention were not observed in the primary occipital cortex, an area affected at the last stage of tau spreading,36 much later than its surrounding neocortex.37,38 Finally, wide heterogeneity in the load and regional rates of change of [18F]THK5317 retention was detected between individual AD patients even within each clinical group. It is therefore likely that patients clinically classified in the same group according to their functional impairment—namely, the prodromal and dementia stages of AD—were actually at different neuropathological stages and exhibited different patterns of tau propagation.

The extent of baseline [18F]THK5317 retention in the parahippocampal and fusiform gyri and the orbitofrontal cortex, areas affected by tau pathology early in the time course of the disease, showed a negative correlation with the annual rates of change of [18F]THK5317 retention. Despite the exploratory nature of this analysis, it is possible that the tau load in these areas could reach a critical value, after which the rate of accumulation would decelerate or even plateau, similarly to the changes in amyloid-β pathology early in the disease course.39

In contrast with the lack of changes in [18F]THK5317 retention in patients with AD at a group level over the 17-month period, there were significant decreases in [18F]FDG uptake in widespread areas of the temporoparietal cortex—a signature area for AD.40 Interestingly, our data indicated the lack of a temporal association between changes in tau pathology and glucose metabolism, especially in areas affected early by tau. Based on the imaging profiles of the individual patients, it is conceivable that there may be a lag phase between the build-up of tau pathology and changes in glucose metabolism; during this phase, the neurons could tolerate a substantial amount of tau pathology,41,42 possibly through endogenous compensatory mechanisms, before exhibiting metabolic changes. The latter would be consistent with cross-sectional findings indicating that the associations between markers of tau and metabolism become extensive in the dementia stage of the disease.12,13,15,43 The existence of a lag phase is further supported by our analyses across the major brain lobes, where the relationship between tau pathology and hypometabolism proved closer with disease progression, when hypometabolism becomes more prevalent. The use of non-linear models in future studies with larger samples might be better suited to investigate the association between the two markers with time. Furthermore, longitudinal studies with adequate power, to study multi-level interactions, are warranted to investigate further the relationship between baseline tau pathology and the rate of change of glucose metabolism, as well as the possible synergistic effect of tau and amyloid-beta burden, as previously suggested.14,44

Global cognition (measured with MMSE) declined significantly in patients with AD; this decline correlated with decreased [18F]FDG uptake but not with changes in [18F]THK5317 retention over time. Changes in episodic memory correlated with changes in both longitudinal [18F]FDG uptake (positively) and longitudinal [18F]THK5317 retention (negatively); the correlations with episodic memory, however, did not survive correction for multiple comparisons and can be considered preliminary due to the small number of patients completing the episodic memory testing. The current data are in line with earlier cross-sectional findings45 and support the notion that while tau deposition in the temporal cortex may be related to the earliest memory impairment, global cognitive changes are more closely related with hypometabolism in the relevant areas than with measures of tau propagation. Taken together, these findings hint that, in the symptomatic stages of AD, [18F]FDG uptake and [18F]THK5317 retention are not tightly associated, and that [18F]FDG PET might be a better tool for monitoring symptom progression.

At baseline, the patients with a clinical diagnosis of CBS showed a pattern of high [18F]THK5317 retention and low [18F]FDG uptake that was distinct from that in patients with AD, as previously described.12,16,17,18 [18F]THK5317 retention increased in both these patients over time, predominantly in the basal ganglia and fronto-temporal areas, and [18F]FDG uptake decreased isocortically. Nonetheless, there were differences between the two cases that might indicate that they were at different stages in the progression of the pathology. Notwithstanding the limitations of case studies, our findings illustrate the different regional distribution and progression patterns of tau pathology in patients with CBS, in contrast to patients with AD. More research in large groups of patients with CBS is required to determine: (1) whether the high binding of [18F]THK5317 in the basal ganglia, which is dominant in both CBS and AD, would allow discrimination between the two, or if the assessment of AD-vulnerable regions would serve better this purpose; and (2) whether tau pathology monitors better disease progression in CBS in comparison to AD, as preliminarily observed in the profiles of our individual patients.

The results of this work are in contrast to the, so far, only longitudinal tau PET study published that reported changes in tracer retention over time.26 We believe that this discrepancy could be attributed to the much smaller sample size of the latter study—five patients with AD dementia—as well as to the methodology used. In our cohort, the greater sample size increases the power of our calculations, and therefore increases the probability (positive predictive value) that our significant observed effects would reflect true effects, and not type I errors. Moreover, the use of dynamic data in a kinetic model-based approach to quantify [18F]THK5317 retention (DVR) is another advantage of the current work, as DVR provide more accurate quantification of the tracer retention over time than late frame summation images (SUVR).10

The strength of this work lies in the longitudinal, multimodal design, which allowed the comparison between different clinical and imaging markers of the underlying pathology. The small sample size, however, is a clear limitation of the study, since it did not allow for separate inferential statistical analyses for the prodromal and AD dementia groups. Such analyses would probably provide additional information on the inter-relationships between the markers. The processing of the follow-up data based on the baseline MRI, similarly to previous work,46 represents another limitation, although the reported low rates of atrophy in isocortical areas at similar time intervals47 would not be expected to have more than minimal effects given the spatial resolution of PET.

The interpretation of the observed relative stability in [18F]THK5317 retention, at a group level, over time, however, is limited by our understanding of the tracer’s binding. More specifically, the stability in retention could not exclude changes in the maturation status of the underlying tau deposits-preliminary evidence from in vitro studies suggests that the tracer binds to different types of tau deposits.48

In summary, this study suggests that patients in the symptomatic stages of AD exhibit heterogeneous regional load and propagation patterns of tau pathology over a 17-month period, despite being classified into similar groups of functional impairment. In direct contrast, glucose metabolism showed a more homogeneous pattern of decrease, which better tracked clinical progression. Our findings provide evidence of the lack of a linear relationship between tau pathology and more downstream markers, namely hypometabolism and measures of global cognitive impairment, during the course of the disease. Considerably more work is warranted to investigate tau propagation in larger sample sizes, at different stages of disease and using longer follow-up intervals, with special care around the possible uncoupling of clinical and neuropathological staging of AD.