Participants
We studied 220 participants (mean age = 61.4y, range: 46.9-76.8y, 151 female) selected from the WRAP cohort based on the availability of at least one 11C-PiB-PET scan. Participants underwent the first 11C-PiB-PET scan on average 6.4 years after the first neuropsychological assessment. Among participants, 159 (72%) had at least one biological parent diagnosed with AD dementia. A subset of 157 participants had available imaging data from a second 11C-PiB-PET scan (average time interval of 2.5 years since the first 11C-PiB-PET scan), and 60 had additional imaging data from a third 11C-PiB-PET scan (average time interval of 6.1 years since the first scan). Participants were healthy and unimpaired at baseline. Eight participants had a diagnosis of clinical MCI at a WRAP study visit prior to the first 11C-PiB-PET scan and one had a dementia diagnosis shortly after the first 11C-PiB-PET scan. We excluded data from these participants from the longitudinal analysis of cognitive performance. Out of 189 participants with available APOE-ε4 data, 34% were APOE-ε4 positive.
Imaging data
Acquisition of 11C-PiB-PET and MRI imaging data in the WRAP cohort has been described in detail previously (23). Briefly, 11C-PiB-PET scans were acquired in a 3-D mode with a dynamic 70-minute acquisition protocol after an injection of a 15 mCi target dose of 11C-PiB bolus. Dynamic acquisition frames consisted of 17 time frames, including 5 × 2 minutes and 12 × 5 minutes frames. A filtered back-projection algorithm was used for reconstructing the data. For anatomical reference, a high-resolution T1-weighted MRI scan was acquired using a 3.0 Tesla GE MR750 scanner with an 8 or 32 channel head coil. The 3-D inversion recovery prepared fast spoiled gradient-echo sequence had following parameters: inversion time (TI) = 450 ms; repetition acquisition matrix = 256 x 256 x 156 mm, field of view (FOV) = 256 mm; slice thickness = 1.0 mm. The reconstructed time series of 11C-PiB-PET data were realigned, corrected for motion, de-noised, and coregistered to the subject’s T1-weighted MRI scan based on co-registration of the time-integrated PET scan utilizing Statistical Parametric Mapping software (SPM12; www.fil.ion.ucl.ac.uk/spm). Parametric distribution volume ratio (DVR) maps were generated using Logan graphical analysis methods (24, 25) with t* = 35 minutes and cerebellar gray matter as a reference region of non-displaceable binding.
Image analysis
The imaging data were further pre-processed for regional staging analysis using previously described procedures (10). MRI images were segmented into different tissue types and spatially normalized to a customized aging/AD-specific reference template space (26) using the high-dimensional spatial registration algorithm DARTEL (27). 11C-PiB-PET DVR maps were corrected for PVE using the 3-compartment “Müller-Gärtner” method in subject’s native space (28, 29), and then spatially normalized to the reference template space using transformation parameters from the corresponding MRI. Regional DVR values were then extracted from 52 regions of interest within the reference template defined using the Harvard-Oxford atlas, which included 48 cortical regions, as well as the hippocampus, amygdala, striatum, and thalamus. We also extracted the average global non-PVE-corrected and PVE-corrected 11C-PiB-PET DVR signal within a cortical composite mask (30).
We used a two-dimensional Gaussian Mixture Model (GMM) approach utilizing regional and global mean PVE-corrected DVR values to establish region-specific thresholds for amyloid-positivity. Analogous to previous studies using one-dimensional GMM (8, 31, 32), we fit low and high amyloid distributions for each region. The two-dimensional GMM approach is different in that it estimates distribution of two variables at once so that the contribution of each regional DVR value to the low or high amyloid distribution is estimated in conjunction with the global amyloid signal of each participant. This approach was intended to decrease the susceptibility of the procedure to potential noisiness of regional signal resulting in more robust and biologically plausible regional estimates. Regional thresholds were defined as 1.65 standard deviations above the mean value of the low Aβ distribution corresponding to the 95th percentile (16, 33).
In analogy to neuropathological staging models and our previous PET-based staging study, we determined a regional amyloid progression model based on the frequency of regional amyloid-positivity across individuals as an indicator of progressive temporal involvement (4, 9, 10, 34). Regional frequencies of amyloid-positivity were calculated from the baseline PET data using 10,000 bootstrap resamples, and the obtained range of frequencies was split into four equal parts to obtain a discrete stage model of amyloid progression across four larger anatomical divisions (10, 17).
In sensitivity analyses, we additionally assessed the effect of alternative strategies for estimating regional positivity thresholds, including a 1-dimensional GMM approach based on regional values only, and a regional resampling approach in a subsample of the 20 youngest, APOE-ε4 negative subjects without familial history of AD (mean age = 59.8y, 16 female). For both of these methods, the thresholds were analogously estimated as 1.65 standard deviations above the mean value, and regional frequencies of amyloid-positivity were calculated using 10,000 bootstrap resamples. The correspondence between the regional amyloid-positivity frequencies derived from the different cut-off methods was assessed using pair-wise Spearman rank correlations.
Individual amyloid deposition profiles were staged according to the regional hierarchy indicated by the estimated amyloid progression model. For that, each of the four larger anatomical divisions defined by the 4-stage model was considered amyloid-positive if at least half of the included regions displayed a suprathreshold signal (10, 17). The individual stage was then determined based on the corresponding amyloid-positive anatomical divisions. For example, a classification of stage III requires positivity in the anatomical divisions 1, 2 and 3, but not 4. Participants whose regional amyloid-positivity profile did not adhere to the expected regional hierarchy (e.g., positivity in anatomical division 2, but not in 1) were classified as non-stageable. For comparison, we dichotomized the 11C-PiB-PET scans into standard amyloid-positive/-negative categories based on a previously established threshold of 1.08 applied to the global composite DVR value in non-PVE-corrected data (8).
Longitudinal imaging analysis
The longitudinal validity of the cross-sectionally estimated regional amyloid staging model was assessed in two complementary analyses. First, we assessed individual longitudinal changes in amyloid stages from baseline to the furthest available follow-up PET scan. Among stageable participants at baseline, 155 had a follow-up 11C-PiB-PET scan with an average time delay of 4 years (range: 1.7-7.7). In a complementary analysis independent from the estimated staging model, we assessed the first longitudinal appearance of regional amyloid-positivity in subjects who had no suprathreshold signal in any of the 52 brain regions at baseline (n = 64) by recording the regional amyloid-positivity occurring at the follow-up 11C-PiB-PET scans, on average 3.9 years later (range: 1.8-7.6 years).
Neuropsychological testing
Finally, to examine the potential clinical relevance of the amyloid staging approach we analyzed longitudinal cognitive trajectories of participants at different amyloid stages using previously developed domain-specific and global cognitive composite scores (35). These scores included a delayed recall composite (THEO-DEL-REC), an executive function composite (THEO-EXEC-FN), an immediate learning composite (THEO-IMM-LRN), as well as a global cognitive composite score – a three test version of the preclinical Alzheimer’s cognitive composite (PACC3).
In order to assess differences in future cognitive trajectories across in-vivo amyloid stages, in the regression analysis we selected neuropsychological scores obtained at visits taking place at earliest three months before the first 11C-PiB-PET measurement and later. Five participants did not have available neuropsychological test scores after that time point and were excluded. The composite scores were only available from WRAP study visit 2 and onwards, because the more extensive cognitive testing required for the composite score calculation was not yet introduced at the first WRAP study visits. As a result, the closest WRAP visit with neuropsychological testing was on average 1.1 years after the first 11C-PiB-PET scan. We included cognitive data from a median of 3 WRAP study visits per participant conducted on average at 2.5-year intervals. The mean duration of the total follow-up was 6.5 years from the first 11C-PiB-PET scan until the last available cognitive assessment, with a maximum of 8.7 years of follow-up. Longitudinal trajectories of the four cognitive composite scores were analyzed using linear mixed-effects regression models implemented in R 3.6.0 (36). The effect of amyloid stage on longitudinal cognitive decline was assessed by the time*amyloid stage interaction, controlled for age at the first analyzed WRAP visit, sex and years of education.