Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Intermediate Phenotypes Identify Divergent Pathways to Alzheimer's Disease

  • Joshua M. Shulman,

    Affiliations Program in Translational NeuroPsychiatric Genomics, Department of Neurology, Brigham and Women's Hospital, Boston, Massachusetts, United States of America, Harvard Medical School, Boston, Massachusetts, United States of America, Program in Medical and Population Genetics, Broad Institute, Cambridge, Massachusetts, United States of America

  • Lori B. Chibnik,

    Affiliations Program in Translational NeuroPsychiatric Genomics, Department of Neurology, Brigham and Women's Hospital, Boston, Massachusetts, United States of America, Harvard Medical School, Boston, Massachusetts, United States of America, Program in Medical and Population Genetics, Broad Institute, Cambridge, Massachusetts, United States of America

  • Cristin Aubin,

    Affiliations Program in Translational NeuroPsychiatric Genomics, Department of Neurology, Brigham and Women's Hospital, Boston, Massachusetts, United States of America, Harvard Medical School, Boston, Massachusetts, United States of America, Program in Medical and Population Genetics, Broad Institute, Cambridge, Massachusetts, United States of America

  • Julie A. Schneider,

    Affiliations Department of Pathology, Rush University Medical Center, Chicago, Illinois, United States of America, Rush Alzheimer's Disease Center, Department of Neurological Sciences, Rush University Medical Center, Chicago, Illinois, United States of America

  • David A. Bennett,

    Affiliation Rush Alzheimer's Disease Center, Department of Neurological Sciences, Rush University Medical Center, Chicago, Illinois, United States of America

  • Philip L. De Jager

    pdejager@rics.bwh.harvard.edu

    Affiliations Program in Translational NeuroPsychiatric Genomics, Department of Neurology, Brigham and Women's Hospital, Boston, Massachusetts, United States of America, Harvard Medical School, Boston, Massachusetts, United States of America, Program in Medical and Population Genetics, Broad Institute, Cambridge, Massachusetts, United States of America

Abstract

Background

Recent genetic studies have identified a growing number of loci with suggestive evidence of association with susceptibility to Alzheimer's disease (AD). However, little is known of the role of these candidate genes in influencing intermediate phenotypes associated with a diagnosis of AD, including cognitive decline or AD neuropathologic burden.

Methods/Principal Findings

Thirty-two single nucleotide polymorphisms (SNPs) previously implicated in AD susceptibility were genotyped in 414 subjects with both annual clinical evaluation and completed brain autopsies from the Religious Orders Study and the Rush Memory and Aging Project. Regression analyses evaluated the relation of SNP genotypes to continuous measures of AD neuropathology and cognitive function proximate to death. A SNP in the zinc finger protein 224 gene (ZNF224, rs3746319) was associated with both global AD neuropathology (p = 0.009) and global cognition (p = 0.002); whereas, a SNP at the phosphoenolpyruvate carboxykinase locus (PCK1, rs8192708) was selectively associated with global cognition (p = 3.57×10−4). The association of ZNF224 with cognitive impairment was mediated by neurofibrillary tangles, whereas PCK1 largely influenced cognition independent of AD pathology, as well as Lewy bodies and infarcts.

Conclusions/Significance

The findings support the association of several loci with AD, and suggest how intermediate phenotypes can enhance analysis of susceptibility loci in this complex genetic disorder.

Introduction

Alzheimer's disease (AD), the most common cause of dementia, leads to progressive loss of memory and other cognitive domains, and is characterized pathologically by the accumulation of extracelluar amyloid plaques and intracellular neurofibrillary tangles. AD likely develops from an interaction of numerous genes along with environmental risk factors, each with modest and incompletely penetrant effects. Linkage studies have identified rare gene mutations as causal in familial, early age-of-onset AD, but these Mendelian variants only explain a small fraction of disease burden in the general population [1]. The identification of susceptibility loci for sporadic, late age-of-onset AD has been more challenging, with numerous reports of candidate gene associations, most of which have not been consistently replicated in follow-up studies [2][4]. One notable exception is the apolipoprotein E locus (APOE): the ε4 allele is common, increases AD susceptibility 3-fold, and is estimated to explain at least 10% of the population-attributable risk of disease [1]. In addition, the APOE ε2 allele is a validated AD protective allele though it is less common, and its effect size is more modest than ε4.

Genome-wide association (GWA) studies have emerged as a promising approach to identify susceptibility loci in common diseases with complex genetic inheritance, but until recently, most GWA scans in AD have been relatively underpowered, and identified loci have not been consistently replicated [5][14]. Increasing sample size is one approach for boosting statistical power, and this strategy has recently led to the identification of several promising new AD susceptibility loci, including CR1, CLU, and PICALM [15], [16]. However, clinical heterogeneity remains a significant confounder of the case/control study design in AD, due to the likely inclusion of dementia cases with multiple pathologies, such as cerebrovascular disease or other neurodegenerative conditions. In addition, since AD develops following a protracted pre-clinical phase consisting of mild symptoms, control groups are susceptible to contamination by latent disease cases. Substantial AD pathology is often present in advanced age, including in those with minimal or no cognitive impairment at death [17]. Subjects with significant pathology but subclinical disease are likely to dilute power in an AD case/control association analysis.

