Skip to main content

Coping with brain amyloid: genetic heterogeneity and cognitive resilience to Alzheimer’s pathophysiology


Although abnormal accumulation of amyloid in the brain is an early biomarker of Alzheimer’s disease (AD), wide variation in cognitive trajectories during life can be seen in the setting of brain amyloidosis, ranging from maintenance of normal function to progression to dementia. It is widely presumed that cognitive resilience (i.e., coping) to amyloidosis may be influenced by environmental, lifestyle, and inherited factors, but relatively little in specifics is known about this architecture. Here, we leveraged multimodal longitudinal data from a large, population-based sample of older adults to discover genetic factors associated with differential cognitive resilience to brain amyloidosis determined by positron emission tomography (PET). Among amyloid-PET positive older adults, the AD risk allele APOE ɛ4 was associated with worse longitudinal memory trajectories as expected, and was thus covaried in the main analyses. Through a genome-wide association study (GWAS), we uncovered a novel association with cognitive resilience on chromosome 8 at the MTMR7/CNOT7/ZDHHC2/VPS37A locus (p = 4.66 × 10–8, β = 0.23), and demonstrated replication in an independent cohort. Post-hoc analyses confirmed this association as specific to the setting of elevated amyloid burden and not explained by differences in tau deposition or cerebrovascular disease. Complementary gene-based analyses and publically available functional data suggested that the causative variant at this locus may tag CNOT7 (CCR4-NOT Transcription Complex Subunit 7), a gene linked to synaptic plasticity and hippocampal-dependent learning and memory. Pathways related to cell adhesion and immune system activation displayed enrichment of association in the GWAS. Our findings, resulting from a unique study design, support the hypothesis that genetic heterogeneity is one of the factors that explains differential cognitive resilience to brain amyloidosis. Further characterization of the underlying biological mechanisms influencing cognitive resilience may facilitate improved prognostic counseling, therapeutic application, and trial enrollment in AD.


Observational studies support a dynamic model of Alzheimer’s disease (AD), in which amyloidosis is an early event that is eventually followed by other biomarker abnormalities and clinical impairment [24, 27]. However, a considerable fraction of older adults display abnormal accumulation of amyloid in the absence of overt cognitive impairment, with the frequency of this finding modified by age, APOE (apolipoprotein E) ɛ4 status, and sex, among other factors [7, 28, 30]. This observation highlights the importance of elucidating the underlying biological mechanisms of resilience to (i.e., coping with) amyloid pathology within frameworks attempting to forecast cognitive outcomes in AD [3].

It is widely presumed that resilience may be influenced by the environment and lifestyle [41] in addition to inherited factors [23]. The recent discovery of the APOE (apolipoprotein E) Christchurch mutation as protective in the face of high amyloid burden and an AD-causative PSEN1 (presenilin 1) mutation represents a germane exemplar for this concept [1]. However, there is a relative paucity of literature on genetic resilience factors specific to AD. Analyses of general cognitive decline in the wider setting of older adults have most consistently implicated variants in the APOE region on chromosome 19 [13,14,15, 70]. A few studies have examined putative downstream effects in the setting of AD pathology for specific candidate loci, including the APOE ɛ4 [46], KL (klotho) VS [6], and BDNF (brain-derived neurotrophic factor) Val66Met [19] alleles, among others. However, much of the presumed genetic architecture underlying resilience to AD pathology is still unaccounted for by known candidate genes.

The validation and expansion of positron emission tomography (PET) imaging biomarkers creates the opportunity to noninvasively assess in vivo AD pathology in large samples conducive for discovery-oriented genetic analyses of resilience. In this study, we hypothesized that novel gene variants and biological pathways would be associated with differential cognitive resilience to amyloidosis. We tested this hypothesis by conducting a genome-wide association study (GWAS) of longitudinal memory performance in a large, population-based sample of amyloid-PET positive older adults.


Selection of participants

The discovery cohort in this study was drawn from the Mayo Clinic Study of Aging (MCSA), a population-based prospective study of older adults residing in Olmsted County, Minnesota [49, 54]. Individuals were identified for recruitment using the Rochester Epidemiology Project (REP) medical records linkage system [56, 61]. Enrollment began in 2004 for individuals 70–89 years of age, and the study was subsequently extended to include those aged 50 and older (2012) and 30 and older (2015). Clinical data through questionnaires and in-person history, neuropsychological assessment, multimodal neuroimaging, and laboratory tests were assessed at selected visits based on study protocols. Clinical diagnoses were made by a consensus panel, incorporating all available information. All study protocols were approved by the Mayo Clinic and Olmsted Medical Center Institutional Review Boards. Written informed consent was obtained from all participants or their surrogates.

The replication cohort analyzed separately in this study was drawn from the Alzheimer’s Disease Neuroimaging Initiative (ADNI), a longitudinal multicenter study launched in 2004 as a public–private partnership [63, 68]. The goal of the ADNI is to facilitate development of clinical, imaging, genetic, and biochemical biomarkers for the early detection and tracking of AD. Individuals 55–90 years of age were recruited from over 50 sites across the United States and Canada, and initially followed at 6–12-month intervals for 2–3 years. Subsequent study phases have extended follow-up for existing participants and have enrolled additional individuals. All participants provided written informed consent and study protocols were approved by each site’s institutional review board. Further information about the ADNI can be found at

Inclusion criteria for this study included the following: age 50 years or older, the presence of amyloid positivity by PET, at least 2 subsequent time points with neuropsychological assessment data, and genome-wide single nucleotide polymorphism (SNP) genotype data. This resulted in 546 individuals in the MCSA discovery cohort and 545 individuals in the ADNI replication cohort. A separate, non-overlapping group of 953 amyloid PET negative individuals from the MCSA was analyzed post-hoc for comparison.

Demographic and clinical data

Age at the time of neuroimaging, sex, and years of education were ascertained for each patient. In the discovery cohort, a measure of cerebrovascular disease risk (CMC) was ascertained from health care records as a summation of the presence or absence of hypertension, hyperlipidemia, cardiac arrhythmias, coronary artery disease, congestive heart failure, diabetes, and stroke [67].

Longitudinal cognitive data

For each applicable participant visit in the MCSA discovery cohort, a composite memory domain z-score was generated as described previously [34, 55], based on the delayed recall tasks of the Wechsler Memory Scale-Revised Logical Memory II, Wechsler Memory Scale-Revised Visual Reproduction II, and Auditory Verbal Learning Test from cognitively unimpaired participants aged 50 and older and weighted back to the Olmsted County population. For each subject, linear regression on these memory z-scores was used to generate a subject-specific intercept and slope. The variability in slope across all of the subjects captures the variability in longitudinal trajectories of memory functioning across all of the subjects, and was used as our primary outcome measure. The intercept for each subject was used as a covariate in analyses, representing an estimate of early-age. For the ADNI replication cohort, the 13-item Alzheimer's Disease Assessment Scale-Cognitive (ADAS-Cog) Subscale was used as the analogous measure of interest based on previous work employing longitudinal cognitive data in ADNI [2, 39]. Linear regression was similarly used to generate a subject-specific intercept and slope for each participant based on longitudinal ADAS-Cog scores for the replication cohort, with the variability in slope across all subjects again used as the primary outcome.

Genetic data

