Abstract
Although several studies have shown decreased cerebral blood flow (CBF) in Alzheimer’s disease (AD), the role of hypoperfusion in the disease pathogenesis remains unclear. Combining arterial spin labeling MRI, positron emission tomography, and biomarkers of cerebrospinal fluid, we investigated the associations between CBF and the key mechanisms in AD including amyloid-β (Aβ) and tau pathology, synaptic dysfunction and axonal degeneration. Further, we applied a disease progression modeling to characterize the temporal sequence of different AD biomarkers. Lower perfusion was observed in the temporo-occipito-parietal regions in the Aβ-positive cognitively impaired compared to both the Aβ-positive and Aβ-negative cognitively unimpaired individuals. In participants along the AD spectrum (those with Aβ pathology regardless of their cognitive status), CBF was inversely associated with tau and synaptic dysfunction, but not Aβ, in similar cortical regions. Moreover, the disease progression modeling revealed that CBF disruption followed the abnormality of biomarkers of Aβ, tau and brain atrophy. These findings indicate that tau tangles and synaptic degeneration are more closely connected with CBF changes rather than Aβ pathology. This supports the notion that hypoperfusion is not an early event associated with the build-up of Aβ during the preclinical phase of AD.
Introduction
Alzheimer’s disease (AD) is a neurodegenerative disorder characterized by the accumulation of amyloid-β (Aβ) plaques and neurofibrillary tangles of hyperphosphorylated tau. These pathological changes precede the onset of dementia by many years, leading to the conceptualization of AD as a continuum evolving from preclinical to prodromal stages and finally to dementia (Aisen et al., 2017; McKhann et al., 2011). Multiple neuroimaging studies have shown alterations in functional connectivity (Berron et al., 2020; Franzmeier et al., 2019; Hahn et al., 2019), brain metabolism (Mosconi et al., 2008; Riederer et al., 2018; Verclytte et al., 2016) and perfusion (Huang et al., 2018; Mattsson et al., 2014) in different cortical regions through the course of AD. The latter can be noninvasively measured with arterial spin labeling (ASL) which utilizes magnetically labeled water in the blood as an endogenous tracer (Williams et al., 1992). Due to the close coupling of glucose metabolism and cerebral blood flow (CBF), the patterns of regional hypometabolism on well-established positron emission tomography with [18F]-fluorodeoxyglucose (FDG-PET) overlap with reduced CBF on ASL scans (Riederer et al., 2018; Musiek et al., 2012). This renders ASL as a promising radiation-free candidate for the diagnosis and monitoring of AD.
Despite several reports on CBF changes in AD (Binnewijzend et al., 2016; Wolk and Detre, 2012; Fazlollahi et al., 2020), their association with the core pathologies including Aβ and tau burden remains elusive. In a previous study, a negative correlation was reported between amyloid-PET and ASL-CBF in temporo-parietal regions across the AD continuum (Mattsson et al., 2014). In contrast, other studies have suggested a link between cerebrospinal fluid (CSF) biomarkers of tau, but not Aβ, with perfusion abnormalities on single photon emission computed tomography (Habert et al., 2010; Stomrud et al., 2012). The few available studies on the relationship between ASL-CBF and tau-PET have yielded mixed findings. While one study has found negative associations in widespread temporo-parietal regions (Albrecht et al., 2020), another one has shown that higher tau-PET up-take is correlated with hypoperfusion exclusively in the entorhinal cortex (Rubinski et al., 2021). Moreover, no studies to date have investigated the associations of ASL-CBF with synaptic and axonal markers in the CSF. There is also the important question whether decreased CBF triggers the chain of pathological events in AD, or if it is a consequence of the downstream effects of Aβ and tau. Thus, the overarching aims of the current study were to examine the multimodal relationships between ASL-CBF with distinct imaging and CSF biomarkers across the AD spectrum and to identify when along the disease evolution ASL-CBF changes occur.
To this end, we first measured ASL-CBF variations in Aβ-positive or negative cognitively unim-paired individuals and cognitively impaired patients with Aβ pathology. We then evaluated the association between CBF and deposition of Aβ and tau measured with Aβ-PET and CSF Aβ42/40, tau-PET and CSF phosopho-tau217 (P-tau217), as well as the CSF concentrations of synaptic (Neuronal Pentraxin2 to total tau ratio; NPTX2/T-tau), and axonal (neurofilament light chain; NfL) biomarkers. Finally, we ran a disease progression modeling (Young et al., 2018) approach to estimate the order of abnormality of these markers. Specifically, we addressed the following questions. i) When does CBF change in the AD continuum? ii) Which brain regions are primarily affected in AD by CBF alterations? iii) Are CBF changes associated with Aβ or tau burden in the AD spectrum? iv) Are CBF variations corelated with markers of synaptic or axonal integrity? Our results demonstrate that CBF alterations occur at later stages of the AD continuum. These changes are most pronounced in the occipito-temporo-parietal areas and are largely associated with tau load and synaptic degeneration.
Results
Participants and methodology overview
The present study comprised 137 cognitively unimpaired (CU) controls, and 119 cognitively impaired (CI) patients diagnosed either with mild cognitive impairment or AD with dementia. The groups were further stratified into Aβ-positive and negative individuals based on the CSF Aβ42/40 ratio. Demographic and clinical characteristics of the participants are summarized in Table 1.
Voxel-wise and region of interest (ROI)-based general linear models were performed to assess between-group differences in CBF and its relationship with the PET and CSF measures. The ROI-based analyses were conducted in six a priori defined ROIs including lateral and medial parts of the parietal and temporal lobes, middle frontal gyrus and superior lateral occipital cortex (Figure 1). All the analyses were adjusted for age, sex and ASL sequence version (see Methods; MRI acquisition and processing). To cover brain regions affected by tau and Aβ pathology, composite-ROIs were created. Tau-PET ROIs (illustrated in Supplementary Figure 1) approximated Braak staging scheme (Braak and Braak, 1991) and encompassed an early (entorhinal) ROI corresponding to Braak stage I-II, a temporal meta-ROI reflecting Braak stage I-IV, and a late neocortical meta-ROI corresponding to Braak stage V-VI (Ossenkoppele et al., 2018). The Aβ-PET neocortical composite ROI covered prefrontal, parietal, temporal lateral, anterior/posterior cingulate, and precuneus (Lundqvist et al., 2013; Landau et al., 2014).
To estimate the timing of CBF disruption along the disease pathogenesis, we applied Subtype and Stage Inference (SuStaIn) approach, which combines disease progression modeling with traditional clustering to achieve probabilistic spatiotemporal classification (Young et al., 2018). CSF levels of Aβ42/40 and P-tau217, tau-PET deposition in the aforementioned three ROIs, a composite volumetric measure in AD signature regions (Jack Jr et al., 2017; Ossenkoppele et al., 2018) and weighted mean CBF in Aβ-negative CU individuals and those on the AD continuum i.e., all Aβ-positive individuals were used as inputs for this modeling framework (see Methods; disease progression modeling).
CBF differences between diagnostic groups
Aβ-negative CU vs. Aβ-positive CU
No significant difference in CBF was observed when comparing Aβ-negative and Aβ-positive CU individuals, neither with the ROI-based nor with the voxel-wise analyses.
Aβ-negative CU vs. Aβ-positive CI
The ROI-based analyses revealed significant CBF reductions in the lateral parietal (β = -0.059, PFDR = 0.0052), superior lateral occipital (β = -0.056, PFDR = 0.0443) and middle frontal (β = -0.074, PFDR = 0.0021) ROIs in Aβ-positive CI individuals compared to Aβ-negative CU participants. These findings were further confirmed by the voxel-wise analyses (Figure 2 A), demonstrating decreased CBF in the bilateral superior lateral occipital, bilateral lateral/medial parietal and the left lateral temporal regions (p < 0.05, FWE-corrected).
Aβ-positive CU vs. Aβ-positive CI
The overall hypoperfusion pattern in Aβ-positive CI with respect to Aβ-positive CU individuals was similar to the pattern observed in the same group when compared to Aβ-negative CU subjects, albeit less pronounced. Specifically, lower CBF was observed in the lateral parietal (β = -0.047, PFDR = 0.0426), superior lateral occipital (β = -0.070, PFDR = 0.0426) and middle frontal (β = -0.050, PFDR = 0.0426) ROIs in Aβ-positive CI individuals. As depicted in Figure 2 B, the voxel-wise analyses showed decreased CBF in the lateral temporal, medial/lateral parietal and superior lateral occipital cortices, largely in line with the ROI-based analyses (p < 0.05, FWE-corrected).
Additional correction for the grey matter (GM) volume did not alter the observed findings between the two CU groups. However, the patterns of decreased CBF in Aβ-positive CI patients compared to both the Aβ-negative and Aβ-positive CU participants were diminished after controlling for the GM volume (see Supplementary Table 1 and Supplementary Figure 2).
Correlations between CBF and Aβ pathology
Next, we studied the associations between CBF and Aβ pathology measured by Aβ-PET SUVR in the neocortical composite ROI and CSF levels of Aβ42/40. In CU individuals, no associations were observed between CBF and markers of Aβ, neither with the ROI-based nor with the voxel-wise analyses. Similarly, CBF was not correlated with Aβ pathology in individuals along the AD continuum.
Correlations between CBF and tau pathology
To study associations between CBF and tau pathology, we initially used the tau-PET data. No significant associations were found between CBF and tau load in any of the early, temporal and late tau-PET ROIs in CU individuals regardless of Aβ status. Additionally, CBF was not associated with CSF levels of P-tau217 in this group. In individuals on the AD spectrum, however, both the ROI-based and voxel-wise analyses demonstrated inverse associations between CBF and tau biomarkers. Higher tau-PET uptake in the early ROI was significantly associated with decreased CBF in the lateral parietal (β = -0.052, PFDR = 0.0134), lateral temporal (β = -0.054, PFDR = 0.0134), and superior lateral occipital (β = -0.081, PFDR = 0.0068) ROIs. Tau burden in the temporal meta-ROI was negatively associated with CBF in the lateral parietal (β = -0.060, PFDR = 0.0001), lateral temporal (β = -0.064, PFDR = 0.0001), the superior lateral occipital (β = -0.088, PFDR = 0.00005) and middle frontal (β = -0.034, PFDR = 0.0499) ROIs. Moreover, tau load in the late meta-ROI had inverse associations with CBF in the same regions [lateral parietal (β = -0.122, PFDR = 0.00005), lateral temporal (β = -0.086, PFDR = 0.0046), superior lateral occipital (β = -0.128, PFDR = 0.0011) and middle frontal (β = -0.087, PFDR = 0.0046) ROIs]. The associations of CBF with tau-PET ROIs were still significant after controlling for CSF concentrations of Aβ42/40 (Supplementary Table 2). Analogous findings were observed after controlling for cognitive status and GM volume (see Supplementary Table- and Figure 3). The overall significant associations between CBF and tau uptake in each of the three ROIs are plotted in Supplementary Figure 4. Further using the voxel-wise analyses, similar patterns of hypoperfusion were observed for all three tau-PET ROIs mainly in the occipito-temporo-parietal regions (Figure 3). Except for the tau-PET SUVRs in the early ROI, additional adjustment for the cognitive status did not alter the observed voxel-wise associations, as shown in Supplementary Figure 5. This also held true when controlling for GM volume (Supplementary Figure 6).
Afterwards, we examined the associations between CBF and CSF P-tau217 which is an earlier marker than tau PET (Mattsson-Carlgren et al., 2020). Elevated CSF levels of P-tau217 were correlated with lower CBF in lateral parietal (β = -9.643e-05, PFDR = 0.0117), lateral temporal (β = -8.518e-05, PFDR = 0.0381), superior lateral occipital (β = -1.708e-04, PFDR = 0.0017) and middle frontal (β = -8.738e-05, PFDR = 0.0381), see Supplementary Figure 7. After adjustment for GM volume, only the associations in the lateral parietal and superior lateral occipital remained significant [β = -8.18e-05, p = 0.0185, and β = -1.65e-04, PFDR = 0.0049, respectively]. Similarly, the voxel-wise analyses showed a negative correlation between higher levels of P-tau217 and hypoperfusion in the occipito-temporo-parietal areas (Figure 4). As depicted in Supplementary Figure 8, following further adjustment for the GM volume, the observed voxel-wise associations were restricted to clusters within the lateral temporal and superior lateral occipital cortex.
Correlations between CBF and CSF biomarkers of synaptic and axonal integrity
To study the associations between CBF and markers of neurodegeneration, we used CSF levels of NfL, which increases in relation to ongoing axonal degeneration (Gaetani et al., 2019), and NPTX2/T-tau ratio which decreases with respect to synaptic degeneration (Galasko et al., 2019). None of the ROI-based and voxel-wise analyses revealed any significant correlation between CBF and CSF levels of NfL or NPTX2/T-tau in CU individuals. In contrast, a positive association was observed between NPTX2/T-tau and CBF in the AD spectrum for lateral parietal (β = 0.062, PFDR = 0.0069), lateral temporal (β = 0.064, PFDR = 0.0072), superior lateral occipital (β = 0.096, PFDR = 0.0028) and middle frontal (β = 0.074, PFDR = 0.0028) ROIs, see Supplementary Figure 9. The observed correlations were largely corroborated by the voxel-wise analyses showing hypoperfusion with decreasing concentrations of NPTX2/T-tau (see Figure 5). Additional correction for the GM volume did not affect the observed ROI-based and voxel-wise associations (Supplementary Figure 10). Conversely, no significant associations were found between NfL and CBF in this group, neither with the ROI-based nor with the voxel-wise analyses.
Staging of CBF with respect to biomarkers of Aβ, tau and neurodegeneration
SuStaIn was applied to probabilistically assign individuals in Aβ-negative CU group and those along the AD continuum to one of 21 progressive stages (see Figure 6). In accordance with the aforementioned results, the model predicted that abnormal CSF concentrations of Aβ42/40 and P-tau217 initiated the progression of the disease. This was followed by pathological deposition of tau-PET in early and temporal/limbic ROIs, as well as neurodegeneration indicated by the atrophy in the AD signature regions and abnormal tau burden in the late meta-ROI. Only subsequent to these changes, did ASL-CBF alterations become apparent.
Discussion
The present study disentangles the multimodal associations of ASL-CBF with different pathological markers in AD. We found that i) CBF did not differ between the Aβ-negative and Aβ-positive CU groups. However in Aβ-positive CI patients, CBF was markedly reduced compared to CU individuals independent of Aβ pathology; ii) the CBF reduction were mainly localized in the occipito-temporoparietal regions; iii) CBF was not associated with markers of Aβ neither in CU nor in individuals on the AD continuum; iv) higher tau-PET uptake and increased concentrations of CSF P-tau217 were correlated with hypoperfusion patterns in the AD spectrum; v) abnormal levels of the NPTX2/T-tau ratio, a marker of synaptic degeneration, were associated with lower CBF in the same group; and vi) pathological reduction of CBF was evident only after the abnormality of Aβ, tau and regional cortical atrophy. Taken together, these findings indicate that ASL-CBF changes are not an early event associated with aggregation of Aβ plaques during the preclinical phase of AD. Instead, these alterations are more tightly linked to tau pathology, and synaptic degeneration.
The observed CBF differences between Aβ-positive CI and the CU groups are consistent with the accumulating evidence demonstrating widespread hypoperfusion across different cortical regions in AD (Dai et al., 2009; Johnson et al., 2005; Binnewijzend et al., 2013). When additional correction for GM volume was applied, the observed effect size was decreased suggesting that between group CBF differences may not be fully independent of GM atrophy. Conflicting results have been reported on CBF changes in CU individuals with elevated risk of AD. While some studies have shown patterns of hyperperfusion in CU APOE ε4 carriers (Wierenga et al., 2014; McKiernan et al., 2020), others have found no CBF changes in CU individuals with Aβ pathology (Binnewijzend et al., 2016). Our findings endorse the latter indicating that ASL-CBF is a marker of the disease severity, not directly related to Aβ deposition.
Lack of an association between CBF and Aβ burden, as suggested by our results, is in line with a previous report on non-demented humans with unilateral occlusion of precerebral arteries, showing that long-lasting cerebral hypoperfusion does not result in aggregation of Aβ (Hansson et al., 2018). Another study, however, has found positive associations between hyperperfusion and increased Aβ-PET load in a small sample (N = 13) of Aβ-positive CU individuals (Fazlollahi et al., 2020). On the other hand, inverse associations between lower perfusion and increased Aβ-PET uptake have been reported across the spectrum of sporadic AD and in asymptomatic and mildly symptomatic subjects with autosomal dominant AD (Mattsson et al., 2014; McDade et al., 2014). While the observed discrepancy suggests that the potential role of tau should be taken into account when assessing the CBF-Aβ relationship, it may also arise from the differences in the sample size or acquisition of the ASL-CBF images. In the current study, we applied the recommended implementation for the clinical use of ASL i.e., pseudo-continuous sequence with background suppression and 3D gradient- and spin-echo readout, whereas in the latter two studies, a pulsed sequence was used, which is less sensitive and more prone to being influenced by the variations in arterial transit time (Alsop et al., 2015). Our findings on the negative CBF-tau associations mainly in the occipito-temporo-parietal cortex in the AD spectrum are in general agreement with previous studies comparing tau-PET uptake with glucose metabolism using FDG-PET (Whitwell et al., 2018), and neuropathological results showing higher tau pathology with increased concentrations of vascular endothelial growth factor, an indicator of hypoxia-induced hypoperfusion (Thomas et al., 2015). A recent study with direct comparison of ASL-CBF and tau-PET has also shown CBF reduction in similar cortical regions (Albrecht et al., 2020). The underlying mechanisms linking tau pathology to CBF decline in AD are yet to be determined. Evidence from animal models has suggested that chronic lower perfusion promotes the accumulation of tau tangles (Park et al., 2019; Shimada et al., 2019). In contrast, it has been proposed that neurofibrillary tau tangles have detrimental effects on the integrity of the cortical microvasculature (Bennett et al., 2018).
In this study, the unique combination of imaging and CSF measures allowed us to explore the associations of ASL-CBF not only with Aβ and tau as the hallmarks of AD, but also with synaptic and axonal function. Importantly, we found that lower CBF was associated with lower levels of the NPTX2/T-tau ratio, indicating that synaptic degeneration is accompanied by CBF reduction. NPTX2 is a marker of glutamatergic synaptic plasticity (Swanson et al., 2016), which has been shown to be down-regulated in AD correlating with cognitive deterioration and hippocampal volume (Xiao et al., 2017). Here, we used NPTX2/T-tau ratio as it has been shown to have a high classification accuracy in discriminating AD from normal controls and a strong correlation with cognition (Galasko et al., 2019). Together, these findings suggest that CBF alterations in the AD spectrum are elicited by the spread of tau tangles and synaptic damage.
Contrary to what we expected, no correlation was observed between NfL and ASL-CBF. Although there are no published reports on such associations, a few studies have found a link between higher CSF or plasma levels of NfL with cortical hypometabolism in AD (Mayeli et al., 2019; Benedet et al., 2019). Further, another study has reported longitudinal associations between higher concentrations of NfL and decline in several neuroimaging measures including FDG-PET (Mielke et al., 2019). However, no cross-sectional correlations were found between either CSF or plasma levels of NfL with any of the neuroimaging measures in that study. This agrees with our observation on the absences of a cross-sectional relationship between ASL-CBF and NfL. Collectively, our findings suggest that variations in CBF are more related to synaptic integrity rather than the axonal damage.
Whether perfusion disruption is a cause or consequence of AD remains unclear (Austin et al., 2011; Mazza et al., 2011). Some studies have reported that CBF deficiency induces AD neuropathology (Iturria-Medina et al., 2016; Zhao and Gong, 2015; Koike et al., 2010) whereas others have shown that cerebral hypoperfusion is secondary to neurodegeneration (Hansson et al., 2018). By modeling the disease progression using SuStaIn, we found that anomalous concentrations of Aβ and tau along with cortical volume loss preceded the CBF disturbances in AD, lending support to the latter view. Nevertheless, it is still of paramount importance to treat cardiovascular risk factors in the elderly population to decrease the risk of vascular dementia and the additive effects of lacunae and white matter lesions as co-pathologies in AD.
The cross-sectional nature of this study did not allow us to determine the temporal relationship of ASL-CBF with each of the imaging and CSF measures. Future studies are, thus, warranted to elucidate the longitudinal associations between CBF and markers of Aβ, tau, synaptic and axonal integrity. Many ASL-studies have used in-house scripts for image processing, which hampers the reproducibility of the results. Here, we employed open source tools aimed at standardization of ASL image processing. However, it should be stressed that due to the limited spatial resolution of ASL, the observed findings may not represent the dynamics of blood flow within the cerebral capillaries.
In conclusion, these results provide in vivo evidence indicating that tau aggregates and neurodegeneration, but not Aβ, are closely connected to ASL-CBF alterations across the AD continuum. This suggests that cerebral hypoperfusion, at least measured by ASL, does not stimulate the development of AD pathology but rather is a late reactive change due to the reduced metabolic demand.
Methods and Materials
Participants
The participants were recruited from the Swedish BioFINDER-2 study (NCT03174938), which has been previously described in detail (Palmqvist et al., 2020; Spotorno et al., 2020). Out of 447 eligible participants, we included 256 individuals, who were at least 50 years of age and whose ASL scans had passed the quality control procedure (see MRI acquisition and processing). All participants gave written informed consent, and the study was approved by the Ethical Review Board in Lund, Sweden.
MRI acquisition and processing
In accordance with the recommendations of the white paper (Alsop et al., 2015), ASL scans were acquired using a prototype 3D pseudo-continuous (pCASL) sequence with background suppression and gradient- and spin-echo (GRASE) readout in a 3T MAGNETOM Prisma scanner with a 64-channel receiver-coil array (Siemens Healthcare, Erlangen, Germany). The labeling plane was planned using an MR-angiogram and the following parameters were used for the pCASL sequence: TR | TE = 4600 | 21.76 ms, flip angle (FA) = 28°, labeling duration | post-labeling delay = 1500 | 2000 ms, reconstructed voxel size = 1.9 × 1.9 × 4 mm3, FoV = 240 mm, GRAPPA factor = 2, Turbo factor = 12, EPI factor = 31, 3 segments, 36 axial slices, 12 pairs of control/label images and total acquisition time of 5:50 min. A proton density weighted (M0) volume with TR = 4000 ms was also acquired to calibrate estimated CBF into absolute units. Halfway through the study, the ASL prototype sequence was updated. The protocol was adjusted to optimally resemble the pre-update version. Additionally, a whole-brain T1-weighted scan (MPRAGE sequence, TR | TE | TI = 1900 | 2.54 | 900 ms, resolution = 1 × 1 × 1 mm3, FA = 9°, FoV = 256 × 256 mm2, 176 slices, bandwidth = 220 Hz/px TA = 5:15 min) and a T2-weighted FLAIR scan (TR | TE | TI = 5000 | 393 | 1800 ms, band-width = 781 Hz/px TA = 4:37 min, same resolution and FoV as for the T1-weighted image) were obtained. ASL and structural scans were processed using the publicly available ExploreASL tool-box (https://github.com/ExploreASL/ExploreASL) developed in MATLAB which uses the Statistical Parametric Mapping (SPM 12) routines. The processing pipeline has been described in detail else-where (Mutsaerts et al., 2020). Briefly, white matter hyperintensities (WMH) were segmented on the FLAIR image and used to fill the corresponding hypointensities on T1-weighted scan. Voxel intensities in these lesion regions were replaced by bias-field corrected values from the surrounding normally appearing WM on T1 image. Subsequently, the anatomical image was segmented into grey matter (GM), WM and CSF using computational anatomy toolbox 12 (CAT12) to obtain partial volume (PV) maps which were normalized to the MNI152 space. ASL scans were motion-corrected and co-registered to the T1 scan. The M0 image was smoothed with a 16 mm FWHM Gaussian kernel. CBF was quantified using a single compartment model (Alsop et al., 2015) as follows: where λ is the brain/blood partition coeffcient (0.9 ml/g), PLD is post-labeling delay (2000 ms), T1a is the longitudinal relaxation time of arterial blood, α is the labeling effciency (85%), M0 is the proton density image intensity and τ is the labeling duration (1500 ms). Subsequently, PV effects were corrected and the CBF images were transformed into MNI152 space. Finally, ASL scans were quality controlled by an automatic method which ranks the processed images based on their spatial coeffcient of variation, reflecting the amount of arterial transit time artefacts (Mutsaerts et al., 2019). This procedure was reviewed by a trained radiologist (S.I., with 5 years of experience) to detect vascular or other sources of artefacts. Region of interest (ROI) masks were obtained by intersecting the Harvard-Oxford (HO) atlas with individual GM and WM masks. Mean bilateral PV-corrected CBF values were calculated for the atlas ROIs in the standard space (Mutsaerts et al., 2020).
Regions of interest definition for CBF
In line with previous studies, regional CBF values for the atlas ROIs were normalized to CBF of the precentral gyrus to adjust for individual variations in the blood flow (Mattsson et al., 2014; Yew et al., 2017). No between group CBF differences were found in this region, prior to the normalization. A set of a priori ROIs were defined based on previous studies showing CBF alterations within these regions in the AD spectrum. These a priori defined ROIs encompassed lateral and medial parts of the parietal and temporal cortex, middle frontal gyrus and superior division of the lateral occipital cortex (Figure 1). The lateral and medial parietal meta-ROIs consisted of anterior/posterior supramarginal gyri, angular gyrus and posterior cingulate plus precuneus, respectively. The lateral temporal meta-ROI included anterior/posterior divisions of the middle and inferior temporal gyri while hippocampus and anterior/posterior divisions of parahippocampal gyrus (anterior parahippocampal gyrus in HO atlas represents the entorhinal cortex) were grouped into the medial temporal meta-ROI. The weighted-mean CBF for a priori ROIs was obtained by the sum of products of the mean CBF in each sub-ROI and its volume, divided by the sum of the volumes of all the sub-ROIs constituting each priori ROI.
PET acquisition and processing
Tau-PET scans using [18F] RO948 were performed on a GE Discovery MI scanner (General Electric Medical Systems) 70 - 90 minutes post-injection. [18F] flutemetamol Aβ-PET scans were acquired on the same scanner 90 - 110 minutes following the injection as previously described (Leuzy et al., 2020). Standardized uptake value ratio (SUVR) images were calculated using the inferior cerebellar GM and pons as the reference regions for tau- and Aβ-PET, respectively. The aforementioned T1-weighted scan was used for PET image co-registration and template normalization. Note that the Aβ-PET data were only available for a subset of 183 participants, including 87 Aβ-negative, 40 Aβ-positive CU individuals and 56 Aβ-positive CI patients. Tau- and Aβ-PET composite ROIs were created using the Freesurfer-based anatomical parcellation (https://surfer.nmr.mgh.harvard.edu/).
CSF biomarkers
Lumbar CSF sampling and analysis were conducted according to the Alzheimer’s Association Flow Chart (Blennow et al., 2010). As previously described, the analysis of P-tau217 assays was performed using the Meso Scale Discovery (MSD) platform (Janelidze et al., 2020). Total tau (T-tau) was quantified using Innotest® immunoassay (Fujirebio; Gent, Belgium). CSF Aβ42 and Aβ40 levels were measured using MSD or Fujirebio Lumipulse assays according to the manufacturer’s instructions, and the ratio values were subsequently calculated. NfL was measured using a Simoa kit (Quanterix; Billerica, MA). NPTX2 was measured in the facilities of ADx (Gent, Belgium) using a research grade sandwich ELISA. The number of available data per group for each of the CSF and PET biomarkers, together with the corresponding mean values are shown in Supplementary Table 4.
Statistical analyses
The ROI-based analyses were performed using R (version 3.5.2). Group differences across demographics were compared using analysis of variance (ANOVA) or chi-squared X2 tests (table 1). Global CBF differences between the diagnostic groups were evaluated by univariate linear regression models with CBF of a priori ROIs as the dependent variable and the groups as the independent variable. To assess the relationships between CBF and each of the PET and CSF markers of Aβ, tau, synaptic dysfunction and axonal integrity, the regression analysis was separately applied to the AD continuum, and CU participants with the CBF in a priori ROIs as the dependent variable and each biomarker as the independent variable. All the analyses were controlled for the confounding effects of age, sex and ASL sequence version. The distribution of NfL was skewed and required a logarithmic transformation. The results were considered statistically significant at p < 0.05. False discovery rate (FDR) correction was used to adjust for multiple comparisons.
Voxel-wise analyses
In addition to the ROI-based analyses, whole-brain voxel-wise analyses were performed to compare CBF variations between groups and to assess the spatial associations of CBF with Aβ, tau and other CSF markers in both CU individuals and those on the AD spectrum. All voxel-wise analyses were controlled for age, sex and ASL sequence version. The CBF images were multiplied by a GM mask and subsequently smoothed with a Gaussian kernel of 4 mm FWHM. The analyses were performed using threshold-free, cluster-enhanced (TFCE) permutation statistics using FSL randomise command (https://fsl.fmrib.ox.ac.uk/fsl/fslwiki/Randomise) with 5000 permutations. Statistical significance was set at the family-wise error (FWE) corrected threshold of p < 0.05.
Disease progression modeling
Subtype and Stage Inference (SuStaIn) was employed to characterize the progression of AD (Young et al., 2018) in terms of different biomarkers via the pySuStaIn package (https://github.com/ucl-pond/pySuStaIn). SuStaIn is an unsupervised machine learning approach that combines clustering and disease progression modeling to unravel the phenotypic and temporal heterogeneity of any progressive disease (Young et al., 2018, 2020). It simultaneously identifies subgroups of individuals with common phenotypes (subtypes) and distinct temporal progression patterns (stages). This data-driven approach requires inputs expressed as z-scores relative to a control population. It builds on event-based modeling (Fonteijn et al., 2012; Young et al., 2014) delineating the disease progression as a series of events each of which corresponds to a biomarker switching from a normal to an abnormal level. SuStaIn replaces the instantaneous normal to abnormal transition of events with a linear z-score model in which the trajectory of each subtype is described by linear accumulation of a biomarker from one z-score to another over different stages. Markov Chain Monte Carlo sampling is then used to quantify the uncertainty of the sequence ordering (Young et al., 2018) upon which a positional variance diagram can be visualized. Seven biomarkers that are most commonly associated with AD were fed into SuStaIn. They comprised CSF levels of Aβ42/40 and P-tau217, tau-PET deposition in early, temporal/limbic and late ROIs, a composite volumetric measure in AD signature regions i.e., bilateral entorhinal, inferior/middle temporal and fusiform cortices (Ossenkoppele et al., 2018; Jack Jr et al., 2017) and weighted mean CBF of regions having voxel-wise associations with tau-PET load in the early ROI (Figure 3 A) in the AD spectrum. Note that in the Aβ-negative CU individuals, the CBF values were extracted from the same voxel-wise correlation maps (see Figure 3 A). The values for each biomarker were represented as a z-score relative to Aβ-negative CU subjects who were younger than 50 years old (control population; N = 44), and thresholds were set at 1, 2 and 3 standard deviations. Except for being a reference for SuStaIn, the data from the control population were not included in any other analyses. Given that SuStaIn cannot handle missing data, participants with at least one missing value in any of the aforementioned markers were excluded. As such, the model was run on a sample of 203 individuals.
Data Availability
Anonymized data will be shared by request from a qualified academic investigator for the sole purpose of replicating procedures and results presented in the article as long as data transfer is in agreement with EU legislation on the general data protection regulation and decisions by the Ethical Review Board of Sweden and Region Skåne, which should be regulated in a material transfer agreement.
Data availability
Anonymized data will be shared by request from a qualified academic investigator for the sole purpose of replicating procedures and results presented in the article as long as data transfer is in agreement with EU legislation on the general data protection regulation and decisions by the Ethical Review Board of Sweden and Region Skåne, which should be regulated in a material transfer agreement.
Funding
The present study was supported by the European Research Council (grant no. 311292), the Swedish Research Council (grants no. 2016-00906 and 2018-02052), the Knut and Alice Wallenberg Foundation (grant no. 2017-0383), the Marianne and Marcus Wallenberg Foundation (grant no. 2015.0125), the Swedish Alzheimer Foundation (grant no. AF-745911), the Swedish Brain Foundation (grant no. FO2019-0326), the Swedish federal government under the ALF agreement (grants no. 2018-Projekt0279 and 2018-Projekt0226) and the Strategic Research Area MultiPark (Multidisciplinary Research in Parkinson’s disease) at Lund University. F.B. was supported by the NIHR biomedical research Centre at UCLH.
Competing interests
O.H. has acquired research support (for the institution) from AVID Radiopharmaceuticals, Biogen, Eli Lilly, Eisai, GE Healthcare, Pfizer, and Roche. In the past 2 years, he has received consultancy/speaker fees from AC Immune, Alzpath, Biogen, Cerveau and Roche. J.P. is an employee of Siemens Healthcare. All other authors report no competing interests.
Supplementary materials
Acknowledgments
The authors thank the participants of the Swedish BioFINDER-2 study for their enthusiastic commitment.