One approach to overcoming these obstacles is to study quantitative intermediate phenotypes. The manifestation of the AD clinical syndrome is the final culmination of a sequence of events beginning with genetic and environmental risk factors that trigger intermediate pathological changes, synapse loss and cell death, and ultimately cognitive decline and dementia. Outcome measures selected more proximally along this causal chain are expected be less confounded and more strongly associated with susceptibility loci. In addition, compared to the dichotomous clinical diagnosis, quantitative intermediate phenotypes can capture more of the underlying heritable trait variation, further enhancing statistical power. Based on this promise, a number of studies have begun to take advantage of intermediate phenotypes for genetic association analysis in AD, including neuropsychiatric test measures [18], MRI imaging data [19], [20], biomarkers from blood and CSF [21], [22], and direct measurements of AD pathology [23]. The latter approach requires access to large study populations with detailed clinical and neuropathologic characterization. The Religious Orders Study and Rush Memory and Aging Project are prospectively following more than 2,300 older persons, all of whom have agreed to annual clinical evaluation and brain donation at death. More than 800 autopsies have been completed to date, and quantitative analyses of amyloid and tangle burden has been performed on nearly 600. In a recent study of APOE in this cohort, we found that intermediate cognitive and pathological phenotypes substantially increase power for genetic association analysis [23]. In addition, using neuropathologic phenotypes, the association between APOE and cognitive impairment was previously shown to be mediated by a sequential cascade of amyloid plaque formation and subsequent development of neurofibrillary tangle pathology [24], [25]. Therefore, beyond enhancing power for association analysis, intermediate phenotypes hold the additional promise of testing mechanistic hypotheses of gene action.

In this study, we extend our previous work to evaluate several candidate AD susceptibility loci for associations with intermediate phenotypes relevant to AD. Thirty-two candidate SNPs were selected based on their discovery in AD GWA studies and/or evidence from the AlzGene online meta-analyses [2], [26]. SNPs were genotyped in more than 400 subjects with detailed cognitive and pathological data, allowing assessment of genotype relations to quantitative AD pathology and cognitive function proximate to death. We subsequently leveraged the detailed phenotypes available in our cohorts to dissect the functional pathways that link genetic variants to cognitive impairment.

Methods

Ethics Statement

Written informed consent and an Anatomic Gift Act were signed by all Religious Orders Study and Rush Memory and Aging project participants after the procedures were fully explained, and both studies were approved by the Institutional Review Board of Rush University Medical Center. The work described in this report was additionally approved by the Institutional Review Boards of the Brigham and Women's Hospital and Massachusetts Institute of Technology.

Subjects

Clinical and post-mortem data came from participants in the Religious Orders Study and Rush Memory and Aging Project, two longitudinal, epidemiologic clinical-pathologic studies of aging and AD [17]. In both studies, participants without known dementia at baseline agreed to annual detailed clinical evaluation and brain donation at the time of death. Participants in the Religious Orders Study were older Cathololic nuns, priests and brothers from about 40 groups in 12 states across the United States. Subjects in the Rush Memory and Aging Project were older, community-dwelling persons from about 40 retirement communities and subsidized senior housing facilities across northeastern Illinois. Since 1993, more than 2,300 persons agreed to participate in these studies. The overall follow-up rate exceeds 90% of survivors and the overall autopsy rate exceeds 90% of decedents. Of those subjects with completed neuropathologic analyses, and following genotyping quality control filters, 414 persons with genotyping data were available for analysis in February of 2009 when this study was initiated (250 from the Religious Orders Study and 164 from the Rush Memory and Aging Project).

Clinical evaluation

The clinical diagnoses of dementia and AD were made each year following the recommendations of the joint working group of the National Institute of Neurologic and Communicative Disorders and Stroke and the AD and Related Disorders Association [27], as previously described in detail [28]. Probable AD refers to persons with clinical AD and no other clinical condition contributing to cognitive impairment and possible AD refers to persons meeting inclusion criteria for AD who are thought to have another condition (e.g., stroke) contributing to cognitive impairment. MCI referred to those individuals rated as cognitively impaired by the neuropsychologist but not demented by the examining physician, as previously described [29]. At the time of death, clinical data were reviewed by a neurologist without access to post-mortem data and a summary diagnostic opinion was rendered regarding the most likely clinical diagnosis at the time of death. Level of cognition was based on cognitive testing performed proximate to death. The Religious Orders Study and Rush Memory and Aging Project have 19 cognitive performance tests in common, and use identical analytic procedures to develop summary statistics. Mini-Mental State Examination [30] was used to describe the cohort and one test was used for diagnostic classification purposes only. The remaining 17 tests have been previously described [17]. Tests were converted to z scores, using the mean and SD from the baseline evaluation of all participants, and averaged to yield summary measures of global cognition and five cognitive domains: episodic memory, semantic memory, working memory, perceptual speed, and visuospatial ability. Summary measures minimize floor and ceiling effects and other sources of random variability. For the mediation analyses incorporating diagnosis of diabetes, annual clinical evaluations allowed documentation of history of diabetes and use of medications to treat diabetes. Diabetes was determined to be present if the participant was ever taking a medication, such as insulin or an oral hypoglycemic, to treat diabetes, as determined by direct inspection of medication containers, or ever reported a history of diagnosis of diabetes, or both, as previously described [31].

Neuropathological evaluation

Brain autopsies were performed across the US as previously described [17]. Bielschowsky silver stain was used to visualize neuritic plaques, diffuse plaques, and neurofibrillary tangles in tissue sections from the midfrontal, middle temporal, inferior parietal, and entorhinal cortices, and the hippocampal CA1 sector. The neuropathologic diagnosis of AD was made by a board-certified neuropathologist without access to any clinical data as previously reported [17], [28]. We classified persons as having pathologic AD based on intermediate or high likelihood of AD by National Institute on Aging (NIA)-Reagan criteria using CERAD estimates of neuritic plaque density and Braak staging of neurofibrillary pathology [32][34], as previously described [17]. The quantitative composite AD pathology score was based on counts of neuritic plaques, diffuse plaques and neurofibrillary tangles as previously described [35], [36]. Because the means, standard deviations, and ranges of the data varied widely for the pathologic indices, we converted the raw counts to a standard distribution by dividing each person's count by the standard deviation for that particular count and formed a summary measure by averaging the scaled scores. Because the data were skewed, square root of the scaled score was used in analyses. Separate summary measures of neurofibrillary tangles and neuritic and diffuse plaques were also made. Chronic macroscopic cerebral infarctions and alpha-synuclein immunoreactive Lewy bodies were determined as previously described and considered present or absent for analyses [17].