Peripheral blood samples were acquired at the baseline visit for 1783 MCSA participants. Genomic DNA extracted from these samples was used for genotyping of 658,805 SNPs via the Illumina Infinium Global Screening Array-24 v2.0. Standard SNP-level quality control (QC) filters were applied using PLINK version 1.9 [9, 51], including call rate ≥ 95%, Hardy–Weinberg Equilibrium p ≥ 1 × 10–5, and minor allele frequency (MAF) ≥ 1%. Subject-level QC filters included call rate ≥ 98%, sex checks versus clinical data, Caucasian ancestry determined through STRUCTURE version 2.3.4, and ensuring no cryptic first- or second-degree relatedness (PLINK identity by descent PI_HAT < 0.25). Following these procedures, GWAS array data passing QC was available for 506,136 SNPs and 1727 MCSA participants. APOE ɛ2/ɛ3/ɛ4 allele status determined through the GWAS array (via genotyping of rs429358 and rs7412) displayed 100% concordance with results from standard restriction digest methods [20]. As a conservative measure to account for any potential confounding effects of population stratification, principal component analysis of the genotype data with SNPRelate [71] was used to generate eigenvectors for use as covariates. Genome-wide imputation was performed with the Michigan Imputation Server [12] using Minimac version 4-–1.0.2 and the Haplotype Reference Consortium reference panel [43]. Following additional post-imputation QC filters including SNP call rate ≥ 95%, sample call rate ≥ 98%, Hardy–Weinberg Equilibrium p ≥ 1 × 10–6, MAF ≥ 1%, stringent imputation quality measure (r2) ≥ 0.8, and removal of SNPs with no or duplicate identifying rs number, data was available for 6,153,814 SNPs and 1727 MCSA participants. For this study, to ensure appropriate power for analyzed variants given the sample size in the discovery cohort, an additional MAF filter of 5% was applied, leaving 4,456,454 SNPs for analysis in the GWAS.

