Association of Alzheimer’s disease polygenic risk scores with amyloid accumulation in cognitively intact older adults
Alzheimer's Research & Therapy volume 14, Article number: 138 (2022)
Early detection of individuals at risk for Alzheimer’s disease (AD) is highly important. Amyloid accumulation is an early pathological AD event, but the genetic association with known AD risk variants beyond the APOE4 effect is largely unknown. We investigated the association between different AD polygenic risk scores (PRS) and amyloid accumulation in the Flemish Prevent AD Cohort KU Leuven (F-PACK).
We calculated PRS with and without the APOE region in 90 cognitively healthy F-PACK participants (baseline age 67.8 (52–80) years, 41 APOE4 carriers), with baseline and follow-up amyloid-PET (time interval 6.1 (3.4–10.9) years). Individuals were genotyped using Illumina GSA and imputed. PRS were calculated using three p-value thresholds (pT) for variant inclusion: 5 × 10−8, 1 × 10−5, and 0.1, based on the stage 1 summary statistics from Kunkle et al. (Nat Genet 51:414–30, 2019). Linear regression models determined if these PRS predicted amyloid accumulation.
A score based on PRS excluding the APOE region at pT = 5 × 10−8 plus the weighted sum of the two major APOE variants (rs429358 and rs7412) was significantly associated with amyloid accumulation (p = 0.0126). The two major APOE variants were also significantly associated with amyloid accumulation (p = 0.0496). The other PRS were not significant.
Specific PRS are associated with amyloid accumulation in the asymptomatic phase of AD.
Alzheimer’s disease (AD) is the most common form of dementia in the global population, with a complex interplay of both genetic and environmental factors . Changes observed in AD, such as the accumulation of brain amyloid-β protein, can occur more than a decade prior to symptom onset (the preclinical or asymptomatic phase) . Age is one of the largest risk factors for AD risk, where many of the observed pathological changes occur as age increases . It is of high importance to develop methods that detect those individuals at risk of developing the disease prior to symptom onset, as this phase presents a window of opportunity for early disease diagnosis, treatment administration, and determining individuals suitable for prevention trials.
Apolipoprotein E ε4 (APOE4) is the largest genetic risk factor for development of sporadic AD  and has been shown to increase amyloid deposition in a dose-dependent manner, in the asymptomatic, the mild cognitive impairment and the dementia stages of AD [5,6,7]. Accumulation of brain amyloid in the asymptomatic phase occurs in a sigmoidal fashion, prior to reaching a plateau . However, genetic risk for amyloid accumulation is largely unknown above the known effect of APOE4.
Beyond APOE, recent large-scale genome-wide association studies (GWAS) have highlighted low effect size variants in non-APOE genes that modify AD dementia risk [9,10,11,12,13]. In isolation, the identified risk SNPs have limited use in predicting AD risk. Combinations of selected AD risk SNP scores have previously been implemented in case-control studies, and have shown an increased risk of AD with a higher score (e.g. ). A polygenic risk score (PRS) approach to disease risk prediction is able to take into consideration all SNPs that show an association with AD, despite their small effect sizes, by capturing the overall genetic risk for an individual into a single score. This has already been performed in AD case-control studies (e.g. [14,15,16]), but there is a requirement for investigation in the asymptomatic phase, as well as determining the genetic risk for amyloid accumulation.
There are still many factors unknown about how to best model the PRS, for example the optimal threshold for SNP inclusion (pT) and how to represent the effect of APOE (e.g. for AD risk). A seminal study from Escott-Price et al.  suggested a PRS with the highest prediction accuracy (AUC = 78.2%) in predicting AD cases from controls is built using a liberal threshold for SNP inclusion (pT < 0.5), suggesting a polygenic architecture to AD. However, more recent studies have implied a more oligogenic architecture to AD, in which there are fewer SNPs associated with disease risk . This suggests the use of more stringent (closer to genome-wide significant) pTs for the calculation of PRS for AD.
A more recent study from Leonenko et al.  further evaluated how best to model the APOE4 effect in the context of PRS and found that the best prediction accuracy to predict AD cases from controls (AUC = 74.1%) came from a combination of a PRS excluding the APOE region with the addition of the weighted sum of two directly genotyped APOE ε2 and ε4 alleles (rs429358, rs7412, PRSnoAPOE+APOEε2+ε4) and a threshold pT < 0.1. This was when comparing PRS built using different SNP combinations, including PRSnoAPOE and PRS including all possible SNPs, at varying thresholds.
We aimed to evaluate PRS built as PRSnoAPOE+APOEε2+ε4, using different thresholds for SNP inclusion to determine the association with amyloid accumulation. We included cognitively intact older individuals participating in the Flemish Prevent AD Cohort KU Leuven (F-PACK), some of whom are in the asymptomatic stage of AD. We compared the performance of PRSnoAPOE+APOEε2+ε4 to other PRS.
The Laboratory for Cognitive Neurology follows a cohort of 180 community-recruited deeply phenotyped older adults, known as F-PACK. F-PACK individuals were recruited between 2009 and 2015 in three waves of 60 individuals. At recruitment, individuals had to be aged between 50 and 80 years old, score ≥ 27 on the mini mental state examination (MMSE) and score zero on the clinical dementia rating (CDR) scale. Furthermore, individuals had to score within published norms on an extensive neuropsychological test battery [19, 20]. Exclusion criteria included a history of neurological or psychiatric illness, contraindication for magnetic resonance imaging (MRI), focal brain lesions on MRI, history of cancer, or exposure to radiation one year preceding the baseline amyloid positron emission tomography (PET) scan. Recruitment was stratified for two genetic factors: APOE4 allele (present or absent) and BDNF status (66 met allele present or absent). This was carried out such that per 5-year age bin each factorial cell contained the same number of individuals matched for age, sex, and education. All recruited individuals received structural MRI and an 18F-Flutemetamol amyloid-PET scan at baseline. Participants are invited for 2-yearly neuropsychological assessments over a 10-year period.
F-PACK individuals received blood collection for genotyping at recruitment. A subset of 90 participants have additionally received a follow-up amyloid-PET scan, on average 6.1 (3.4–10.9) years after baseline.
The protocol was approved by the Ethics Committee University Hospitals Leuven. All participants provided written informed consent in accordance with the declaration of Helsinki.
Structural MRI acquisition
At baseline and follow-up, each participant received a high resolution T1-weighted structural MRI scan for the PET processing procedure described below. Scans were performed using a 3T Philips Achieva dstream 32-channel headcoil MRI scanner (Philips, Best, The Netherlands). All baseline and 69 follow-up scans were acquired using the same 3D turbo field echo sequence: repetition time = 9.6 ms; echo time = 4.6 ms; flip angle = 8°; field of view = 250 × 250 mm; 182 slices; voxel size 0.98 × 0.98 × 1.2 mm3. Eighteen follow-up scans were acquired using a three-dimensional magnetisation-prepared rapid gradient-echo sequence, due to being acquired as part of The Amyloid Imaging to Prevent Alzheimer’s Disease study (AMYPAD): repetition time = 6.6 ms; echo time = 3.1 ms; flip angle = 9°; field of view = 270 × 252 mm; 170 slices; voxel size 1.05 × 1.05 × 1.2 mm3. Three individuals refused a follow-up scan.
18F-Flutemetamol PET acquisition and pre-processing
As previously described, 18F-Flutemetamol PET scans were acquired on a 16-slice Biograph PET/CT scanner (Siemens, Erlangen, Germany) at baseline and follow-up (6.1 (3.4–10.9) years), with a net injected intravenous dose of 149 MBq (127–162 MBq) and 156 MBq (77–198 MBq), respectively, with an acquisition window of 90–120 min post-injection [19, 21,22,23,24]. Four follow-up PET scans had an acquisition window of 90–110 min due to acquisition prior to a protocol amendment. A low-dose CT scan was acquired prior to all scans for attenuation correction. Random and scatter corrections were applied. Scans were reconstructed as frames of five minutes. All scans were reconstructed using ordered subsets expectation maximisation. All baseline scans and 73 follow-up scans were reconstructed as five iterations in eight subsets. Sixteen follow-up scans were reconstructed as four iterations in 21 subsets, due to being reconstructed as part of AMYPAD. The spatial resolution of the scanner is 4.6 mm full width at half maximum 1cm off centre measured with the NEMA protocol, and all scans were smoothed with a 5-mm full width at half maximum Gaussian filter.
18F-Flutemetamol image analysis
Statistical Parametric Mapping version 12 (SPM12, Wellcome Trust Centre for Neuroimaging, London, UK, http://www.fil.ionucl.ac.uk/spm) running on MATLAB R2018b (Mathworks, Natick, MA, USA) was used to process the images, as described previously [19, 21,22,23,24].
We calculated (mean) standardised uptake value ratios (SUVRs) from the spatially normalised images (voxel size: 2 × 2 × 2 mm3) in a composite cortical volume of interest (SUVRcomp) based on the automated anatomical labelling atlas (AAL). This composite volume of interest included the following bilateral regions: frontal (AAL areas 3–10, 13–16, 23–28), parietal (AAL 57–70), anterior cingulate (AAL 31–32), posterior cingulate (AAL 35–36), and lateral temporal (AAL 81–82, 85–90) and was masked with the participant-specific grey matter (GM) segmentation map (threshold = 0.3) [19, 25]. Cerebellar grey matter was used as the reference region to calculate SUVRcomp, defined as AAL areas 91–108 and was masked by the participant-specific GM map (intensity threshold = 0.3) . SUVRcomp were then converted to Centiloids (CL) using the formula CL = 127.6 × SUVRcomp–149 [26, 27]. Individuals were deemed amyloid positive when their CL was ≥ 23.5, a pathologically confirmed cut-off for amyloid positivity . Furthermore, to model amyloid change longitudinally, we calculated rate of change as:
For use in secondary analyses, we repeated the PET image processing pipeline including a partial volume correction (PVC, 5 mm full-width half-maximum) step, using a modified Müller-Gärtner procedure . Rate of change was derived, as per the above equation using SUVRs, of these corrected images for use in the analyses described below.
Structural MRI longitudinal change
In a secondary analysis, we also examined whether any PRS with significant predictive value for amyloid accumulation also had any predictive value for longitudinal MRI volume loss. As a measure of atrophy, we used the GM maps at baseline and follow-up, produced from the PET processing procedure, to calculate the percental whole brain GM volume change between the two time points as:
Furthermore, for a voxelwise analysis we calculated GM change images, as follow-up GM map minus baseline GM map, using SPM12, for use in voxelwise regression analyses described below.
Modelling of longitudinal cognitive change
For the PRS with predictive value for amyloid accumulation, we also examined whether this predictive power extended to cognitive change. To model cognitive change, based on a prior study , the mean Buschke Selective Reminding Test Total Retention (BSRT TR) score  was selected and slopes were calculated using latent growth curve analysis, using the R package Lavaan . Individuals are being followed for a period of 10 years with two-yearly neuropsychological evaluations, so we used mean BSRT TR scores from all available neuropsychological testing time points (FUY) up until the closest to the follow-up amyloid-PET scan (number of individuals at each FUY: 4 FUY2, 37 FUY4, 32 FUY6, 10 FUY8, 7 FUY10; mean interval between follow-up amyloid-PET and closest neuropsychological assessment is 13.06 ± 12.04 months). Missing values were imputed in R using the CART imputation method in the package mice . The calculated mean BSRT TR slopes were then used further, where a more negative slope represents steeper decline.
Genetic data acquisition and processing
DNA was further available for 177 F-PACK participants, which was subsequently genotyped using the Illumina Global Screening Array (GSA, coverage of 657,598 SNPS) in collaboration with the Institute of Clinical and Medical Biology (University Hospital Schleswig-Holstein, Germany) . Standard quality control (QC) was performed using PLINK (Version 1.9, www.cog-genomics.org/plink/1.9) and included SNP call rate ≥ 0.95, minor allele frequency (MAF) ≥ 0.01, and outlying heterozygosity (± 5 standard deviations) . Hardy-Weinberg equilibrium threshold = 1 × 10−6 was also applied. Ethnic outliers were detected using the Phase 3 1000 Genomes (1KG) dataset (N = 2504 ). Imputation was performed using the Michigan Imputation Server (https://imputationserver.sph.umich.edu)  and Haplotype Reference Consortium reference panel (http://www.haplotype-reference-consortium.org), resulting in 39,131,578 SNPs. Data were filtered with imputation information score > 0.7 and MAF ≥ 0.01. After imputation and QC 7,466,483 SNPs remained for further analysis. Three individuals were removed during the QC process (one due to sample duplication, one due to relatedness (pi-hat > 0.2), and one due to outlying heterozygosity).
Polygenic risk score calculation
PRS calculations were performed using PRSice-2 . The stage 1 summary statistics from the GWAS performed by Kunkle et al.  were used as the base file, and the European individuals from 1KG (N = 503 ) were used as an external reference panel for clumping to remove SNPs in high linkage disequilibrium (clumping window = 250 kilobases, r2 = 0.1). Post-clumping, there were 335,326 SNPs remaining for PRS calculations. We first calculated PRS using the PRSnoAPOE+APOEε2+ε4 approach . We calculated PRS excluding the APOE region (chromosome 19: 45-48.8 Mb, PRSnoAPOE) at three thresholds for SNP inclusion (pT): pT = 0.1; 1 × 10−5; 5 × 10−8, and then added the weighted sum of the two major APOE SNPs (rs429358 and rs7412, APOEε2+ε4), using the effect sizes from Kunkle et al. , to each score (N1 * β1 + N2 * β2, where N1 and N2 are the number of alleles for each APOE SNP, respectively, and β1 and β2 are the corresponding effect sizes).
For comparison, we also calculated PRSAD, in which all available SNPs were included at each pT, as well as PRSAPOEonly, in which PRS were calculated at each pT using only the APOE region specified above. All PRS were z-score normalised prior to further analyses.
Statistical analyses were performed in R version 4.1.0 (2021-05-18; The R Foundation for Statistical Computing; https://cran.r-project.org/). Prior to analyses, Shapiro-Wilk tests were used to determine data normality.
Cohort characteristics were assessed between APOE4 carriers and non-carriers using Wilcoxon rank sum tests with continuity correction or Welch two-sample t-tests for continuous data, depending on normality, and χ2 tests for categorical data.
For our primary analysis, linear regressions were performed for each PRSnoAPOE+APOEε2+ε4 (i.e. at 3 different SNP inclusion thresholds) as predictor, and amyloid rate of change as outcome variable. Baseline age, sex, and the first three principal components were included as covariates. Inference was based on an uncorrected p < 0.05 threshold divided by the number of SNP inclusion thresholds (N = 3).
We investigated how PRSnoAPOE+APOEε2+ε4 performed compared to other PRS, thus linear regressions were also performed as described above with PRSAD, APOEε2+ε4, PRSnoAPOE, or with PRSAPOEonly. Per type of PRS, inference was based on an uncorrected p < 0.05 threshold divided by the number of SNP inclusion thresholds (N = 3), except for APOEε2+ε4 given that this score does not rely on p-value thresholds for SNP inclusion as it is built using only the weighted sum of two APOE SNPs. We also investigated whether all calculated PRS were associated with baseline amyloid burden using the same statistical approach.
In a further secondary analysis, we determined the ability for any significant PRS to discriminate individuals who were amyloid negative at both time points, on the one hand, from individuals who were amyloid negative at baseline and positive at follow-up or amyloid positive at both time points, on the other hand. We performed a set of receiver operating characteristic (ROC) analyses using the R package pROC based on a logistic regression classifier . We included two demographic models (one model with age + sex; one model with age + sex + APOE4 status (yes/no)) and models including any significant PRS determined from the previous analyses, with age and sex. Measures of performance are given as areas under the curve (AUCs) and 95% confidence intervals (CI). Pairwise comparisons of model AUCs were performed using the DeLong method .
To evaluate the effect of PVC, we also carried out the regression analyses with the longitudinal change in PVC SUVRs as predictor.
Finally, we performed the regression analyses for cognitive decline or GM atrophy, with any PRS that had a significant association with amyloid rate of change. In addition, we performed whole-brain voxelwise t-tests using the GM change maps, with the significant PRS and time interval as covariates. We used a FWE corrected cluster-level threshold of p < 0.05 with the voxel-level set at uncorrected p < 0.001 .
From the 90 F-PACK individuals with follow-up amyloid-PET included in this study, nine were amyloid positive at baseline (10%) and 21 at follow-up (23%, Table 1, Fig. 1). APOE4 carriers had a significantly higher amyloid rate of change (rate of change for APOE4 carriers: 1.06 (range: − 2.85–13.92); APOE4 non-carriers: − 0.03 (range: − 3.03–4.22) p = 0.01, outlier excluded p = 0.02) and a significantly lower time interval between amyloid-PET scans than non-carriers (time interval for APOE4 carriers: 5.4 (range: 3.4–10.0) years; APOE4 non-carriers 6.1 (range: 4.0–10.9) years, p = 0.008). There was also a significantly higher number of APOE4 carriers classified as amyloid positive at follow-up (34%) compared to non-carriers (14%, p = 0.049). Two individuals had a CDR that had evolved to 0.5 at the time of follow-up amyloid-PET.
There was an outlier with amyloid rate of change (as determined by Grubb’s test in R, using the package Outliers). Therefore, the amyloid analyses were repeated excluding this individual.
Primary analyses: PRSnoAPOE+APOE ε2+ε4
PRSnoAPOE+APOEε2+ε4 had a significant association with amyloid rate of change when the SNP inclusion threshold (pT) = 5 × 10−8 (p = 0.0126, β = 0.68 (95% CI 0.15, 1.20), R2 = 0.12, Figs. 2 and 3A, Table 2). When the outlier was removed, significance remained (p = 0.004, β = 0.65 (95% CI 0.21, 1.09), R2 = 0.11, Fig. 3B). Figure 2 shows that the highest R2 and lowest p-value occur when the pT is more stringent.
Both regression plots in Fig. 3A and B highlight increasing amyloid rate of change with increasing PRS. One can also appreciate the lower PRS are for those participants with ε2ε3 or ε3ε3 APOE genotypes, whereas, in general, those participants with higher PRS carry at least one ε4 allele. No significance was found at the other pTs.
Secondary analyses: other PRS models and phenotypes
From our secondary analyses, there was a significant association between amyloid rate of change and APOEε2+ε4 (p = 0.0496, β = 0.56 (95% CI 0.0009, 1.09), R2 = 0.09, Fig. 2, Table 2). When the outlier was removed, the significance remained (p = 0.0201, β = 0.56 (95% CI 0.09-1.00), R2 = 0.08). None of the other PRS were significant with the total group. The results from the total cohort can be observed in Fig. 2, where red points represent significant p-values, and blue points are non-significant. From the figure, one can also appreciate the higher variance explained by the PRS on amyloid accumulation in the regressions when using the more stringent SNP inclusion thresholds. Of note is the higher R2 with PRSnoAPOE+APOEε2+ε4 when pT = 5 × 10−8 (R2 = 0.12) compared to when using APOEe2+e4 alone (R2 = 0.09). No PRS were associated with baseline amyloid load (Supplementary Table 1).
After performing PVC, the p-value for PRSnoAPOE+APOEε2+ε4 at pT = 5 × 10−8 was comparable to that obtained without PVC (values for PVC SUVRs: p = 0.0196 (β = 0.008 (95% CI 0.001, 0.014), R2 = 0.13). Note that we used SUVRs for the PVC analysis; hence, this has to be compared with the same analysis for the primary analysis using SUVRs as well. At pT = 5 × 10−8, the regression output for the primary analysis with non-PVC SUVRs is p = 0.0131 (β = 0.005 (95% CI 0.001, 0.009), R2 = 0.11). Significance was lost for APOEε2+ε4 (p = 0.30, β = 0.005 (95% CI − 0.002, 0.011), R2 = 0.10) after PVC correction. For comparison, the regression output with SUVRs for APOEε2+ε4 correspond to: p = 0.051, β = 0.004 (95% CI − 0.00002, 0.009), R2 = 0.08).
Figure 4 shows the results from the ROC analysis. We determined the ability of PRSnoAPOE+APOEε2+ε4 and of APOEε2+ε4 to discriminate amyloid negative individuals at both time points, on the one hand, from individuals who were amyloid negative at baseline and positive at follow-up or amyloid positive at both time points, on the other hand. The best performing model was age + sex + PRSnoAPOE+APOEε2+ε4 with pT = 5 × 10−8 (AUC = 0.74, 95% CI = 0.62–0.86). Numerically, the second best performing model was the demographic model consisting of age + sex + APOE4 status (yes/no) (AUC = 0.71, 95% CI = 0.58–0.84). This had a similar performance to the model with age + sex + APOEε2+ε4 (AUC = 0.70, 95% CI = 0.57–0.83). The demographic model including age + sex (AUC = 0.63, 95% CI = 0.51–0.75) performed the worst. None of the models were significantly different from each other when performing pairwise model comparisons (Supplementary Table 2).
PRSnoAPOE+APOEε2+ε4 and APOEε2+ε4 were not associated with cognitive decline over the time course measured. Neither were these PRS associated with longitudinal whole-brain grey matter atrophy nor were there significant effects based on the whole-brain voxelwise analyses (Supplementary Table 3).
Our study has highlighted that specific AD PRS are associated with amyloid accumulation in the asymptomatic phase of the disease, when built using a stringent SNP inclusion threshold.
The present study primarily investigated the association of amyloid change and PRSnoAPOE+APOEε2+ε4, given PRS built in this way have recently been shown to give the best prediction accuracy to predict AD cases from controls (AUC = 74.1%) with pT < 0.1 . In this previous paper, this PRS was superior at predicting AD cases from controls, compared to PRS built using different SNP combinations, including PRSnoAPOE and PRSAD, at varying thresholds. In that study, the best prediction accuracy for PRSAD reached an AUC of 69.8% (pT = 5 × 10−8) and PRSnoAPOE AUC = 61.3% (pT = 0.1). Furthermore, these were not better performing than the two major APOE SNPs (rs429358 and rs7412) alone (AUC = 70.0%) . The data presented in the present study also suggest that PRSnoAPOE+APOEε2+ε4 with pT = 5 × 10−8 is superior than the other PRS at having an association with amyloid accumulation in the asymptomatic phase of AD. The effect sizes observed for this score are higher than those observed in the other models. This score had a more significant association with amyloid accumulation than the major APOE SNPs, in the linear regressions, and produced, numerically, the highest AUC (AUC = 74%) compared to the demographic models (AUC for age + sex = 63% and AUC for age + sex + APOE4 status (yes/no) = 71%) and the model including the weighted sum of the two major APOE SNPs (AUC = 70%). The AUC of 74% is comparable to those reported in the literature in AD case-control prediction studies. This score also had a borderline significant association with amyloid rate of change after PVC was performed. The slight loss in significance is likely due to the introduction of additional noise from the PVC. This can occur when PVC is performed on PET scans of cognitively intact older adults, which pertains to the F-PACK cohort . However, the results further strengthen the hypothesis that there are other variants above APOE that are important for amyloidogenic processes. The increasing levels of amyloid in this asymptomatic phase may be associated with these lower effect size variants, whereas in the latter disease stages these associations are absent due to amyloid accumulation reaching a plateau .
This study investigated whether the thresholds for SNP inclusion that were found to be optimal in AD case-control studies were also applicable in the asymptomatic phase to detect amyloid changes. The results show stringent thresholds for SNP inclusion have a significant association with amyloid accumulation compared to the more liberal threshold of pT = 0.1. More liberal thresholds may be optimal for predicting cases from controls, but the present study suggests when trying to detect amyloid accumulation in asymptomatic AD, stringent thresholds for SNP inclusion, that reduce noise, are more optimal.
Previous studies have been unable to demonstrate an association between PRS and amyloid accumulation (e.g. [42, 43]). It is consistently found in many studies that APOE4 is associated with higher amyloid load and amyloid accumulation (in the asymptomatic phase [6, 44]), but the PRS-based studies were unable to provide evidence for an association between PRS and baseline amyloid nor PRS and amyloid accumulation. Consistent with previous studies, the data from the present study replicate the effect of APOE4 on amyloid accumulation, as well as providing a lack of association between PRS and baseline amyloid load. However, the present study did find an association between a specific PRS and amyloid accumulation over time. It can be observed that some individuals in the F-PACK cohort exhibit negative rates of change, along with the individuals that increase. Figure 1 shows the change in amyloid over time for the F-PACK cohort. The range in amyloid load below the threshold is narrow, whereas the change in amyloid above threshold is much larger, highlighting that the most change is driven by changes above the threshold of amyloid positivity. This change is not occurring below the threshold and change is thus minimal, even with those participants with a negative amyloid rate of change. Furthermore, a sensitivity analysis was conducted, whereby we removed the individuals in whom amyloid increased but did not surpass the threshold for positivity, and the regression still produced significant results with pT = 5 × 10−8 (data not shown due to the introduction of potential bias from removing data a posteriori). However, this further strengthens the primary analysis results that the observed effect is not driven by spurious accumulations in amyloid. Thus, the use of PRS to predict amyloid accumulation in the asymptomatic phase allows for testing genetic determinants of the amyloidogenic processes in sporadic AD, in the absence of downstream secondary effects in this early disease stage.
The effect of the APOE variants is clearly high given the results presented. This highlights that in the asymptomatic stage APOE4 has a large influence on the association of a PRS with amyloid change. When the threshold for SNP inclusion is more than genome-wide significant (> 5 × 10−8), the variance explained decreases, which is coupled with a loss in significance of the PRS. Thus, lower effect size variants play a role in AD risk, but increasing the threshold for SNP inclusion beyond genome-wide significance results in the addition of SNPs that create noise, thus suggesting an oligogenic architecture to (asymptomatic) AD. Nevertheless, the variance explained by the PRS on amyloid accumulation is low, thus there may be a large influence from other factors, such as gene-environment interactions.
Some limitations need to be considered. The amyloid scans were acquired between 90 and 120 min post injection, and modelling of cerebral blood flow (CBF) changes over time was therefore not possible. According to a cross-sectional study in the Baltimore Longitudinal Study of Aging, amyloid SUVR can be influenced by CBF and mostly so in individuals with amyloid PET values in the upper tertile . The effect of CBF on amyloid PET SUVR has been further quantified in a simulation study . The increases in CBF that are needed to account for a 1% change in SUVR is of the order of 5–15% increase in CBF, and higher CBF changes are needed to affect SUVR when amyloid load is lower (, Fig. 3). In a cognitively normal longitudinal cohort of individuals who remain relatively stable on a cognitive level and who do not manifest a neurological disease, a longitudinal CBF change of that order is implausible, also given the strict regulation of CBF under physiological circumstances. Note that CBF would need to increase over time with higher PRS scores. This can only be empirically excluded by concomitant longitudinal blood flow studies. The time interval between the baseline and follow-up amyloid-PET scans was variable, some individuals with a shorter interval than others. Amyloid accumulation is not a linear process; therefore, a difference in the onset of the rising phase of amyloid accumulation as well as the slope may be driving the differences we observe. The PRS were built using GWAS data from Kunkle et al.  that are not necessarily transferable to other ethnicities; thus, results should be carefully considered when inferring associations.
To conclude, a PRS built as PRSnoAPOE+APOEε2+ε4 with a stringent threshold for SNP inclusion had a more significant association with amyloid accumulation than the major APOE variants alone or than PRS built with other SNP combinations. This suggests an oligogenic, rather than polygenic, architecture to (asymptomatic) AD, in line with recent publications. The results may aid in participant recruitment and stratification for clinical trials, by identifying those individuals who are more susceptible to early brain amyloid changes, and thus more at risk to developing AD. According to the current dataset, PRSnoAPOE+APOEε2+ε4 outperforms the simple use of the APOE4 polymorphism alone for this purpose.
Availability of data and materials
Data are available upon reasonable request.
Automated anatomical labelling atlas
Apolipoprotein E ε4
Area under the curve
Animal Verbal Fluency
Rey’s Auditory Verbal Learning Test
Brain-derived neurotrophic factor
Boston Naming Test
- BSRT TR:
Buschke Selective Reminding Test Total Retention
Cerebral blood flow
Clinical Dementia Rating
Cerebellar grey matter
Flemish Prevent AD Cohort KU Leuven
Follow-up year neuropsychological testing time point
Global Screening Array
Genome-Wide Association Study
Letter Verbal Fluency
Minor allele frequency
Mini Mental State Examination
Magnetic resonance imaging
Psycholinguistic Assessment of Language Processing in Aphasia item 49
Positron emission tomography
Polygenic risk score
p-value threshold for single nucleotide polymorphism inclusion
Partial volume correction
Receiver operating characteristic
Raven’s Progressive Matrices
Single nucleotide polymorphism
Standardised uptake value ratio
- SUVRcomp :
Standardised uptake value ratio in the composite cortical volume of interest
Trial Making Test
Gatz M, Reynolds CA, Fratiglioni L, Johansson B, Mortimer JA, Berg S, et al. Role of genes and environments for explaining Alzheimer disease. Arch Gen Psychiatry. 2006;63:168–74 American Medical Association.
Villemagne VL, Burnham S, Bourgeat P, Brown B, Ellis KA, Salvado O, et al. Amyloid β deposition, neurodegeneration, and cognitive decline in sporadic Alzheimer’s disease: A prospective cohort study. Lancet Neurol. 2013;12:357–67.
Jack CR Jr, Knopman DS, Jagust WJ, et al. Hypothetical model of dynamic biomarkers of the Alzheimer's pathological cascade. Lancet Neurol. 2010;9(1):119–28. https://doi.org/10.1016/S1474-4422(09)70299-6.
Corder EH, Saunders AM, Strittmatter WJ, Schmechel DE, Gaskell PC, Small GW, et al. Gene dose of apolipoprotein E type 4 allele and the risk of Alzheimer’s disease in late onset families. Science (1979). 1993;261:921–3 American Association for the Advancement of Science.
Grimmer T, Tholen S, Yousefi BH, Alexopoulos P, Förschler A, Förstl H, et al. Progression of cerebral amyloid load is associated with the apolipoprotein e ε4 genotype in Alzheimer’s disease. Biol Psychiatry. 2010;68:879–84 NIH Public Access.
Villemagne VL, Pike KE, Chételat G, Ellis KA, Mulligan RS, Bourgeat P, et al. Longitudinal assessment of Aβ and cognition in aging and Alzheimer disease. Ann Neurol. 2011;69:181–92 NIH Public Access.
Yamazaki Y, Zhao N, Caulfield TR, Liu CC, Bu G. Apolipoprotein E and Alzheimer disease: pathobiology and targeting strategies. Nat Rev Neurol. 2019;15(9):501–18. https://doi.org/10.1038/s41582-019-0228-7.
Jack CR, Wiste HJ, Lesnick TG, Weigand SD, Knopman DS, Vemuri P, et al. Brain β-amyloid load approaches a plateau. Neurology. 2013;80:890–6 American Academy of Neurology.
Lambert JC, Ibrahim-Verbaas CA, Harold D, Naj AC, Sims R, Bellenguez C, et al. Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat Genet. 2013;45:1452 NIH Public Access.
Jansen IE, Savage JE, Watanabe K, Bryois J, Williams DM, Steinberg S, et al. Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer’s disease risk. Nat Genet. 2019;51:404–13 Nature Publishing Group.
Kunkle BW, Grenier-Boley B, Sims R, Bis JC, Damotte V, Naj AC, et al. Genetic meta-analysis of diagnosed Alzheimer’s disease identifies new risk loci and implicates Aβ, tau, immunity and lipid processing. Nat Genet. 2019;51(3):414–30 [cited 2021 Dec 20], Nature Publishing Group; Available from: https://www.nature.com/articles/s41588-019-0358-2.
Wightman DP, Jansen IE, Savage JE, Shadrin AA, Bahrami S, Holland D, et al. A genome-wide association study with 1,126,563 individuals identifies new risk loci for Alzheimer’s disease. Nat Genet. 2021;53:1276–82.
Bellenguez C, Küçükali F, Jansen IE, Kleineidam L, Moreno-Grau S, Amin N, et al. New insights into the genetic etiology of Alzheimer’s disease and related dementias. Nat Genet. 2022;54:1–25 Nature Publishing Group; [cited 2022 Apr 13]. Available from: https://www.nature.com/articles/s41588-022-01024-z.
Sleegers K, Bettens K, De Roeck A, Van Cauwenberghe C, Cuyvers E, Verheijen J, et al. A 22-single nucleotide polymorphism Alzheimer’s disease risk score correlates with family history, onset age, and cerebrospinal fluid Aβ42. Alzheimer’s Dement. 2015;11:1452–60 Elsevier Inc.
Escott-Price V, Sims R, Bannister C, Harold D, Vronskaya M, Majounie E, et al. Common polygenic variation enhances risk prediction for Alzheimer’s disease. Brain. 2015;138:3673–84 [cited 2019 Nov 29]. Available from: https://academic.oup.com/brain/article-lookup/doi/10.1093/brain/awv268.
Escott-Price V, Myers A, Huentelman M, Shoai M, Hardy J. Polygenic risk score analysis of Alzheimer’s disease in cases without APOE4 or APOE2 alleles. J Prev Alzheimers Dis NLM (Medline). 2019;6:16–9.
Zhang Q, Sidorenko J, Couvy-Duchesne B, Marioni RE, Wright MJ, Goate AM, et al. Risk prediction of late-onset Alzheimer’s disease implies an oligogenic architecture. Nat Commun. 2020;11:1–11 Nature Research.
Leonenko G, Baker E, Stevenson-Hoare J, Sierksma A, Fiers M, Williams J, et al. Identifying individuals with high risk of Alzheimer’s disease using polygenic risk scores. Nat Commun. 2021;12:1–10 Nature Publishing Group; [cited 2021 Aug 5]. Available from: https://www.nature.com/articles/s41467-021-24082-z.
Adamczuk K, De Weer AS, Nelissen N, Chen K, Sleegers K, Bettens K, et al. Polymorphism of brain derived neurotrophic factor influences β amyloid load in cognitively intact apolipoprotein ε4 carriers. NeuroImage: Clin. 2013;2:512–20 Elsevier.
Schaeverbeke JM, Gabel S, Meersmans K, et al. Baseline cognition is the best predictor of 4-year cognitive change in cognitively intact older adults. Alz Res Therapy. 2021;13:75. https://doi.org/10.1186/s13195-021-00798-4.
Koole M, Lewis DM, Buckley C, Nelissen N, Vandenbulcke M, Brooks DJ, et al. Whole-body biodistribution and radiation dosimetry of 18F-GE067: a radioligand for in vivo brain amyloid imaging. J Nucl Med. 2009;50:818–22 Society of Nuclear Medicine.
Vandenberghe R, Van Laere K, Ivanoiu A, Salmon E, Bastin C, Triau E, et al. 18F-flutemetamol amyloid imaging in Alzheimer disease and mild cognitive impairment a phase 2 trial. Ann Neurol. 2010;68:319–29.
Vandenberghe R, Nelissen N, Salmon E, Ivanoiu A, Hasselbalch S, Andersen A, et al. Binary classification of 18F-flutemetamol PET using machine learning: Comparison with visual reads and structural MRI. Neuroimage. 2013;64:517–25.
Adamczuk K, Schaeverbeke J, Nelissen N, Neyens V, Vandenbulcke M, Goffin K, et al. Amyloid imaging in cognitively normal older adults: comparison between 18F-flutemetamol and 11C-Pittsburgh compound B. Eur J Nucl Med Mol Imaging. 2016;43:142–51 springer berlin.
Tzourio-Mazoyer N, Landeau B, Papathanassiou D, Crivello F, Etard O, Delcroix N, et al. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. Neuroimage. 2002;15:273–89.
Klunk WE, Koeppe RA, Price JC, Benzinger TL, Devous MD, Jagust WJ, et al. The Centiloid project: standardizing quantitative amyloid plaque estimation by PET. Alzheimer’s Dement. 2015;11:1–15.e4 Wiley.
De Meyer S, Schaeverbeke JM, Verberk IMW, Gille B, De Schaepdryver M, Luckett ES, et al. Comparison of ELISA- and SIMOA-based quantification of plasma Aβ ratios for early detection of cerebral amyloidosis. Alzheimer’s Res Ther. 2020;12:162 BioMed Central Ltd.
la Joie R, Ayakta N, Seeley WW, Borys E, Boxer AL, DeCarli C, et al. Multisite study of the relationships between antemortem [ 11 C]PIB-PET Centiloid values and postmortem measures of Alzheimer’s disease neuropathology. Alzheimer’s Dement. 2019;15:205–16 NIH Public Access; [cited 2021 Nov 5]. Available from: /pmc/articles/PMC6368897/.
Müller-Gärtner HW, Links JM, Prince JL, Bryan RN, McVeigh E, Leal JP, et al. Measurement of radiotracer concentration in brain gray matter using positron emission tomography: MRI-based correction for partial volume effects. J Cereb Blood Flow Metab. 1992;12:571–83 SAGE Publications Sage UK: London, England.
Buschke H. Selective reminding for analysis of memory and learning. J Verbal Learn Verbal Behav. 1973;12:543–50 Academic Press.
Rosseel Y. Lavaan: An R package for structural equation modeling. J Stat Softw. 2012;48:1–36 American Statistical Association.
van Buuren S, Groothuis-Oudshoorn K. Mice: multivariate imputation by chained equations in R. J Stat Softw. 2011;45:1–67.
Hong S, Prokopenko D, Dobricic V, Kilpert F, Bos I, Vos SJB, et al. Genome-wide association study of Alzheimer’s disease CSF biomarkers in the EMIF-AD Multimodal Biomarker Discovery dataset. Transl Psychiatry. 2020;10:1–12 Nature Publishing Group.
Marees AT, de Kluiver H, Stringer S, Vorspan F, Curis E, Marie-Claire C, et al. A tutorial on conducting genome-wide association studies: quality control and statistical analysis. Int J Methods Psychiatr Res. 2018;27:e1608 Wiley.
1000 Genomes Project Consortium, Auton A, Brooks LD, Durbin RM, Garrison EP, Kang HM, Korbel JO, Marchini JL, McCarthy S, McVean GA, Abecasis GR. A global reference for human genetic variation. Nature. 2015;526(7571):68–74. https://doi.org/10.1038/nature15393.
Das S, Forer L, Schönherr S, Sidore C, Locke AE, Kwong A, et al. Next-generation genotype imputation service and methods. Nat Genet. 2016;48:1284–7 NIH Public Access.
Choi SW, O’Reilly PF. PRSice-2: Polygenic Risk Score software for biobank-scale data. Gigascience. 2019;8:1–6 Oxford University Press.
Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez JC, et al. pROC: An open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinformatics. 2011;12:1–8 BioMed Central.
DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics. 1988;44:837–45.
Poline JB, Worsley KJ, Evans AC, Friston KJ. Combining spatial extent and peak intensity to test for activations in functional imaging. Neuroimage. 1997;5:83–96.
Schwarz CG, Gunter JL, Lowe VJ, Weigand S, Vemuri P, Senjem ML, et al. A comparison of partial volume correction techniques for measuring change in serial amyloid pet SUVR. J Alzheimers Dis. 2019;67:181–95 IOS Press;[cited 2022 Apr 8]. Available from: /pmc/articles/PMC6398556/.
Ge T, Sabuncu MR, Smoller JW, Sperling RA, Mormino EC. Dissociable influences of APOE e4 and polygenic risk of AD dementia on amyloid and cognition. Neurology. 2018;90:E1605–12 American Academy of Neurology.
Leonenko G, Shoai M, Bellou E, Sims R, Williams J, Hardy J, et al. Genetic risk for Alzheimer disease is distinct from genetic risk for amyloid deposition. Ann Neurol. 2019;86:427–35 Wiley-Blackwell.
Mishra S, Blazey TM, Holtzman DM, Cruchaga C, Su Y, Morris JC, et al. Longitudinal brain imaging in preclinical Alzheimer disease: impact of APOE ϵ4 genotype. Brain. 2018;141:1828–39 Oxford University Press.
Sojkova J, Goh J, Bilgel M, Landman B, Yang X, Zhou Y, et al. Voxel-wise relationships between distribution volume ratio and cerebral blood flow: implications for analysis of β-amyloid images. J Nucl Med. 2015;56:1042 NIH Public Access; [cited 2022 Jul 11]. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5367770/.
Heeman F, Yaqub M, Lopes Alves I, Heurling K, Bullich S, Gispert JD, et al. Simulating the effect of cerebral blood flow changes on regional quantification of [18F]flutemetamol and [18F]florbetaben studies. J Cereb Blood Flow Metab. 2021;41:579 SAGE Publications; [cited 2022 Jul 11]. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7907983/.
We would like to thank Carine Schildermans and the staff of Nuclear Medicine and Neurology at the University Hospitals Leuven, in particular Kwinten Porters and Jef Van Loock, Karen Meersmans, Valerie Goovaerts, Astrid Hofkens, and Eva Dries.
This work was supported by the Foundation for Alzheimer Research SAO-FRMA (09013, 11020, 13007); Research Foundation Flanders (G094418N, G0F8516N and JPND-EraNet Triage G0G1519N); KU Leuven (C14/17/108 and C14/21/109); and Vlaams Agentschap voor Innovatie en Onderzoek HBC.2019.2523; 18F-Flutemetamol for the baseline scans was provided by GE Healthcare free of charge for this academic investigator-driven trials.
The project leading to this application has received funding for the Innovative Medicines Initiative 2 Joint Undertaking under grant agreement No 115952 (http://www.imi.europa.eu). This joint undertaking receives the support from the European Union’s Horizon 2020 research and innovation programme and EFPIA. This communication reflects the views of the authors and neither IMI nor the European Union and EFPIA are liable for any use that may be made of the information contained herein.
JS (12Y1620N) is a junior postdoctoral fellow and SDM is a doctoral fellow (11M0522N) of the Flanders Research Foundation (FWO, Fonds voor Wetenschappelijk Onderzoek, Belgium).
Ethics approval and consent to participate
The protocol was approved by the Ethics Committee University Hospitals Leuven. All participants provided written informed consent in accordance with the declaration of Helsinki.
Consent for publication
RV has received research grants from Research Foundation—Flanders (FWO) and KU Leuven, has had a clinical trial agreement for phase 1 and 2 study between University Hospitals Leuven and GE-HC, has received a non-financial support from GE-HC (provision of 18F-Flutemetamol for conduct of investigator-driven trial free of cost), and has a clinical trial agreement (local principal investigator) between University Hospitals Leuven and Biogen, J&J, NovoNordisk, Roche, and UCB.
KVL has received research grants from FWO and KU Leuven, has clinical trial agreements through University Hospitals Leuven/KU Leuven and J&J, Merck, BMS, GE Healthcare, CHDI, Cerveau, Eikonizo, Lundbeck, UCB and Curasen.
ESL, YA, MR, KA, JS, SV, SDM, PD, and IC have no disclosers.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Regression results for PRS with baseline amyloid load. APOE SNPs: rs429358, rs7412. Baseline age, sex and PCs 1-3 included as covariates. N = 90. Supplementary Table 2. Model comparisons between Receiver Operating Characteristic Areas Under the Curve using the DeLong method. Supplementary Table 3. Regression results for the significant PRS with cognitive decline or grey matter atrophy. Raw values are shown. APOE SNPs: rs429358, rs7412. Baseline age, sex and PCs 1-3 included as covariates. N = 90.
About this article
Cite this article
Luckett, E.S., Abakkouy, Y., Reinartz, M. et al. Association of Alzheimer’s disease polygenic risk scores with amyloid accumulation in cognitively intact older adults. Alz Res Therapy 14, 138 (2022). https://doi.org/10.1186/s13195-022-01079-4