Genotyping

DNA was extracted from lymphocytes or frozen post-mortem brain tissue. APOE genotyping was performed by Agencourt Bioscience Corporation (Beverly, MA) utilizing high throughput sequencing of codon 112 (position 3937) and codon 158 (position 4075) of exon 4 of the APOE gene on chromosome 19. In addition to the APOE ε4 and ε2 alleles, 32 SNPs were selected for genotyping in our cohort, based on prior evidence from the literature, as of February, 2009. Thus, the more recently discovered CR1, CLU, and PICALM [15], [16] loci were not included in this study, but are the focus of a separate study (Chibnik et al., submitted). The selected SNPs were equally divided between the top results of AD case/control GWA studies [5][7], [9], [10][14] (16 SNPs) and candidate gene association studies (16 SNPs), which were chosen based on their top ranking in AlzGene meta-analyses [2]. The 32 candidate SNPs were genotyped using matrix-assisted laser desorption-ionization time-of-flight mass spectrometry on a MassARRAY platform (Sequenom). After excluding subjects for failed genotyping exceeding the 10% threshold, 414 individuals remained for subsequent analysis (genotyping rate in these subjects was >99%). All SNP allele frequencies satisfied Hardy-Weinberg equilibrium (p>0.001). Allele frequencies were not significantly different between the Religious Orders Study and Rush Memory and Aging Project subjects, supporting the validity of pooled analyses.

Statistical Analysis

Given the complementary study designs and similar procedures for data collection and generation of the cognitive and neuropathologic outcome traits, we pooled data from the Religious Orders Study and Rush Memory and Aging Project for our analyses, consistent with numerous prior studies [17], [23]. Genetic association was performed using the PLINK analysis software toolkit [37]. Linear regression was used to evaluate the association of allele genotypes with level of cognition proximate to death in a 2-degree-of-freedom, genotypic test of association, with covariates included for age, gender, and years of education. In order to refine the genetic model, we additionally tested selected SNPs using a 1-degree-of-freedom test to examine for additive, dominant, or recessive allelic effects. These studies were performed using PLINK as well as the R statistical computing program (www.r-project.org). Linear regression modeling in R was used to calculate residual quantitative trait variance explained, and to perform statistical mediation analyses. For the case/control association analysis based on AD clinical diagnosis, logistic regression was performed in PLINK under both additive and genotypic models, and again including covariates for age, gender, and education. All p-values reported are unadjusted for multiple hypothesis testing. A Bonferroni-corrected significance threshold of p<0.001 was calculated for the 34 SNPs tested for associations with our two primary outcomes, a quantitative measure of global cognitive performance and global AD pathology. Given the high correlation between the pathologic and cognitive traits, applying an adjustment for 68 tests would be overly conservative. Otherwise, the threshold of p<0.01 was selected to indicate suggestive statistical evidence of association. All other evaluated phenotypes, including those for AD clinical diagnosis, pathology sub-types, and cognitive domains were considered secondary analyses.

Results

Associations with global AD pathology and global cognition

Subject demographics, clinical and neuropathologic diagnoses, cognitive status, and APOE genotypes for the cohort analyzed in this study are presented in Table 1. In clinical evaluations proximate to death, of the 414 subjects in our study cohort, 131 (31.6%) had normal cognition, 98 (23.7%) had mild cognitive impairment, and 185 (44.7%) were demented (173 met criteria for possible or probable AD). As expected, a significant proportion (41.5%) of individuals without dementia satisfied NIA-Reagan pathological criteria for intermediate or high likelihood AD.

thumbnail
Table 1. Demographic, clinical and pathological characteristics of the study cohort.

https://doi.org/10.1371/journal.pone.0011244.t001

We initially tested for associations between each of the 34 polymorphisms and our two primary outcomes, intermediate phenotypes representing a measure of global AD pathologic burden on autopsy and a measure of global cognitive function proximate to death (Table 2). Linear regression models were used to examine the relation of SNP genotypes to the quantitative neuropathologic and cognitive traits in a 2 degree-of-freedom statistical test, adjusting for the effects of age at death, gender, and years of education. As expected, APOE ε4 was significantly associated with both cognition (p = 3.4×10−10) and AD pathology (p = 1.6×10−24) in our cohort, whereas an association with APOE ε2 was only seen for the pathological phenotype (p = 9.1×10−4). In addition, we found associations with AD intermediate phenotypes for two SNPs, within the zinc finger protein 224 (ZNF224) and phosphoenolpyruvate carboxykinase 1 (PCK1) genes, both of which were selected for genotyping based on their identification in AD case/control GWA studies [6], [10]. The ZNF224 SNP (rs3746319) was associated with both global cognition (p = 0.009) and global AD pathology (p = 0.004). In contrast, the PCK1 SNP (rs8192708) was significantly associated with global cognition (p = 3.57×10−4) but not global AD pathology (p = 0.056), suggesting that this locus may influence cognitive impairment through mechanisms other than AD pathology. Besides APOE ε4, none of the SNP associations surpass the currently accepted threshold for genome-wide significance (p<5.0×10−8); however, the association between PCK1 and global cognition exceeds a Bonferroni-corrected significance threshold of p<0.001 for 34 independent tests. Given the high correlation between the pathologic and cognitive traits, applying an adjustment for 68 tests would be overly conservative; however, the PCK1 association still exceeds that standard (p<7×10−4).