ADNI-1, ADNI-GO, and ADNI-2 participants were genotyped on one of three Illumina GWAS arrays as described previously [58]. Processed, post-QC genotype data files were downloaded from the ADNI database ( Imputation was performed within groups based on the genotyping array utilized, and then the independently imputed datasets were merged using PLINK. Imputation methods and post-imputation QC were performed as in the MCSA dataset and resulted in 5,599,642 SNPs for 1662 unique ADNI participants. In the merged sample, a total of 18 individuals (representing 6 pairs and 2 trios) were found to have cryptic relatedness (PI_HAT ≥ 0.25) based on identity by descent analysis using common (MAF ≥ 5%) SNPs roughly pruned for LD (r2 < 0.6). After one individual from each pair or trio was randomly selected for retention, data was available for analysis for 1652 unique ADNI participants. As with the MCSA dataset, the first 5 principal component eigenvectors were used as covariates in genetic analyses.

Neuroimaging data

In the MCSA, amyloid PET and tau PET scans were acquired and analyzed using an in-house fully-automated image processing pipeline as described in detail elsewhere [26]. Amyloid PET imaging was performed with Pittsburgh compound B (PiB) [33] and tau PET was performed with AV-1451 (18F-flortaucipir), synthesized on site with precursor supplied by Avid Radiopharmaceuticals [42]. Standardized uptake value ratio (SUVR) measures for amyloid and tau PET were generated by normalizing median tracer uptake in target regions of interest (ROIs) to the cerebellar crus grey matter. The target amyloid PET measure was global cortical amyloid load, computed from the prefrontal, orbitofrontal, parietal, temporal, anterior cingulate, and posterior cingulate/precuneus ROIs. Amyloid PET positivity was defined by global cortical SUVR ≥ 1.48 as previously described [42]. The target tau PET measure was a meta-ROI computed from the entorhinal, amygdala, parahippocampal, fusiform, and inferior and middle temporal ROIs [26]. Tau PET positivity was defined by meta-ROI SUVR ≥ 1.25 as previously described [25]. Entorhinal cortex tau PET burden was also assessed.

MRI for MCSA participants was performed on 3T MRI systems (General Electric Healthcare, Waukesha, WI). Target MRI biomarkers of neurodegeneration were derived using FreeSurfer version 5.3 and included hippocampal volume and the cortical thickness in an AD-signature meta-ROI comprised of the entorhinal cortex and inferior and middle temporal and fusiform gyri [59]. Target MRI biomarkers of cerebrovascular disease included white matter hyperintensities (WMHs) from FLAIR MRI [17] and fractional anisotropy (FA) of the genu of the corpus callosum from diffusion tensor imaging (DTI) MRI [66].

In the ADNI, amyloid PET was performed with 18F-florbetapir using acquisition and processing protocols as described at Amyloid PET summary measures generated at the University of California, Berkeley [29] were downloaded from the ADNI database ( The target amyloid PET measure was global cortical amyloid load, assessed from FreeSurfer-defined regions of interest including the anterior and posterior cingulate, frontal, lateral parietal, and lateral temporal cortices, and normalized to a whole cerebellum reference region [39]. Amyloid PET positivity was defined by global SUVR ≥ 1.11, as previously described [38].

Statistical analysis

Genetic association tests with the memory slope phenotype were performed using PLINK version 1.9. Prior to the genome-wide scan in the discovery dataset, the APOE ɛ4 and ɛ2 alleles were analyzed for association signal given their well-validated relationship with AD and age-related cognitive trajectories. These candidate analyses and the GWAS were conducted on the discovery sample using linear regression under an additive genetic model and including age at the time of neuroimaging, sex, years of education, memory domain intercept (as a measure of premorbid performance), and the first 5 genetic principal component eigenvectors as covariates. For the GWAS, APOE ɛ4 status (presence versus absence) was also included as a covariate to account for its demonstrated association. The standard conservative threshold for genome-wide significance (p < 5 × 10–8) was used in the GWAS [48], and SNPs exceeding this association threshold in the discovery cohort were analyzed for confirmatory signal (p < 0.05) in the replication cohort. Effect sizes were denoted by standardized beta coefficients. Manhattan and Q-Q plots were generated using Haploview version 4.2 [4], and regional Manhattan plots were generated using the web-based tool LocusZoom [50].

To accompany the SNP-level GWAS results, gene-based analyses were performed using H-MAGMA [60], which utilizes functional genomics evidence via tissue-specific Hi-C chromatin interactions to map SNPs to genes. We applied the purified human astrocyte Hi-C data for mapping of SNPs to genes in H-MAGMA, given the high expression of degenerative-disorder-associated genes in these cells [60]. Summary p values were generated for 50,777 Ensembl transcript IDs [22], representing 18,038 unique HGNC gene IDs, and accounting for SNP-level association statistics and gene size and density. Complementary pathway analyses were performed using the GWAS SNP-based summary statistics, GSA-SNP2 [69], and the Canonical Pathways collection version 7.2 from the Molecular Signatures Database [62]. For the pathway analyses, set size was restricted to 5–200 genes to limit the potential for size-influenced spurious associations, and the false discovery rate (FDR) was used to account for multiple comparisons [52].

Additional analyses to extend and complement the GWAS were performed on the MCSA discovery cohort with SPSS Statistics version 22.0 (IBM Corp., Armonk, NY) and SAS version 9.4 (SAS Institute Inc., Cary, NC). Post-hoc models including cycle number to account for potential learning effects were substantially unchanged from the primary results. Within the discovery sample, gene variant associations with PET and MRI biomarkers were assessed, including age (at scan), sex, APOE ɛ4 status, and genetic principal components as covariates. For volume and thickness measures from MRI (hippocampal volume, AD meta-ROI cortical thickness, and FLAIR WMH), total intracranial volume was also included as a covariate. Post-hoc analyses of SNP associations with the primary memory phenotype utilized a supramaximal threshold of global PiB SUVR ≥ 2.0 (based on previous work [35]) to sub-stratify the sample into individuals with extremely high amyloid burden versus those with abnormal amyloidosis below this threshold (1.48 ≤ global PiB SUVR < 2.0). Separate post-hoc analyses also stratified the sample based on tau PET positivity.


Sample characteristics

The discovery cohort (Table 1) included 546 individuals aged 50 years or older and amyloid PET positive from the population-based MCSA. All subjects had a baseline visit and at least one clinical follow-up with complete neuropsychological assessment data. Median follow-up time was 4.2 years (range 1–13 years) and the median number of longitudinal visits was 4 (range 2–11). Most of the sample (85%) had a diagnosis of cognitively unimpaired at the initial imaging visit, and just over half of the sample (58%) did not carry the APOE ɛ4 allele. The replication cohort, used for targeted validation of top SNPs from the discovery sample GWAS, included 545 amyloid PET positive adults aged 55 years or older from the ADNI. An additional sample of 953 amyloid PET negative older adults from the MCSA was analyzed post-hoc for comparison.

Table 1 Sample characteristics

Hypothesis-driven analyses of APOE alleles with cognitive resilience to amyloidosis

We first analyzed the APOE ɛ4 and ɛ2 alleles given the body of literature relating these to AD and cognitive decline. Age at the imaging visit, sex, years of education, regression-derived memory domain intercept (as a measure of early-life performance), and genetic principal component eigenvectors were included as covariates. As expected, APOE ɛ4 was associated with worse memory trajectory in the setting of amyloid PET positivity under both dominant (p = 0.005, standardized β = -0.12) and additive (p = 0.003, β = -0.13) genetic models. APOE ɛ2 was present in 11% of the sample and was not associated with memory trajectory under additive (p = 0.31, β = 0.04) or dominant (p = 0.28, β = 0.05) models. Based on the validated association of the APOE ɛ4 allele with the primary outcome and to prioritize discovery of novel and independent SNP associations, APOE ɛ4 status (present versus absent) was included as an additional covariate in the GWAS.

GWAS identifies novel association on chromosome 8 with cognitive resilience to amyloidosis

There was no evidence of spurious inflation of association statistics in the GWAS (λ = 1.01). A novel genome-wide significant association (p < 5 × 10–8) with cognitive resilience in the setting of amyloid PET positivity was identified on chromosome 8 (Fig. 1). The top associated SNP in this region (Fig. 2) was rs12056505 (p = 4.66 × 10–8, β = 0.23), an intronic variant in MTMR7 (myotubularin-related protein 7) that overlaps with the 3′ untranslated region of VPS37A (vacuolar protein sorting-associated protein 37A) and is nearby CNOT7 (CCR4-NOT transcription complex subunit 7). Suggestive associations (p < 5 × 10–6) were identified at additional loci, including CSMD3 (CUB and sushi multiple domains 3) which has been proposed as a regulator of dendrite development in hippocampal neurons [45], and OXCT1 (3-oxoacid CoA-transferase 1) which resides in a region found to have differential DNA methylation in dementia cases versus controls [18] (Table 2).

Fig. 1
figure 1

Manhattan plot for the GWAS of cognitive resilience to amyloidosis. The Manhattan plot displays observed −log10 p values (y-axis) for all single nucleotide polymorphisms (SNPs) tested in the GWAS of longitudinal cognitive trajectory in the setting of amyloid PET positivity. Age, sex, years of education, memory domain intercept (as a measure of premorbid performance), genetic principal components, and APOE ɛ4 status (presence versus absence) were included as covariates. An additive genetic model was utilized, along with a minor allele frequency filter of 5%. A genome-wide significant association (p < 5 × 10–8; red line) with cognitive resilience to amyloidosis was identified on chromosome 8. Suggestive associations (p < 5 × 10–6; blue line) were identified on additional chromosomes. Plot created using Haploview

Fig. 2
figure 2

Regional Manhattan plot for the GWAS hit on chromosome 8. Regional association data are displayed for the genome-wide significant hit on chromosome 8. The top association signal in this region was for rs12056505 (denoted by the purple diamond), an intronic variant in the vicinity of MTMR7 (myotubularin-related protein 7), VPS37A (vacuolar protein sorting-associated protein 37A), and CNOT7 (CCR4-NOT transcription complex subunit 7). All variants within a 500 kb region surrounding the index SNP are plotted based on their association − log10 p values, NCBI build 37 genomic position, and recombination rates calculated from the 1000 Genomes Project reference data. The color scale of r2 values is used to label SNPs based on their degree of linkage disequilibrium with the index SNP. Genes in the region are labeled below the plots, with arrows denoting 5′-to-3′ orientation. Plot created using the LocusZoom software suite

Table 2 Top independent loci from the genome-wide association study

Association of Chromosome 8 rs12056505 with resilience is specific to the amyloid-positive setting and is validated in an independent sample

The GWAS hit for rs12056505 implicated a common variant, with MAF = 12% for the C allele. This association was dose-dependent, with more positive memory trajectories observed for rs12056505-CC (mean slope = 0.18, standard error = 0.09, n = 8) and rs12056505-CT (mean slope = 0.02, standard error = 0.02, n = 131) individuals than for rs12056505-TT individuals (mean slope = -0.09, standard error = 0.02, n = 407). The association of rs12056505-C with cognitive resilience was specific to the amyloid-positive setting, as no corresponding signal was observed in a comparison cohort of amyloid-negative individuals from the MCSA (p = 0.78, β = 0.01). In analyses of an independent cohort of 545 amyloid-positive individuals from the ADNI, we identified modest evidence of replication of the association of rs12056505 with cognitive resilience to amyloidosis (p = 0.04, β = 0.09).

Post-hoc analyses in the MCSA provide context for the effect of rs12056505 on resilience

For deeper characterization of the GWAS hit, we further examined rs12056505 against relevant imaging biomarkers (Table 3). In the MCSA discovery sample, rs12056505-C was weakly associated with lower global amyloid PET burden, but its association with cognitive trajectory remained robust after additionally covarying for amyloid burden in the main model (p = 9.95 × 10–8, β = 0.22). In addition, when the sample was stratified into “extreme amyloid” (global PiB SUVR ≥ 2.0, n = 163) versus “high amyloid” (1.48 ≤ global PiB SUVR < 2.0, n = 383) groups based on prior work [35], the effect sizes for the association of rs12056505 with resilience were not substantially different (β = 0.26 in extreme amyloid group versus β = 0.22 in high amyloid group). Collectively, these results suggest that the association of rs12056505 with resilience was not purely mediated by differences in amyloid load.

Table 3 Associations of the top GWAS Hit with neuroimaging biomarkers in the MCSA

In addition, rs12056505 was not significantly associated with MRI WMH burden or DTI-based FA of the corpus callosum genu, both measures of vascular brain injury. Its association with cognitive resilience was also unchanged after additionally covarying for an index score of cerebrovascular disease risk factors [67]. Altogether, these results indicate that the association of rs12056505 with resilience was not mediated by better vascular brain health.

This SNP was also nominally associated with higher AD-signature cortical thickness and hippocampal volume, and among the small proportion of the sample with corresponding tau PET (106/546 = 19%), the association of rs12056505 with lower entorhinal cortex tau PET burden was marginally nonsignificant. However, the effect size for rs12056505 remained strong (β = 0.35) even after additionally covarying for tau PET burden. Among the 34 individuals who were A + /T + (abnormally elevated amyloid and tau by PET), the effect size for rs12056505 (β = 0.26) was comparable to that of the full sample, indicating that its protective association was not restricted to individuals without substantial AD biomarker abnormalities. In addition, for rs12056505-C carriers compared to non-carriers there were no differences in sex (49% men versus 55% men, p = 0.23) or years of education (14.5 years versus 14.7 years, p = 0.46), indicating that rs12056505 genotype was not a proxy for these factors.

Gene- and pathway-based associations with resilience to amyloidosis

To augment the SNP-level GWAS results, we used gene- and pathway-based analyses to assess for broader patterns of association signal across biologically-relevant combinations of variants. Utilizing H-MAGMA [60], the top gene-based association was for CNOT7 (p = 1.65 × 10–7), which was genome-wide significant based on a stringent Bonferroni correction (p < 0.05/50,777 = 9.85 × 10–7). Additional strong gene-based associations are listed in Table 4. Complementary pathway analysis using GSA-SNP2 [69] revealed enrichment of association with resilience for pathways related to integrin-related cell adhesion, TCF (T-cell factor) signaling related to the Wnt/β-catenin pathway, and IFN-γ (interferon gamma) signaling for immune activation, including two pathway-level associations significant after FDR correction (Table 5).

Table 4 Top HMAGMA gene-based associations in the GWAS of resilience to amyloidosis
Table 5 Top GSA-SNP2 pathway associations in the GWAS of resilience to amyloidosis

Candidate SNP associations with resilience in the GWAS

Within the GWAS, we additionally looked closely at a small set of a priori SNPs of interest, comprising AD risk loci and variants previously related to resilience phenotypes. Among 40 SNPs which previously displayed genome-wide significant associations with clinical AD diagnosis in large case–control studies through the IGAP consortium [36, 37], nominal associations (p < 0.05) were identified for three variants: the AD risk allele BIN1 rs6733839-T was associated with worse memory trajectory (p = 0.01, β = -0.11); the AD protective allele SORL1 rs11218343-C was associated with better memory trajectory (p = 0.04, β = 0.09); and the AD protective allele MEF2C rs190982-G was associated with worse memory trajectory (p = 0.03, β = -0.09). Although TOMM40 rs2075650 displayed nominal association within a basic model not accounting for APOE ɛ4 (p = 0.001, β = -0.14), this association was attenuated in the GWAS which included APOE ɛ4 as a covariate (p = 0.07, β = -0.11). Replication signal was not identified for SNPs reported in prior work to have associations with resilience proxy measures, including RAB10 (RAS oncogene family Rab10) rs142787485 (p = 0.50, β = 0.03) [53], BDNF rs6265 (p = 0.59, β = 0.02) [19], KL rs9536314 (p = 0.44, β = 0.03) [6], and ATP8B1 (ATPase phospholipid transporting 8B1) rs2571244 (p = 0.28, β = 0.05) [16].


This study leveraged a large, population-based sample of amyloid PET positive older adults who had genome-wide SNP data and serial longitudinal cognitive assessments to understand genetic factors that contribute to cognitive resilience. Through this design, we discovered a novel association with cognitive resilience to amyloidosis for a locus on chromosome 8 which was specific to the amyloid positive setting and displayed replication in an independent cohort. We also identified biological pathways with enrichment of association, including a preponderance related to immune system activation. Our data support the hypothesis that genetic heterogeneity is one of the factors that explains differential cognitive resilience to brain amyloidosis (Fig. 3).

Fig. 3
figure 3

adapted from Jack et al. [27] with the author’s permission

Conceptual model of genetic heterogeneity influencing cognitive resilience. A conceptual model displays the role of genetic heterogeneity in cognitive resilience to amyloidosis. Time is shown on the x-axis and points of maximum outcome abnormality (PET or cognitive functioning) are indicated by higher values along the y-axis. Sample cognitive trajectories, depicted by the blue curves, are shown in relation to the sigmoidal red curve depicting increasing brain PET amyloidosis over time. Genetic heterogeneity is shown as a modifier of cognitive resilience to amyloidosis, with more severe impairment (dark blue) related to a genetic risk profile and less impairment (light blue) related to a genetic protective profile. Figure

Resilience to brain amyloidosis is likely complex. It is well-understood that environmental/lifestyle factors, such as early- and late-life intellectual enrichment and cerebrovascular disease, will influence cognitive trajectories in older adults [41, 47, 64, 65]. Genetic factors are widely presumed to impact the degree of AD pathology, but relatively little of this architecture has been studied in the context of individuals with cognitive resilience in the setting of high levels of AD pathology. Further, the work so far at a population level has focused mainly on APOE. Case–control GWAS designs have been fruitful in identifying non-APOE risk variants for AD dementia, but may not be ideal to discover disease-specific resilience factors given the discrepancies between clinically diagnosed probable AD dementia versus biologically defined AD [5] as well as the heterogeneity amongst controls, some of whom may have extant AD pathophysiology while remaining non-demented at the time of study inclusion. Our approach using an easily interpretable and generalizable setting (amyloid positivity in a population-based sample of older adults) and integrating imaging biomarkers, genetics, and longitudinal clinical follow-up expands knowledge about heritable resilience factors in AD.

The top hit in our GWAS resided in MTMR7, which is ubiquitously expressed in the brain and contains a different SNP previously associated with susceptibility to Creutzfeldt-Jakob disease [57]. However, gene-based analyses of this region utilizing brain Hi-C for SNP-to-gene mapping suggested that the causative variant at this locus may instead tag CNOT7, a gene linked to synaptic plasticity and hippocampal-dependent learning and memory in model systems [44]. In addition, rs12056505 is a splicing quantitative trait locus (sQTL) for CNOT7 in cultured fibroblasts (p = 1.5 × 10–7) [10] and may be an expression quantitative trait locus (eQTL) for CNOT7 in brain cortex (p = 5.82 × 10–3) [22]. The nearby gene VPS37A has also been linked to hereditary spastic paraplegia [72] and belongs to a family of sorting proteins which may be important for tau clearance [40]. Our new findings call attention to this gene-rich region on chromosome 8 for further functional characterization, particularly given that top associated SNPs may not always represent the true functional variant at a locus.

It is still an open question whether the mechanisms underlying resilience are fundamentally mediated through differences in brain structure, metabolic maintenance, functional network compensation, or a combination of these and other avenues [3]. Emerging evidence suggests that global brain amyloid PET burden is independently (from entorhinal cortex tau PET burden and cortical thickness) associated with memory decline, but only in the setting of very high amyloid levels present for many years [35]. The association of chromosome 8 rs12056505 with cognitive protection even in the setting of extremely high amyloid burden (global PiB SUVR ≥ 2.0) suggests that its role as a resilience factor cannot be purely ascribed to upstream processes yielding a lesser burden of amyloidosis. We considered other potential explanations for the cognitive resilience observed in rs12056505-C carriers. We did not find evidence that this protective variant was associated with better vascular brain health or with significantly lower tau PET burden in key AD regions. However, the subset of our sample with tau PET data was modest, and it is possible that lower tau accumulation in the setting of amyloidosis could be a mechanism for resilience in this variant that our subsample was underpowered to detect. Non-AD pathophysiology related to TDP-43 [32] or other concomitant degenerative factors can impact cognitive trajectories and could not be directly accounted for in this study in the absence of validated in vivo biomarkers. Alternatively, other mechanisms outside of typical neurodegenerative and vascular disease pathways may be contributing to resilience, such as the immune system and inflammation-related pathways identified in our poly-SNP analyses. A more comprehensive modeling of these putative resilience mechanisms in the context of genetic, sex-related, and lifestyle mediators may facilitate enhanced preventive and therapeutic targeting in AD and related disorders. Specific to this work, Fig. 4 summarizes our overarching experimental design and the discovery of genetic protective factors (and negative testing for other related factors) that were found in the context of resilience to amyloidosis.

Fig. 4
figure 4

Experimental design for the discovery of genetic factors influencing cognitive resilience to amyloidosis. This study analyzed older individuals with significant brain amyloidosis and found evidence for genetic factors associated with cognitive resilience in that setting (green check mark). Alternative potential mediators for resilience were assessed, including intellectual enrichment, sex-related factors, resistance to tau accumulation, and cerebrovascular disease burden, and not explanatory for the genetic associations (red backslash icon). Our findings nominate new targets which warrant further study of the underlying molecular processes which impact brain structure and metabolism, functional network connectivity, and neuronal and synaptic health, to ultimately account for differential coping with amyloid pathology

Despite these strengths, there are reasons why our main finding should be interpreted with caution. The sample size of our discovery cohort was small in comparison to other GWAS of AD-relevant outcomes, and although the confirmatory association of APOE Ɛ4 with resilience in our dataset is reassuring for its broader interpretation, the possibility of a winner’s curse phenomenon regarding our GWAS hit must be acknowledged. Although similar in size and data scope to the MCSA discovery cohort, the ADNI replication sample was also not a perfect match in some characteristics, including differences in selection framework (population-based versus clinical trial sample), amyloid PET tracer, and specific cognitive outcome. In this context, the presence of a validation signal for rs12056505 is encouraging, but is still modest overall and does not rule out the possibility of a false positive. In examining published results from the largest available AD case–control GWAS [36], we cannot find evidence supporting an association of rs12056505 with lower risk of clinically diagnosed probable AD dementia (p = 0.90, β = 0.997). Clinical diagnosis of AD dementia is by no means a proxy for resilience to AD pathology, and our data more broadly mirrors conclusions from a recent study on resilience in suggesting that the genetic architectures underlying these outcomes are likely to be meaningfully different [16]. Nevertheless, the lack of clear protective relationship of rs12056505 on clinical AD dementia diagnosis in published data is a limitation. Further, cognitive performance is impacted by numerous factors and as such may not be as specific as fluid- or imaging-based AD-relevant quantitative endophenotypes. In summary, our top discovery is promising but additional validation studies in other cohorts and molecular and functional characterization are needed.

This work has several other limitations. For discovery we leveraged a rich dataset from a population-based sample, which offers the benefit of generalizability to the broader setting of older adults. However, this setting may meaningfully differ from clinical trial samples (including our replication cohort) where a higher proportion of participants have AD biomarker abnormalities and/or cognitive impairment at baseline. An advantage of our study design was that it utilized a well-validated AD biomarker (amyloid PET) and linear regression to model cognitive trajectories, approaches that can be straightforwardly applied to other datasets. However, the optimal approach to modeling resilience is still an active question, and alternative methods based on latent variables (estimating cognitive dysfunction in cross-sectional data over and above that expected from pathology) [16], gene expression data [21], and longitudinal linear mixed models [8] may provide complementary information. In addition, we assessed differential cognitive resilience through the lens of a continuous measure of rate of decline, and acknowledge that our findings may not be applicable to a paradigm whereby cognitive resilience is operationalized as a marker of clinical status (i.e., coping without impairment). We also focused on memory performance as an outcome given its relevance to the setting of clinically typical AD, but recognize that cognitive-genetic associations may well be domain-specific [31] related to the involvement of diverse brain regions and functional networks, and that there may be weaknesses in the use of z-scores for longitudinal studies of cognitive trajectory across individuals who may start from different levels of cognitive ability [11]. Further, we acknowledge that the respective outcome measures for cognitive resilience in the discovery and replication cohorts were not identical, and that future efforts incorporating our study design to additional large samples may benefit from phenotype harmonization to facilitate meta-analysis.

In summary, this study of a population-based sample discovered a novel gene variant and biological pathways associated with cognitive resilience to brain amyloidosis. Additional study of these putative resilience mechanisms may be vital for improved risk stratification and drug targeting for individuals with biomarker-confirmed AD pathophysiology.

Availability of data and materials

Data from this study is available from the authors upon reasonable request.


  1. Arboleda-Velasquez JF, Lopera F, O’Hare M, Delgado-Tirado S, Marino C, Chmielewska N, Saez-Torres KL, Amarnani D, Schultz AP, Sperling RA et al (2019) Resistance to autosomal dominant Alzheimer’s disease in an APOE3 Christchurch homozygote: a case report. Nat Med 25:1680–1683.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Arenaza-Urquijo EM, Przybelski SA, Lesnick TL, Graff-Radford J, Machulda MM, Knopman DS, Schwarz CG, Lowe VJ, Mielke MM, Petersen RC et al (2019) The metabolic brain signature of cognitive resilience in the 80+: beyond Alzheimer pathologies. Brain 142:1134–1147.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Arenaza-Urquijo EM, Vemuri P (2018) Resistance vs resilience to Alzheimer disease: clarifying terminology for preclinical studies. Neurology 90:695–703.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Barrett JC, Fry B, Maller J, Daly MJ (2005) Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 21:263–265.

    Article  CAS  PubMed  Google Scholar 

  5. Beach TG, Monsell SE, Phillips LE, Kukull W (2012) Accuracy of the clinical diagnosis of Alzheimer disease at National Institute on Aging Alzheimer Disease Centers, 2005–2010. J Neuropathol Exp Neurol 71:266–273.

    Article  PubMed  Google Scholar 

  6. Belloy ME, Napolioni V, Han SS, Le Guen Y, Greicius MD, Alzheimer’s Disease Neuroimaging I (2020) Association of klotho-VS heterozygosity with risk of Alzheimer disease in individuals who carry APOE4. JAMA Neurol 77:849–862.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Brookmeyer R, Abdalla N (2018) Estimation of lifetime risks of Alzheimer’s disease dementia using biomarkers for preclinical disease. Alzheimers Dement 14:981–988.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Buchman AS, Yu L, Boyle PA, Schneider JA, De Jager PL, Bennett DA (2016) Higher brain BDNF gene expression is associated with slower cognitive decline in older adults. Neurology 86:735–741.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ (2015) Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience 4:7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Consortium GT (2013) The Genotype-Tissue Expression (GTEx) project. Nat Genet 45:580–585.

    Article  CAS  Google Scholar 

  11. Crane PK, Narasimhalu K, Gibbons LE, Mungas DM, Haneuse S, Larson EB, Kuller L, Hall K, van Belle G (2008) Item response theory facilitated cocalibrating cognitive tests and reduced bias in estimated rates of decline. J Clin Epidemiol 61(1018–1027):e1019.

    Article  Google Scholar 

  12. Das S, Forer L, Schonherr S, Sidore C, Locke AE, Kwong A, Vrieze SI, Chew EY, Levy S, McGue M et al (2016) Next-generation genotype imputation service and methods. Nat Genet 48:1284–1287.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Davies G, Harris SE, Reynolds CA, Payton A, Knight HM, Liewald DC, Lopez LM, Luciano M, Gow AJ, Corley J et al (2014) A genome-wide association study implicates the APOE locus in nonpathological cognitive ageing. Mol Psychiatry 19:76–87.

    Article  CAS  PubMed  Google Scholar 

  14. De Jager PL, Shulman JM, Chibnik LB, Keenan BT, Raj T, Wilson RS, Yu L, Leurgans SE, Tran D, Aubin C et al (2012) A genome-wide scan for common variants affecting the rate of age-related cognitive decline. Neurobiol Aging 33(1017):e1011-1015.

    Article  CAS  Google Scholar 

  15. Debette S, Ibrahim Verbaas CA, Bressler J, Schuur M, Smith A, Bis JC, Davies G, Wolf C, Gudnason V, Chibnik LB et al (2015) Genome-wide studies of verbal declarative memory in nondemented older people: the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium. Biol Psychiatry 77:749–763.

    Article  PubMed  Google Scholar 

  16. Dumitrescu L, Mahoney ER, Mukherjee S, Lee ML, Bush WS, Engelman CD, Lu Q, Fardo DW, Trittschuh EH, Mez J et al (2020) Genetic variants and functional pathways associated with resilience to Alzheimer’s disease. Brain 143:2561–2575.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Fatemi F, Kantarci K, Graff-Radford J, Preboske GM, Weigand SD, Przybelski SA, Knopman DS, Machulda MM, Roberts RO, Mielke MM et al (2018) Sex differences in cerebrovascular pathologies on FLAIR in cognitively unimpaired elderly. Neurology 90:e466–e473.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Fransquet PD, Lacaze P, Saffery R, Phung J, Parker E, Shah R, Murray A, Woods RL, Ryan J (2020) Blood DNA methylation signatures to detect dementia prior to overt clinical symptoms. Alzheimer’s Dement (Amsterdam, Netherlands) 12:e12056.

    Article  Google Scholar 

  19. Franzmeier N, Ren J, Damm A, Monte-Rubio G, Boada M, Ruiz A, Ramirez A, Jessen F, Duzel E, Rodriguez Gomez O et al (2019) The BDNFVal66Met SNP modulates the association between beta-amyloid and hippocampal disconnection in Alzheimer’s disease. Mol Psychiatry.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Hixson JE, Vernier DT (1990) Restriction isotyping of human apolipoprotein E by gene amplification and cleavage with HhaI. J Lipid Res 31:545–548

    Article  CAS  PubMed  Google Scholar 

  21. Hohman TJ, Dumitrescu L, Cox NJ, Jefferson AL, Alzheimer’s Neuroimaging I (2017) Genetic resilience to amyloid related cognitive decline. Brain Imaging Behav 11:401–409.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Hunt SE, McLaren W, Gil L, Thormann A, Schuilenburg H, Sheppard D, Parton A, Armean IM, Trevanion SJ, Flicek P et al (2018) Ensembl variation resources. Database (Oxford).

    Article  PubMed Central  Google Scholar 

  23. Huq AJ, Fransquet P, Laws SM, Ryan J, Sebra R, Masters CL, Winship IM, James PA, Lacaze P (2019) Genetic resilience to Alzheimer’s disease in APOE epsilon4 homozygotes: a systematic review. Alzheimers Dement 15:1612–1623.

    Article  PubMed  Google Scholar 

  24. Jack CR Jr, Knopman DS, Jagust WJ, Petersen RC, Weiner MW, Aisen PS, Shaw LM, Vemuri P, Wiste HJ, Weigand SD et al (2013) Tracking pathophysiological processes in Alzheimer’s disease: an updated hypothetical model of dynamic biomarkers. Lancet Neurol 12:207–216.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Jack CR Jr, Wiste HJ, Therneau TM, Weigand SD, Knopman DS, Mielke MM, Lowe VJ, Vemuri P, Machulda MM, Schwarz CG et al (2019) Associations of amyloid, tau, and neurodegeneration biomarker profiles with rates of memory decline among individuals without dementia. JAMA 321:2316–2325.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Jack CR Jr, Wiste HJ, Weigand SD, Therneau TM, Lowe VJ, Knopman DS, Gunter JL, Senjem ML, Jones DT, Kantarci K et al (2017) Defining imaging biomarker cut points for brain aging and Alzheimer’s disease. Alzheimers Dement 13:205–216.

    Article  PubMed  Google Scholar 

  27. Jack CR, Knopman DS, Jagust WJ, Shaw LM, Aisen PS, Weiner MW, Petersen RC, Trojanowski JQ (2010) Hypothetical model of dynamic biomarkers of the Alzheimer’s pathological cascade. Lancet Neurol 9:119–128.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Jack CR, Wiste HJ, Weigand SD, Rocca WA, Knopman DS, Mielke MM, Lowe VJ, Senjem ML, Gunter JL, Preboske GM et al (2014) Age-specific population frequencies of cerebral beta-amyloidosis and neurodegeneration among people with normal cognitive function aged 50–89 years: a cross-sectional study. Lancet Neurol 13:997–1005

    Article  PubMed  PubMed Central  Google Scholar 

  29. Jagust WJ, Bandy D, Chen K, Foster NL, Landau SM, Mathis CA, Price JC, Reiman EM, Skovronsky D, Koeppe RA et al (2010) The Alzheimer’s disease neuroimaging initiative positron emission tomography core. Alzheimers Dement 6:221–229.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Jansen WJ, Ossenkoppele R, Knol DL, Tijms BM, Scheltens P, Verhey FR, Visser PJ, Amyloid Biomarker Study G, Aalten P, Aarsland D et al (2015) Prevalence of cerebral amyloid pathology in persons without dementia: a meta-analysis. JAMA 313:1924–1938.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Kamboh MI, Fan KH, Yan Q, Beer JC, Snitz BE, Wang X, Chang CH, Demirci FY, Feingold E, Ganguli M (2019) Population-based genome-wide association study of cognitive decline in older adults free of dementia: identification of a novel locus for the attention domain. Neurobiol Aging 84:239 e215-239 e224.

    Article  Google Scholar 

  32. Kapasi A, Yu L, Boyle PA, Barnes LL, Bennett DA, Schneider JA (2020) Limbic-predominant age-related TDP-43 encephalopathy (LATE), ADNC pathology, and cognitive decline in aging. Neurology.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Klunk WE, Engler H, Nordberg A, Wang Y, Blomqvist G, Holt DP, Bergstrom M, Savitcheva I, Huang GF, Estrada S et al (2004) Imaging brain amyloid in Alzheimer’s disease with Pittsburgh Compound-B. Ann Neurol 55:306–319.

    Article  CAS  PubMed  Google Scholar 

  34. Knopman DS, Beiser A, Machulda MM, Fields J, Roberts RO, Pankratz VS, Aakre J, Cha RH, Rocca WA, Mielke MM et al (2015) Spectrum of cognition short of dementia: Framingham Heart Study and Mayo Clinic Study of Aging. Neurology 85:1712–1721.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Knopman DS, Lundt ES, Therneau TM, Vemuri P, Lowe VJ, Kantarci K, Gunter JL, Senjem ML, Mielke MM, Machulda MM et al (2019) Entorhinal cortex tau, amyloid-beta, cortical thickness and memory performance in non-demented subjects. Brain 142:1148–1160.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Kunkle BW, Grenier-Boley B, Sims R, Bis JC, Damotte V, Naj AC, Boland A, Vronskaya M, van der Lee SJ, Amlie-Wolf A et al (2019) Genetic meta-analysis of diagnosed Alzheimer’s disease identifies new risk loci and implicates Abeta, tau, immunity and lipid processing. Nat Genet 51:414–430.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Lambert JC, Ibrahim-Verbaas CA, Harold D, Naj AC, Sims R, Bellenguez C, DeStafano AL, Bis JC, Beecham GW, Grenier-Boley B et al (2013) Meta-analysis of 74,046 individuals identifies 11 new susceptibility loci for Alzheimer’s disease. Nat Genet 45:1452–1458.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Landau SM, Breault C, Joshi AD, Pontecorvo M, Mathis CA, Jagust WJ, Mintun MA, Alzheimer’s Disease Neuroimaging I (2013) Amyloid-beta imaging with Pittsburgh compound B and florbetapir: comparing radiotracers and quantification methods. J Nucl Med 54:70–77.

    Article  CAS  PubMed  Google Scholar 

  39. Landau SM, Mintun MA, Joshi AD, Koeppe RA, Petersen RC, Aisen PS, Weiner MW, Jagust WJ (2012) Amyloid deposition, hypometabolism, and longitudinal cognitive decline. Ann Neurol 72:578–586.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Li JG, Chiu J, Pratico D (2020) Full recovery of the Alzheimer’s disease phenotype by gain of function of vacuolar protein sorting 35. Mol Psychiatry 25:2630–2640.

    Article  CAS  PubMed  Google Scholar 

  41. Livingston G, Huntley J, Sommerlad A, Ames D, Ballard C, Banerjee S, Brayne C, Burns A, Cohen-Mansfield J, Cooper C et al (2020) Dementia prevention, intervention, and care: 2020 report of the Lancet Commission. Lancet 396:413–446.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Lowe VJ, Bruinsma TJ, Wiste HJ, Min HK, Weigand SD, Fang P, Senjem ML, Therneau TM, Boeve BF, Josephs KA et al (2019) Cross-sectional associations of tau-PET signal with cognition in cognitively unimpaired adults. Neurology 93:e29–e39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. McCarthy S, Das S, Kretzschmar W, Delaneau O, Wood AR, Teumer A, Kang HM, Fuchsberger C, Danecek P, Sharp K et al (2016) A reference panel of 64,976 haplotypes for genotype imputation. Nat Genet 48:1279–1283.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. McFleder RL, Mansur F, Richter JD (2017) Dynamic control of dendritic mRNA expression by CNOT7 regulates synaptic efficacy and higher cognitive function. Cell Rep 20:683–696.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Mizukami T, Kohno T, Hattori M (2016) CUB and Sushi multiple domains 3 regulates dendrite development. Neurosci Res 110:11–17.

    Article  CAS  PubMed  Google Scholar 

  46. Ossenkoppele R, Lyoo CH, Jester-Broms J, Sudre CH, Cho H, Ryu YH, Choi JY, Smith R, Strandberg O, Palmqvist S et al (2020) Assessment of demographic, genetic, and imaging variables associated with brain resilience and cognitive resilience to pathological tau in patients with Alzheimer disease. JAMA Neurol 77:632–642.

    Article  PubMed  Google Scholar 

  47. Oveisgharan S, Wilson RS, Yu L, Schneider JA, Bennett DA (2020) Association of early-life cognitive enrichment with Alzheimer disease pathological changes and cognitive decline. JAMA Neurol.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Pe’er I, Yelensky R, Altshuler D, Daly MJ (2008) Estimation of the multiple testing burden for genomewide association studies of nearly all common variants. Genet Epidemiol 32:381–385.

    Article  PubMed  Google Scholar 

  49. Petersen RC, Roberts RO, Knopman DS, Geda YE, Cha RH, Pankratz VS, Boeve BF, Tangalos EG, Ivnik RJ, Rocca WA (2010) Prevalence of mild cognitive impairment is higher in men. The Mayo Clinic Study of Aging. Neurology 75:889–897.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, Boehnke M, Abecasis GR, Willer CJ (2010) LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 26:2336–2337.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ et al (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81:559–575.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Ramanan VK, Shen L, Moore JH, Saykin AJ (2012) Pathway analysis of genomic data: concepts, methods, and prospects for future development. Trends Genet: TIG 28:323–332.

    Article  CAS  PubMed  Google Scholar 

  53. Ridge PG, Karch CM, Hsu S, Arano I, Teerlink CC, Ebbert MTW, Gonzalez Murcia JD, Farnham JM, Damato AR, Allen M et al (2017) Linkage, whole genome sequence, and biological data implicate variants in RAB10 in Alzheimer’s disease resilience. Genome Med 9:100.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Roberts RO, Geda YE, Knopman DS, Cha RH, Pankratz VS, Boeve BF, Ivnik RJ, Tangalos EG, Petersen RC, Rocca WA (2008) The Mayo Clinic Study of Aging: design and sampling, participation, baseline measures and sample characteristics. Neuroepidemiology 30:58–69.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Roberts RO, Knopman DS, Syrjanen JA, Aakre JA, Vassilaki M, Kremers WK, Mielke MM, Machulda MM, Graff-Radford J, Geda YE et al (2017) Weighting and standardization of frequencies to determine prevalence of AD imaging biomarkers. Neurology 89:2039–2048.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Rocca WA, Yawn BP, St. Sauver JL, Grossardt BR, Melton LJ, (2012) History of the Rochester Epidemiology Project: half a century of medical records linkage in a US population. Mayo Clin Proc 87:1202–1213

    Article  PubMed  PubMed Central  Google Scholar 

  57. Sanchez-Juan P, Bishop MT, Aulchenko YS, Brandel JP, Rivadeneira F, Struchalin M, Lambert JC, Amouyel P, Combarros O, Sainz J et al (2012) Genome-wide study links MTMR7 gene to variant Creutzfeldt-Jakob risk. Neurobiol Aging 33(1487):e1421-1488.

    Article  CAS  Google Scholar 

  58. Saykin AJ, Shen L, Xiaohui Y, Kim S, Nho K, Risacher SL, Ramanan VK, Foroud TM, Faber KM, Sarwar N et al (2015) Genetic studies of quantitative MCI and AD phenotypes in ADNI: progress, opportunities, and plan. Alzheimer’s Dement 11(7):792–814

    Article  Google Scholar 

  59. Schwarz CG, Gunter JL, Wiste HJ, Przybelski SA, Weigand SD, Ward CP, Senjem ML, Vemuri P, Murray ME, Dickson DW et al (2016) A large-scale comparison of cortical thickness and volume methods for measuring Alzheimer’s disease severity. NeuroImage Clin 11:802–812.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Sey NYA, Hu B, Mah W, Fauni H, McAfee JC, Rajarajan P, Brennand KJ, Akbarian S, Won H (2020) A computational tool (H-MAGMA) for improved prediction of brain-disorder risk genes by incorporating brain chromatin interaction profiles. Nat Neurosci 23:583–593.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. St Sauver JL, Grossardt BR, Yawn BP, Melton LJ 3rd, Pankratz JJ, Brue SM, Rocca WA (2012) Data resource profile: the Rochester Epidemiology Project (REP) medical records-linkage system. Int J Epidemiol 41:1614–1624.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES et al (2005) Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA 102:15545–15550.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Veitch DP, Weiner MW, Aisen PS, Beckett LA, Cairns NJ, Green RC, Harvey D, Jack CR Jr, Jagust W, Morris JC et al (2019) Understanding disease progression and improving Alzheimer’s disease clinical trials: Recent highlights from the Alzheimer’s Disease Neuroimaging Initiative. Alzheimers Dement 15:106–152.

    Article  PubMed  Google Scholar 

  64. Vemuri P (2018) “Exceptional brain aging” without Alzheimer’s disease: triggers, accelerators, and the net sum game. Alzheimer’s Res Ther 10:53–53.

    Article  CAS  Google Scholar 

  65. Vemuri P, Lesnick TG, Knopman DS, Przybelski SA, Reid RI, Mielke MM, Graff-Radford J, Lowe VJ, Machulda MM, Petersen RC et al (2019) Amyloid, vascular, and resilience pathways associated with cognitive aging. Ann Neurol.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Vemuri P, Lesnick TG, Przybelski SA, Graff-Radford J, Reid RI, Lowe VJ, Zuk SM, Senjem ML, Schwarz CG, Gunter JL et al (2018) Development of a cerebrovascular magnetic resonance imaging biomarker for cognitive aging. Ann Neurol 84:705–716.

    Article  PubMed  PubMed Central  Google Scholar 

  67. Vemuri P, Lesnick TG, Przybelski SA, Knopman DS, Lowe VJ, Graff-Radford J, Roberts RO, Mielke MM, Machulda MM, Petersen RC et al (2017) Age, vascular health, and Alzheimer disease biomarkers in an elderly sample. Ann Neurol 82:706–718.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Weiner MW, Aisen PS, Jack CR Jr, Jagust WJ, Trojanowski JQ, Shaw L, Saykin AJ, Morris JC, Cairns N, Beckett LA et al (2010) The Alzheimer’s disease neuroimaging initiative: progress report and future plans. Alzheimers Dement 6(202–211):e207.

    Article  Google Scholar 

  69. Yoon S, Nguyen HCT, Yoo YJ, Kim J, Baik B, Kim S, Kim J, Kim S, Nam D (2018) Efficient pathway enrichment and network analysis of GWAS summary data using GSA-SNP2. Nucleic Acids Res 46:e60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Zhang C, Pierce BL (2014) Genetic susceptibility to accelerated cognitive decline in the US Health and Retirement Study. Neurobiol Aging 35(1512):e1511-1518.

    Article  Google Scholar 

  71. Zheng X, Levine D, Shen J, Gogarten SM, Laurie C, Weir BS (2012) A high-performance computing toolset for relatedness and principal component analysis of SNP data. Bioinformatics 28:3326–3328.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Zivony-Elboum Y, Westbroek W, Kfir N, Savitzki D, Shoval Y, Bloom A, Rod R, Khayat M, Gross B, Samri W et al (2012) A founder mutation in Vps37A causes autosomal recessive complex hereditary spastic paraparesis. J Med Genet 49:462–472.

    Article  CAS  PubMed  Google Scholar 

Download references


The authors thank all of the study participants and staff in the Mayo Clinic Study of Aging, Mayo Alzheimer’s Disease Research Center, and Mayo Clinic Aging and Dementia Imaging Research laboratory. Figures 3 and 4 created with


This work was supported by NIH grants U01 AG006786 (PI: Petersen/Mielke/Jack), R01 NS097495 (PI: Vemuri), R01 AG56366 (PI: Vemuri), P50 AG016574 (PI: Petersen), P30 AG062677 (PI: Petersen), R37 AG011378 (PI: Jack), R01 AG041851 (PIs: Jack and Knopman), RF1 AG55151 (PI: Mielke), U54 NS100693 (PI: Ross), and R01 AG034676 (PI: Rocca); the GHR Foundation, the Alexander Family Alzheimer’s Disease Research Professorship of the Mayo Clinic, the Mayo Foundation for Medical Education and Research, the Liston Award, the Elsie and Marvin Dekelboum Family Foundation, the Schuler Foundation, and Opus Building NIH grant C06 RR018898. Data collection and sharing for the ADNI data utilized in this project was funded by the ADNI NIH grant U01 AG024904, other funding through the the National Institute of Biomedical Imaging and Bioengineering, and private sector contributions from the following (facilitated by the Foundation for the National Institutes of Health with the grantee organization as the Northern California Institute for Research and Education): AbbVie, Alzheimer’s Association; Alzheimer’s Drug Discovery Foundation; Araclon Biotech; BioClinica, Inc.; Biogen; Bristol-Myers Squibb Company; CereSpir, Inc.; Cogstate; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; EuroImmun; F. Hoffmann-La Roche Ltd and its affiliated company Genentech, Inc.; Fujirebio; GE Healthcare; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research & Development, LLC.; Johnson & Johnson Pharmaceutical Research & Development LLC.; Lumosity; Lundbeck; Merck & Co., Inc.; Meso Scale Diagnostics, LLC.; NeuroRx Research; Neurotrack Technologies; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Piramal Imaging; Servier; Takeda Pharmaceutical Company; and Transition Therapeutics. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. The study is coordinated by the Alzheimer’s Therapeutic Research Institute at the University of Southern California. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of Southern California. We would like to greatly thank AVID Radiopharmaceuticals, Inc., for their support in supplying AV-1451 precursor, chemistry production advice and oversight, and FDA regulatory cross-filing permission and documentation needed for this work. The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Author information

Authors and Affiliations




Study conception: VKR, OAR, and PV. Study design: VKR, TGL, SAP, MGH, OAR, and PV. Data collection: SAP, DSK, JGR, VJL, MMM1, MMM2, CRJ, RCP, OAR, and PV. Data analysis: VKR, TGL, and SAP. Drafting of the manuscript: VKR. Critical revision of the manuscript: TGL, SAP, MGH, DSK, JGR, VJL, MMM1, MMM2, CRJ, RCP, OAR, and PV. All authors read and approved the final manuscript

Corresponding authors

Correspondence to Vijay K. Ramanan or Prashanthi Vemuri.

Ethics declarations

Ethics approval and consent to participate

All study protocols were approved by each participating site’s Institutional Review Board, and written informed consent was obtained from all study participants or their surrogates.

Consent for publication

No individual patient-identifying data is presented in this manuscript.

Competing interests

Dr. Knopman reported serving on a data safety monitoring board for the DIAN study, serving as an investigator in clinical trials sponsored by Biogen, Lilly Pharmaceuticals, and the University of Southern California, and receiving research support from the National Institutes of Health (NIH) outside of the submitted work. Dr. Graff-Radford reported receiving research support from the National Institute on Aging outside the submitted work. Dr. Lowe reported receiving consulting fees from Bayer Schering Pharma, Piramal Life Sciences, and Merck Research and grants from GE Healthcare, Siemens Molecular Imaging, AVID Radiopharmaceuticals, and the NIH outside the submitted work. Dr. Mielke reported receiving research support from the NIH, Department of Defense, and unrestricted research grants from Biogen outside the submitted work. Dr. Jack reported consulting for Eli Lilly, serving on an independent data monitoring board for Roche, and serving as a speaker for Eisai but receives no personal compensation from any commercial entity; he also reported receiving research support from the NIH and the Alexander Family Alzheimer's Disease Research Professorship of the Mayo Clinic. Dr. Petersen reported receiving consulting fees from Hoffman-La Roche Inc, Merck Inc, Genentech Inc, Biogen Inc, GE Healthcare, and Eisai Inc. outside the submitted work. Dr. Vemuri reported receiving grants from the NIH during the conduct of the study. No other disclosures were reported.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.



Data used in preparation of this article were obtained from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) database ( As such, the investigators within the ADNI contributed to the design and implementation of ADNI and/or provided data but did not participate in analysis or writing of this report. A complete listing of ADNI investigators can be found at:

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ramanan, V.K., Lesnick, T.G., Przybelski, S.A. et al. Coping with brain amyloid: genetic heterogeneity and cognitive resilience to Alzheimer’s pathophysiology. acta neuropathol commun 9, 48 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: