Abstract
Alzheimer’s disease typically progresses in stages, which have been defined by the presence of disease-specific biomarkers: Amyloid (A), Tau (T) and neurodegeneration (N). This progression of biomarkers has been condensed into the ATN framework, where each of the biomarkers can be either positive (+) or negative (-). Over the past decades genome wide association studies have implicated about 90 different loci involved with the development of late onset Alzheimer’s disease. Here we investigate whether genetic risk for Alzheimer’s disease contributes equally to the progression in different disease stages or whether it exhibits a stage-dependent effect.
Amyloid (A) and tau (T) status was defined using a combination of available PET and CSF biomarkers in the Alzheimer’s Disease Neuroimaging Initiative cohort. In 312 participants with biomarker-confirmed A-T- status, we employed Cox proportional hazards models to estimate the contribution of APOE and polygenic risk scores (beyond APOE) to convert to A+T- status (65 conversions). Furthermore, we repeated the analysis in 290 participants with A+T- status and investigated the genetic contribution to conversion to A+T+ (45 conversions). Both survival analyses were adjusted for age, sex, and years of education.
For progression from A-T- to A+T-, APOE-e4 burden showed significant effect (HR=2.88; 95% CI: 1.70-4.89; P<0.001), while polygenic risk did not (HR=1.09; 95% CI: 0.84-1.42; P=0.53). Conversely, for the transition from A+T- to A+T+, the APOE-e4 burden contribution was reduced (HR=1.62 95% CI: 1.05-2.51; P=0.031), while the polygenic risk showed an increased contribution (HR=1.73; 95% CI:1.27-2.36; P<0.001). The marginal APOE effect was driven by e4 homozygotes (HR=2.58; 95% CI: 1.05-6.35; P=0.039) as opposed to e4 heterozygotes (HR=1.74; 95% CI: 0.87-3.49; P=0.12).
The genetic risk for late-onset Alzheimer’s disease unfolds in a disease stage-dependent fashion. A better understanding of the interplay between disease stage and genetic risk can lead to a more mechanistic understanding of transition between ATN stages, a better understanding of the molecular processes leading to Alzheimer’s disease as well as opening therapeutic windows for targeted interventions.
Introduction
Alzheimer’s disease is characterized, on the neuropathological level, by the buildup of two proteins: amyloid plaques and neurofibrillary tangles of phosphorylated tau1. Both these pathological features can be observed long before the memory loss and decline in executive function that is characteristic of patients with Alzheimer’s disease. Accumulation of the amyloid plaques in the brain predates the clinical symptoms of Alzheimer’s disease by two decades2, while the spatial distribution of tau tangles reflects more closely the reported cognitive deficits and neurodegeneration3.
The amyloid cascade hypothesis postulates that the deposition of the amyloid-β protein – the main component of the amyloid plaques – is the cause of Alzheimer’s disease, and that neurofibrillary tangles, cell loss, vascular damage and dementia are a direct consequence4. In keeping with the amyloid cascade hypothesis, a theoretical framework for the progression of biomarkers during the course of Alzheimer’s disease has been developed5. Here, amyloid pathology is the first to appear, followed by tau pathology, neurodegeneration and finally cognitive decline. Support for this theoretical framework comes from a number of lines of evidence, including a variety of data-driven modeling approaches based on biomarker data 6–10. In an attempt to operationalize this, a simplified ATN model has been proposed11. The components of the ATN model refer to the status of three different key biomarkers in Alzheimer’s disease: amyloid (A), tau (T) and neurodegeneration (N). In this approach, each of the three biomarkers can either be positive or negative. That is, exceeding the centiloid threshold in amyloid PET imaging would place a participant into the A positive (A+) group, while a scan just below the threshold would be considered amyloid negative (A-). One practical advantage (but also a major source of criticism) of this model is that the biomarker status can be assessed using a variety of methods: wet biomarkers (CSF or plasma) or brain imaging (centiloids or visual reads)12. Although the progression from A-T-N- to A+T-N- to A+T+N- to A+T+N+ would be the most typical progression in Alzheimer’s disease and in keeping with the theoretical framework of biomarker progression, all combinations of biomarker statuses emerge in observational cohorts13,14.
Genome-wide association studies (GWAS) have broadened the understanding of the genetic basis of late-onset Alzheimer’s disease (LOAD) over the last decades15. Initially, these studies were restricted to cases with a clinical diagnosis of Alzheimer’s disease and healthy controls16–18. Recently, these GWAS have been enriched with participants with a family history of Alzheimer’s disease (diagnosis-by-proxy)19–21, leading to a drastic increase in sample size and expanding the set of genetic risk loci for Alzheimer’s disease to 9015. These loci have been linked to various molecular processes such as immunity, cholesterol processing and endocytosis22. Further studies investigated the genetic effects on Alzheimer’s disease-related biomarkers ranging from tau and amyloid levels in CSF23,24, in the brain25,26 to MRI-based measures such as hippocampal volume27 and phenotypes derived from disease progression modeling28. The strongest common genetic risk factor for Alzheimer’s disease is the e4 allele of the APOE gene: carriers of the e4 allele have a 2-4-fold increased risk to develop Alzheimer’s disease, while e4 homozygous subjects have an 8-12-fold increased risk29. The genetic risk outside the APOE region is often summarized using polygenic scores, which have been shown to improve predictions of clinical diagnosis30 and pathology-confirmed cases31. Likewise, the effect of APOE as well as the polygenic risk on various imaging and non-imaging biomarkers have been investigated32–35, with ongoing work suggesting that risk accumulated along different molecular pathways exerts differential effects on different biomarkers in Alzheimer’s disease35–37.
For Alzheimer’s disease and for other disorders, genetic risk is often considered as a time-invariant constant. That is, genetic risk identified through case-controls studies is assumed to affect both onset and progression. However, as Alzheimer’s disease is now understood to unfold in stages, we hypothesized that in line with the amyloid cascade hypothesis and the A/T/N framework that genetic risk in Alzheimer’s disease is disease-stage dependent. That is, some genetic risk factors will aid the transition from A- to A+, while other, distinct genetic risk factors will increase the risk to transition from T- to T+.
In this work we explore whether genetic vulnerability to Alzheimer’s disease varies with disease stage. Using longitudinal data from the Alzheimer’s disease neuroimaging imitative (ADNI) and survival analysis we show that APOE contributes to progression from A-T- to A+T-, but only marginally from A+T- to A+T+. Conversely, polygenic risk contributes to the progression from A+T- to A+T+, but not from A-T- to A+T-.
Materials and methods
Data
Data used in the preparation of this article were obtained from the ADNI database (http://adni.loni.usc.edu). The ADNI was launched in 2003 as a public–private partnership, led by Principal Investigator Michael W. Weiner, MD. The primary goal of the ADNI has been to test whether serial MRI, PET, other biological markers, and clinical and neuropsychological assessment can be combined to measure the progression of mild cognitive impairment (MCI) and early Alzheimer’s disease. For up-to-date information, see www.adni-info.org. ADNI study data were accessed through the R-package ADNIMERGE (accessed: 20th July 2023).
Preparation of genetic data
The genetic data preparation followed the procedure described in Altmann et al.32. The additional genetic data contributed by the ADNI-3 cohort was integrated with the existing data using the same processing pipeline. Briefly, single nucleotide polymorphism (SNP) genotyping data were available for n = 2001 subjects across ADNI phases 1,2, GO and 3. Genotyping was conducted using four different platforms: Human610-Quad, HumanOmniExpress, Omni 2.5 M and Illumina Infinium Global Screening Array v2 (Illumina)38. Prior to imputation, we applied subject-level quality control (QC) steps based on call rate (10% cutoff) and concordance between chip-inferred sex and self-reported sex separately for the four genotyping arrays, all subjects were retained. On SNP-level, we conducted standard QC steps ensuring compatibility with the reference panel used for imputation [strand consistency, allele names, position, Ref/Alt assignments and minor allele frequency (MAF) discrepancy (0.2 cutoff)]. Imputation was carried out using the Sanger Imputation Server (https://imputation.sanger.ac.uk/) with SHAPEIT for phasing39, Positional Burrows-Wheeler Transform40 for imputation and the entire Haplotype Reference Consortium (release 1.1) reference panel41. Data from the four different genotyping platforms were imputed separately. As part of post-imputation QC, multi- allelic variants and SNPs with imputation INFO score <0.3 were removed and genotype calls with posterior probability <0.9 were set to missing (i.e., hard called). Following the initial QC, genotypes from the four platforms were merged. Additional information on the imputation and QC process is detailed in Scelsi et al.28 and https://rpubs.com/maffleur/452627. Using the merged data, we retained SNPs with MAF ≥1% and genotyping rate >0.9.
SNPweights42 was used to infer genetic ancestry from genotyped SNPs. The reference panel comprised Central European, Yoruba Africans and East Asian samples from HapMap 343 and native Americans from Reich et al.44. For this study, only participants with predicted central European ancestry of 80% or more were retained (n=1851). Next, using the imputed and merged data genetic relatedness between central European participants was computed. First, the SNP content was restricted to SNPs with MAF ≥5% and linkage disequilibrium (LD)- pruning was carried out in PLINK v1.9 (–indep-pairwise 1000 50 0.1). The genetic relatedness matrix was computed using the remaining autosomal SNPs and the dataset was trimmed to remove subjects with relatedness >0.1 (–rel-cutoff 0.1), leading to n=1833 unrelated participants.
Definition of genetic risk
In this study, we focused on two sources of genetic risk: (1) the risk conferred through the APOE gene based on the genetic markers for APOE-ε2 and APOE-ε4 and (2) the polygenic risk conferred by the remaining genome. As described previously32, PRSs were computed using the software PRSice v2.1.945. As base GWAS the stage 1 results of the Alzheimer’s disease GWAS featuring a clinically defined Alzheimer’s disease phenotype was used18. For PRS computation, SNPs with MAF ≥5% were considered and SNPs were selected using LD- clumping (1000 KB, R2 of 0.1 and P-value threshold of 1.0) within the ADNI cohort, missing SNPs were simply ignored on subject level (using the setting–missing SET_ZERO) and the APOE region was excluded (hg19 coordinates chr19 from 44,400,000 to 46,500,000). For this study, we used only the P-value cutoff of 1.0e-8 to build the PRS (Supplementary Table 1). PRSs were computed for all ADNI participants with genome-wide genotyping data. Of the remaining subjects, n=417 ADNI-1 participants who contributed to the Alzheimer’s disease GWAS18 were excluded from the analysis to ensure independence between training and application dataset for PRS. Thus, n=1416 unrelated participants with central European ancestry were eligible for inclusion into the study.
Definition of amyloid status
For this project we relied on two modalities to define amyloid status: Amyloid-β PET using the 18F-florbetapir (FBP) and the 18F-florbetaben (FBB) PET tracer, and CSF measures of amyloid-β (1-42) using the Roche Elecsys® immunoassay. We used processed data by ADNI for both modalities. Detailed information on the PET processing is available elsewhere (http://adni-info.org)46; and information on CSF amyloid-β(1-42) processing is detailed in47,48. For CSF amyloid-β(1-42) we used the cutoff of 880 pg/mL48 and for amyloid-β PET we used the tracer-specific SUVR (whole cerebellum reference) cutoffs of 1.11 and 1.08 for FBP and FBB49, respectively.
A participant’s visit was labelled as A+ if either the PET result or the CSF result indicated a positive amyloid finding (i.e., in cases where the results were discordant, the visit would be labelled as A+). Visits with only negative amyloid findings were labeled as A- and visits without any information on amyloid (i.e., neither a PET nor a CSF result) were labelled as ‘Amyloid missing’.
Definition of tau status
In keeping with our definition of amyloid positivity we used available data from CSF and PET imaging to define the tau status. More precisely, we used the Phospho-Tau(181P) CSF Roche Elecsys® immunoassay and PET imaging using the 18F-Flortaucipir tracer. Details on the processing are available elsewhere (http://adni-info.org)47,48,50. For CSF we used phosphorylated 181P tau (pTau) with a cutoff of 34.61 pg/mL and for tau PET we used the cutoff of 1.42 in the meta temporal ROI comprising amygdala, entorhinal cortex, fusiform gyrus, inferior and middle temporal gyri51 when normalized to the inferior cerebellar gray matter52. Both cut-offs were data-driven: (1) the tau PET cutoff corresponds to a z-score of 2.0 in the cognitively normal participants in ADNI (n=506) and is close to the ‘high tau’ cutoff of 1.43 defined in by Jack Jr et al. 53; (2) the pTau cutoff was set to maximize the Youden’s index between CSF pTau and tau PET (at the 1.42 cutoff) in (n=502) ADNI participants with concurrent CSF and PET measurements. The same labeling scheme as for amyloid was applied: a visit was labelled as T+ if either the CSF or the PET indicated a positive finding, T- if there were only negative tau findings and ‘Tau missing’ if neither data on CSF pTau nor on tau PET were available.
Statistical modeling
We used cox-proportional hazards models to investigate the genetic contribution of progressing from (1) A-T- to A+T- and from (2) A+T- to A+T+. For (1) we included every eligible participant with genetic data and who was A-T- based on their biomarker results as described above. This earliest A-T- visit was considered the ‘start’ visit (i.e., the status of previous visits was ignored). A subject was considered a converter when the available biomarkers indicated A+T- at a later visit. The time of the first A+T- biomarker finding after their initial A-T- visit was used as the conversion time. For non-converters we recorded the last visit where both amyloid and tau biomarker information was available to define the maximum follow-up time. Similarly, for (2) we included every participant with genetic data and who was A+T-. This first A+T- visit was considered the ‘start’ visit. A subject was considered a converter when the available biomarkers indicated A+T+ at a later visit. The time of the first A+T+ biomarker finding after their initial A+T- visit was used as the conversion time. As before, for non- converters we recorded the last visit where amyloid and tau biomarker information was available as their last point of contact.
For the analyses, the cox-proportional hazards model included age at the inclusion visit (i.e., either the A-T- or the A+T- visit), sex, and education. As variables of interest, we further included genetic variables for PRS as well as allele counts of APOE-e2 and APOE-e4. The proportional hazards assumption was tested for each covariate and for the global model using statistical tests and graphical diagnostics based on scaled Schoenfeld residuals. As a measure of overall model performance, the concordance index (C-index) was computed for the (i) full models (ii) models without APOE, (iii) models without PRS, and (iv) models without any genetics. Statistical tests were carried out in R (v 4.1.0) using the survival (v 3.5.5), survminer (v 0.4.9), and rms (v 6.8.0) packages.
Sensitivity analyses
In addition to the two main analyses, we conducted a series of sensitivity analyses addressing the conversion definition, biomarker cutoffs, biomarker source, and polygenic score source.
Conversion definition
To maximize the available data, we relaxed the requirement for both amyloid and tau biomarker to be available at the same visit to define conversion. As before, both biomarkers at the same time were required to define the ‘start’ visit as either A-T- or A+T-. However, for defining progression from A-T- or A+T-, a single A+ or T+ visit was sufficient, respectively.
Biomarker cutoffs
We varied the tau PET cutoff from 1.0 to 3.0 standard deviations above the mean in the cognitively normal ADNI participants. Notably, the lowest cutoff resulted in 1.31 and was close to the neuropathologically defined cutoff of 1.29 by Lowe et al.54. The pTau cutoffs were adjusted accordingly to maximize Youden’s index between CSF pTau and tau PET.
Biomarker source
To maximize data and follow-up time, the main analyses combined data from two biomarker sources: CSF and PET. Additional sensitivity analyses relied exclusively on either CSF biomarkers or PET biomarkers. For this analysis pTau cutoffs were varied in the range from 22 to 31, covering the values of 24.25 and 29.19, which were found to indicate tau PET positivity in Braak III/IV and Braak V/VI regions, respectively55. Tau PET cutoffs were varied, as before, from 1.0 to 3.0 standard deviations above the mean in cognitively normal participants (i.e., in the same range as the main analysis).
Polygenic Score source
To include more SNPs into the PRS, we explored the summary statistics on Alzheimer’s Disease and related dementias by Bellenguez et al.21. We followed the same PRS pipeline as above and applied a p-value threshold of 5.0e-08 (i.e., genome-wide significant) leading to 77 included SNPs (Supplementary Table 2).
Results
Out of 16,401 visits recorded in ADNI, 3789 (23.1%) visits had both amyloid and tau biomarker data available. Concordance between PET- and CSF-based assessment was 80% and 81% for amyloid and tau, respectively. For both survival models we identified around 300 subjects in the ADNI database with both biomarkers available (Table 1): 312 individuals were A-T- of whom 65 converted to A+T-. The mean age of the A-T- group was 71.3 (6.65) years and there was an almost equal number of males and females (49.4% females). n=290 individuals were A+T- at any stage, of whom 45 converted to A+T+. The mean age of the A+T- group was 73.2 (6.9), significantly older than the A-T- group (unpaired t-test: T=3.50, P<0.001). The fraction of females in that cohort was lower compared to the A-T- group (43.8% females), but not at a significant level (chi-square-test; chi-square=1.65; df=1; P=0.19). The distribution of APOE-e4 alleles (chi-square test; chi-square= 81.1; df=2; P=2.4e-18) and APOE-e2 alleles (Fisher’s exact test; P= 0.0004) differed significantly between A-T- and A+T. There were more APOE-e2 carriers and fewer APOE-e4 carriers in the A-T- group than in the A+T- group. The PRS did not differ between the A-T- and A+T- groups (two-sided t-test; t=0.84; df=600; P=0.39). The cohort with the relaxed conversion criterion showed comparable characteristics (Supplementary Table 3).
APOE-e4 influences progression from A-T- to A+T-
The survival analysis showed a significant contribution by APOE-e4 allele count (HR=2.88; 1.70-4.89; P=8.7e-05) but no significant contribution by PRS (HR=1.09; 0.84-1.42; P=0.53) (Fig. 1). The APOE-e2 allele count directionally favoured a protective effect, but this was not significant (HR=0.73; 0.28-1.86; P=0.51). The Cox proportional hazards assumption held for this model (Global Schoenfeld test P=0.24) (Supplementary Fig. 1). The C-indices align with this pattern: full model (C=0.612), no APOE (C=0.525), no PRS (C=0.611), no genetics (C=0.531). This pattern of associations of APOE and PRS with progression was largely independent of the tau PET and pTau thresholds (Supplementary Table 4). Furthermore, using the more relaxed conversion criteria (i.e., confirmed A+ status was sufficient instead of a confirmation of A+T-) yielded more conversions (85 instead of 65) but qualitatively the same result (Supplementary Fig. 2), i.e., a significant contribution by APOE-e4 allele count (HR=3.34; 2.14-5.22; P=1.0e=07) but not by PRS (HR=1.06; 0.84-1.34; P=0.61).
Polygenic risk affects the progression from A+T- to A+T+
The survival analysis showed a marginally significant contribution by APOE-e4 burden (HR=1.62; 1.05-2.51; P=0.031), which was mainly driven by APOE-e4 homozygotes (HR=2.58; 1.05-6.35; P=0.039) rather than APOE-e4 heterozygotes (HR=1.74; 0.87-3.49; P=0.12) (Fig. 2). Furthermore, there was a significant contribution by PRS (HR=1.72; 1.27- 2.36; P=0.00057). The APOE-e2 allele count again was directionally consistent with a protective effect, but this was not significant (HR=0.41; 0.09-1.78; P=0.23). The Cox proportional hazards assumption held for this model (Supplementary Fig. 3). The C-indices drop marginally when either APOE or PRS are removed from the model: full model (C=0.657), no APOE (C=0.634), no PRS (C=0.615), no genetics (C=0.549). Moreover, the association pattern of APOE and PRS with progression was largely independent of the tau PET and pTau cutoffs (Supplementary Table 4). In addition, education showed a marginally protective association with conversion to A+T+ (HR=0.89; 0.80-0.99; P=0.039). Applying the more relaxed conversion criteria (i.e., confirmed T+ status was sufficient instead of a confirmation of A+T+) yielded more conversions (51 instead of 45) but qualitatively the same result (Supplementary Figure 5): a significant contribution by PRS (HR=1.62; 1.22-2.17; P=0.001) and a marginal contribution by APOE-e4 allele burden (HR=1.56; 1.04-2.35; P=0.031).
Results are independent of PRS source
Using an the alternative PRS with 77 SNPs led to the same observation of a significant effect by APOE-e4 burden on A-T- to A+T- conversion (HR=2.84; 1.68-4.82; P=0.0001) and a lack of contribution by the PRS (HR=0.97; 0.77-1.2; P=0.80; Supplementary Fig. 5A). Conversely, for A+T- to A+T+ conversion, the contribution by APOE-e4 was reduced (HR=1.71; 1.09-2.68; P=0.019), where the PRS exhibited a strong contribution (HR=1.59; 1.19-2.13; P=0.00163; Supplementary Fig. 5B).
Results are independent of biomarker source
The main analysis combined different biomarker sources to maximize the available data and the observation time for the conversion analysis. Just relying on PET biomarkers alone, only 7.5% of visits (1237 of 16,401) had both biomarkers and it resulted in shorter observation times for A-T- to A+T- conversion (3.62 (sd=1.07) years) and A+T- to A+T+ conversion (2.94 (sd=1.13) years) compared to the main analysis. Relying solely on CSF biomarkers, 19.2% of visits (3155 of 16,401) had both biomarkers. Overall, this led to shorter observation time with respect to the main analysis for A-T- to A+T- conversion (3.94 (sd=2.7) years) and A+T- to A+T+ conversions (3.6 (sd=2.38) years). Furthermore, relying on a single source of biomarkers led to reduced sample sizes (from around 300 in the main analyses to 70-250 in the sensitivity analyses) and observed conversions (Supplementary Table 5 and Supplementary Table 6). However, despite the reduced statistical power these sensitivity analyses confirmed the pattern observed in the main analysis: APOE-e4 contributed mainly to the A-T- to A+T- conversion, while PRS contributed to the A+T- to A+T+ conversion (Fig. 3, Supplementary Table 5, Supplementary Table 6).
Discussion
This longitudinal survival analysis demonstrated that APOE-e4 plays an important role for the progression from A-T- to A+T-, but APOE is only of marginal importance in A+T- participants who progress to A+T+. Conversely, polygenic risk for AD exhibited the inverse pattern: there was no contribution to the progression from A-T- to A+T-, but a significant contribution to faster progression from A+T- to A+T+. This held true for an alternative PRS defined using a different genetic study and involving a larger number of genetic loci. Notably, when assessing covariates, a differential effect of years of education was observed: higher education had no effect (Fig. 1) or was marginally harmful (Supplementary Fig. 2) for the progression from A- T- to A+T-, but was protective for the conversion from A+T- to A+T+ (Fig. 2). The bisection of the genetic risk by disease stage was largely independent of the applied biomarker cutoffs. Moreover, relying on only a single source of biomarkers for defining stage and conversion confirmed the findings of the main analysis despite reduced sample size and observation time.
The finding of a stronger effect of APOE-e4 earlier in the disease progress may explain the observation of stronger genetic effect of APOE-e4 on Alzheimer’s disease in the group of 60-80 year-old people compared to people 80 years and older56–58. Amyloid deposition occurs 10+ years before other Alzheimer’s disease processes2 and APOE-e4 is the strongest common genetic risk factor for amyloid deposition, then it would be expected for APOE-e4 to exert its maximum effect in younger people. Still, APOE remains the strongest risk factor in individuals 80 years and older57. Thus, the age-dependent heterogeneity of APOE is likely compounded by a survivor bias: individuals with a very late onset despite carrying APOE-e4 may harbor protective variants59 such as KLOTHO-VS, where a protective effect on amyloid deposition and Alzheimer’s disease was only observed in 60-80 year-olds, but not in the 80+ group56,60.
The findings from this longitudinal analysis are also in line with previous reports in the ADNI cohort of APOE and polygenic risk on amyloid and tau. For instance, cross-sectional amyloid biomarkers in the CSF and in the brain were mainly driven by APOE, while cross- sectional CSF tTau and pTau measurements were associated with PRS beyond the APOE locus32. Moreover, APOE was found to predict amyloid status while polygenic risk for Alzheimer’s disease improved predictions of diagnosis as well as clinical progression from MCI to Alzheimer’s disease above APOE alone35. These observations extend to plasma markers of tau pathology: PRS (that excluded the APOE region) were found to be only associated with plasma p-tau181 in A+ participants61. This association between polygenic risk (beyond APOE) and CSF tau biomarkers rather than with amyloid and neurodegeneration was also observed outside of the ADNI study62. Polygenic risk (beyond APOE) was associated with non-amyloid endophenotypes in a large cohort of people with MCI. This suggests that these variants are more closely linked with neuronal degeneration than with Alzheimer’s disease- related amyloid pathology63. All these previous studies made the connection between existing amyloid pathology and correlations between polygenic risk and tau pathology using cross- sectional study designs. In a recent longitudinal study of tau PET in the ADNI cohort, higher polygenic risk for Alzheimer’s disease was associated with accelerated increase in tau signal in brain and this effect was modulated by amyloid pathology: A+ participants showed a stronger effect of PRS on tau accumulation64. Our longitudinal analysis, which combined CSF and PET data to maximize the sample size, confirms these observations and indicates that polygenic Alzheimer’s disease risk (outside the APOE region) contributes to tau pathology in A+ participants, but has no meaningful contribution in A- participants. The genetic architecture of Alzheimer’s disease, as captured by the PRS, involves multiple different pathways, mainly β-amyloid processing, tau, immunity, and lipid processing18,21. The employed PRS covers three genes that have been associated with tau binding: BIN1, CLU and PICALM. BIN1 mediates Alzheimer’s disease risk by modulating tau pathology65, and BIN1 risk variants increase tau PET (but not amyloid PET)66 in an amyloid-dependent fashion67. Our observations of an amyloid-dependent effect of the PRS align with these earlier single gene studies. While there are currently no mechanistic analyses that explain the amyloid-dependent effect of BIN1 on tau pathology, recent data from animal models in Alzheimer’s disease suggest a state-dependent effect of genetic risk factors related to microglia: Deletion of Trem2 in mouse models exacerbated tau accumulation and spreading leading to brain atrophy, but only in the presence of existing β-amyloid pathology68. Along the same line, physical contact between microglia and plaques as well as a functioning TREM2 gene are necessary for the appropriate microglia response to amyloid pathology69. Thus, defects in TREM2 can only contribute to neurodegeneration once amyloid pathology has been established. Consequently, other genes contributing to the PRS may also exert their effect in an amyloid-dependent fashion. Our longitudinal analysis presented here is the first to support such a state-dependent genetic risk model in humans, and further fine-grained examination of how the pathways involved in the PRS contribute to sequential disease progression are needed.
The partition of Alzheimer’s disease genetic risk into APOE-related and polygenic risk beyond APOE is a simplification in this analysis. Recent works have linked established Alzheimer’s disease risk loci outside APOE, such as CR1, to amyloid biomarker levels23. Conversely, studies of biomarker levels of tau repeatedly highlight the APOE locus23,24. However, if being A+ were a pre-requisite to exhibit pathological accumulation of tau, then the strong genetic association with APOE in these GWASs would merely reflect the necessary condition rather than a genuine direct molecular process that affects tau levels. The strong dependency of tau levels on established amyloid pathology is supported by mediation analyses in recent cross-sectional62 and longitudinal64 studies. Moreover, the known Alzheimer’s disease genetic risk variants are contributing differently to molecular pathways18,21,22, where each pathway in turn will exercise differential effects on the AD biomarkers including markers for vascular pathology37. Therefore, a pathway PRS only comprising genes associated with the regulation of the amyloid precursor protein catabolic process (e.g., Gene Ontology term GO:1902991) may contribute significantly to the conversion from A-T- to A+T-. Likewise, a pathway PRS using only genes known to bind the tau protein (GO:0048156) may exhibit an even stronger association with A+T- to A+T+ conversion than the general PRS employed here.
In addition to the cascading effect of genetic risk in Alzheimer’s disease, we also observed a stage-dependent effect of non-genetic risk factors. Education has been shown to have a protective effect against dementia70: here we show that higher rates of education do not influence transition to amyloid positivity, but do slow progression from A+T- to A+T+. Consequently, other non-genetic risk factors may show a similarly state-dependent effect on the pathological pathway from A-T- to A+T+ and further neurodegeneration.
The study has several limitations. Firstly, the current analysis was limited to just two biomarkers in Alzheimer’s disease: β-amyloid and tau. It would be desirable to include neurodegeneration (N; of the ATN framework) or potentially more fine-grade staging from advanced data-driven disease progression modeling7. However, at this point adding further stages would reduce the available sample size. Secondly, we partitioned the genetic risk in Alzheimer’s disease into two components: APOE and other top variants combined into a single polygenic risk score. Further work should explore a more fine-grained partition of the polygenic risk into individual SNPs or into pathway-PRS37. Thirdly, the study population was of central European ancestry, therefore it is unclear whether the findings would generalize to other genetic backgrounds. Finally, although the ADNI cohort is a large cohort, the number of subjects who were eligible for our analysis was reduced due to the requirement of concordant and longitudinal recordings of multiple biomarkers as well as genetics. The available sample size may have limited statistical power to render the estimated hazards ratios significant in some settings. However, the two conversion analyses were based on similar sample sizes (around 300 participants), thus allowing us to make a relative comparison between the genetic effect (of APOE or PRS) under two different, biomarker defined, disease stages. Moreover, uncertainty of the estimated effects may also be increased due to disease heterogeneity in Alzheimer’s disease71–74, which is likely underpinned by differences in genetic architecture. Thus, analyses in further large longitudinal cohorts are required to confirm the observation of stage-dependent genetic vulnerability in Alzheimer’s disease and to uncover more fine-grained associations with Alzheimer’s disease subtypes.
In this work we demonstrated, in a simplified setting, that genetic risk for late onset Alzheimer’s disease unfolds in a disease stage-dependent fashion. A better understanding of the interplay between disease stage and genetic risk can lead to a better understanding of the molecular processes leading to Alzheimer’s disease as well as opening therapeutic windows for targeted interventions and personalized approaches to dementia prevention.
Data Availability Statement
Data used in the preparation of this article were obtained from the ADNI database (http://adni.loni.usc.edu) and are freely available after registration. Analysis scripts are available at https://github.com/andrealtmann/cascading_genetic_risk.
Funding
This study was supported by the Early Detection of Alzheimer’s Disease Subtypes (E-DADS) project, an EU Joint Programme - Neurodegenerative Disease Research (JPND) project (see www.jpnd.eu). The project is supported under the aegis of JPND through the following funding organizations: United Kingdom, Medical Research Council (MR/T046422/1); Netherlands, ZonMW (733051106); France, Agence Nationale de la Recherche (ANR-19- JPW2–000); Italy, Italian Ministry of Health (MoH); Australia, National Health & Medical Research Council (1191535); Hungary, National Research, Development and Innovation Office (2019–2.1.7-ERA-NET-2020–00008). N.P.O. is a UKRI Future Leaders Fellow (MR/S03546X/1). JMS acknowledges the support of the National Institute for Health Research University College London Hospitals Biomedical Research Centre, Wolfson Foundation, Alzheimer’s Research UK, Brain Research UK, Weston Brain Institute, Medical Research Council, British Heart Foundation, and Alzheimer’s Association. This research was funded in whole, or in part, by the Wellcome Trust [227341/Z/23/Z]. For the purpose of open access, the author has applied a CC BY-ND public copyright license to any Author Accepted Manuscript version arising from this submission.
Competing interests
The authors report no competing interests.
Supplementary material
Supplementary material is available at Brain online.
Data Availability
Data used in the preparation of this article were obtained from the ADNI database (http://adni.loni.usc.edu) and are freely available after registration.
Supplementary Material
Footnotes
↵* Data used in preparation of this article were obtained from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database (adni.loni.usc.edu). As such, the investigators within the ADNI contributed to the design and implementation of ADNI and/or provided data but did not participate in the analysis or the writing of this report. A complete list of ADNI investigators can be found at: http://adni.loni.usc.edu/wpcontent/uploads/how_to_apply/ADNI_Acknowledgement_List.pdf.
(1) added sensitivity analysis using a PRS with more SNPs (2) extended discussion