thumbnail
Table 2. Relation of candidate AD polymorphisms to intermediate cognitive and pathologic phenotypes.

https://doi.org/10.1371/journal.pone.0011244.t002

Although the risk alleles for the associations of both the ZNF224 and PCK1 loci with the intermediate phenotypes in our cohort also increase risk of AD diagnosis (Table 3), their effects are opposite to that reported in the original GWA studies [6], [10]. In the case of ZNF224, we find that the minor allele, rs3746319A, is associated with both increased AD pathologic burden and decreased cognitive performance; whereas this variant was protective against AD in the GWA study (G. Beecham and M. Pericak-Vance, personal communication). Similarly, for PCK1, the minor allele, rs8192708G, significantly protected against cognitive decline in our cohort but was in fact the AD risk allele in the original GWA study [2]. Interestingly, two subsequent replication analyses of rs8192708 documented associations of decreased AD risk with the minor allele, consistent with our findings [4], [38]. Therefore, while the effects of the ZNF224 and PCK1 loci on a diagnosis of AD and on intermediate phenotypes are consistent within our study (and in two other PCK1 replication studies); they are not consistent with the original GWA analyses. In the discussion section, we further address possible explanations for these discrepancies.

thumbnail
Table 3. Relation of candidate AD polymorphisms to clinical AD diagnosis.

https://doi.org/10.1371/journal.pone.0011244.t003

Associations with neuritic and diffuse plaques, neurofibrillary tangles, and cognitive subdomains

The global AD pathology score averages the post-mortem density of neuritic and diffuse plaques and neurofibrillary tangles in multiple brain regions; however, we hypothesized that certain AD susceptibility loci might selectively promote one type of pathology, in which case the composite pathologic outcome might dilute statistical power to detect associations. We therefore performed secondary analyses to determine whether any of the candidate SNPs tested demonstrate selective or more robust association signals with separate quantitative measures of plaques or tangle pathology (Table 4). All analyses were again performed using linear regression models to test for associations with SNP genotypes, adjusted for the effects of age, gender, and education. A SNP at the GALP locus (rs3745833) showed suggestive evidence for association with diffuse plaques (p = 0.003), but not with neurofibrillary tangles (p = 0.373). In contrast, the ZNF224 SNP (rs3746319) was strongly associated with neurofibrillary tangle burden (p = 1.49×10−4), whereas no significant association was seen with either neuritic plaque (p = 0.018) or diffuse plaque (p = 0.290) pathology. Therefore, the association with the tangle subscore is likely the primary driver for the ZNF224 locus association with global AD pathology (p = 0.009), and the composite score appears to dilute statistical power. Interestingly, the PCK1 SNP (rs8192708), which was not associated with the global pathology measure, did show suggestive evidence for association with neuritic plaque pathology (p = 0.007); however, this did not appear to explain the strong association with global cognition (p = 3.57×10−4), as investigated further below.

thumbnail
Table 4. Relation of polymorphisms to amyloid plaques and neurofibrillary tangles.

https://doi.org/10.1371/journal.pone.0011244.t004

Similar to the approach taken with the pathological phenotypes, we performed secondary analyses to assess whether any of the SNPs showed more robust association with the five cognitive subdomains that comprise the global cognition score. Linear regression was again used to test for association of each SNP with separate quantitative trait outcomes representing episodic memory, semantic memory, working memory, perceptual speed, and visuospatial ability (Table 5). Episodic memory impairment, the most characteristic cognitive deficit of AD, was associated with both the ZNF224 locus (p = 0.003) and the PCK1 locus (p = 3.69×10−4). ZNF224 was additionally associated with decline in visuospatial function (p = 0.007), and PCK1 showed evidence for association with semantic memory impairment (p = 0.001).

thumbnail
Table 5. Relation of polymorphisms to measures of cognitive performance.

https://doi.org/10.1371/journal.pone.0011244.t005

Divergent pathways from genes to cognitive impairment

For the SNPs at the ZNF224 and PCK1 loci, we performed additional linear regression analyses to refine the genetic model for the relation with our intermediate phenotypes (additive, dominant, or recessive), better characterize the strength of the observed effects, and develop statistical models to test hypotheses about mechanistic pathways. Our core regression model, consisting of age at death, gender, and years of education, explained 3% and 7% of the variation in our pathological and cognitive traits, respectively. Using the optimal dominant model of inheritance, the ZNF224 SNP (rs3746319) explained an additional 2% (Beta = 0.13, p = 0.003) of the residual variance in global AD pathology and 2.1% (Beta = −0.39, p = 0.002) of the variance in global cognition (Table 6). We next explored whether the effect of this locus on AD pathology might mediate its association with cognition (Table 7). When a term for global AD pathology was incorporated in our linear regression model, the magnitude of the association between ZNF224 and global cognition was attenuated by 44% (Beta = −0.22, p = 0.05). In our analyses of the neuropathologic subtypes, we found that the association of ZNF224 on AD pathology appeared to be due to a predominant effect on neurofibrillary tangles. Indeed, when we substituted a term for neurofibrillary tangles instead of the global pathology variable in our regression model, the effect of the ZNF224 variant on global cognition was reduced by 64%, and was no longer significant (Beta = −0.14, p = 0.21), whereas tangles showed a robust association with cognitive impairment (p<2×10−16). These results are consistent with a sequence of events whereby an effect on the formation of neurofibrillary tangles accounts for the association of the ZNF224 allele with cognitive function.

thumbnail
Table 6. Detailed genotype-phenotype data and statistical modeling.

https://doi.org/10.1371/journal.pone.0011244.t006

thumbnail
Table 7. Distinct pathways of ZNF224 and PCK1 association with cognition.

https://doi.org/10.1371/journal.pone.0011244.t007

In contrast to ZNF224, the PCK1 locus showed a relatively selective association with global cognition, but not with global pathology. For this locus, an additive model of inheritance was a best fit for our data, and PCK1 explained 3.4% (Beta = 0.49, p = 1.02×10−4) of the variance in global cognition proximate to death. We next used multiple linear regression to test whether the association of this SNP on cognitive impairment is predominantly independent of AD pathology. Indeed, after inclusion of a model term for global AD pathology, the PCK1 SNP (rs8192708) remained associated with global cognition (Beta = 0.37, p = 0.002), despite the strong, independent association between pathology and cognition (Beta = −1.22, p<2×10−16; Table 7). Given our finding of an association with neuritic plaques, we substituted a model term for neuritic pathology instead of the global pathology variable; however, despite a modest reduction in the effect size, the relation between PCK1 and cognitive impairment remained significant (Beta = 0.31, p = 0.005; Table 7). Besides AD-related pathology, Lewy bodies and infarcts are the two additional brain pathologies most commonly seen in association with age-related cognitive decline [39]. We therefore investigated whether the association of PCK1 on global cognition might be mediated by either Lewy bodies or infarcts, by including relevant terms into our regression model (Table 8). Again, the association between PCK1 and global cognition remained significant (Beta = 0.39, p = 2.04×10−4), and PCK1 continued to explain 3% of the residual variance in global cognition in our cohort after adjusting for the three most common brain pathologies associated with dementia. Common variation at the PCK1 locus has also been associated with type 2 diabetes in a number of independent studies [40][42]. Since diabetes has also been implicated as a risk factor for age-related cognitive decline [31], we examined whether diabetes mediates the association of PCK1 with cognitive impairment. However, adjusting for diabetes diagnosis in our linear regression model did not significantly attenuate the association of PCK1 with global cognition (Beta = 0.370, p = 0.001).

thumbnail
Table 8. PCK1 association with cognition is largely independent of AD pathology, infarcts, and Lewy bodies.

https://doi.org/10.1371/journal.pone.0011244.t008

Discussion

In this study, by genotyping a panel of loci within two cohorts of subjects with detailed cognitive and neuropathological characterization, we evaluate intermediate phenotypes as a tool for the functional dissection of candidate AD susceptibility loci. Using the extensively validated APOE locus, we previously demonstrated that intermediate traits enhance statistical power to detect associations, even in studies of modest sample size. Here, using the same strategy, we present evidence supporting the possible role of two additional loci in influencing age-related cognitive decline and AD neuropathology. Specifically, the ZNF224 locus is associated with a quantitative measure of global AD pathology, and both ZNF224 and PCK1 are associated with a summary measure of global cognition proximate to death. Using separate quantitative traits for each of the predominant AD pathological features, we document associations between GALP and PCK1 and diffuse and neuritic plaque pathology, respectively, whereas ZNF224 showed a relatively selective association with neurofibrillary tangle pathology. Finally, in a series of statistical mediation analyses, we tested hypotheses about the causal chain of events linking genetic variation in the ZNF224 and PCK1 loci with cognitive decline, with strikingly different outcomes. In the case of ZNF224, we find that AD pathology, and more specifically, neurofibrillary tangles mediate an association with cognitive impairment. In contrast, we find that the association between PCK1 and cognition is largely independent of not only AD pathology, but also Lewy bodies, and infarcts, which together comprise the three most common known brain pathologies associated with dementia [39], [43].

Both ZNF224 and PCK1 were initially implicated by AD GWA studies; however, neither locus has yet been consistently replicated in subsequent genetic studies, and little is known about their potential mechanism of action in disease pathogenesis. The ZNF224 locus encodes a Kruppel-associated box-containing zinc-finger protein that is widely expressed, including in the adult brain, and likely functions as a transcriptional repressor [44], [45]. The SNP evaluated in this study, rs3746319, encodes a missense mutation causing a Lys to Glu change at position 640, which falls near the C-terminus within one of 19 zinc-finger repeat motifs. However, we do not yet know enough about ZNF224 protein structure and function to speculate further on how this variant might promote neurofibrillary tangle formation and subsequent cognitive impairment, and further investigation will be required to determine if rs3746319 is the causal variant and whether ZNF224 is indeed the causal gene. The PCK1 gene encodes phosphoenolpyruvate carboxykinase 1, which catalyzes the rate-limiting step of gluconeogenesis [46]. The SNP genotyped in our study, rs8192708, is also a missense mutation, causing an Ile to Val change at position 267; however, the functional consequences of this change, if any, are not known. PCK1 variants have also been suggested to be associated with diabetes [40][42], and independently, diabetes has been identified as a risk factor for the development of dementia [31]. In our mediation analysis, adjusting for the effect of diabetes diagnosis did not account for the association of PCK1 and cognitive impairment; however, it is possible that an appropriate intermediate phenotype, such as direct measurements of blood glucose or hemoglobin A1c, might allow detection of mediation. In another study performed in the same cohort, a relation was found between diabetes and infarcts [47]; however, we were also unable to mediate the PCK1 association by including a model term for cerebral infarctions. Our finding that the PCK1 association with cognitive decline is not explained by AD pathology, Lewy bodies, or infarcts suggests that this locus might influence additional, unmeasured pathologies. For example, whereas our analyses adjusted for macroscopic infarcts, PCK1 may instead primarily influence microscopic forms of cerebrovascular injury. Further, while our intermediate pathologic phenotype accounts for amyloid plaques and neurofibrillary tangles, it does not capture levels of soluble, but potentially still neurotoxic, forms of amyloid or tau pathology [48], [49]. Alternatively, variation at PCK1 might influence one or multiple steps in the cascade of events predicted to occur downstream of amyloid, tangles and other pathologies, such as synapse loss, inflammation, and/or cell death pathways.

Unexpectedly, the variants in ZNF224 and PCK1 show opposite directions of allelic effects for association with AD intermediate phenotypes in our cohort compared to their association with AD diagnosis in the initial GWA studies. In other words, the alleles associated with increased AD risk in the initial reports (rs3746319G and rs8192708G) are actually protective against cognitive decline in our cohort. Importantly, this discrepancy is not accounted for by our use of intermediate phenotypes, as the ZNF224 and PCK1 SNPs show consistent direction of affect on AD diagnosis in our study population (Table 3). Such “flip-flop” associations have been reported with increasing frequency as GWA scans are completed for many common diseases, and replication efforts are subsequently undertaken [50]. Indeed, in the case of PCK1, two prior replication studies found evidence that the major allele, rs8192708A, may increase risk for dementia, consistent with our results suggesting an association between this allele and both cognitive decline and AD [4], [38].

The interpretation of reversals in the direction of variant associations between different study cohorts remains controversial [50]. The most common explanation for such observations are that they are in fact spurious and representative of chance fluctuations around the null hypothesis. However, in our study, the strongly suggestive statistical evidence for the associations between ZNF224 and PCK1 with AD intermediate phenotypes makes their arising by chance less likely; and additionally, the reversals of allelic effect are seen with both loci in our analysis. Instead, we propose that differences in subject ascertainment and recruitment are more likely to be responsible for our observations. The Religious Orders Study (ROS) and Rush Memory and Aging Project (MAP), from which our study cohort is based, are prospective, longitudinal studies in which subjects from the community are recruited non-demented at baseline (mean age  = 75 and 79, for ROS and MAP respectively). All cases of clinical AD are therefore incident within our cohort. In contrast, similar to nearly all AD GWA studies performed to date, the initial reports of association with the PCK1 and ZNF224 loci come from AD cases recruited from a neurology clinic population with prevalent dementia. In addition, whereas subjects in our study were recruited at approximately similar ages to the GWA cohorts, they were significantly older at the time of last clinical evaluation and autopsy (mean age of death  = 87). Studies with different designs (cross-sectional vs. prospective) and varying methods of subject ascertainment can generate contradictory epidemiological findings, for example due to survival bias. If an AD risk allele is associated with earlier age of dementia onset; it might be under-represented in the prospective cohort, which requires subjects to be non-demented at enrollment; and therefore, might subsequently appear to be a protective allele. “Flip-flop” associations might additionally arise from variation in linkage disequilibrium structure in the genomic region of interest between the cohorts in different studies. In fact, both the ZNF224 and PCK1 SNPs fall under modest recombination peaks, based on HapMap data [51]. Although both our study and the GWA analyses were conducted in subjects of European ancestry, it remains possible that sampling variation between two populations of similar ethnicity might lead to the association reversal that we have observed, as recombination could distribute our tag SNP onto haplotypes that are different from that harboring the causal variant [50]. Ultimately, further analysis of both SNPs and fine mapping of each locus in larger study samples will be required to validate both PCK1and ZNF224 as AD susceptibility loci, and resolve which allele may increase risk for disease.

Of the thirty-four SNPs evaluated in our study, both of the loci that we found to be associated with AD intermediate phenotypes were initially identified by GWA studies, suggesting the power of this unbiased approach to identify genes that might be overlooked by prevailing hypotheses of disease biology. Our study was initiated prior to the recent report of two large AD case/control GWA studies which independently identified three new susceptibility genes, CLU, CR1, and PICALM [15], [16]. In a parallel effort, we recently found that CR1 is associated with age-related cognitive decline in our study cohorts; and further, that this association was mediated by an effect on amyloid pathology (Chibnik et al., submitted). The power of a GWA study design and the types of genes one expects to discover are tightly linked to the selected phenotypic outcome. To the extent possible, the chosen outcome measure should be closely matched to the underlying biology responsible for the heritable trait variation of interest. In autopsy cohort studies of aged individuals in the community setting, most subjects with probable AD demonstrate multiple brain pathologies [52]. Based on our results, we believe that intermediate pathological and cognitive traits have great promise to enhance gene discovery and for functional characterization of loci that emerge from current AD GWA studies.

Author Contributions

Conceived and designed the experiments: JMS JAS DB PLDJ. Performed the experiments: JMS CA. Analyzed the data: JMS LBC. Contributed reagents/materials/analysis tools: CA JAS DB. Wrote the paper: JMS LBC JAS DB PLDJ.

References

  1. 1. Ertekin-Taner N (2007) Genetics of Alzheimer's disease: a centennial review. Neurol Clin 25: 611–667.
  2. 2. Bertram L, McQueen M, Mullin K, Blacker D, Tanzi R (2007) Systematic meta-analyses of Alzheimer disease genetic association studies: the AlzGene database. Nat Genet 39: 17–23.
  3. 3. Cousin E, Macé S, Rocher C, Dib C, Muzard G, et al. (2009) No replication of genetic association between candidate polymorphisms and Alzheimer's disease. Neurobiology of Aging. in press.
  4. 4. Feulner T, Laws S, Friedrich P, Wagenpfeil S, Wurst S, et al. (2009) Examination of the current top candidate genes for AD in a genome-wide association study. Mol Psychiatry. in press.
  5. 5. Abraham R, Moskvina V, Sims R, Hollingworth P, Morgan A, et al. (2008) A genome-wide association study for late-onset Alzheimer's disease using DNA pooling. BMC medical genomics 1: 44.
  6. 6. Beecham GW, Martin ER, Li Y-J, Slifer MA, Gilbert JR, et al. (2009) Genome-wide association study implicates a chromosome 12 risk locus for late-onset Alzheimer disease. Am J Hum Genet 84: 35–43.
  7. 7. Bertram L, Lange C, Mullin K, Parkinson M, Hsiao M, et al. (2008) Genome-wide association analysis reveals putative Alzheimer's disease susceptibility loci in addition to APOE. Am J Hum Genet 83: 623–632.
  8. 8. Bertram L, Tanzi RE (2009) Genome-wide association studies in Alzheimer's disease. Human Molecular Genetics 18: R137–145.
  9. 9. Carrasquillo MM, Zou F, Pankratz VS, Wilcox SL, Ma L, et al. (2009) Genetic variation in PCDH11X is associated with susceptibility to late-onset Alzheimer's disease. Nat Genet 41: 192–198.
  10. 10. Grupe A, Abraham R, Li Y, Rowland C, Hollingworth P, et al. (2007) Evidence for novel susceptibility genes for late-onset Alzheimer's disease from a genome-wide association study of putative functional variants. Hum Mol Genet 16: 865–873.
  11. 11. Li H, Wetten S, Li L St, Jean PL, Upmanyu R, et al. (2008) Candidate Single-Nucleotide Polymorphisms From a Genomewide Association Study of Alzheimer Disease. Archives of Neurology 65: 45–53.
  12. 12. Li Y, Grupe A, Rowland C, Nowotny P, Kauwe J, et al. (2006) DAPK1 variants are associated with Alzheimer's disease and allele-specific expression. Hum Mol Genet 15: 2560–2568.
  13. 13. Liu F, Arias-Vásquez A, Sleegers K, Aulchenko YS, Kayser M, et al. (2007) A genomewide screen for late-onset Alzheimer disease in a genetically isolated Dutch population. Am J Hum Genet 81: 17–31.
  14. 14. Reiman E, Webster J, Myers A, Hardy J, Dunckley T, et al. (2007) GAB2 alleles modify Alzheimer's risk in APOE epsilon4 carriers. Neuron 54: 713–720.
  15. 15. Harold D, Abraham R, Hollingworth P, Sims R, Gerrish A, et al. (2009) Genome-wide association study identifies variants at CLU and PICALM associated with Alzheimer's disease. Nat Genet 41: 1088–1093.
  16. 16. Lambert J-C, Heath S, Even G, Campion D, Sleegers K, et al. (2009) Genome-wide association study identifies variants at CLU and CR1 associated with Alzheimer's disease. Nat Genet 41: 1094–1099.
  17. 17. Bennett D, Schneider J, Arvanitakis Z, Kelly J, Aggarwal N, et al. (2006) Neuropathology of older persons without cognitive impairment from two community-based studies. Neurology 66: 1837–1844.
  18. 18. Mcqueen MB, Bertram L, Lange C, Becker KD, Albert MS, et al. (2007) Exploring candidate gene associations with neuropsychological performance. Am J Med Genet 144B: 987–991.
  19. 19. Potkin SG, Guffanti G, Lakatos A, Turner JA, Kruggel F, et al. (2009) Hippocampal atrophy as a quantitative trait in a genome-wide association study identifying novel susceptibility genes for Alzheimer's disease. PLoS ONE 4: e6501.
  20. 20. Seshadri S, DeStefano A, Au R, Massaro J, Beiser A, et al. (2007) Genetic correlates of brain aging on MRI and cognitive test measures: a genome-wide association and linkage analysis in the Framingham Study. BMC Med Genet 8: Suppl 1S15.
  21. 21. Papassotiropoulos A, Streffer JR, Tsolaki M, Schmid S, Thal D, et al. (2003) Increased brain beta-amyloid load, phosphorylated tau, and risk of Alzheimer disease associated with an intronic CYP46 polymorphism. Archives of Neurology 60: 29–35.
  22. 22. Peskind E, Li G, Shofer J, Quinn J, Kaye J, et al. (2006) Age and apolipoprotein E*4 allele effects on cerebrospinal fluid beta-amyloid 42 in adults with normal cognition. Arch Neurol 63: 936–939.
  23. 23. Bennett DA, De Jager PL, Leurgans SE, Schneider JA (2009) Neuropathologic intermediate phenotypes enhance association to Alzheimer susceptibility alleles. Neurology 72: 1495–1503.
  24. 24. Bennett D, Schneider J, Wilson R, Bienias J, Berry-Kravis E, et al. (2005) Amyloid mediates the association of apolipoprotein E e4 allele to cognitive function in older people. J Neurol Neurosurg Psychiatry 76: 1194–1199.
  25. 25. Mortimer JA, Snowdon DA, Markesbery WR (2009) The effect of APOE-epsilon4 on dementia is mediated by Alzheimer neuropathology. Alzheimer Dis Assoc Disord 23: 152–157.
  26. 26. Bertram L, Tanzi RE (2008) Thirty years of Alzheimer's disease genetics: the implications of systematic meta-analyses. Nat Rev Neurosci 9: 768–778.
  27. 27. McKhann G, Drachman D, Folstein M, Katzman R, Price D, et al. (1984) Clinical diagnosis of Alzheimer's disease: report of the NINCDS-ADRDA Work Group under the auspices of Department of Health and Human Services Task Force on Alzheimer's Disease. Neurology 34: 939–944.
  28. 28. Bennett DA, Schneider JA, Aggarwal NT, Arvanitakis Z, Shah RC, et al. (2006) Decision rules guiding the clinical diagnosis of Alzheimer's disease in two community-based cohort studies compared to standard practice in a clinic-based cohort study. Neuroepidemiology 27: 169–176.
  29. 29. Bennett DA, Wilson RS, Schneider JA, Evans DA, Beckett LA, et al. (2002) Natural history of mild cognitive impairment in older persons. Neurology 59: 198–205.
  30. 30. Folstein MF, Folstein SE, McHugh PR (1975) “Mini-mental state”. A practical method for grading the cognitive state of patients for the clinician. J Psychiatr Res 12: 189–198.
  31. 31. Arvanitakis Z, Wilson RS, Bienias JL, Evans DA, Bennett DA (2004) Diabetes mellitus and risk of Alzheimer disease and decline in cognitive function. Archives of Neurology 61: 661–666.
  32. 32. The National Institute on Aging and Reagan Institute Working Group on Diagnostic Criteria for the Neuropathological Assessment of Alzheimer's Disease (1997) Consensus recommendations for the postmortem diagnosis of Alzheimer's disease. Neurobiol Aging 18: S1–2.
  33. 33. Braak H, Braak E (1991) Neuropathological stageing of Alzheimer-related changes. Acta Neuropathol 82: 239–259.
  34. 34. Mirra SS, Heyman A, McKeel D, Sumi SM, Crain BJ, et al. (1991) The Consortium to Establish a Registry for Alzheimer's Disease (CERAD). Part II. Standardization of the neuropathologic assessment of Alzheimer's disease. Neurology 41: 479–486.
  35. 35. Bennett DA, Schneider JA, Tang Y, Arnold SE, Wilson RS (2006) The effect of social networks on the relation between Alzheimer's disease pathology and level of cognitive function in old people: a longitudinal cohort study. Lancet Neurol 5: 406–412.
  36. 36. Bennett DA, Wilson RS, Schneider JA, Evans DA, Aggarwal NT, et al. (2003) Apolipoprotein E epsilon4 allele, AD pathology, and the clinical expression of Alzheimer's disease. Neurology 60: 246–252.
  37. 37. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira M, et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81: 559–575.
  38. 38. Figgins JA, Minster RL, Demirci FY, Dekosky ST, Kamboh MI (2009) Association studies of 22 candidate SNPs with late-onset Alzheimer's disease. Am J Med Genet B Neuropsychiatr Genet 150B: 520–526.
  39. 39. Sonnen JA, Larson EB, Crane PK, Haneuse S, Li G, et al. (2007) Pathological correlates of dementia in a longitudinal, population-based sample of aging. Ann Neurol 62: 406–413.
  40. 40. Cao H, van der Veer E, Ban MR, Hanley AJG, Zinman B, et al. (2004) Promoter polymorphism in PCK1 (phosphoenolpyruvate carboxykinase gene) associated with type 2 diabetes mellitus. J Clin Endocrinol Metab 89: 898–903.
  41. 41. Hamilton G, Proitsi P, Jehu L, Morgan A, Williams J, et al. (2007) Candidate gene association study of insulin signaling genes and Alzheimer's disease: Evidence for SOS2, PCK1, and PPARγas susceptibility loci. Am J Med Genet 144B: 508–516.
  42. 42. Willer CJ, Bonnycastle LL, Conneely KN, Duren WL, Jackson AU, et al. (2007) Screening of 134 single nucleotide polymorphisms (SNPs) previously associated with type 2 diabetes replicates association with 12 SNPs in nine genes. Diabetes 56: 256–264.
  43. 43. Schneider J, Arvanitakis Z, Bang W, Bennett D (2007) Mixed brain pathologies account for most dementia cases in community-dwelling older persons. Neurology 69: 2197–2204.
  44. 44. Medugno L, Florio F, Cesaro E, Grosso M, Lupo A, et al. (2007) Differential expression and cellular localization of ZNF224 and ZNF255, two isoforms of the Krüppel-like zinc-finger protein family. Gene 403: 125–131.
  45. 45. Medugno L, Florio F, De Cegli R, Grosso M, Lupo A, et al. (2005) The Krüppel-like zinc-finger protein ZNF224 represses aldolase A gene transcription by interacting with the KAP-1 co-repressor protein. Gene 359: 35–43.
  46. 46. Yang J, Kalhan SC, Hanson RW (2009) What is the metabolic role of phosphoenolpyruvate carboxykinase? J Biol Chem 284: 27025–27029.
  47. 47. Arvanitakis Z, Schneider JA, Wilson RS, Li Y, Arnold SE, et al. (2006) Diabetes is related to cerebral infarction but not to AD pathology in older persons. Neurology 67: 1960–1965.
  48. 48. McLean CA, Cherny RA, Fraser FW, Fuller SJ, Smith MJ, et al. (1999) Soluble pool of Abeta amyloid as a determinant of severity of neurodegeneration in Alzheimer's disease. Ann Neurol 46: 860–866.
  49. 49. Santacruz K, Lewis J, Spires T, Paulson J, Kotilinek L, et al. (2005) Tau suppression in a neurodegenerative mouse model improves memory function. Science 309: 476–481.
  50. 50. Lin P-I, Vance JM, Pericak-Vance MA, Martin ER (2007) No gene is an island: the flip-flop phenomenon. Am J Hum Genet 80: 531–538.
  51. 51. Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, et al. (2007) A second generation human haplotype map of over 3.1 million SNPs. Nature 449: 851–861.
  52. 52. Schneider JA, Arvanitakis Z, Leurgans SE, Bennett DA (2009) The neuropathology of probable Alzheimer disease and mild cognitive impairment. Ann Neurol 66: 200–208.